Quantum Chemical Studies of Proton-Coupled Electron Transfer in Metalloenzymes

Origin and Evolution of Water Oxidation before the Last Common Ancestor of the Cyanobacteria

Photosystem II, the water oxidizing enzyme, altered the course of evolution by filling the atmosphere with oxygen. Here, we reconstruct the origin and evolution of water oxidation at an unprecedented level of detail by studying the phylogeny of all D1 subunits, the main protein coordinating the water oxidizing cluster (Mn₄CaO₅) of Photosystem II. We show that D1 exists in several forms making well-defined clades, some of which could have evolved before the origin of water oxidation and presenting many atypical characteristics. The most ancient form is found in the genome of Gloeobacter kilaueensis JS-1 and this has a C-terminus with a higher sequence identity to D2 than to any other D1. Two other groups of early evolving D1 correspond to those expressed under prolonged far-red illumination and in darkness. These atypical D1 forms are characterized by a dramatically different Mn₄CaO₅ binding site and a Photosystem II containing such a site may assemble an unconventional metal cluster. The first D1 forms with a full set of ligands to the Mn₄CaO₅ cluster are grouped with D1 proteins expressed only under low oxygen concentrations and the latest evolving form is the dominant type of D1 found in all cyanobacteria and plastids. In addition, we show that the plastid ancestor had a D1 more similar to those in early branching Synechococcus. We suggest each one of these forms of D1 originated from transitional forms at different stages toward the innovation and optimization of water oxidation before the last common ancestor of all known cyanobacteria.

Computational organic chemistry: bridging theory and experiment in establishing the mechanisms of chemical reactions

Understanding the mechanisms of chemical reactions, especially catalysis, has been an important and active area of computational organic chemistry, and close collaborations between experimentalists and theorists represent a growing trend. This Perspective provides examples of such productive collaborations. The understanding of various reaction mechanisms and the insight gained from these studies are emphasized. The applications of various experimental techniques in elucidation of reaction details as well as the development of various computational techniques to meet the demand of emerging synthetic methods, e.g., C-H activation, organocatalysis, and single electron transfer, are presented along with some conventional developments of mechanistic aspects. Examples of applications are selected to demonstrate the advantages and limitations of these techniques. Some challenges in the mechanistic studies and predictions of reactions are also analyzed.

The role of molecular oxygen in the iron(iii)-promoted oxidative dehydrogenation of amines

A mechanistic study is presented of the oxidative dehydrogenation of the iron(iii) complex [Fe^IIIL³]³⁺, 1, (L³ = 1,9-bis(2′-pyridyl)-5-[(ethoxy-2′′-pyridyl)methyl]-2,5,8-triazanonane) in ethanol in the presence of molecular oxygen.

The pathway for serial proton supply to the active site of nitrogenase: enhanced density functional modeling of the Grotthuss mechanism

Proton translocation along a chain of eight waters to the active site of nitrogenase is described in detail, using density functional simulations with a 269 atom system that includes surrounding amino acids.

A quantum chemical study of the ω-transaminase reaction mechanism

The detailed half-transamination mechanism of Chromobacterium violaceum ω-transaminase is investigated by means of density functional theory calculations.

How cytochrome c oxidase can pump four protons per oxygen molecule at high electrochemical gradient

Experiments have shown that the A-family cytochrome c oxidases pump four protons per oxygen molecule, also at a high electrochemical gradient. This has been considered a puzzle, since two of the reduction potentials involved, Cu(II) and Fe(III), were estimated from experiments to be too low to afford proton pumping at a high gradient. The present quantum mechanical study (using hybrid density functional theory) suggests a solution to this puzzle. First, the calculations show that the charge compensated Cu(II) potential for CuB is actually much higher than estimated from experiment, of the same order as the reduction potentials for the tyrosyl radical and the ferryl group, which are also involved in the catalytic cycle. The reason for the discrepancy between theory and experiment is the very large uncertainty in the experimental observations used to estimate the equilibrium potentials, mainly caused by the lack of methods for direct determination of reduced CuB. Second, the calculations show that a high energy metastable state, labeled EH, is involved during catalytic turnover. The EH state mixes the low reduction potential of Fe(III) in heme a3 with another, higher potential, here suggested to be that of the tyrosyl radical, resulting in enough exergonicity to allow proton pumping at a high gradient. In contrast, the corresponding metastable oxidized state, OH, is not significantly higher in energy than the resting state, O. Finally, to secure the involvement of the high energy EH state it is suggested that only one proton is taken up via the K-channel during catalytic turnover.

Mechanism for O-O bond formation in a biomimetic tetranuclear manganese cluster--A density functional theory study

Density functional theory calculations have been used to study the reaction mechanism of water oxidation catalyzed by a tetranuclear Mn-oxo cluster Mn4O4L6 (L=(C6H4)2PO4(-)). It is proposed that the O-O bond formation mechanism is different in the gas phase and in a water solution. In the gas phase, upon phosphate ligand dissociation triggered by light absorption, the O-O bond formation starting with both the Mn4(III,III,IV,IV) and Mn4(III,IV,IV,IV) oxidation states has to take place via direct coupling of two bridging oxo groups. The calculated barriers are 42.3 and 37.1 kcal/mol, respectively, and there is an endergonicity of more than 10 kcal/mol. Additional photons are needed to overcome these large barriers. In water solution, water binding to the two vacant sites of the Mn ions, again after phosphate dissociation triggered by light absorption, is thermodynamically and kinetically very favorable. The catalytic cycle is suggested to start from the Mn4(III,III,III,IV) oxidation state. The removal of three electrons and three protons leads to the formation of a Mn4(III,IV,IV,IV)-oxyl radical complex. The O-O bond formation then proceeds via a nucleophilic attack of water on the Mn(IV)-oxyl radical assisted by a Mn-bound hydroxide that abstracts a proton during the attack. This step was calculated to be rate-limiting with a total barrier of 29.2 kcal/mol. This is followed by proton-coupled electron transfer, O2 release, and water binding to start the next catalytic cycle.

With or without light: comparing the reaction mechanism of dark-operative protochlorophyllide oxidoreductase with the energetic requirements of the light-dependent protochlorophyllide oxidoreductase

The addition of two electrons and two protons to the C₁₇=C₁₈ bond in protochlorophyllide is catalyzed by a light-dependent enzyme relying on NADPH as electron donor, and by a light-independent enzyme bearing a (Cys)₃Asp-ligated [4Fe–4S] cluster which is reduced by cytoplasmic electron donors in an ATP-dependent manner and then functions as electron donor to protochlorophyllide. The precise sequence of events occurring at the C₁₇=C₁₈ bond has not, however, been determined experimentally in the dark-operating enzyme. In this paper, we present the computational investigation of the reaction mechanism of this enzyme at the B3LYP/6-311+G(d,p)//B3LYP/6-31G(d) level of theory. The reaction mechanism begins with single-electron reduction of the substrate by the (Cys)₃Asp-ligated [4Fe–4S], yielding a negatively-charged intermediate. Depending on the rate of Fe–S cluster re-reduction, the reaction either proceeds through double protonation of the single-electron-reduced substrate, or by alternating proton/electron transfer. The computed reaction barriers suggest that Fe–S cluster re-reduction should be the rate-limiting stage of the process. Poisson–Boltzmann computations on the full enzyme–substrate complex, followed by Monte Carlo simulations of redox and protonation titrations revealed a hitherto unsuspected pH-dependence of the reaction potential of the Fe–S cluster. Furthermore, the computed distributions of protonation states of the His, Asp and Glu residues were used in conjuntion with single-point ONIOM computations to obtain, for the first time, the influence of all protonation states of an enzyme on the reaction it catalyzes. Despite exaggerating the ease of reduction of the substrate, these computations confirmed the broad features of the reaction mechanism obtained with the medium-sized models, and afforded valuable insights on the influence of the titratable amino acids on each reaction step. Additional comparisons of the energetic features of the reaction intermediates with those of common biochemical redox intermediates suggest a surprisingly simple explanation for the mechanistic differences between the dark-catalyzed and light-dependent enzyme reaction mechanisms.

Fourier transform infrared detection of a polarizable proton trapped between photooxidized tyrosine YZ and a coupled histidine in photosystem II: relevance to the proton transfer mechanism of water oxidation

The redox-active tyrosine YZ (D1-Tyr161) in photosystem II (PSII) functions as an immediate electron acceptor of the Mn4Ca cluster, which is the catalytic center of photosynthetic water oxidation. YZ is also located in the hydrogen bond network that connects the Mn4Ca cluster to the lumen and hence is possibly related to the proton transfer process during water oxidation. To understand the role of YZ in the water oxidation mechanism, we have studied the hydrogen bonding interactions of YZ and its photooxidized neutral radical YZ(•) together with the interaction of the coupled His residue, D1-His190, using light-induced Fourier transform infrared (FTIR) difference spectroscopy. The YZ(•)-minus-YZ FTIR difference spectrum of Mn-depleted PSII core complexes exhibited a broad positive feature around 2800 cm(-1), which was absent in the corresponding spectrum of another redox-active tyrosine YD (D2-Tyr160). Analyses by (15)N and H/D substitutions, examination of the pH dependence, and density functional theory and quantum mechanics/molecular mechanics (QM/MM) calculations showed that this band arises from the N-H stretching vibration of the protonated cation of D1-His190 forming a charge-assisted strong hydrogen bond with YZ(•). This result provides strong evidence that the proton released from YZ upon its oxidation is trapped in D1-His190 and a positive charge remains on this His. The broad feature of the ~2800 cm(-1) band reflects a large proton polarizability in the hydrogen bond between YZ(•) and HisH(+). QM/MM calculations further showed that upon YZ oxidation the hydrogen bond network is rearranged and one water molecule moves toward D1-His190. From these data, a novel proton transfer mechanism via YZ(•)-HisH(+) is proposed, in which hopping of the polarizable proton of HisH(+) to this water triggers the transfer of the proton from substrate water to the luminal side. This proton transfer mechanism could be functional in the S2 → S3 transition, which requires proton release before electron transfer because of an excess positive charge on the Mn4Ca cluster.

Applications of density functional theory to iron-containing molecules of bioinorganic interest

The past decades have seen an explosive growth in the application of density functional theory (DFT) methods to molecular systems that are of interest in a variety of scientific fields. Owing to its balanced accuracy and efficiency, DFT plays particularly useful roles in the theoretical investigation of large molecules. Even for biological molecules such as proteins, DFT finds application in the form of, e.g., hybrid quantum mechanics and molecular mechanics (QM/MM), in which DFT may be used as a QM method to describe a higher prioritized region in the system, while a MM force field may be used to describe remaining atoms. Iron-containing molecules are particularly important targets of DFT calculations. From the viewpoint of chemistry, this is mainly because iron is abundant on earth, iron plays powerful (and often enigmatic) roles in enzyme catalysis, and iron thus has the great potential for biomimetic catalysis of chemically difficult transformations. In this paper, we present a brief overview of several recent applications of DFT to iron-containing non-heme synthetic complexes, heme-type cytochrome P450 enzymes, and non-heme iron enzymes, all of which are of particular interest in the field of bioinorganic chemistry. Emphasis will be placed on our own work.

Proton pumping in cytochrome c oxidase: energetic requirements and the role of two proton channels

Cytochrome c oxidase is a superfamily of membrane bound enzymes catalyzing the exergonic reduction of molecular oxygen to water, producing an electrochemical gradient across the membrane. The gradient is formed both by the electrogenic chemistry, taking electrons and protons from opposite sides of the membrane, and by proton pumping across the entire membrane. In the most efficient subfamily, the A-family of oxidases, one proton is pumped in each reduction step, which is surprising considering the fact that two of the reduction steps most likely are only weakly exergonic. Based on a combination of quantum chemical calculations and experimental information, it is here shown that from both a thermodynamic and a kinetic point of view, it should be possible to pump one proton per electron also with such an uneven distribution of the free energy release over the reduction steps, at least up to half the maximum gradient. A previously suggested pumping mechanism is developed further to suggest a reason for the use of two proton transfer channels in the A-family. Since the rate of proton transfer to the binuclear center through the D-channel is redox dependent, it might become too slow for the steps with low exergonicity. Therefore, a second channel, the K-channel, where the rate is redox-independent is needed. A redox-dependent leakage possibility is also suggested, which might be important for efficient energy conservation at a high gradient. A mechanism for the variation in proton pumping stoichiometry over the different subfamilies of cytochrome oxidase is also suggested. This article is part of a Special Issue entitled: 18th European Bioenergetic Conference.

Discovery of potent inhibitor for matrix metalloproteinase-9 by pharmacophore based modeling and dynamics simulation studies

Matrix metalloproteinase-9 (MMP-9) is an attractive target for anticancer therapy. In the present study ligand based pharmacophore modeling was performed to elucidate the structural elements for a diverse class of MMP-9 inhibitors. The pharmacophore model was validated through Güner-Henry (GH) scoring method. The final pharmacophore model consisted of three hydrogen bond acceptors (HBA), and two ring aromatic regions (RA). This model was utilized to screen the natural compound database to seek novel compounds as MMP-9 inhibitors. The identified hits were validated using molecular docking and molecular dynamics simulation studies. Finally, one compound named Hinokiflavone from Juniperus communis had high binding free energy of -26.54kJ/mol compared with the known inhibitors of MMP-9. Cytotoxicity for hinokiflavone was evaluated by MTT assay. Inhibition of MMP-9 in the presence of hinokiflavone was detected by gelatin zymography and gelatinolytic inhibition assay. Results revealed that the natural compounds derived based on the developed pharmacophore model would be useful for further design and development of MMP-9 inhibitors.

A dominant homolytic O-Cl bond cleavage with low-spin triplet-state Fe(IV)=O formed is revealed in the mechanism of heme-dependent chlorite dismutase

Chlorite dismutase (Cld) is a heme-dependent enzyme that catalyzes the decomposition of toxic chlorite (ClO2(-)) into innocuous chloride and O2. In this paper, using the hybrid B3LYP density functional theory (DFT) method including dispersion interactions, the Cld reaction mechanism has been studied with a chemical model constructed on the X-ray crystal structure. The calculations indicate that the reaction proceeds along a stepwise pathway in the doublet state, i.e. a homolytic O-Cl bond cleavage of the substrate leading to an O-Fe(heme) species and a ClO˙ radical, followed by a rebinding O-O bond formation between them. The O-Fe(heme) species is demonstrated to be a low-spin triplet-state Fe(IV)=O diradicaloid. A low-spin singlet-state Fe(IV)=O is much less stable than the former, with an energy difference of 9.2 kcal mol(-1). The O-Cl bond cleavage is rate-limiting with a barrier of 10.6 kcal mol(-1), in good agreement with the experimental reaction rate of 2.0 × 10(5) s(-1). Furthermore, a heterolytic O-Cl bond dissociation in the initial step is shown to be unreachable, which ensures the high efficiency of the Cld enzyme by avoiding the generation of chlorate byproduct observed in the reactions of synthetic Fe porphyrins. Also, the pathways in the quartet and sextet states are unfavorable for the Cld reaction. The present results reveal a detailed mechanism III (defined in the text) including an interesting di-radical intermediate composed of a low-spin triplet-state Fe(IV)=O and a ClO˙ radical. Compared to a competitive heterolytic Cl-O cleavage in synthetic Fe porphyrins, the revelation of the domination of homolysis in Cld indicates not only the high efficiency of enzyme, but also the sensitivity of a heme and the significance of the enzymatic active-site surroundings (the His170 and Arg183 residues in the present case), which gives more insights into heme chemistry.

Molecular mechanisms for generating transmembrane proton gradients

Membrane proteins use the energy of light or high energy substrates to build a transmembrane proton gradient through a series of reactions leading to proton release into the lower pH compartment (P-side) and proton uptake from the higher pH compartment (N-side). This review considers how the proton affinity of the substrates, cofactors and amino acids are modified in four proteins to drive proton transfers. Bacterial reaction centers (RCs) and photosystem II (PSII) carry out redox chemistry with the species to be oxidized on the P-side while reduction occurs on the N-side of the membrane. Terminal redox cofactors are used which have pKas that are strongly dependent on their redox state, so that protons are lost on oxidation and gained on reduction. Bacteriorhodopsin is a true proton pump. Light activation triggers trans to cis isomerization of a bound retinal. Strong electrostatic interactions within clusters of amino acids are modified by the conformational changes initiated by retinal motion leading to changes in proton affinity, driving transmembrane proton transfer. Cytochrome c oxidase (CcO) catalyzes the reduction of O2 to water. The protons needed for chemistry are bound from the N-side. The reduction chemistry also drives proton pumping from N- to P-side. Overall, in CcO the uptake of 4 electrons to reduce O2 transports 8 charges across the membrane, with each reduction fully coupled to removal of two protons from the N-side, the delivery of one for chemistry and transport of the other to the P-side.

Why is the molybdenum-substituted tungsten-dependent formaldehyde ferredoxin oxidoreductase not active? A quantum chemical study

Formaldehyde ferredoxin oxidoreductase is a tungsten-dependent enzyme that catalyzes the oxidative degradation of formaldehyde to formic acid. The molybdenum ion can be incorporated into the active site to displace the tungsten ion, but is without activity. Density functional calculations have been employed to understand the incapacitation of the enzyme caused by molybdenum substitution. The calculations show that the enzyme with molybdenum (Mo-FOR) has higher redox potential than that with tungsten, which makes the formation of the Mo(VI)=O complex endothermic by 14 kcal/mol. Following our previously suggested mechanism for this enzyme, the formaldehyde substrate oxidation was also investigated for Mo-FOR using the same quantum-mechanics-only model, except for the displacement of tungsten by molybdenum. The calculations demonstrate that formaldehyde oxidation occurs via a sequential two-step mechanism. Similarly to the tungsten-catalyzed reaction, the Mo(VI)=O species performs the nucleophilic attack on the formaldehyde carbon, followed by proton transfer in concert with two-electron reduction of the metal center. The first step is rate-limiting, with a total barrier of 28.2 kcal/mol. The higher barrier is mainly due to the large energy penalty for the formation of the Mo(VI)=O species.

Computational insight into the initial steps of the Mars-van Krevelen mechanism: electron transfer and surface defects in the reduction of polyoxometalates

Metal oxides as a rule oxidize and oxygenate substrates via the Mars-van Krevelen mechanism. A well-defined α-Keggin polyoxometalate, H(5)PV(2)Mo(10)O(40), can be viewed as an analogue of discrete structure that reacts via the Mars-van Krevelen mechanism both in solution and in the gas phase. Guided by previous experimental observations, we have studied the key intermediates on the reaction pathways of its reduction by various compounds using high-level DFT calculations. These redox reactions of polyoxometalates require protons, and thus such complexes were explicitly considered. First, the energetics of outer-sphere proton and electron transfer as well as coupled proton and electron transfer were calculated for seven substrates. This was followed by identification of possible key intermediates on the subsequent reaction pathways that feature displacement of the metal atom from the Keggin structure and coordinatively unsaturated sites on the H(5)PV(2)Mo(10)O(40) surface. Such metal defects are favored at vanadium sites. For strong reducing agents the initial outer-sphere electron transfer, alone or possibly coupled with proton transfer, facilitates formation of metal defects. Subsequent coordination allows for formation of reactive ensembles on the catalyst surface, for which the selective oxygen-transfer step becomes feasible. Weak reducing agents do not facilitate defect formation by outer-sphere electron and/or proton transfers, and thus formation of metal defect structures prior to the substrate activation is suggested as an initial step. Calculated geometries and energies of metal defect structures support experimentally observed intermediates and demonstrate the complex nature of the Mars-van Krevelen mechanism.

The Third Dimension of a More O'Ferrall-Jencks Diagram for Hydrogen Atom Transfer in the Isoelectronic Hydrogen Exchange Reactions of (PhX)(2)H(•) with X = O, NH, and CH(2)

A critical element in theoretical characterization of the mechanism of proton-coupled electron transfer (PCET) reactions, including hydrogen atom transfer (HAT), is the formulation of the electron and proton localized diabatic states, based on which a More O'Ferrall-Jencks diagram can be represented to determine the step-wise and concerted nature of the reaction. Although the More O'Ferrall-Jencks diabatic states have often been used empirically to develop theoretical models for PCET reactions, the potential energy surfaces for these states have never been determined directly based on first principles calculations using electronic structure theory. The difficulty is due to a lack of practical method to constrain electron and proton localized diabatic states in wave function or density functional theory calculations. Employing a multistate density functional theory (MSDFT), in which the electron and proton localized diabatic configurations are constructed through block-localization of Kohn-Sham orbitals, we show that distinction between concerted proton-electron transfer (CPET) and HAT, which are not distinguishable experimentally from phenomenological kinetic data, can be made by examining the third dimension of a More O'Ferrall-Jencks diagram that includes both the ground and excited state potential surfaces. In addition, we formulate a pair of effective two-state valence bond models to represent the CPET and HAT mechanisms. We found that the lower energy of the CPET and HAT effective diabatic states at the intersection point can be used as an energetic criterion to distinguish the two mechanisms. In the isoelectronic series of hydrogen exchange reaction in (PhX)(2)H(•), where X = O, NH, and CH(2), there is a continuous transition from a CPET mechanism for the phenoxy radical-phenol pair to a HAT process for benzyl radical and toluene, while the reaction between PhNH(2) and PhNH(•) has a mechanism intermediate of CPET and HAT. The electronically nonadiabatic nature of the CPET mechanism in the phenol system can be attributed to the overlap interactions between the ground and excited state surfaces, resulting in roughly orthogonal minimum energy paths on the adiabatic ground and excited state potential energy surfaces. On the other hand, the minimum energy path on the adiabatic ground state for the HAT mechanism coincides with that on the excited state, producing a large electronic coupling that separates the two surfaces by more than 120 kcal/mol.

Multiple-site concerted proton-electron transfer reactions of hydrogen-bonded phenols are nonadiabatic and well described by semiclassical Marcus theory

Photo-oxidations of hydrogen-bonded phenols using excited-state polyarenes are described to derive fundamental understanding of multiple-site concerted proton-electron transfer reactions (MS-CPET). Experiments have examined phenol bases having -CPh(2)NH(2), -Py, and -CH(2)Py groups ortho to the phenol hydroxyl group and tert-butyl groups in the 4,6-positions for stability (HOAr-NH(2), HOAr-Py, and HOAr-CH(2)Py, respectively; Py = pyridyl; Ph = phenyl). The photo-oxidations proceed by intramolecular proton transfer from the phenol to the pendent base concerted with electron transfer to the excited polyarene. For comparison, 2,4,6-(t)Bu(3)C(6)H(2)OH, a phenol without a pendent base and tert-butyl groups in the 2,4,6-positions, has also been examined. Many of these bimolecular reactions are fast, with rate constants near the diffusion limit. Combining the photochemical k(CPET) values with those from prior thermal stopped-flow kinetic studies gives data sets for the oxidations of HOAr-NH(2) and HOAr-CH(2)Py that span over 10(7) in k(CPET) and nearly 0.9 eV in driving force (ΔG(o)'). Plots of log(k(CPET)) vs ΔG(o)', including both excited-state anthracenes and ground state aminium radical cations, define a single Marcus parabola in each case. These two data sets are thus well described by semiclassical Marcus theory, providing a strong validation of the use of this theory for MS-CPET. The parabolas give λ(CPET) ≅ 1.15-1.2 eV and H(ab) ≅ 20-30 cm(-1). These experiments represent the most direct measurements of H(ab) for MS-CPET reactions to date. Although rate constants are available only up to the diffusion limit, the parabolas clearly peak well below the adiabatic limit of ca. 6 × 10(12) s(-1). Thus, this is a very clear demonstration that the reactions are nonadiabatic. The nonadiabatic character slows the reactions by a factor of ~45. Results for the oxidation of HOAr-Py, in which the phenol and base are conjugated, and for oxidation of 2,4,6-(t)Bu(3)C(6)H(2)OH, which lacks a base, show that both have substantially lower λ and larger pre-exponential terms. The implications of these results for MS-CPET reactions are discussed.

Molecular basis of intramolecular electron transfer in proteins during radical-mediated oxidations: computer simulation studies in model tyrosine-cysteine peptides in solution

Experimental studies in hemeproteins and model Tyr/Cys-containing peptides exposed to oxidizing and nitrating species suggest that intramolecular electron transfer (IET) between tyrosyl radicals (Tyr-O(·)) and Cys residues controls oxidative modification yields. The molecular basis of this IET process is not sufficiently understood with structural atomic detail. Herein, we analyzed using molecular dynamics and quantum mechanics-based computational calculations, mechanistic possibilities for the radical transfer reaction in Tyr/Cys-containing peptides in solution and correlated them with existing experimental data. Our results support that Tyr-O(·) to Cys radical transfer is mediated by an acid/base equilibrium that involves deprotonation of Cys to form the thiolate, followed by a likely rate-limiting transfer process to yield cysteinyl radical and a Tyr phenolate; proton uptake by Tyr completes the reaction. Both, the pKa values of the Tyr phenol and Cys thiol groups and the energetic and kinetics of the reversible IET are revealed as key physico-chemical factors. The proposed mechanism constitutes a case of sequential, acid/base equilibrium-dependent and solvent-mediated, proton-coupled electron transfer and explains the dependency of oxidative yields in Tyr/Cys peptides as a function of the number of alanine spacers. These findings contribute to explain oxidative modifications in proteins that contain sequence and/or spatially close Tyr-Cys residues.

Electron transfer triggered by optical excitation of phenothiazine-tris(meta-phenylene-ethynylene)-(tricarbonyl)(bpy)(py)rhenium(I)

We have investigated excited-state electron transfer in a donor-bridge-acceptor complex containing phenothiazine (PTZ) linked via tris(meta-phenylene-ethynylene) to a tricarbonyl(bipyridine)(pyridine)Re(I) unit. Time-resolved luminescence experiments reveal two excited-state (*Re) decay regimes, a multiexponential component with a mean lifetime of 2.7 ns and a longer monoexponential component of 530 ns in dichloromethane solution. The faster decay is attributed to PTZ → *Re electron transfer in a C-shaped PTZ-bridge-Re conformer (PTZ-Re ≈ 7.5 Å). We assign the longer lifetime, which is virtually identical to that of free *Re, to an extended conformer (PTZ-Re > 20 Å). The observed biexponential *Re decay requires that interconversion of PTZ-bridge-Re conformers be slower than 10(6) s(-1).

Coordination of Bi3+ to metal-free metallothionein: spectroscopy and density functional calculation of structure, coordination, and electronic excitations

Understanding the structure of mammal Bi-containing metallothionein-2 (Bi-MT2) is of great physiological significance due to the importance of Bi-MT2 in alleviating adverse effect of anti-cancer drugs. A unique feature of rabbit liver Bi-MT2 is the metal-oxygen bond (BiO), which is absent in well-characterized Zn-MT2 and Cd-MT2. However, the ligand contributing to the BiO bonding in Bi-MT2 remains unidentified. In this study, the coordination of Bi(3+) to rabbit liver metal-free metallothionein was investigated using both experimental and theoretical methods. UV-visible and circular dichroism spectra indicate that Bi-MT2 has a different secondary structure from those of Zn-MT2 and Cd-MT2. Three possible Bi(3+) coordination structures in Bi(7)-MT2 and relative binding free energies were calculated using the density functional theory. Absorption spectra corresponding to these coordination structures were evaluated by time-dependent density functional theory. Our computation results are consistent with the UV-vis spectroscopic data and strongly suggest that the carboxyl group in the aspartic acid residues contributes to the BiO bond formation.

Probing quantum and dynamic effects in concerted proton-electron transfer reactions of phenol-base compounds

The oxidation of three phenols, which contain an intramolecular hydrogen bond to a pendent pyridine or amine group, has been shown, in a previous experimental study, to undergo concerted proton-electron transfer (CPET). In this reaction, the electron is transferred to an outer-sphere oxidant, and the proton is transferred from the oxygen to nitrogen atom. In the present study, this reaction is studied computationally using a version of Hammes-Schiffer's multistate continuum theory where CPET is formulated as a transmission frequency between neutral and cation vibrational-electronic states. The neutral and cation proton vibrational wave functions are computed from one-dimensional potential energy surfaces (PESs) for the transferring proton in a fixed heavy atom framework. The overlap integrals for these neutral/cation wave functions, considering several initial (i.e., neutral) and final (i.e., cation) vibrational states, are used to evaluate the relative rates of oxidation. The analysis is extended to heavy atom configurations with various proton donor-acceptor (i.e., O-N) distances to assess the importance of heavy atom "gating". Such changes in d(ON) dramatically affect the nature of the proton PESs and wave functions. Surprisingly, the most reactive configurations have similar donor-acceptor distances despite the large (~0.2 Å) differences in the optimized structures. These theoretical results qualitatively reproduce the experimental faster reactivity of the reaction of the pyridyl derivative 1 versus the CH(2)-pyridyl 2, but the computed factor of 5 is smaller than the experimental 10(2). The amine derivative is calculated to react similarly to 1, which does not agree with the experiments, likely due to some of the simplifying assumptions made in applying the theory. The computed kinetic isotope effects (KIEs) and their temperature dependence are in agreement with experimental results.

Dynamic mechanism of proton transfer in mannitol 2-dehydrogenase from Pseudomonas fluorescens: mobile GLU292 controls proton relay through a water channel that connects the active site with bulk solvent

The active site of mannitol 2-dehydrogenase from Pseudomonas fluorescens (PfM2DH) is connected with bulk solvent through a narrow protein channel that shows structural resemblance to proton channels utilized by redox-driven proton pumps. A key element of the PfM2DH channel is the "mobile" Glu(292), which was seen crystallographically to adopt distinct positions up and down the channel. It was suggested that the "down → up" conformational change of Glu(292) could play a proton relay function in enzymatic catalysis, through direct proton shuttling by the Glu or because the channel is opened for water molecules forming a chain along which the protons flow. We report evidence from site-directed mutagenesis (Glu(292) → Ala) substantiated by data from molecular dynamics simulations that support a role for Glu(292) as a gate in a water chain (von Grotthuss-type) mechanism of proton translocation. Occupancy of the up and down position of Glu(292) is influenced by the bonding and charge state of the catalytic acid base Lys(295), suggesting that channel opening/closing motions of the Glu are synchronized to the reaction progress. Removal of gatekeeper control in the E292A mutant resulted in a selective, up to 120-fold slowing down of microscopic steps immediately preceding catalytic oxidation of mannitol, consistent with the notion that formation of the productive enzyme-NAD(+)-mannitol complex is promoted by a corresponding position change of Glu(292), which at physiological pH is associated with obligatory deprotonation of Lys(295) to solvent. These results underscore the important role of conformational dynamics in the proton transfer steps of alcohol dehydrogenase catalysis.