Partially folded state of the disulfide-reduced N-terminal half-molecule of ovotransferrin as a renaturation intermediate

A previous report (Hirose, M., Akuta, T., and Takahashi, N. (1989) J. Biol. Chem. 264, 16867-16872) has shown that for the efficient oxidative refolding of disulfide-reduced ovotransferrin, a preincubation under reduced conditions at a low temperature is essential. To study the renaturation pathway, the disulfide-reduced N-terminal half-molecule of ovotransferrin was analyzed by CD spectrum. The reduced protein was found to take, at low temperatures, a partially folded conformation that can be distinguished from both the native and denatured states. The folded protein was in a metastable state with delta GD value of 2.2-2.8 kcal/mol at 6 degrees C. The conformation was variable depending on temperature conditions; its stability was decreased at a lower temperature (1.0-1.2 kcal/mol at 0 degrees C). Subsequent reoxidation at 6 degrees C by oxidized glutathione led efficiently the reduced protein to the correctly renatured form having the iron-binding capacity, indicating that the partially folded state is the immediate precursor to subsequent oxidative refolding.     

Escherichia coli thioredoxin folds into two compact forms of different stability to urea denaturation

The urea-induced denaturation of Escherichia coli thioredoxin and thioredoxin variants has been examined by electrophoresis on urea gradient slab gels by the method of Creighton [Creighton, T. (1986) Methods Enzymol. 131, 156-172]. Thioredoxin has only two cysteine residues, and these form a redox-active disulfide at the active site. Oxidized thioredoxin-S2 and reduced thioredoxin-(SH)2 each show two folded isomers with a large difference in stability to urea denaturation. The difference in stability is greater for the isomers of oxidized than for the isomers of reduced thioredoxin. At 2 degrees C, the urea concentrations at the denaturation midpoint are approximately 8 and 4.3 M for the oxidized isomers and 4.8 and 3.7 M for the reduced isomers. The difference between the gel patterns of samples applied in native versus denaturing buffer, and at 2 and 25 degrees C, is characteristic for the involvement of a cis-proline-trans-proline isomerization. The data very strongly suggest that the two folded forms of different stabilities correspond to the cis and trans isomers of the highly conserved Pro 76 peptide bond, which is cis in the crystal structure of oxidized thioredoxin. Urea gel experiments with the mutant thioredoxin P76A, with alanine substituted for proline at position 76, corroborate this interpretation. The electrophoretic banding pattern diagnostic for an involvement of proline isomerization in urea denaturation is not observed for oxidized P76A. In broad estimates of delta G degree for the native-denatured transition, the difference in delta G degree (no urea) between the putative cis and trans isomers of the Ile 75-Pro 76 peptide bond is approximately 3 kcal/mol for oxidized thioredoxin and approximately 1.5 kcal/mol for reduced thioredoxin. Since cis oxidized thioredoxin is much more stable than trans, folded oxidized thioredoxin is essentially all cis. In folded reduced thioredoxin, cis and trans interconvert slowly, on the minute time scale at 2 and 25 degrees C. In the absence of urea, the folded reduced thioredoxin is less than a few percent trans. Three additional mutants with additions or substitutions at the active site also show electrophoresis banding patterns consistent with a difference in stability between cis and trans isomers.     

QM/MM study of the active site of free papain and of the NMA-papain complex.

Hybrid quantum mechanical/molecular mechanical (QM/MM) calculations using restricted and unrestricted Hartree-Fock and B3LYP ab initio (QM) and Amber force field (MM), respectively, have been applied to study the catalytic site of papain in both free and substrate bonded forms. Ab initio geometry optimizations have been performed for the active site of papain and the N-methyl-acetamide (NMA)-papain complex within the molecular mechanical treatment of the protein environment. A covalent tetrahedral intermediate structure could be obtained only when the amide N atom of the substrate molecule was protonated through a proton transfer from the His-159 in the catalytic site. Our results support the previous assumption that a proton transfer from His-159 to the amide N atom of the substrate occurs prior to or concerted with the nucleophilic attack of the Cys-25 sulfur atom to the carbonyl group of the substrate. The electron correlation effect will reduce the proton transfer barrier. Therefore, this proton transfer can be easily observed in the B3LYP/6-31G* calculations. The HF/6-31G* method overestimates the reaction barrier against this proton transfer. The sulfur atom of Cys-25 and the imidazole ring of His-159 are found to be coplanar in the free form of the enzyme. However, the rotation of the imidazole ring of His-159 was observed during the formation of the tetrahedral intermediate. Without the papain environment, the coplanar thiolate-imidazolium ion pair RS-...ImH+ is much less stable than the neutral form of RSH....Im. Within the protein environment, however, the thiolate-imidazolium ion pair becomes more stable than its neutral form by 4.1 and 0.4 kcal/mol in HF/6-31G* and B3LYP/6-31G* calculations, respectively. The barrier of proton transfer from S-H group of Cys-25 to the imidazole ring of His-159 was reduced from 22.0 kcal/mol to 15.2 kcal/mol by the protein environment in HF/6-31G* calculations. This barrier is found to be much smaller (2.5 kcal/mol) in B3LYP/6-31G* calculations.     

Proteolysis as a measure of the free energy difference between cytochrome c and its derivatives.

Limited cleavage of oxidized and reduced horse heart cytochrome c (Cyt c) and the azide complex of Cyt c by proteinase K at room temperature yields a single cut within the central loop (36-60 in the sequence). Using an assay that allows spectroscopic evaluation of the fraction of intact protein as a function of time, together with a simple kinetic model for proteolysis, fluctuation opening of the loop can be related to the free energy of the corresponding protein. This allows us to estimate quantitatively the free energy difference between the oxidized form of Cyt c and other states using proteolysis as a probe. The results we obtain indicate that oxidized Cyt c is 2.0 kcal mol(-1) less stable than the reduced form, and 0.07 kcal mol(-1) is more stable than the Cyt c: azide complex at 25 degrees C. These values agree in magnitude with results from hydrogen exchange and unfolding studies, suggesting that the stability of a protein can be directly related to its structural dynamics.     

Replacement of a conserved proline eliminates the absorbance-detected slow folding phase of iso-2-cytochrome c

As a test of the proline isomerization model, we have used oligonucleotide site-directed mutagenesis to construct a mutant form of iso-2-cytochrome c in which proline-76 is replaced by glycine [Wood, L. C., Muthukrishnan, K., White, T. B., Ramdas, L., & Nall, B. T. (1988) Biochemistry (preceding paper in this issue)]. For the oxidized form of Gly-76 iso-2, an estimate of stability by guanidine hydrochloride induced unfolding indicates that the mutation destabilizes the protein by 1.2 kcal/mol under standard conditions of neutral pH and 20 degrees C (delta G degrees u = 3.8 kcal/mol for normal Pro-76 iso-2 versus 2.6 kcal/mol for Gly-76 iso-2). The kinetics of folding/unfolding have been monitored by fluorescence changes throughout the transition region using stopped-flow mixing. The rates for fast and slow fluorescence-detected refolding are unchanged, while fast unfolding is increased in rate 3-fold in the mutant protein compared to normal iso-2. A new kinetic phase in the 1-s time range is observed in fluorescence-detected unfolding of the mutant protein. The presence of the new phase is correlated with the presence of species with an altered folded conformation in the initial conditions, suggesting assignment of the phase to unfolding of this species. The fluorescence-detected and absorbance-detected slow folding phases have been monitored as a function of final pH by manual mixing between pH 5.5 and 8 (0.3 M guanidine hydrochloride, 20 degrees C).(ABSTRACT TRUNCATED AT 250 WORDS)     

Analysis of the ionization constants and heats of ionization of reduced and oxidized horse heart cytochrome c

Simultaneous curve fitting for the ionization parameters of oxidized and reduced horse heart cytochrome c in 0.15M KCl and 20°C yields values for the ionization constants (as pK′) and the heats of ionization (ΔH_i) which can reconstruct either the potentiometric or thermal titration curves. Reduced cytochrome c requires 8 sets of groups, whereas oxidized cytochrome c requires 10 sets of groups. The additional groups in the oxidized preparation appear to involve the ferriheme (pK′, 9.25; ΔH_i, 13.7 kcal/mol) and a tyrosine (pK′ ? 10.24) that is not present in the reduced form. The potentiometric and thermal difference curves (reduced – oxidized) involve the appearance of 17 kcal/mol centered at pH 9.7 and 5.8 kcal/mol centered at pH 4.9. The carboxyl groups in both species appear to be normal for the hydrogen-bonded form. Only one histidine has normal ionization properties (pK′, 6.7; ΔH_i, 7.5 kcal/mol), as do 17 of the lysine residues (pK′, 10.8; ΔH_i, 11.5 kcal/mol).     

Conformational stability and activity of ribonuclease T1 with zero, one, and two intact disulfide bonds

Ribonuclease T1 has two disulfide bonds linking cysteine residues 2-10 and 6-103. We have prepared a derivative of ribonuclease T1 in which one disulfide bond is broken and the cysteine residues carboxymethylated, (2-10)-RCM-T1, and three derivatives in which both disulfides are broken and the cysteine residues reduced, R-T1, carboxamidomethylated, RCAM-T1, or carboxymethylated, RCM-T1. The RNA hydrolyzing activity of these proteins has been measured, and urea and thermal denaturation studies have been used to determine conformational stability. The activity, melting temperature, and conformational stability of the proteins are: ribonuclease T1 (100%, 59.3 degrees C, 10.2 kcal/mol), (2-10)-RCM-T1 (86%, 53.3 degrees C, 6.8 kcal/mol), R-T1 (53%, 27.2 degrees C, 3.0 kcal/mol), RCAM-T1 (43%, 21.2 degrees C, 1.5 kcal/mol), and RCM-T1 (35%, 16.6 degrees C, 0.9 kcal/mol). Thus, the conformational stability is decreased by 3.4 kcal/mol when one disulfide bond is broken and by 7.2-9.3 kcal/mol when both disulfide bonds are broken. It is quite remarkable that RNase T1 can fold and function with both disulfide bonds broken and the cysteine residues carboxymethylated. The large decrease in the stability is due mainly to an increase in the conformational entropy of the unfolded protein which results when the constraints of the disulfide bonds on the flexibility are removed. We propose a new equation for predicting the effect of a cross-link on the conformational entropy of a protein: delta Sconf = -2.1 - (3/2)R 1n n, where n is the number of residues between the side chains which are cross-linked. This equation gives much better agreement with experimental results than other forms of this equation which have been used previously.     

Alteration of thermal stability of glucose oxidase associated with the redox states

By the method of differential scanning calorimetry, it was found that thermal stability of glucose oxidase was dependent on its redox states. The oxidized form showed an apparent denaturation temperature at 76°C and the denaturation enthalpy was approximately 865 kcal/mol. On reduction of the enzyme, the denaturation temperature increased by about 10°, but no significant change was seen in the denaturation enthalpy. The activation energies of the denaturation of the oxidized and the reduced enzymes were about 89 and 103 kcal/mol, respectively. These results may imply conformational changes in the catalytic turnover of this enzyme.     

Structures of nitroreductase in three states: effects of inhibitor binding and reduction.

The crystal structure of the nitroreductase enzyme from Enterobacter cloacae has been determined for the oxidized form in separate complexes with benzoate and acetate inhibitors and for the two-electron reduced form. Nitroreductase is a member of a group of enzymes that reduce a broad range of nitroaromatic compounds and has potential uses in chemotherapy and bioremediation. The monomers of the nitroreductase dimer adopt an alpha+beta fold and together bind two flavin mononucleotide prosthetic groups at the dimer interface. In the oxidized enzyme, the flavin ring system adopts a strongly bent (16 degrees ) conformation, and the bend increases (25 degrees ) in the reduced form of the enzyme, roughly the conformation predicted for reduced flavin free in solution. Because free oxidized flavin is planar, the induced bend in the oxidized enzyme may favor reduction, and it may also account for the characteristic inability of the enzyme to stabilize the one electron-reduced semiquinone flavin, which is also planar. Both inhibitors bind over the pyrimidine and central rings of the flavin in partially overlapping sites. Comparison of the two inhibitor complexes shows that a portion of helix H6 can flex to accommodate the differently sized inhibitors suggesting a mechanism for accommodating varied substrates.     

An ab initio study of valence isomerization in the HOCl-HClO system

Hartree-Fock level and post Hartree-Fock level molecular orbital calculations have been completed for HOCl and its valence isomer, HClO. Ground state geometries have been determined for each molecule. The energy change for the reaction HOCl→HClO is estimated to be 67±5 kcal/mol endothermic and the activation energy for the process is 74±5 kcal/ mol. The transition state for the reaction is identified and discussed. A vibrational analysis has been carried out for both HOCl and HClO. The calculated photoelectron spectrum is discussed for each of the molecules.     

Molecular dynamics simulations of rubredoxin from Clostridium pasteurianum: Changes in structure and electrostatic potential duringredox reactions

Molecular dynamics simulations of Clostridium pasteurianum rubredoxin in the oxidized and reduced forms have been performed. Good agreement between both forms and crystal data has been obtained (rms deviation of backbone atoms of 1.06 and 1.42 Å, respectively), which was due in part to the use of explicit solvent and counterions. The reduced form exhibits an unexpected structural change: the redox site becomes much more solvent-accessible, so that water enters a channel between the surface and the site, but with little actual structural rearrangement (the rms deviation of backbone atoms between the oxidized and reduced is 0.77 Å). The increase in solvent accessibility is also seen, although to a much lesser extent, between the oxidized and reduced crystal structures of Pyrococcus furiosus rubredoxin, but no high resolution crystal or nuclear magnetic resonance solution data exist for reduced C. pasteurianum rubredoxin. The electrostatic potential at the iron site and fluctuations in the potential, which contribute to both the redox and electron transfer properties, have also been evaluated for both the oxidized and the reduced simulations. These results show that the backbone plays a significant role (62–70 kcall/mol/e) and the polar sidechains contribute relatively little (0–4 kcal/mol/e) to the absolute electrostatic potential at the iron of rubredoxin for both forms. However, both groups contribute significantly to the change in redox state by becoming more polarized and more densely packed around the redox site upon reduction. Furthermore, these results show that the solvent becomes much more polarized in the reduced form than in the oxidized form, even excluding the penetrating water. Finally, the simulation indicates that the contribution of the charged side chains to the electrostatic potential is largely canceled by that of the counterions. © 1995 Wiley-Liss, Inc.     

Conformational and thermodynamic characterization of the molten globule state occurring during unfolding of cytochromes-c by weak salt denaturants

The denaturation of bovine and horse cytochromes-c by weak salt denaturants (LiCl and CaCl(2)) was measured at 25 degrees C by observing changes in molar absorbance at 400 nm (Delta epsilon(400)) and circular dichroism (CD) at 222 and 409 nm. Measurements of Delta epsilon(400) and mean residue ellipticity at 409 nm ([theta](409)) gave a biphasic transition for both modes of denaturation of cytochromes-c. It has been observed that the first denaturation phase, N (native) conformation <--> X (intermediate) conformation and the second denaturation phase, X conformation <--> D (denatured) conformation are reversible. Conformational characterization of the X state by the far-UV CD, 8-anilino-1-naphthalene sulfonic acid (ANS) binding, and intrinsic viscosity measurements led us to conclude that the X state is a molten globule state. Analysis of denaturation transition curves for the stability of different states in terms of Gibbs energy change at pH 6.0 and 25 degrees C led us to conclude that the N state is more stable than the X state by 9.55 +/- 0.32 kcal mol(-1), whereas the X state is more stable than the D state by only 1.40 +/- 0.25 kcal mol(-1). We have also studied the effect of temperature on the equilibria, N conformation <--> X conformation and X conformation <--> D conformation in the presence of different denaturant concentrations using two different optical probes, namely, [theta](222) and Delta epsilon(400). These measurements yielded T(m), (midpoint of denaturation) and Delta H(m) (enthalpy change) at T(m) as a function of denaturant concentration. A plot of Delta H(m) versus corresponding T(m) was used to determine the constant-pressure heat capacity change, Delta C(p) (= ( partial differential Delta H(m)/ partial differential T(m))(p)). Values of Delta C(p) for N conformation <--> X conformation and X conformation <--> D conformation is 0.92 +/- 0.02 kcal mol(-1) K(-1) and 0.41 +/- 0.01 kcal mol(-1) K(-1), respectively. These measurements suggested that about 30% of the hydrophobic groups in the molten globule state are not accessible to the water.     

A test of the model to predict unusually stable RNA hairpin loop stability

To investigate the accuracy of a model [Giese et al., 1998, Biochemistry37:1094-1100 and Mathews et al., 1999, JMol Biol 288:911-940] that predicts the stability of RNA hairpin loops, optical melting studies were conducted on sets of hairpins previously determined to have unusually stable thermodynamic parameters. Included were the tetraloops GNRA and UNCG (where N is any nucleotide and R is a purine), hexaloops with UU first mismatches, and the hairpin loop of the iron responsive element, CAGUGC. The experimental values for the GNRA loops are in excellent agreement (deltaG degrees 37 within 0.2 kcal/mol and melting temperature (TM) within 4 degrees C) with the values predicted by the model. When the UNCG hairpin loops are treated as tetraloops, and a bonus of 0.8 kcal/mol included in the prediction to account for the extra stable first mismatch (UG), the measured and predicted values are also in good agreement (deltaG degrees 37 within 0.7 kcal/mol and TM within 3 degrees C). Six hairpins with unusually stable UU first mismatches also gave good agreement with the predictions (deltaG degrees 37 within 0.5 kcal/mol and TM within 8 degrees C), except for hairpins closed by wobble base pairs. For these hairpins, exclusion of the additional stabilization term for UU first mismatches improved the prediction (AG degrees 37 within 0.1 kcal/mol and TM within 3 degrees C). Hairpins with the iron-responsive element loop were not predicted well by the model, as measured deltaG degrees 37 values were at least 1 kcal/mol greater than predicted.     

Effect of polyanions on the refolding of human acidic fibroblast growth factor.

Acidic fibroblast growth factor (aFGF) is unstable at physiological temperatures in the absence of polyanions such as heparin. Therefore, the effect of temperature on the kinetics of refolding of aFGF has been examined in the presence and absence of several polyanions. The protein folds into its native state at temperatures up to 30 degrees C without polyanions with an activation energy of approximately 14 kcal/mol, but does not acquire native structure above this temperature. When heparin, inositol hexasulfate, or sulfate ion are present, aFGF refolds below 30 degrees C with a slightly reduced activation energy (10-11 kcal/mol). In addition, the protein now also renatures between 30 and 50 degrees C with activation energies of 1-2 (heparin), 16 (inositol hexasulfate), and 7 (sulfate) kcal/mol. Trace heavy metals appear to inhibit the refolding process, but a molecular chaperone (bovine 70-kDa heat shock cognate protein) and a peptidylprolyl isomerase (the FK506-binding protein) have no effect. It is concluded that the rate of refolding of aFGF at physiological temperatures is probably controlled by the interaction of a native-like state of the protein with an unknown polyanionic species.     

Intramolecular electron transfer in trimethylamine dehydrogenase from bacterium W3A1.

Reductive optical/EPR titrations of trimethylamine dehydrogenase with sodium dithionite have been performed, indicating that the equilibrium distribution of reducing equivalents between the covalently bound FMN and 4Fe/4S centers in partially reduced trimethylamine dehydrogenase is pH-dependent. In the case of two-electron reduced enzyme, formation of fully reduced flavin with oxidized iron-sulfur is favored below pH 7.5, whereas above pH 8 formation of flavin semiquinone with reduced iron-sulfur is preferred. The rates of electron transfer between the sites have been measured with the stopped-flow rapid mixing technique using a pH jump. The observed rate constants fall in the range of 200 s-1 to 1000 s-1 at 25 degrees C with the larger values occurring at higher values of final pH. The values of the rate constants depend on the final pH and are independent of observation wave-length. The temperature dependencies of these reactions give linear Arrhenius plots with activation energies in the range of 12 to 16 kcal/mol, consistent with prototropic equilibria being associated with electron transfer. The pH dependence of EPR spectral line widths for the flavin semiquinone and static optical spectra suggest that the semiquinone form of flavin present at pH 10 is anionic, whereas the neutral form is present at pH 7. The observed rate constants at 25 degrees C are greater than or equal to 100-fold larger than kcat for this enzyme and indicate that intramolecular electron transfer is not intrinsically rate-limiting in overall catalysis.     

Equilibrium and kinetic measurements of the conformational transition of reduced thioredoxin

The single disulfide bond in Escherichia coli thioredoxin was reduced by reaction with a 20-fold excess of reduced dithiothreitol at neutral pH and 25 degrees C. For some measurements, reduced thioredoxin was further reacted with iodoacetamide to alkylate the cysteinyl residues. The denaturation transitions of oxidized, reduced, and reduced alkylated thioredoxin were observed by using far-ultraviolet circular dichroic and exclusion chromatographic measurements. Cleavage of the disulfide bond lowers the stability of the native thioredoxin to denaturation by about 2.4 kcal/mol, and subsequent alkylation lowers the stability by a further 1.6 kcal/mol. The kinetics of the conformational change of reduced thioredoxin in guanidine hydrochloride were observed by using exclusion chromatography at moderate pressure and 2 degrees C. Analyses of single and multimixing protocols are consistent with a predominant nonnative configuration in the denatured state and the transient accumulation of a compact nativelike intermediate during refolding. The intermediate can incorporate the nonnative configuration and can accommodate its isomerization. No compelling chromatographic evidence was found for a conformation having an elution time different from that characteristic for either the native or the denatured protein.     

An electron spin resonance investigation and molecular orbital calculation of the anion radical intermediate in the enzymatic cis-trans isomerization of furylfuramide, a nitrofuran derivative of ethylene

The enzymatic cis-trans isomerization of nitrofuran derivatives has been proposed to occur via the formation of a radical anion intermediate. ESR investigations, in conjunction with intermediate neglect of differential overlap (INDO) molecular orbital calculations, support this concept by demonstrating the enzymatic generation of cis and trans radical anions of 3-(5-nitro-2-furyl)-2-(2-furyl) acrylamide. The INDO calculations further indicate that the rotational barrier between the cis and trans anion radicals of this compound is only 5--10 kcal/mol, whereas a 70 kcal/mol barrier exists for the parent geometric isomers. Hyperfine splitting constants for the cis-trans conformers have been assigned on the basis of INDO calculations. Surprisingly, only the nitrogen hyperfine splitting of the nitro group is distinguishably different in the two conformers, a result which is not inconsistent with the INDO calculations.     

The severity of osteogenesis imperfecta: a comparison to the relative free energy differences of collagen model peptides

Molecular dynamics simulations were carried out to calculate free energy differences between the folded and unfolded states of wild type and mutant collagen model peptides. The calculated stability of the collagen models was compared with the severity of osteogenesis imperfecta. Free energy differences of Gly → Xaa (Xaa: Ser, Cys, Glu, and Asp) mutations between the wild type and the mutants at position 15 of the model peptide were 3.8, 4.2, 5.6, and 8.8 kcal/mol, respectively. The corresponding free energy differences of a second Gly mutation at the same position in different chains were, on average, 1.3, 1.5, 2.9, and 5.4 kcal/mol, respectively. Free energy simulations were also performed to estimate the relative stability between an oxidized form and a reduced form of the mutants containing two Cys residues, which indicated that the mutant of the collagen-like peptide containing an intramolecular disulfide bond was more stable than the mutant containing one Cys residue but less stable than the wild type. The calculated free energy differences between an oxidized and a reduced form of the mutants containing two Cys residues are 0.8 and 2.6 kcal/mol for the disulfide bonds between Chains A and B and between Chains A and C, respectively.     

Effects of temperature on the mitochondrial malate dehydrogenase of adult muscle of Toxocara canis

Purified mitochondrial malate dehydrogenase isoenzyme (m-MDH) of Toxocara canis muscle presented maximum activity at 48 degrees C. A clear change in slope of the Arrhenius plot was observed. The energy of activation calculated for the catalytic process showed values of 3.2 kcal/mol and 10.5 kcal/mol. Thermal inactivation of m-MDH showed that it is more thermolabile than the s-isoenzyme. The inactivation of the enzyme by heat could be reduced at least in part by the addition of 0.1 mM NADH. The heat denaturation showed to be a first-order process. The rate constant (k) was calculated as being of the order of 5.28 X 10(-4) s-1 at 40 degrees C. The activation energy for the heat inactivation process was 16.45 kcal/mol between 30 degrees C and 40 degrees C and 13.79 kcal/mol between 40 degrees C and 48 degrees C.     

A theoretical study on nitric oxide reductase activity in a ba(3)-type heme-copper oxidase

The mechanism of nitric oxide reduction in a ba(3)-type heme-copper oxidase has been investigated using density functional theory (B3LYP). Four possible mechanisms have been studied and free energy surfaces for the whole catalytic cycle including proton and electron transfers have been constructed by comparison to experimental data. The first nitric oxide coordinates to heme a(3) and is partly reduced having some nitroxyl anion character ((3)NO(-)), and it is thus activated toward the attack by the second N-O. In this reaction step a cyclic hyponitrous acid anhydride intermediate with the two oxygens coordinating to Cu(B) is formed. The cyclic hyponitrous acid anhydride is quite stable in a local minimum with high barriers for both the backward and forward reactions and should thus be observable experimentally. To break the N-O bond and form nitrous oxide, the hyponitrous acid anhydride must be protonated, the latter appearing to be an endergonic process. The endergonicity of the proton transfer makes the barrier of breaking the N-O bond directly after the protonation too high. It is suggested that an electron should enter the catalytic cycle at this stage in order to break the N-O bond and form N(2)O at a feasible rate. The cleavage of the N-O bond is the rate limiting step in the reaction mechanism and it has a barrier of 17.3 kcal/mol, close to the experimental value of 19.5 kcal/mol. The overall exergonicity is fitted to experimental data and is 45.6 kcal/mol.