Jumpstarting the cytochrome P450 catalytic cycle with a hydrated electron

Cytochrome P450cam (CYP101Fe3+) regioselectively hydroxylates camphor. Possible hydroxylating intermediates in the catalytic cycle of this well-characterized enzyme have been proposed on the basis of experiments carried out at very low temperatures and shunt reactions, but their presence has not yet been validated at temperatures above 0 °C during a normal catalytic cycle. Here, we demonstrate that it is possible to mimic the natural catalytic cycle of CYP101Fe3+ by using pulse radiolysis to rapidly supply the second electron of the catalytic cycle to camphor-bound CYP101[FeO2]2+. Judging by the appearance of an absorbance maximum at 440 nm, we conclude that CYP101[FeOOH]2+ (compound 0) accumulates within 5 μs and decays rapidly to CYP101Fe3+, with a k440 nm of 9.6 × 104 s−1. All processes are complete within 40 μs at 4 °C. Importantly, no transient absorbance bands could be assigned to CYP101[FeO2+por•+] (compound 1) or CYP101[FeO2+] (compound 2). However, indirect evidence for the involvement of compound 1 was obtained from the kinetics of formation and decay of a tyrosyl radical. 5-Hydroxycamphor was formed quantitatively, and the catalytic activity of the enzyme was not impaired by exposure to radiation during the pulse radiolysis experiment. The rapid decay of compound 0 enabled calculation of the limits for the Gibbs activation energies for the conversions of compound 0 → compound 1 → compound 2 → CYP101Fe3+, yielding a ΔG‡ of 45, 39, and 39 kJ/mol, respectively. At 37 °C, the steps from compound 0 to the iron(III) state would take only 4 μs. Our kinetics studies at 4 °C complement the canonical mechanism by adding the dimension of time.


REACTION 1
The electrode potential of this half-reaction is estimated to be ϩ1.5 V at pH 7 (3).
For the best characterized cytochrome P450 isozyme, CYP101, it has been shown that the supply of electrons from putidaredoxin to cytochrome P450 is rate-limiting (15). Hence, it is not possible to observe with standard spectroscopic methods reactive intermediates that occur after the formation of the [FeO 2 ] 2ϩ intermediate. However, the hydrated electron (e aq Ϫ ), which can be generated by the pulse radiolysis technique, reduces proteins at diffusion-controlled rates on the microsecond time scale. Under pseudo-first order conditions, with CYP101 in excess, CYP101Fe 3ϩ reacts with e aq Ϫ with a rate constant of 3 ϫ 10 10 M Ϫ1 s Ϫ1 ; the yield of reduced heme is ϳ10% relative to e aq Ϫ (16,17). Kobayashi et al. (18) used pulse radiolysis to add an electron to the oxy-form of 2,4-diacetyldeutero-cytochrome P450cam, and the difference absorption spectrum observed was assigned to Cpd 2; however, these authors were unable to analyze formation of hydroxylated camphor; thus, no conclusions could be drawn about the catalytic competence of Cpd 2.
We set out to study the kinetics of the reaction of e aq Ϫ with the Fe 3ϩ and [FeO 2 ] 2ϩ forms of camphor-hydroxylating cytochrome P450, CYP101, by UV-visible spectroscopy, with the expectation to observe reactions 2 and 3. We show that the pulse radiolysis technique can be used to "jump-start" the catalytic cycle of cytochrome P450 by providing the second electron to reduce the oxy-form of the enzyme. Although we find no spectral evidence for accumulation of Cpd 1, our evidence indicates that it must be formed.

Half-life of the hydrated electron
The half-life, t1 ⁄ 2 , of e aq Ϫ decay, measured after irradiation (100 Gy) of a corresponding N 2 -saturated 10 mM potassium phosphate buffer at pH 7.4 in the presence of 10 M camphor but without CYP101, is 0.65 s ( Fig. 2A).

Reduction of CYP101Fe 3؉ by hydrated electrons
The rate constant for the reaction of the e aq Ϫ with the heme of CYP101Fe 3ϩ was determined by competition kinetics, i.e. via the yield of this reaction. The yield of reduction was determined by irradiating a N 2 -saturated solution of 2.0 M CYP101Fe 3ϩ containing 2 mM t-BuOH at a dose of 80 Gy, which results in the formation of 22 M e aq Ϫ . The observed difference absorbance spectrum (Fig. 2B) is consistent with conversion of CYP101Fe 3ϩ to CYP101Fe 2ϩ . The absorbance at 390 nm, which was used to follow the reduction of CYP101Fe 3ϩ (Fig. 2C), decreased by a total of ϳ65 mAbs in 50 s. However, over the course of 5 half-lives of e aq Ϫ decay (ϳ3 s), the corresponding loss of absorbance at 390 nm was ϳ40 mAbs, which, given the extinction coefficients of the oxidized and reduced proteins at that wavelength (2) and the 6-cm path length, corresponds to a concentration of CYP101Fe 2ϩ of 0.18 M. This first phase of absorbance change reflects reduction of ϳ9% of the heme iron(III) relative to protein and a yield of reduction of the heme relative to e aq Ϫ of 0.8%. The same result was found when traces recorded at 410, 415, 420, and 425 nm were analyzed. The second stage of heme reduction, the loss of the remaining 25 mAbs at 390 nm, occurred after essentially all of the e aq Ϫ decayed, proceeding for ϳ50 s, with k obs_reduction ϭ 5.4 ϫ 10 4 s Ϫ1 (Fig. 2C). This stage of heme reduction is likely caused by adducts of e aq Ϫ and H ⅐ with amino acids (21,22) formed during the first 3 s, with reducing equivalents then transferred to the heme by a tunneling or hopping mechanism.

Rate constant for the formation of CYP101[FeO 2 ] 2؉
To validate pulse radiolysis as the method of choice to provide reducing equivalents, we followed the formation of CYP101[FeO 2 ] 2ϩ , as shown in the following reaction, by reducing CYP101Fe 3ϩ with e aq Ϫ in the presence of excess O 2 . The rate constant measured, k 4 ϭ 2.1 ϫ 10 5 M Ϫ1 s Ϫ1 (Fig. S1a), agrees with that determined by stopped flow, 2.4 ϫ 10 5 M Ϫ1 s Ϫ1 (Fig. S1b), and is the same order of magnitude as the literature value, 7.7 ϫ 10 5 M Ϫ1 s Ϫ1 (19). Moreover, these experiments demonstrate that, at 4°C, CYP101[FeO 2 2ϩ ] does not autoxidize significantly for at least 1 min (Fig. S1c).

The reaction of hydrated electrons with CYP101[FeO 2 ] 2؉
To observe intermediates of the catalytic cycle, we formed substrate-bound CYP101[FeO 2 ] 2ϩ at 4°C in the irradiation cell by mixing an anaerobic solution of CYP101Fe 2ϩ , produced by reducing CYP101Fe 3ϩ in the presence of camphor with a slight excess of S 2 O 4 2Ϫ inside a N 2 -filled glove box, with a buffered solution containing O 2 . The CYP101[FeO 2 ] 2ϩ formed can autoxidize; thus, the mixture was irradiated as soon as possible, generally within 30 s. Two distinct absorption bands with maxima at 410 and 440 nm were detected upon exposure to e aq Ϫ decay at different rates (Fig. 3A), which indicates that two processes take place simultaneously. Both appear within 5 s, i.e. "instantaneously." The absorbance at 440 nm decays by a firstorder reaction, k 440 nm ϭ 9.6 ϫ 10 4 s Ϫ1 (Fig. 3B), whereas that at 410 nm decays by two consecutive processes with rate constants k 410 nm-I ϭ 4.2 ϫ 10 4 s Ϫ1 and k 410 nm-II ϭ 2.5 ϫ 10 3 s Ϫ1 (Fig. 3C). All processes are complete within Ͻ1 ms. All rate data have been collected in Table 1.

Jumpstarting cytochrome P450
Importantly, no transient absorbance bands that could correspond to Cpd 1 are detected, and the spectrum recorded after irradiation indicates that all iron has returned to the Fe 3ϩ state (Fig. 3D), which leads us to conclude that the heme did not suffer radiation damage. We also determined that the rate of NADH oxidation, which has been shown to reflect the rate of substrate hydroxylation by CYP101 (20), was similar before and after irradiation; the observed rate constant of k obs_activity after PR ϭ (3.7 Ϯ 0.3) ϫ 10 2 min Ϫ1 demonstrates that the enzyme activity is not significantly impaired by radiolysis (Fig. S3).
Analysis of kinetics experiments carried out in D 2 O followed at 440 nm give a result for k 440 nm (D 2 O) ϭ 3.9 ϫ 10 4 s Ϫ1 (Fig.  S4a). The difference relative to the rate constant for the reaction in water (Fig. 3B) amounts to a kinetic isotope effect of ϳ2, In contrast, the kinetics at 410 nm are not affected (Fig. S4b).

Product analysis after pulse radiolysis
To confirm that the transient species observed occur along the pathway of camphor hydroxylation, two samples of ϳ1 ml each were collected for product analysis. These samples, containing 50 M of the oxy-enzyme, had each been exposed to a dose of 200 Gy. The substrate and products were extracted with chloroform and analyzed by GC-MS. The GC-MS elution profile shown in Fig. 4A shows, in addition to the peak corresponding to the camphor substrate at 4.9 min, two additional peaks with elution times of 5.9 and 6.9 min. The peak at 6.9 min corresponds to 5-hydroxycamphor (167 Da) with m/z ϭ 168 (after single protonation by electrospray ionization; Fig. 4B), as follows from a comparison with the mass spectrum of authentic 5-hydroxycamphor (Fig. 4C). The yields, quantified by GC with flame ionization detection (FID), were, from two separate experiments, 2.7 and 3.7 M. The elution profile for m/z ϭ 166 with the peak at 5.9 min shown in Fig. 4D, together with the mass spectrum ( Fig. 4E), indicates that a small amount (ϳ0.9 Ϯ 0.3 M, quantified by GC-FID) of 5-ketocamphor, a secondary oxidation product, was also formed. No oxidation products of camphor were observed in control experiments performed without CYP101 present or in the presence of lysozyme, a non-heme protein, in place of CYP101.

Simulation of pulse radiolysis reaction conditions and theoretical product yield
The theoretical yield of 5-hydroxycamphor was calculated with a chemical kinetics simulator (IBM Research) program by simulating the reactions shown in Table S1. To provide a value for k 3 for use in this simulation, we assume it to be equal to k 2 , which we calculate according to Equation 1, where k control ϭ 1.0 ϫ 10 6 s Ϫ1 is the rate at which e aq Ϫ disappears in the absence of protein ( Fig. 2A). We found the yield of reduc- Ϫ with the heme has a t1 ⁄2 shorter than 0.65 s, which is the lifetime in the absence of the protein. Therefore reaction 2 is completed within 3 s or 5 half-lives. The absorbance change between 0 and 3 s, 40 mAbs, was used to calculate the rate constant of reaction 2. Reduction of the heme after 3 s is caused by highly reducing species generated by the reaction of H ⅐ and e aq Ϫ with the amino acids of the protein. In the absence of O 2 , this process contributes to roughly half of the total yield of reduction.

Jumpstarting cytochrome P450
tion (reaction 2), i.e. of the heme relative to e aq Ϫ ([CYP101Fe 2ϩ ]/ [e aq Ϫ ]), to be 0.8%, from which a value of k 3 (ϭ k 2 ) of (3.8 Ϯ 0.7) ϫ 10 9 M Ϫ1 s Ϫ1 is derived. This finding is consistent with literature reports of a rate constant for the overall reaction of CYP 101Fe 3ϩ with e aq Ϫ of 3 ϫ 10 10 M Ϫ1 s Ϫ1 and that ϳ9% of the e aq Ϫ reaches the heme (16,17). The concentration of O 2 available to CYP101Fe 2ϩ was determined by a kinetics simulation (Fig. S2), which shows that the original concentration of 50 M O 2 decreases by 18 M in 1 s. From a simulation of substrate oxidation by the primary oxidant generated via reduction of the heme of 50 M CYP101[FeO 2 ] 2ϩ at 200 Gy, we estimate the product yield to be 3.1 Ϯ 0.6 M (Fig. S5). We found experimental yields of 2.7 and 3.7 M hydroxycamphor product, from which we conclude that the reaction is, within the error, quantitative. Because there is no indication that pulse radiolysis impairs enzyme function or its spectral properties, we assume that reaction of e aq Ϫ with the protein causes no damage to the protein, which allows us to ignore protein damaging reactions in the simulations.
In addition to 5-hydroxycamphor, formation of a further oxidation product, 5-ketocamphor, was observed. To test the hypothesis that the alcohol is converted to the ketone by reaction with H 2 O 2 in the presence of Fe 3ϩ , we incubated 5-hydroxycamphor with CYP101Fe 3ϩ in the presence of H 2 O 2 at concentrations of up to 1 mM, a concentration more than 1 order of magnitude higher than can be achieved by the radia-

Jumpstarting cytochrome P450
tion dose used, and found no evidence for the formation of 5-ketocamphor, even after overnight incubation. Hydrogen peroxide clearly does not interfere with our measurements. Formation of 5-ketocamphor was observed only with H 2 O 2 present at concentrations of Ն2 mM (Fig. S6a). To test the hypothesis that 5-hydroxycamphor remaining in the active site could be oxidized during a second reaction cycle, we irradiated a solution of 50 M CYP101Fe 3ϩ with 100 M of 5-hydroxycamphor as substrate in the presence of 50 M O 2 ; analysis of a sample by GC-MS shows formation of 0.7 M 5-ketocamphor (Fig. S6c). As a negative control, a solution of 75 M of CYP101Fe 2ϩ was irradiated in the presence of 100 M O 2 without substrate present, after which the sample was mixed with a solution of 100 M 5-hydroxycamphor before analysis by GC-MS. Surprisingly, formation of 0.7 M 5-ketocamphor was found as well in this negative control (Fig. S6b).

Discussion
Our main findings are that the pulse radiolysis technique can be used to jump-start the catalytic cycle of cytochrome P450, providing the electrons to reduce the enzyme, and that Cpd 0 so formed is converted to Cpd 1 very rapidly. The latter reacts either with bound substrate or with a protein moiety even faster, such that Cpd 1 does not accumulate. Additionally, we demonstrate that reduction of CYP101 by pulse radiolysis is productive, i.e. that 5-hydroxycamphor is formed.

Validation of the pulse radiolysis technique
We demonstrate that it is possible to use the pulse radiolysis technique, which allows electrons to be added to the heme one at a time and permits reactions that occur within a few microseconds to be followed spectroscopically, to probe the natural catalytic cycle of CYP101. Pulse radiolysis is an inherently quantitative technique in which the dose is used to calculate [e aq Ϫ ], and the yields of all reaction products can be estimated relative to [e aq Ϫ ]. Although reactions after formation of Cpd 0 are too rapid to be accessible by stopped-flow methods, the third reaction in the CYP101cycle-formation of CYP101[FeO 2 ] 2ϩ from CYP101Fe 2ϩ and O 2 , Fig. 1-is considerably slower. We could analyze this relatively slow reaction by both the stopped-flow and pulse radiolysis techniques to confirm that pulse radiolysis is suitable

Jumpstarting cytochrome P450
for analysis of the more rapid reactions in the cycle. The rate constant measured by pulse radiolysis, k 4 ϭ 2.1 ϫ 10 5 M Ϫ1 s Ϫ1 , is essentially identical to that measured by stopped flow, 2.4 ϫ 10 5 M Ϫ1 s Ϫ1 , and is in agreement with the literature (19).

Intermediates in the CYP101 reaction cycle identified by oneelectron reduction steps
The one-electron reduction of CYP101[FeO 2 ] 2ϩ initially gives rise to an absorbance spectrum that closely resembles that published for Cpd 0. There are minor differences at higher wavelengths that may originate from the 100-fold higher dose of radiation used by Denisov et al. (23). We assign the maximum at 440 nm (Fig. 3A) to Cpd 0. The decay kinetics at 440 nm are monophasic, with k 440 nm ϭ 9.6 ϫ 10 4 s Ϫ1 (Fig. 3B). The initially increased absorbance observed at 410-nm decays in a biphasic first-order fashion, but with different rate constants, k 410 nm-I ϭ 4.2 ϫ 10 4 s Ϫ1 and k 410 nm-II ϭ 2.5 ϫ 10 3 s Ϫ1 (Fig. 3C). Given that the absorbance bands at both 410 and 440 nm are present within s after the pulse, we conclude that the absorbance decays at these wavelengths reflect not serial, but parallel processes. The transient absorbance at 410 nm may indicate the presence of a tyrosyl radical (TyrO ⅐ ), which is unlikely to result from reaction with O 2 . , present in high concentrations directly after the pulse, because this radical does not react with amino acids (24). Instead, the absorbance at 410 nm is probably due to camphor-free CYP101, which has been shown to form Cpd ES upon reaction with a peracid (25,26 ] 2ϩ without camphor bound allows us to appreciate the time scale of the events after formation of Cpd 0. The formation and decay of the camphor radical, which has no absorbance maxima in the UVvisible range, cannot be observed, but if the absorbance at 410 nm is indeed caused by formation of TyrO ⅐ , this may be taken as a "reporter" for formation of Cpd 1. Thus, although we do not directly observe formation of Cpd 1, we conclude that any Cpd 1 formed immediately oxidizes camphor, or, in the absence of camphor, oxidizes Tyr and does not accumulate. In agreement with this interpretation, we observe that the decay of TyrO ⅐ "trails" that of Cpd 0; the first rate constant for the decay of TyrO ⅐ is about half that of Cpd 0 (Table 1).
Within the error, the yield of product formation approaches 100%, that is, for every molecule of Cpd 0 produced, one of 5-hydroxycamphor was formed. This finding agrees with the cryoreduction study (7), with one difference. By pulse radiolysis, we find quantitative conversion of CYP101 to Cpd 0, whereas the reported yield of Cpd 0 by cryoreduction was ϳ50% (7), which may reflect changes in conformation of CYP101[FeO 2 ] 2ϩ upon reduction at 4°C that would not be pos-sible at the much lower temperatures used in the cryoreduction studies.

Isotope effect
The observed kinetic deuterium isotope effect of 2 corroborates our assignment of the observed absorption band of 440 nm to Cpd 0, which becomes protonated after reduction according to the canonical cycle shown in Fig. 1. In contrast, the rate of decay of absorbance at 410 nm is essentially the same in D 2 O and H 2 O (Fig. 3C and S4b); the absence of an isotope effect fits with the hypothesis that any TyrO ⅐ formed would be reduced by other amino acids, with no isotope effect.

Cpd 1 does not accumulate
It is intriguing that, with camphor bound to CYP101, we find no indication for formation of Cpd 1, i.e. no absorbance increase near 370 or 690 nm, whereas Rittle and Green (12) found convincing support for production of such a species in the absence of substrate, after the reaction of CYP119 with chloroperbenzoate. Given the small absorptivity, 4 ϫ 10 3 M Ϫ1 cm Ϫ1 at 690 nm (29), it is unlikely that under our experimental conditions, we would observe any increase in absorbance at that wavelength, but an intermediate giving rise to a band at a wavelength lower than 400 nm should have been observable. Formation of Cpd 1 was also not observed in cryostudies (7). Without substrate bound, CYP101[FeO 2 ] 2ϩ is too unstable to be effectively reduced by pulse radiolysis.
Although we did not observe any accumulation of Cpd 1, we must infer from isotope studies at low temperature that Cpd 1 is involved in the mechanism (30): the hydrogen attached to camphor is retained in the hydroxylated form when the reaction is carried out with deuterated camphor (31). Furthermore, it would seem unlikely that Cpd 0, with electrode potential E°Ј ϭ ϩ0.93 V (3), could rapidly oxidize Tyr, which would require a redox partner with E°Ј ϭ ϩ0.97 V (32). Indeed, Cpd 1, with E°Ј ϭ ϩ1.35 V (3), is a far more likely candidate.
Our results suggest that hydroxylation is rather fast, as inferred from the observation that at 4°C, the conversion of Cpd 0 to CYP101Fe 3ϩ takes 37 s, calculated from the half-life of the decay at 440 nm, 7.4 s, based on k 440 nm ϭ 9.6 ϫ 10 4 s Ϫ1 , and multiplying the half-life by 5. Given this rate constant for the disappearance of Cpd 0 and under the assumption that Cpd 1 is formed, the Gibbs activation energy is about 45 kJ/mol (Fig.  5), as calculated via Equation 2, in which k is the rate constant, k B is Boltzmann's constant, h is Planck's constant, T is the temperature, ⌬G ‡ is the Gibbs energy of activation, and R is the gas constant. Because the next two steps, abstraction of a hydrogen from camphor and the rebound (33) step, are not observable, these reactions must have rates of Ն10 6 s Ϫ1 and thus activation energies of Յ39 kJ/mol. With radical clock experiments, it has been established that the rebound step is extremely fast, with rate constants on the order of 10 10 -10 11 s Ϫ1 (34), values that imply activation energies much smaller than 39 kJ/mol, indeed, close to 0 kJ/mol. Taking

Jumpstarting cytochrome P450
into consideration that the experiments were carried out at 4°C and using the value of 45 kJ/mol in Equation 2, we calculate that, at a physiological temperature of 37°C, conversion of Cpd 0 to Cpd 1 is faster by a factor of 11. This factor implies that, at physiological temperatures, this process requires only 3.5 s.

Further oxidation of 5-hydroxycamphor
In addition to the 5-hydroxycamphor, we observed formation of the further oxidized product 5-ketocamphor (Fig. 4A), production of which has also been reported from multiple turnovers of the enzyme (35,36). The ketone product could also be formed by reaction with H 2 O 2 , which is a direct product of pulse radiolysis that can oxidize alcohols in the presence of Fe 3ϩ (37). We concluded that the Յ20 M of H 2 O 2 generated during pulse radiolysis is unlikely to cause further oxidation of 5-hydroxycamphor, because even bulk addition of 1 mM H 2 O 2 to a mixture of CYP101-Fe ϩ3 and hydroxycamphor did not produce any 5-ketocamphor (Fig. S6a). The organic hydroperoxides (CH 3 ) 2 C(OH)CH 2 OOH formed in small amounts from the scavenger molecule t-BuOH directly after irradiation could possibly react with 5-hydroxycamphor to form 5-ketocamphor via a yet undocumented reaction mechanism. Conversion of the primary alcohol product to a ketone may also be catalyzed during multiple turnovers of CYP101 (36); 5-hydroxycamphor remaining in the active site could be oxidized during a second reaction cycle initiated by transfer of electrons to the heme from aromatic residues, followed by binding of O 2 , etc. We observed the formation of very small amounts of 5-ketocamphor when 5-hydroxycamphor was added as a substrate to CYP101Fe 2ϩ , believed to have been stripped of camphor, either before irradiation in the presence of O 2 or after irradiation as a negative control.

Energetics
The abstraction of H from camphor by Cpd 1 is, within the error of the estimates, thermoneutral (reaction S1), whereas the rebound step is very favorable, Ϫ241 kJ/mol, (Reaction S4). The energetics of the hydroxylation of camphor by CYP101 are shown in Fig. 5.
The value of E°(C ⅐ , H ϩ /CH) ϭ ϩ1.7 V (see supporting materials), or ϩ1.3 V at pH 7, in which CH represents camphor and C ⅐ represents the camphor radical, corresponds to a bond dissociation energy ⌬ BDE H°ϭ 397 kJ/mol, which agrees very well with that calculated by in silico methods (38). Oxidation of methane or other small alkanes would require dissociation of a stronger C-H bond. For CH 4 , ⌬ BDE H°ϭ 435 kJ/mol, which corresponds to an electrode potential of ϩ2.0 V at pH 0, and ϩ1.6 V at pH 7; slightly smaller numbers are obtained for other small alkanes. Compared with the E°Ј of ϩ1.3 V for the oneelectron oxidation of camphor, the value of ϩ1.6 V required for CH 4 implies that the activation barrier for hydrogen abstraction is higher by ϳ24 kJ/mol, making the reaction much slower, with an important consequence. A small alkane bound to the active site of cytochrome P450 would diffuse readily away, and hydroxylation is possible only when the substrate cannot escape, as has been achieved by the addition of perfluorocarboxylic acids that block the entrance of the active site of cytochrome P450 (39).

Tethering in silico to reality
In silico calculations yield energies larger than those calculated by means of thermodynamic cycles. For example, the reaction Cpd 0 ϩ H ϩ 3 Cpd 1 ϩ H 2 O, is "either almost thermoneutral or slightly exothermic, with a 'best' estimated reaction energy of Ϫ8 kcal/mol" (40). However, the value of Ϫ8 Figure 5. Energetics of the hydroxylation of camphor. Cpd 0 ([FeOOH] 2ϩ ), Cpd 1 (FeO 2ϩ , por ⅐ϩ ), and Cpd 2 (FeO 2ϩ ) with the camphor radical are isoenergetic, as explained in the supporting materials. The reaction is "pulled" by the very favorable final step, the formation of 5-hydroxycamphor and Fe 3ϩ CYP101. Given the rate constant of the decay of Cpd 0, 10 5 s Ϫ1 , and that Cpds 1 and 2 are not observed, we assume that the latter are formed and decay with rate constants of 10 6 s Ϫ1 or higher. These rate constants give rise to the Gibbs energies of the transition states of 45 and 39 kJ/mol; 39 kJ/mol is the upper limit for the transition states between Cpds 1 and 2 and Cpd 2 and Fe 3ϩ .

Jumpstarting cytochrome P450
kcal/mol cannot be correct, because either Cpd 1 would not be sufficiently oxidizing to abstract H, or H 2 O 2 would dissociate from Cpd 0, or both (3). Thus, "almost thermoneutral" must be correct. The activation enthalpy of the transition between Cpd 0 and Cpd 1 has been reported as ϳ15 kcal/mol or 63 kJ/mol (38); the Gibbs energy would be higher. For this reaction, we estimate a Gibbs activation energy of only 45 kJ/mol (Fig. 5). If, at a minimum, the barrier were 63 kJ/mol, t1 ⁄ 2 would be 9 ms, and Cpd 0 should have been observed in stopped-flow experiments but was not. Thus, as no formation of Cpd 0 was detected, we question that the enthalpy could be as high as 63 kJ/mol. The third example is the energy for the release of O 2 . from CYP101[FeO 2 ] 2ϩ , which has been reported as 19.5 kcal/mol (41). The value from a simple thermodynamic cycle, for which experimental values are used, is 34 kJ/mol (3) or 8.1 kcal/mol. In silico calculations applied to complicated reactions have been questioned (42). It is possible and essential to make in silico calculations more quantitative by "tethering" them to semiempirical values and considerations of relatively simple thermodynamic cycles.

Conclusions
We demonstrate that it is possible to mimic the natural catalytic cycle of CYP101 by using pulse radiolysis to rapidly supply electrons to the protein. This rapid electron delivery allowed us to visualize the steps following the reduction of substrate bound CYP101[FeO 2 ] 2ϩ on the microsecond timescale. We found that, as soon as Cpd 0 is formed, the rest of the cycle is completed within 40 s at 4°C. We found no spectral evidence that Cpd 1 accumulates. We demonstrate that product formation is quantitative once CYP101[FeO 2 ] 2ϩ is reduced by one electron. Further, reduction of the enzyme by this method does not damage the enzyme, as is demonstrated by the complete recovery of enzyme activity and of Fe 3ϩ after pulse radiolysis.

Protein preparation
CYP101 was expressed and purified as previously described, with slight modifications (41) (see supporting text and Fig. S7 for details). The ratio A 392 /A 280 was 1.62.

Pulse radiolysis
Pulse radiolysis produces ionizing radiation that leads to the formation of primary species from water with the following yields: G(HO ⅐ ) ϭ 2.65, G(e aq Ϫ ) ϭ 2.65, G(H ⅐ ) ϭ 0.54, G(H 2 O 2 ) ϭ 0.54, where G ϭ 1.0 refers to 0.1036 M species generated per 1 J/kg energy absorbed. Experiments were performed with a 2 MeV Febetron 705 accelerator (Titan Systems Corp., presently L-3 Communications, San Leandro, CA), with Ͻ50-ns pulses of 80 Gy and a custom-built 6-cm quartz irradiation cell (Hellma, Mülhausen, Germany) with a volume of 180 l. Dosimetry was carried out with KSCN (43). Absorbance maxima of 380 -450 nm were measured by means of a 300 mm Czerny-Turner type monochromator equipped with a grating of 150 grooves/mm with 300 nm blazing and a photomultiplier for detection; for wavelengths Ն 450 nm, a grating of 150 grooves/mm with 500 nm blazing was used. The monochromator slit width was set to 0.4 mm.

Sample preparation for reaction of hydrated electrons with [FeO 2 ] 2؉ CYP101
All experiments were carried out at 4°C. CYP101Fe 3ϩ in the presence of camphor was reduced with a slight excess of S 2 O 4 2Ϫ inside a N 2 -filled glove box and loaded into a 10-ml gas-tight syringe (Hamilton Sample Lock). A second gas-tight syringe contained a solution of 100 M O 2 in 10 mM pH 7.4 potassium phosphate buffer, 10 M camphor with 2 mM t-BuOH. t-BuOH was used to scavenge HO ⅐ . The two air-tight syringes were placed in a syringe pump prior to pulse irradiation. To maintain the temperature at 4°C before mixing, the capillaries from the syringes to the mixer were kept on ice. Irradiation took place within 30 s after mixing, such that CYP101[FeO 2 ] 2ϩ did not decay significantly before irradiation. The buffer used in pulse radiolysis, UV-visible, GC-MS, and stopped-flow experiments was 10 mM potassium phosphate buffer, pH 7.4, containing 2 equivalents of camphor relative to CYP101, or as otherwise specified.

Product analysis
A solution of 100 M reduced CYP101 containing 100 M of camphor was mixed with air-saturated buffer prior to irradiation. The applied dose was 200 Gy in a single pulse. A total of ϳ1 ml of irradiated solution was collected. Substrate and products were extracted with chloroform (1:1 v/v) containing 1 mM decane as an external standard; the organic layer was concentrated by evaporation under a stream of N 2 and analyzed by GC-MS or GC-FID. Peaks were assigned by GC-MS and quantified by peak integration of the GC-FID traces and comparison with standard curves of the substrate and product. The product 5-hydroxycamphor was synthesized enzymatically as described previously (41).
Author contributions-H. E. performed the majority of the experiments. T. N. guided and assisted with the pulse radiolysis experiments. H. E., A. V., and P. L. B. were involved in the preparation and purification of CYP 101. H. E. and A. V. performed the product analysis and quantification studies. W. H. K. directed the experiments and introduced the energetics in the discussion. All authors contributed to the writing.