Quantum chemical calculations of the NHA bound nitric oxide synthase active site: O2 binding and implications for the catalytic mechanism

Density functional theory has been employed to model the binding of the intermediate substrate NHA, by nitric oxide synthases. In particular, the orientation and interactions of possibly catalytically important substrate hydrogens, with and without molecular oxygen bound to the active site heme group, are considered. Without O(2), three possible conformers have been found, with the energetically most favored structure being that in which both protons of the -NHOH moiety of NHA are directed toward the heme group. With oxygen bound, four different structures were found. The energetically lowest structure is again found to have both hydrogens of the -NHOH group pointing toward the heme group, thus forming hydrogen bonds between -NH- and the terminal oxygen, and between -OH and the inner oxygen of the heme-O(2) group. In addition, unprotonated structures of the substrate bound to the active site are considered and the proton affinity calculated.     

A density functional theory investigation on the mechanism of the second half-reaction of nitric oxide synthase

Density functional theory methods have been employed to systematically investigate the overall mechanism of the second half-reaction of nitric oxide synthases. The initial heme-bound hydrogen peroxide intermediate previously identified is found to first undergo a simple rotation about its O-O peroxide bond. Then, via a "ping-pong" peroxidase-like mechanism the -O(in)H- proton is transferred back onto the substrate's -NO oxygen then subsequently onto the outer oxygen of the resulting Fe(heme)-OOH species. As a result, O(out) is released as H2O with concomitant formation of a compound I-type (Fe(heme)-O) species. Formation of the final citrulline and NO products can then be achieved in one step via a tetrahedral transition structure resulting from direct attack of the Fe(heme)-O moiety at the substrate's guanidinium carbon center. The possible role of alternative mechanisms involving a protonated compound II-type species or an initial transfer of only the -NH- hydrogen of the =NHOH+ group of N(omega)-hydroxy-L-arginine is also discussed.     

Proton-coupled electron transfer reactions at a heme-propionate in an iron-protoporphyrin-IX model compound

A heme model system has been developed in which the heme-propionate is the only proton donating/accepting site, using protoporphyrin IX-monomethyl esters (PPIX(MME)) and N-methylimidazole (MeIm). Proton-coupled electron transfer (PCET) reactions of these model compounds have been examined in acetonitrile solvent. (PPIX(MME))Fe(III)(MeIm)(2)-propionate (Fe(III)~CO(2)) is readily reduced by the ascorbate derivative 5,6-isopropylidine ascorbate to give (PPIX(MME))Fe(II)(MeIm)(2)-propionic acid (Fe(II)~CO(2)H). An excess of the hydroxylamine TEMPOH or of hydroquinone similarly reduces Fe(III)~CO(2), and TEMPO and benzoquinone oxidize Fe(II)~CO(2)H to return to Fe(III)~CO(2). The measured equilibrium constants, and the determined pK(a) and E(1/2) values, indicate that Fe(II)~CO(2)H has an effective bond dissociation free energy (BDFE) of 67.8 ± 0.6 kcal mol(-1). In these PPIX models, electron transfer occurs at the iron center and proton transfer occurs at the remote heme propionate. According to thermochemical and other arguments, the TEMPOH reaction occurs by concerted proton-electron transfer (CPET), and a similar pathway is indicated for the ascorbate derivative. Based on these results, heme propionates should be considered as potential key components of PCET/CPET active sites in heme proteins.     

A DFT study of the mechanism of Ni superoxide dismutase (NiSOD): role of the active site cysteine-6 residue in the oxidative half-reaction

In the present DFT study, the catalytic mechanism of H2O2 formation in the oxidative half-reaction of NiSOD, E-Ni(II) + O2- + 2H+ --> E-Ni(III) + H2O2, has been investigated. The main objective of this study is to investigate the source of two protons required in this half-reaction. The proposed mechanism consists of two steps: superoxide coordination and H2O2 formation. The effect of protonation of Cys6 and the proton donating roles of side chains (S) and backbones (B) of His1, Asp3, Cys6, and Tyr9 residues in these two steps have been studied in detail. For protonated Cys6, superoxide binding generates a Ni(III)-O2H species in a process that is exothermic by 17.4 kcal/mol (in protein environment using the continuum model). From the Ni(III)-O2H species, H2O2 formation occurs through a proton donation by His1 via Tyr9, which relative to the resting position of the enzyme is exothermic by 4.9 kcal/mol. In this pathway, a proton donating role of His1 residue is proposed. However, for unprotonated Cys6, a Ni(II)-O2- species is generated in a process that is exothermic by 11.3 kcal/mol. From the Ni(II)-O2- species, the only feasible pathway for H2O2 formation is through donation of protons by the Tyr9(S)-Asp3(S) pair. The results discussed in this study elucidate the role of the active site residues in the catalytic cycle and provide intricate details of the complex functioning of this enzyme.     

Fundamental reaction pathway and free energy profile for inhibition of proteasome by Epoxomicin

First-principles quantum mechanical/molecular mechanical free energy calculations have been performed to provide the first detailed computational study on the possible mechanisms for reaction of proteasome with a representative peptide inhibitor, Epoxomicin (EPX). The calculated results reveal that the most favorable reaction pathway consists of five steps. The first is a proton transfer process, activating Thr1-O(γ) directly by Thr1-N(z) to form a zwitterionic intermediate. The next step is nucleophilic attack on the carbonyl carbon of EPX by the negatively charged Thr1-O(γ) atom, followed by a proton transfer from Thr1-N(z) to the carbonyl oxygen of EPX (third step). Then, Thr1-N(z) attacks on the carbon of the epoxide group of EPX, accompanied by the epoxide ring-opening (S(N)2 nucleophilic substitution) such that a zwitterionic morpholino ring is formed between residue Thr1 and EPX. Finally, the product of morpholino ring is generated via another proton transfer. Noteworthy, Thr1-O(γ) can be activated directly by Thr1-N(z) to form the zwitterionic intermediate (with a free energy barrier of only 9.9 kcal/mol), and water cannot assist the rate-determining step, which is remarkably different from the previous perception that a water molecule should mediate the activation process. The fourth reaction step has the highest free energy barrier (23.6 kcal/mol) which is reasonably close to the activation free energy (～21-22 kcal/mol) derived from experimental kinetic data. The obtained novel mechanistic insights should be valuable for not only future rational design of more efficient proteasome inhibitors but also understanding the general reaction mechanism of proteasome with a peptide or protein.     

QM/MM study of mechanisms for compound I formation in the catalytic cycle of cytochrome P450cam

In the catalytic cycle of cytochrome P450cam, after molecular oxygen binds as a ligand to the heme iron atom to yield a ferrous dioxygen complex, there are fast proton transfers that lead to the formation of the active species, Compound I (Cpd I), which are not well understood because they occur so rapidly. In the present work, the conversion of the ferric hydroperoxo complex (Cpd 0) to Cpd I has been investigated by combined quantum-mechanical/molecular-mechanical (QM/MM) calculations. The residues Asp(251) and Glu(366) are considered as proton sources. In mechanism I, a proton is transported to the distal oxygen atom of the hydroperoxo group via a hydrogen bonding network to form protonated Cpd 0 (prot-Cpd0: FeOOH(2)), followed by heterolytic O-O bond cleavage that generates Cpd I and water. Although a local minimum is found for prot-Cpd0 in the Glu(366) channel, it is very high in energy (more than 20 kcal/mol above Cpd 0) and the barriers for its decay are only 3-4 kcal/mol (both toward Cpd 0 and Cpd I). In mechanism II, an initial O-O bond cleavage followed by a concomitant proton and electron transfer yields Cpd I and water. The rate-limiting step in mechanism II is O-O cleavage with a barrier of about 13-14 kcal/mol. According to the QM/MM calculations, the favored low-energy pathway to Cpd I is provided by mechanism II in the Asp(251) channel. Cpd 0 and Cpd I are of similar energies, with a slight preference for Cpd I.     

Intrinsic acidity and redox properties of the adenine radical cation

The intrinsic chemical properties of the gaseous adenine radical cation were examined by using dual cell Fourier transform ion cyclotron resonance mass spectrometry. The adiabatic recombination energy of the radical cation (ionization energy of neutral adenine) was found by bracketing experiments to be 8.55 ± 0.1 eV (at 298 K; earlier literature values range from 8.3 to 8.9 eV). Based on this value, the heat of formation (ΔH_f²⁹⁸) of the adenine radical cation is estimated to be 246 ± 3 kcal/mol. The acidity (ΔH_acid²⁹⁸) of the adenine radical cation was bracketed to be 221 ± 2 kcal/mol. These thermochemical values suggest that the adenine radical cation reacts with neutral guanine by electron abstraction or proton transfer, with neutral cytosine by proton transfer, and via neither pathway with neutral thymine, molecular water or a sugar moiety of DNA (modeled by tetrahydrofuran). Experimental examination of the gas-phase reactivity of the adenine radical cation revealed a slow deuterium atom abstraction from perdeuterated tetrahydrofuran. Hence, in the absence of a nearby guanine or cytosine, the adenine radical cation may be able to abstract a hydrogen atom from a sugar moiety of DNA.     

A theoretical study of the iminoxyl/oxime self-exchange reaction. A five-center, cyclic proton-coupled electron transfer

In solution, the self-exchange reactions for oxygen-centered pi-radicals, e.g., PhO. + PhOH <==>PhOH + PhO., are known to occur with low activation enthalpies (E(a) approximately equal to 2 kcal/mol). For the PhO./PhOH couple and, we conclude, for other O-centered pi-radicals, exchange occurs by proton-coupled electron transfer (PCET) with the proton transferred between oxygen electron pairs while the electron migrates between oxygen orbitals orthogonal to the -O- - -H- - -O- transition state plane (Mayer et al. J. Am. Chem. Soc. 2002, 123, 11142). Iminoxyls, R(2)C=NO., are sigma-radicals with substantial spin density on the nitrogen. The R(2)C=NO./R(2)C=NOH self-exchange has a significant E(a) (Mendenhall et al. J. Am. Chem. Soc. 1973, 95, 627). For this exchange, DFT calculations have revealed a counterintuitive cisoid transition state in which the seven atoms, >C=NO- - -H- - -ON=C<, lie in a plane (R = H, Me) or, for steric reasons, two planes twisted at 45.2 degrees (R = Me(3)C). The planar transition state has the two N-O dipoles close to each other and pointing in the same direction and an O- - -H- - -O angle of 165.4 degrees . A transoid transition state for R = H lies 3.4 kcal/mol higher in energy than the cisoid despite a more favorable arrangement of the dipoles and a near linear O- - -H- - -O. It is concluded that iminoxyl/oxime self-exchange reactions occur by a five-center, cyclic PCET mechanism with the proton being transferred between electron pairs on the oxygens and the electron migrating between in-plane orbitals on the two nitrogens (R(N-N) = 2.65 A). The calculated E(a) values (8.8-9.9 kcal/mol) are in satisfactory agreement with the limited experimental data.     

NMR study of ligand exchange and electron self-exchange between oxo-centered trinuclear clusters [Fe3(μ3-O)(μ-O2CR)6(4-R'py)3](+/0)

The syntheses, single crystal X-ray structures, and magnetic properties of the homometallic μ?-oxo trinuclear clusters [Fe?(μ?-O)(μ-O?CCH?)?(4-Phpy)?](ClO?) (1) and [Fe?(μ?-O)(μ-O?CAd)?(4-Mepy)?](NO?) (2) are reported (Ad = adamantane). The persistence of the trinuclear structure within 1 and 2 in CD?Cl? and C?D?Cl? solutions in the temperature range 190-390 K is demonstrated by 1H NMR. An equilibrium between the mixed pyridine clusters [Fe?(μ?-O)(μ-O?CAd)?(4-Mepy)(3-x)(4-Phpy)(x)](NO?) (x = 0, 1, 2, 3) with a close to statistical distribution of these species is observed in CD?Cl? solutions. Variable-temperature NMR line-broadening made it possible to quantify the coordinated/free 4-Rpy exchanges at the iron centers of 1 and 2: k(ex)2?? = 6.5 ± 1.3 × 10?1 s?1, ΔH(?) = 89.47 ± 2 kJ mol?1, and ΔS(?) = +51.8 ± 6 J K?1 mol?1 for 1 and k(ex)2?? = 3.4 ± 0.5 × 10?1 s?1, ΔH(?) = 91.13 ± 2 kJ mol?1, and ΔS(?) = +51.9 ± 5 J K?1 mol?1 for 2. A limiting D mechanism is assigned for these ligand exchange reactions on the basis of first-order rate laws and positive and large entropies of activation. The exchange rates are 4 orders of magnitude slower than those observed for the ligand exchange on the reduced heterovalent cluster [Fe(III)?Fe(II)(μ?-O)(μ-O?CCH?)?(4-Phpy)?] (3). In 3, the intramolecular Fe(III)/Fe(II) electron exchange is too fast to be observed. At low temperatures, the 1/3 intermolecular second-order electron self-exchange reaction is faster than the 4-Phpy ligand exchange reactions on these two clusters, suggesting an outer-sphere mechanism: k?2?? = 72.4 ± 1.0 × 103 M?1 s?1, ΔH(?) = 18.18 ± 0.3 kJ mol?1, and ΔS(?) = -90.88 ± 1.0 J K?1 mol?1. The [Fe?(μ?-O)(μ-O?CCH?)?(4-Phpy)?](+/0) electron self-exchange reaction is compared with the more than 3 orders of magnitude faster [Ru?(μ?-O)(μ-O?CCH?)?(py)?](+/0) self-exchange reaction (ΔΔG(exptl)(?298) = 18.2 kJ mol?1). The theoretical estimated self-exchange rate constants for both processes compare reasonably well with the experimental values. The equilibrium constant for the formation of the precursor to the electron-transfer and the free energy of activation contribution for the solvent reorganization to reach the electron transfer step are taken to be the same for both redox couples. The larger ΔG(exptl)(?298) for the 1/3 iron self-exchange is attributed to the larger (11.1 kJ mol?1) inner-sphere reorganization energy of the 1 and 3 iron clusters in addition to a supplementary energy (6.1 kJ mol?1) which arises as a result of the fact that each encounter is not electron-transfer spin-allowed for the iron redox couple.     

The "somersault" mechanism for the p-450 hydroxylation of hydrocarbons. The intervention of transient inverted metastable hydroperoxides

A series of model theoretical calculations are described that suggest a new mechanism for the oxidation step in enzymatic cytochrome P450 hydroxylation of saturated hydrocarbons. A new class of metastable metal hydroperoxides is described that involves the rearrangement of the ground-state metal hydroperoxide to its inverted isomeric form with a hydroxyl radical hydrogen bonded to the metal oxide (MO-OH --> MO....HO). The activation energy for this somersault motion of the FeO-OH group is 20.3 kcal/mol for the P450 model porphyrin iron(III) hydroperoxide [Por(SH)Fe(III)-OOH(-)] to produce the isomeric ferryl oxygen hydrogen bonded to an *OH radical [Por(SH)Fe(III)-O....HO(-)]. This isomeric metastable hydroperoxide, the proposed primary oxidant in the P450 hydroxylation reaction, is calculated to be 17.8 kcal/mol higher in energy than the ground-state iron(III) hydroperoxide Cpd 0. The first step of the proposed mechanism for isobutane oxidation is abstraction of a hydrogen atom from the C-H bond of isobutane by the hydrogen-bonded hydroxyl radical to produce a water molecule strongly hydrogen bonded to anionic Cpd II. The hydroxylation step involves a concerted but nonsynchronous transfer of a hydrogen atom from this newly formed, bound, water molecule to the ferryl oxygen with a concomitant rebound of the incipient *OH radical to the carbon radical of isobutane to produce the C-O bond of the final product, tert-butyl alcohol. The TS for the oxygen rebound step is 2 kcal/mol lower in energy than the hydrogen abstraction TS (DeltaE() = 19.5 kcal/mol). The overall proposed new mechanism is consistent with a lot of the ancillary experimental data for this enzymatic hydroxylation reaction.     

Quantum mechanics/molecular mechanics investigation of the chemical reaction in Dpo4 reveals water-dependent pathways and requirements for active site reorganization

The nucleotidyl-transfer reaction coupled with the conformational transitions in DNA polymerases is critical for maintaining the fidelity and efficiency of DNA synthesis. We examine here the possible reaction pathways of a Y-family DNA polymerase, Sulfolobus solfataricus DNA polymerase IV (Dpo4), for the correct insertion of dCTP opposite 8-oxoguanine using the quantum mechanics/molecular mechanics (QM/MM) approach, both from a chemistry-competent state and a crystal closed state. The latter examination is important for understanding pre-chemistry barriers to interpret the entire enzyme mechanism, since the crystal closed state is not an ideal state for initiating the chemical reaction. The most favorable reaction path involves initial deprotonation of O3'H via two bridging water molecules to O1A, overcoming an overall potential energy barrier of approximately 20.0 kcal/mol. The proton on O1A-P(alpha) then migrates to the gamma-phosphate oxygen of the incoming nucleotide as O3' attacks P(alpha), and the P(alpha)-O3A bond breaks. The other possible pathway in which the O3'H proton is transferred directly to O1A on P(alpha) has an overall energy barrier of 25.0 kcal/mol. In both reaction paths, the rate-limiting step is the initial deprotonation, and the trigonal-bipyramidal configuration for P(alpha) occurs during the concerted bond formation (O3'-P(alpha)) and breaking (P(alpha)-O3A), indicating the associative nature of the chemical reaction. In contrast, the Dpo4/DNA complex with an imperfect active-site geometry corresponding to the crystal state must overcome a much higher activation energy barrier (29.0 kcal/mol) to achieve a tightly organized site due to hindered O3'H deprotonation stemming from larger distances and distorted conformation of the proton acceptors. This significant difference demonstrates that the pre-chemistry reorganization in Dpo4 costs approximately 4.0 to 9.0 kcal/mol depending on the primer terminus environment. Compared to the higher fidelity DNA polymerase beta from the X-family, Dpo4 has a higher chemical reaction barrier (20.0 vs 15.0 kcal/mol) due to the more solvent-exposed active site.     

Evaluation of electron-withdrawing group effects on heme binding in designed proteins: implications for heme a in cytochrome c oxidase

Heme a, the metalloporphyrin cofactor unique to cytochrome c oxidases, differs from the more common heme b by two chemical modifications, a C-2 hydroxyethylfarnesyl group and a C-8 formyl group. To elucidate a role of the C-8 formyl group, we compare the heme affinity, spectroscopy, and electrochemistry of a heme a mimic, Fe(diacetyldeuterioporphyrin IX) or Fe(DADPIX), with heme b, Fe(protoporphryrin IX) or Fe(PPIX), incorporated into a designed heme protein. The [Delta7-H3m]2 protein ligand, or maquette, selected for this study contains two equivalent bis-(3-methyl-L-histidine) heme binding sites within a four-alpha-helix bundle scaffold. The spectroscopic data on Fe(PPIX) and Fe(DADPIX) bound to [Delta7-H3m]2 demonstrate that these complexes are excellent synthetic analogues for natural cytochromes b and a, respectively. Comparison of the spectroscopic, electrochemical, and equilibrium thermodynamic data measured for the Fe(PPIX)-[Delta7-H3m]2 maquette with the previously reported Fe(PPIX)-[Delta7-His]2 complex demonstrates that changing the heme axial ligands to 3-methyl-L-histidine from L-histidine does not alter the resulting heme protein properties significantly in either oxidation state. Heme binding studies demonstrate that [Delta7-H3m]2 binds two ferrous Fe(DADPIX) or Fe(PPIX) moieties with similar dissociation constant values. However, in the ferric state, the data show that [Delta7-H3m]2 only binds a single Fe(DADPIX) and that one 2500-fold weaker than oxidized Fe(PPIX). The data demonstrate that the 4.6 kcal mol(-1) weakened affinity of [Delta7-H3m]2 for oxidized Fe(DADPIX) results in the majority of the 160 mV, 3.7 kcal mol(-1), positive shift in the heme reduction potential relative to Fe(PPIX). These data indicate that a role of the formyl group on heme a is to raise the iron reduction potential, thus making it a better electron acceptor, but that it does so by destabilizing the affinity of bis-imidazole sites for the ferric state.     

Experimental and computational studies of hydrogen bonding and proton transfer to [Cp*Fe(dppe)H

The present contribution reports experimental and computational investigations of the interaction between [Cp*Fe(dppe)H] and different proton donors (HA). The focus is on the structure of the proton transfer intermediates and on the potential energy surface of the proton transfer leading to the dihydrogen complex [Cp*Fe(dppe)(H2)]+. With p-nitrophenol (PNP) a UV/Visible study provides evidence of the formation of the ion-pair stabilized by a hydrogen bond between the nonclassical cation [Cp*Fe(dppe)(H2)]+ and the homoconjugated anion ([AHA]-). With trifluoroacetic acid (TFA), the hydrogen-bonded ion pair containing the simple conjugate base (A-) in equilibrium with the free ions is observed by IR spectroscopy when using a deficit of the proton donor. An excess leads to the formation of the homoconjugated anion. The interaction with hexafluoroisopropanol (HFIP) was investigated quantitatively by IR spectroscopy and by 1H and 31P NMR spectroscopy at low temperatures (200-260 K) and by stopped-flow kinetics at about room temperature (288-308 K). The hydrogen bond formation to give [Cp*Fe(dppe)H]HA is characterized by DeltaH degrees =-6.5+/-0.4 kcal mol(-1) and DeltaS degrees = -18.6+/-1.7 cal mol(-1) K(-1). The activation barrier for the proton transfer step, which occurs only upon intervention of a second HFIP molecule, is DeltaH(not equal) = 2.6+/-0.3 kcal mol(-1) and DeltaS(not equal) = -44.5+/-1.1 cal mol(-1) K(-1). The computational investigation (at the DFT/B3 LYP level with inclusion of solvent effects by the polarizable continuum model) reproduces all the qualitative findings, provided the correct number of proton donor molecules are used in the model. The proton transfer process is, however, computed to be less exothermic than observed in the experiment.     

Elimination of water from the carboxyl group of GlyGlyH+

The elimination of water from the carboxyl group of protonated diglycine has been investigated by density functional theory calculations. The resulting structure is identical to the b(2) ion formed in the mass spectrometric fragmentation of protonated peptides (therefore named "b2" in this study). The most stable geometry of the fragment ion ("b2") is an O-protonated diketopiperazine. However, its formation is kinetically disfavored as it requires a free energy of 58.2 kcal/mol. The experimentally observed N-protonated oxazolone is 3.0 kcal/mol less stable. The lowest energy pathway for the formation of the "b2" ion requires a free energy of 37.5 kcal/mol and involves the proton transfer from the amide oxygen of protonated diglycine to the hydroxyl oxygen. Fragmentation initiated by proton transfer from the terminal nitrogen has also a comparable free energy of activation (39.4 kcal/mol). Proton transfer initiating the fragmentation, from the highly basic terminal nitrogen or amide oxygen to the less basic hydroxyl oxygen is feasible at energies reached in usual mass spectrometric experiments. Amide N-protonated diglycine structures are precursors of mainly y(1) ions rather than "b2" ions. In the lowest energy fragmentation channels, proton transfer to the hydroxylic oxygen, bond breaking and formation of an oxazolone ring occur concertedly but asynchronously. Proton transfer to hydroxyl oxygen and cleavage of the corresponding C-O bond take place at the early stages of the fragmentation step, while ring closure to form an oxazolone geometry occurs at the later stages of the transition. The experimentally observed low kinetic energy release is expected to be due to the existence of a strongly hydrogen bonded protonated oxazolone-water complex in the exit channel. Whereas the threshold energy for "b2" ion formation (37.1 kcal/mol) is lower than for the y(1) ion (38.4 kcal/mol), the former requires a tight transition state with an activation entropy, DeltaS++ = -1.2 cal/mol.K and the latter has a loose transition state with DeltaS++ = +8.8 cal/mol.K. This leads to y(1) being the major fragment ion over a wide energy range.     

Stabilization of very rare tautomers of 1-methylcytosine by an excess electron

We characterized valence anionic states of 1-methylcytosine using various electronic structure methods. We found that the most stable valence anion is related to neither the canonical amino-oxo nor a rare imino-oxo tautomer, in which a proton is transferred from the N4 to N3 atom. Instead, it is related to an imino-oxo tautomer, in which the C5 atom is protonated. This anion is characterized by an electron vertical detachment energy (VDE) of 2.12 eV and it is more stable than the anion based on the canonical tautomer by 1.0 kcal/mol. The latter is characterized by a VDE of 0.31 eV. Another unusual low-lying imino-oxo tautomer with a VDE of 3.60 eV has the C6 atom protonated and is 3.6 kcal/mol less stable than the anion of the canonical tautomer. All these anionic states are adiabatically unbound with respect to the canonical amino-oxo neutral, with the instability of 5.8 kcal/mol for the most stable valence anion. The mechanism of formation of anionic tautomers with carbon atoms protonated may involve intermolecular proton transfer or dissociative electron attachment to the canonical neutral tautomer followed by a barrier-free attachment of a hydrogen atom to the C5 or C6 atom. The six-member ring structure of anionic tautomers with carbon atoms protonated is unstable upon an excess electron detachment. Indeed the neutral systems collapse without a barrier to a linear or a bicyclo structure, which might be viewed as lesions to DNA or RNA. Within the PCM hydration model, the anions become adiabatically bound with respect to the corresponding neutrals, and the two most stable tautomers have a carbon atom protonated.     

Distinct reaction pathways followed upon reduction of oxy-heme oxygenase and oxy-myoglobin as characterized by Mössbauer spectroscopy

Activation of O(2) by heme-containing monooxygenases generally commences with the common initial steps of reduction to the ferrous heme and binding of O(2) followed by a one-electron reduction of the O(2)-bound heme. Subsequent steps that generate reactive oxygen intermediates diverge and reflect the effects of protein control on the reaction pathway. In this study, M?ssbauer and EPR spectroscopies were used to characterize the electronic states and reaction pathways of reactive oxygen intermediates generated by 77 K radiolytic cryoreduction and subsequent annealing of oxy-heme oxygenase (HO) and oxy-myoglobin (Mb). The results confirm that one-electron reduction of (Fe(II)-O(2))HO is accompanied by protonation of the bound O(2) to generate a low-spin (Fe(III)-O(2)H(-))HO that undergoes self-hydroxylation to form the alpha-meso-hydroxyhemin-HO product. In contrast, one-electron reduction of (Fe(II)-O(2))Mb yields a low-spin (Fe(III)-O(2)(2-))Mb. Protonation of this intermediate generates (Fe(III)-O(2)H(-))Mb, which then decays to a ferryl complex, (Fe(IV)=O(2-))Mb, that exhibits magnetic properties characteristic of the compound II species generated in the reactions of peroxide with heme peroxidases and with Mb. Generation of reactive high-valent states with ferryl species via hydroperoxo intermediates is believed to be the key oxygen-activation steps involved in the catalytic cycles of P450-type monooxygenases. The M?ssbauer data presented here provide direct spectroscopic evidence supporting the idea that ferric-hydroperoxo hemes are indeed the precursors of the reactive ferryl intermediates. The fact that a ferryl intermediate does not accumulate in HO underscores the determining role played by protein structure in controlling the reactivity of reaction intermediates.     

Photodissociation of diiodomethane in acetonitrile solution and fragment recombination into iso-diiodomethane studied with ab initio molecular dynamics simulations

Series of hydrates of the most stable glycine-H+/2H2+ in the gas phase are presented at the B3LYP level. The results show that only the amino hydrogens and hydroxyl hydrogens can be monohydrated for the glycine-H+, and the amino hydrogens are preferred. The H6(O4) of glycine-2H2+ is the best site for a water molecule to attach, i.e., the corresponding hydrate is the most stable one among its isomers. Calculations reveal that the binding energies of hydrated hydrogens decrease relative to their counterparts in the isolated glycine-H+/2H2+ complexes and they are positive values and without proton transfer except those of monohydrated glycine-2H2+ complexes with the combination modes of H3O+...(glycine-H+). The complex H3O+...(glycine-H+) is formed by the combination of a H2O molecule and one hydroxyl-site proton of glycine-2H2+, and with the proton transfer to H2O. Here the interaction between the proton of H3O+ and the glycine-H+ mainly depends on an electronic one instead of an initial covalent one of the isolated glycine-2H2+. The generation of the bond between the H3O+ and the glycine-H+ makes the energy of the complex higher than the energy sum of its two separated species (or two reactants of the complex), just like the case of M+...(glycine-H+) bond (M = Li,Na). The observation can explain satisfactorily why the combinations of both a proton and an alkali ion or two alkali ions to a glycine molecule can make the corresponding complex hold reservation energy bond(s), while the combination of two protons and a glycine in our previous work cannot [H. Ai et al., J. Chem. Phys. 117, 7593 (2002)]. For the glycine-2H2+, monohydration at the any site of its amino hydrogens can make the binding strength of any other neighboring proton (hydrogens) stronger relative to its counterpart in the isolated glycine-2H2+. Further hydration, especially at the site of either of hydroxyl hydrogens, would disfavor the reservation energy of the system.     

Kinetics and mechanism of ketone enolization mediated by magnesium bis(hexamethyldisilazide)

Magnesium bis(hexamethyldisilazide), Mg(HMDS)(2), reacts with substoichiometric amounts of propiophenone in toluene solution at ambient temperature to form a 74:26 mixture of the enolates (E)- and (Z)-[(HMDS)(2)Mg(2)(mu-HMDS){mu-OC(Ph)=CHCH(3)}], (E)-1 and (Z)-1, which contain a pair of three-coordinate metal centers bridged by an amide and an enolate group. The compositions of (E)-1 and (Z)-1 were confirmed by solution NMR studies and also by crystallographic characterization in the solid state. Rate studies using UV-vis spectroscopy reveal the rapid and complete formation of a reaction intermediate, 2, between the ketone and magnesium, which undergoes first-order decay with rate constants independent of the concentration of excess Mg(HMDS)(2) (DeltaH++ = 17.2 +/- 0.8 kcal/mol, DeltaS++ = -11 +/- 3 cal/mol.K). The intermediate 2 has been characterized by low-temperature (1)H NMR, diffusion-ordered NMR, and IR spectroscopy and investigated by computational studies, all of which are consistent with the formulation of 2 as a three-coordinate monomer, (HMDS)(2)Mg{eta(1)-O=C(Ph)CH(2)CH(3)}. Further support for this structure is provided by the synthesis and structural characterization of two model ketone complexes, (HMDS)(2)Mg(eta(1)-O=C(t)Bu(2)) (3) and (HMDS)(2)Mg{eta(1)-O=C((t)Bu)Ph} (4). A large primary deuterium isotope effect (k(H)/k(D) = 18.9 at 295 K) indicates that proton transfer is the rate-limiting step of the reaction. The isotope effect displays a strong temperature dependence, indicative of tunneling. In combination, these data support the mechanism of enolization proceeding through the single intermediate 2 via intramolecular proton transfer from the alpha carbon of the bound ketone to the nitrogen of a bound hexamethyldisilazide.     

The possibility of the decomposition of 2'-deoxyribose moiety of thymidine induced by the low energy electron attachment

To evaluate the possibility of the decomposition of 2-deoxyribose moiety of thymidine induced by low energy electrons (LEE) attachment, the transition states and the energy barriers of the bond breaking processes of the ribose of the nucleoside have been studied theoretically by applying the density functional theory with the double zeta basis sets (DZP++). The energy barriers for the breakage of the C-C bonds (C(1')-C(2'), C(2')-C(3'), C(3')-C(4'), and C(4')-C(5')) of the ribose group of the radical anion of thymidine are found to be high (ca. 42-57 kcal/mol). The total energies of the C-C bond-broken products are significantly higher than that of the radical anion dT(*-). The decomposition of dT(*-) through the C-C bond rupture is unlikely to take place. The rupture of the C(1')-O(4') bond of dT(*-) needs an activation energy as low as 10.4 kcal/mol. However, the reversed reaction (C(1')-O(4') bond formation) needs the activation energy low as 0.3 kcal/mol. Therefore, the intermediate product LM1(C1')-(O4') is unlikely to be stable and the C(1')-O(4') bond-broken is not favored. The activation energy of the C(4')-O(4') bond rupture process amounts to 20.5 kcal/mol. The total energy of the C(4')-O(4') bond broken product is about 6.5 kcal/mol lower than that of the reactant dT(*-). The subsequent N1-glycosidic bond breaking process is found to have a very low energy barrier. Therefore, the LEE-induced base release through the C(4')-O(4') bond rupture might be a possible pathway.