Long Distance Electron-Transfer Mechanism in Peptidylglycine α-Hydroxylating Monooxygenase: A Perfect Fitting for a Water Bridge

Antioxidant Activity of MgSO4 Ion Pairs by Spin-Electron Stabilization of Hydroxyl Radicals through DFT Calculations: Biological Relevance

Magnesium sulfate has been of great interest as an antioxidant for its ability to decrease the oxidizing capacity of the hydroxyl radical. Previously, it was shown that the contact ion pair of this salt could stabilize •OH by coordinating with Mg and delocalizing the unpaired electron over sulfate. The present study explores in detail the MgSO4 antioxidant properties, considering all its ion pairs with •OH in different conformations. The analyses were based on structural, spin, and energetic properties using the DFT approach. As a result, the high antioxidant potential of MgSO4 is related to the spin-electron transfer from SO4 -2 to •OH causing electron spin delocalization and electrostatic stabilization. This transfer occurs for all ion pairs when •OH approaches the Mg first solvation shell, without being coordinated to Mg. The direct Mg-•OH interaction further stabilizes the radical system. These results show that spin-electron transfers are feasible in all hydrated ion pairs MgSO4-•OH, even at a •OH-sulfate distance greater than 10 Å.

High-Power-Density and Excellent-Flexibility Thermoelectric Generator Based on All-SWCNTs/PVP Composites

Flexible polymer/single-wall carbon nanotube (SWCNT) composites are a vital component for wearable/portable electronics, but the development of their n-type counterpart is laggard. Furthermore, little attention is paid to the interaction between SWCNT and polymers, especially the unconjugated polymers, as well as the conversion mechanism of conduction characteristics. Here, the n-type flexible SWCNTs/Polyvinyl Pyrrolidone (PVP) films are successfully fabricated, where the oxygen atoms in PVP interacted with SWCNT via hydrogen bonds, which can lower the energy barrier of electron tunneling, providing the pathway for the electron transfer. Furthermore, with the increasing synthesis temperature, the hydrogen bonds strengthened and the thermal activation energy further improved, both of which enhanced the electron-donating ability of PVP, resulting in a high-power-factor value of 260 µW m-1 K-2. Based on the optimized SWCNTs/PVP films, a thermoelectric module is assembled, which achieved a power density of 400 µW cm-2 at a temperature difference of 56 K, coupled with excellent flexibility, showing a less than 1% variation of resistance after 5000 bending cycles. It shows the highest output-performance and the best flexibility among the reported SWCNT-based thermoelectric modules. This work provides significant insights into the interaction mechanism and performance optimization of hybrid thermoelectric composites, based on SWCNTs/unconjugated polymers.

The critical role of a conserved lysine residue in periplasmic nitrate reductase catalyzed reactions

Periplasmic nitrate reductase NapA from Campylobacter jejuni (C. jejuni) contains a molybdenum cofactor (Moco) and a 4Fe-4S cluster and catalyzes the reduction of nitrate to nitrite. The reducing equivalent required for the catalysis is transferred from NapC → NapB → NapA. The electron transfer from NapB to NapA occurs through the 4Fe-4S cluster in NapA. C. jejuni NapA has a conserved lysine (K79) between the Mo-cofactor and the 4Fe-4S cluster. K79 forms H-bonding interactions with the 4Fe-4S cluster and connects the latter with the Moco via an H-bonding network. Thus, it is conceivable that K79 could play an important role in the intramolecular electron transfer and the catalytic activity of NapA. In the present study, we show that the mutation of K79 to Ala leads to an almost complete loss of activity, suggesting its role in catalytic activity. The inhibition of C. jejuni NapA by cyanide, thiocyanate, and azide has also been investigated. The inhibition studies indicate that cyanide inhibits NapA in a non-competitive manner, while thiocyanate and azide inhibit NapA in an uncompetitive manner. Neither inhibition mechanism involves direct binding of the inhibitor to the Mo-center. These results have been discussed in the context of the loss of catalytic activity of NapA K79A variant and a possible anion binding site in NapA has been proposed.

Detection of Water Molecules on the Radical Transfer Pathway of Ribonucleotide Reductase by 17O Electron-Nuclear Double Resonance Spectroscopy

The role of water in biological proton-coupled electron transfer (PCET) is emerging as a key for understanding mechanistic details at atomic resolution. Here we demonstrate 17O high-frequency electron-nuclear double resonance (ENDOR) in conjunction with H217O-labeled protein buffer to establish the presence of ordered water molecules at three radical intermediates in an active enzyme complex, the α2β2 E. coli ribonucleotide reductase. Our data give unambiguous evidence that all three, individually trapped, intermediates are hyperfine coupled to one water molecule with Tyr-O···17O distances in the range 2.8-3.1 Å. The availability of this structural information will allow for quantitative models of PCET in this prototype enzyme. The results also provide a spectroscopic signature for water H-bonded to a tyrosyl radical.

Theory Demonstrated a "Coupled" Mechanism for O2 Activation and Substrate Hydroxylation by Binuclear Copper Monooxygenases

Multiscale simulations have been performed to address the longstanding issue of "dioxygen activation" by the binuclear copper monooxygenases (PHM and DβM), which have been traditionally classified as "noncoupled" binuclear copper enzymes. Our QM/MM calculations rule out that Cu_M(II)-O₂^• is an active species for H-abstraction from the substrate. In contrast, Cu_M(II)-O₂^• would abstract an H atom from the cosubstrate ascorbate to form a Cu_M(II)-OOH intermediate in PHM and DβM. Consistent with the recently reported structural features of DβM, the umbrella sampling shows that the "open" conformation of the Cu_M(II)-OOH intermediate could readily transform into the "closed" conformation in PHM, in which we located a mixed-valent μ-hydroperoxodicopper(I,II) intermediate, (μ-OOH)Cu(I)Cu(II). The subsequent O-O cleavage and OH moiety migration to Cu_H generate the unexpected species (μ-O^•)(μ-OH)Cu(II)Cu(II), which is revealed to be the reactive intermediate responsible for substrate hydroxylation. We also demonstrate that the flexible Met ligand is favorable for O-O cleavage reactions, while the replacement of Met with the strongly bound His ligand would inhibit the O-O cleavage reactivity. As such, the study not only demonstrates a "coupled" mechanism for O₂ activation by binuclear copper monooxygenases but also deciphers the full catalytic cycle of PHM and DβM in accord with the available experimental data. These findings of O₂ activation and substrate hydroxylation by binuclear copper monooxygenases could expand our understanding of the reactivities of the synthetic monocopper complexes.

Tuning the reactivity of copper complexes supported by tridentate ligands leading to two-electron reduction of dioxygen

Dinuclear copper complex with tridentate ligand and anthracene linkage catalyses 2-electron reduction of O₂.

Quantum chemical calculations of tryptophan → heme electron and excitation energy transfer rates in myoglobin

The development of optical multidimensional spectroscopic techniques has opened up new possibilities for the study of biological processes. Recently, ultrafast two-dimensional ultraviolet spectroscopy experiments have determined the rates of tryptophan → heme electron transfer and excitation energy transfer for the two tryptophan residues in myoglobin (Consani et al., Science, 2013, 339, 1586). Here, we show that accurate prediction of these rates can be achieved using Marcus theory in conjunction with time-dependent density functional theory. Key intermediate residues between the donor and acceptor are identified, and in particular the residues Val68 and Ile75 play a critical role in calculations of the electron coupling matrix elements. Our calculations demonstrate how small changes in structure can have a large effect on the rates, and show that the different rates of electron transfer are dictated by the distance between the heme and tryptophan residues, while for excitation energy transfer the orientation of the tryptophan residues relative to the heme is important. © 2017 The Authors. Journal of Computational Chemistry Published by Wiley Periodicals, Inc.

Mechanism of O2 activation and substrate hydroxylation in noncoupled binuclear copper monooxygenases

Peptidylglycine α-hydroxylating monooxygenase (PHM) and dopamine β-monooxygenase (DβM) are copper-dependent enzymes that are vital for neurotransmitter regulation and hormone biosynthesis. These enzymes feature a unique active site consisting of two spatially separated (by 11 Å in PHM) and magnetically noncoupled copper centers that enables 1e^- activation of O₂ for hydrogen atom abstraction (HAA) of substrate C-H bonds and subsequent hydroxylation. Although the structures of the resting enzymes are known, details of the hydroxylation mechanism and timing of long-range electron transfer (ET) are not clear. This study presents density-functional calculations of the full reaction coordinate, which demonstrate: (i) the importance of the end-on coordination of superoxide to Cu for HAA along the triplet spin surface; (ii) substrate radical rebound to a Cu^II hydroperoxide favors the proximal, nonprotonated oxygen; and (iii) long-range ET can only occur at a late step with a large driving force, which serves to inhibit deleterious Fenton chemistry. The large inner-sphere reorganization energy at the ET site is used as a control mechanism to arrest premature ET and dictate the correct timing of ET.

Electron transfer, decoherence, and protein dynamics: insights from atomistic simulations

Electron transfer in biological systems drives the processes of life. From cellular respiration to photosynthesis and enzymatic catalysis, electron transfers (ET) are chemical processes on which essential biological functions rely. Over the last 40 years, scientists have sought understanding of how these essential processes function in biology. One important breakthrough was the discovery that Marcus theory (MT) of electron transfer is applicable to biological systems. Chemists have experimentally collected both the reorganization energies (λ) and the driving forces (ΔG°), two parameters of Marcus theory, for a large variety of ET processes in proteins. At the same time, theoretical chemists have developed computational approaches that rely on molecular dynamics and quantum chemistry calculations to access numerical estimates of λ and ΔG°. Yet another crucial piece in determining the rate of an electron transfer is the electronic coupling between the initial and final electronic wave functions. This is an important prefactor in the nonadiabatic rate expression, since it reflects the probability that an electron tunnels from the electron donor to the acceptor through the intervening medium. The fact that a protein matrix supports electron tunneling much more efficiently than vacuum is now well documented, both experimentally and theoretically. Meanwhile, many chemists have provided examples of the rich physical chemistry that can be induced by protein dynamics. This Account describes our studies of the dynamical effects on electron tunneling. We present our analysis of two examples of natural biological systems through MD simulations and tunneling pathway analyses. Through these examples, we show that protein dynamics sustain efficient tunneling. Second, we introduce two time scales: τcoh and τFC. The former characterizes how fast the electronic coupling varies with nuclear vibrations (which cause dephasing). The latter reflects the time taken by the system to leave the crossing region. In the framework of open quantum systems, τFC is a short time approximation of the characteristic decoherence time of the electronic subsystem in interaction with its nuclear environment. The comparison of the respective values of τcoh and τFC allows us to probe the occurrence of non-Condon effects. We use ab initio MD simulations to analyze how decoherence appears in several biological cofactors. We conclude that we cannot account for its order of magnitude by considering only the atoms or bonds directly concerned with the transfer. Decoherence results from contributions from all atoms of the system appearing with a time delay that increases with the distance from the primarily concerned atoms or bonds. The delay and magnitude of the contributions depend on the chemical nature of the system. Finally, we present recent developments based on constrained DFT for efficient and accurate evaluations of the electronic coupling in ab initio MD simulations. These are promising methods to study the subtle fluctuations of the electronic coupling and the mechanisms of electronic decoherence in biological systems.

Charge transfer in dynamical biosystems, or the treachery of (static) images

The image is not the thing. Just as a pipe rendered in an oil painting cannot be smoked, quantum mechanical coupling pathways rendered on LCDs do not convey electrons. The aim of this Account is to examine some of our recent discoveries regarding biological electron transfer (ET) and transport mechanisms that emerge when one moves beyond treacherous static views to dynamical frameworks.

Studies over the last two decades introduced both atomistic detail and macromolecule dynamics to the description of biological ET. The first model to move beyond the structureless square-barrier tunneling description is the Pathway model, which predicts how protein secondary motifs and folding-induced through-bond and through-space tunneling gaps influence kinetics. Explicit electronic structure theory is applied routinely now to elucidate ET mechanisms, to capture pathway interferences, and to treat redox cofactor electronic structure effects. Importantly, structural sampling of proteins provides an understanding of how dynamics may change the mechanisms of biological ET, as ET rates are exponentially sensitive to structure. Does protein motion average out tunneling pathways? Do conformational fluctuations gate biological ET? Are transient multistate resonances produced by energy gap fluctuations? These questions are becoming accessible as the static view of biological ET recedes and dynamical viewpoints take center stage.

This Account introduces ET reactions at the core of bioenergetics, summarizes our team’s progress toward arriving at an atomistic-level description, examines how thermal fluctuations influence ET, presents metrics that characterize dynamical effects on ET, and discusses applications in very long (micrometer scale) bacterial nanowires. The persistence of structural effects on the ET rates in the face of thermal fluctuations is considered. Finally, the flickering resonance (FR) view of charge transfer is presented to examine how fluctuations control low-barrier transport among multiple groups in van der Waals contact. FR produces exponential distance dependence in the absence of tunneling; the exponential character emerges from the probability of matching multiple vibronically broadened electronic energies within a tolerance defined by the rms coupling among interacting groups. FR thus produces band like coherent transport on the nanometer length scale, enabled by conformational fluctuations. Taken as a whole, the emerging context for ET in dynamical biomolecules provides a robust framework to design and interpret the inner workings of bioenergetics from the molecular to the cellular scale and beyond, with applications in biomedicine, biocatalysis, and energy science.

Reaction mechanism of the bicopper enzyme peptidylglycine α-hydroxylating monooxygenase

Peptidylglycine α-hydroxylating monooxygenase is a noninteracting bicopper enzyme that stereospecifically hydroxylates the terminal glycine of small peptides for its later amidation. Neuroendocrine messengers, such as oxytocin, rely on the biological activity of this enzyme. Each catalytic turnover requires one oxygen molecule, two protons from the solvent, and two electrons. Despite this enzyme having been widely studied, a consensus on the reaction mechanism has not yet been found. Experiments and theoretical studies favor a pro-S abstraction of a hydrogen atom followed by the rebinding of an OH group. However, several hydrogen-abstracting species have been postulated; because two protons are consumed during the reaction, several protonation states are available. An electron transfer between the copper atoms could play a crucial role for the catalysis as well. This leads to six possible abstracting species. In this study, we compare them on equal footing. We perform quantum mechanics/molecular mechanics calculations, considering the glycine hydrogen abstraction. Our results suggest that the most likely mechanism is a protonation of the abstracting species before the hydrogen abstraction and another protonation as well as a reduction before OH rebinding.

The key role of water in the dioxygenase function of Escherichia coli flavohemoglobin

Flavohemoglobins (FHbs) are members of the globin superfamily, widely distributed among prokaryotes and eukaryotes that have been shown to carry out nitric oxide dioxygenase (NOD) activity. In prokaryotes, such as Escherichia coli, NOD activity is a defence mechanism against the NO release by the macrophages of the hosts' immune system during infection. Because of that, FHbs have been studied thoroughly and several drugs have been developed in an effort to fight infectious processes. Nevertheless, the protein's structural determinants involved in the NOD activity are still poorly understood. In this context, the aim of the present work is to unravel the molecular basis of FHbs structural dynamics-to-function relationship using state of the art computer simulation tools. In an effort to fulfill this goal, we studied three key processes that determine NOD activity, namely i) ligand migration into the active site ii) stabilization of the coordinated oxygen and iii) intra-protein electron transfer (ET). Our results allowed us to determine key factors related to all three processes like the presence of a long hydrophobic tunnel for ligand migration, the presence of a water mediated hydrogen bond to stabilize the coordinated oxygen and therefore achieve a high affinity, and the best possible ET paths between the FAD and the heme, where water molecules play an important role. Taken together the presented results close an important gap in our understanding of the wide and diverse globin structural-functional relationships.

New insights into the mechanism of electron transfer within flavohemoglobins: tunnelling pathways, packing density, thermodynamic and kinetic analyses

Flavohemoglobins (FlavoHb) are metalloenzymes catalyzing the reaction of nitric oxide dioxygenation. The iron cation of the heme group needs to be preliminarily reduced to the ferrous state to be catalytically competent. This reduction is triggered by a flavin adenine dinucleotide (FAD) prosthetic group which is localized in a distinct domain of the protein. In this paper we obtain new insights into the internal long range electron transfer (over ca. 12 Å) using a combination of experimental and computational approaches. Employing a time-resolved pulse radiolysis technique we report the first direct measurement of the FADH˙→ HemeFe(III) electron transfer rate. A rate constant of (6.8 ± 0.5) × 10(3) s(-1) is found. A large panel of computational approaches are used to provide the first estimation of the thermodynamic characteristics of the internal electron transfer step within flavoHb: both the driving force and the reorganization energy are estimated as a function of the protonated state of the flavin semi-quinone. We also report an analysis of the electron pathways involved in the tunnelling of the electron through the aqueous interface between the globin and the flavin domains.

Exploring biological electron transfer pathway dynamics with the Pathways plugin for VMD

We describe the new Pathways plugin for the molecular visualization program visual molecular dynamics. The plugin identifies and visualizes tunneling pathways and pathway families in biomolecules, and calculates relative electronic couplings. The plugin includes unique features to estimate the importance of individual atoms for mediating the coupling, to analyze the coupling sensitivity to thermal motion, and to visualize pathway fluctuations. The Pathways plugin is open source software distributed under the terms of the GNU's Not Unix (GNU) public license.

Probing the peptidylglycine alpha-hydroxylating monooxygenase active site with novel 4-phenyl-3-butenoic acid based inhibitors

Specific inhibition of the copper-containing peptidylglycine alpha-hydroxylating monooxygenase (PHM), which catalyzes the post-translational modification of peptides involved in carcinogenesis and tumor progression, constitutes a new approach for combating cancer. We carried out a structure-activity study of new compounds derived from a well-known PHM substrate analogue, the olefinic compound 4-phenyl-3-butenoic acid (PBA). We designed, synthesized, and tested various PBA derivatives both in vitro and in silico. We show that it is possible to increase PBA affinity for PHM by appropriate functionalization of its aromatic nucleus. Compound 2 d, for example, bears a meta-benzyloxy substituent, and exhibits better inhibition features (K(i)=3.9 microM, k(inact)/K(i)=427 M(-1) s(-1)) than the parent PBA (K(i)=19 microM, k(inact)/K(i)=82 M(-1) s(-1)). Docking calculations also suggest two different binding modes for PBA derivatives; these results will aid in the development of further PHM inhibitors with improved features.

Molecular dynamics free energy calculations to assess the possibility of water existence in protein nonpolar cavities

Are protein nonpolar cavities filled with water molecules? Although many experimental and theoretical investigations have been done, particularly for the nonpolar cavity of IL-1 beta, the results are still conflicting. To study this problem from the thermodynamic point of view, we calculated hydration free energies of four protein nonpolar cavities by means of the molecular dynamics thermodynamic integration method. In addition to the IL-1 beta cavity (69 A(3)), we selected the three largest nonpolar cavities of AvrPphB (81 A(3)), Trp repressor (87 A(3)), and hemoglobin (108 A(3)) from the structural database, in view of the simulation result from another study that showed larger nonpolar cavities are more likely to be hydrated. The calculations were performed with flexible and rigid protein models. The calculated free energy changes were all positive; hydration of the nonpolar cavities was energetically unfavorable for all four cases. Because hydration of smaller cavities should happen more rarely, we conclude that existing protein nonpolar cavities are not likely to be hydrated. Although a possibility remains for much larger nonpolar cavities, such cases are not found experimentally. We present a hypothesis to explain this: hydrated nonpolar cavities are quite unstable and the conformation could not be maintained.

Surface residues dynamically organize water bridges to enhance electron transfer between proteins

Cellular energy production depends on electron transfer (ET) between proteins. In this theoretical study, we investigate the impact of structural and conformational variations on the electronic coupling between the redox proteins methylamine dehydrogenase and amicyanin from Paracoccus denitrificans. We used molecular dynamics simulations to generate configurations over a duration of 40 ns (sampled at 100-fs intervals) in conjunction with an ET pathway analysis to estimate the ET coupling strength of each configuration. In the wild-type complex, we find that the most frequently occurring molecular configurations afford superior electronic coupling due to the consistent presence of a water molecule hydrogen-bonded between the donor and acceptor sites. We attribute the persistence of this water bridge to a "molecular breakwater" composed of several hydrophobic residues surrounding the acceptor site. The breakwater supports the function of nearby solvent-organizing residues by limiting the exchange of water molecules between the sterically constrained ET region and the more turbulent surrounding bulk. When the breakwater is affected by a mutation, bulk solvent molecules disrupt the water bridge, resulting in reduced electronic coupling that is consistent with recent experimental findings. Our analysis suggests that, in addition to enabling the association and docking of the proteins, surface residues stabilize and control interprotein solvent dynamics in a concerted way.

Understanding selectivity of hard and soft metal cations within biological systems using the subvalence concept. I. Application to blood coagulation: direct cation-protein electronic effects vs. indirect interactions through water networks

Following a previous study by de Courcy et al. ((2009) Interdiscip. Sci. Comput. Life Sci. 1, 55-60), we demonstrate in this contribution, using quantum chemistry, that metal cations exhibit a specific topological signature in the electron localization of their density interacting with ligands according to its "soft" or "hard" character. Introducing the concept of metal cation subvalence, we show that a metal cation can split its outer-shell density (the so-called subvalent domains or basins) according to it capability to form a partly covalent bond involving charge transfer. Such behaviour is investigated by means of several quantum chemical interpretative methods encompasing the topological analysis of the Electron Localization Function (ELF) and Bader's Quantum Theory of Atoms in Molecules (QTAIM) and two energy decomposition analyses (EDA), namely the Restricted Variational Space (RVS) and Constrained Space Orbital Variations (CSOV) approaches. Further rationalization is performed by computing ELF and QTAIM local properties such as electrostatic distributed moments and local chemical descriptors such as condensed Fukui Functions and dual descriptors. These reactivity indexes are computed within the ELF topological analysis in addition to QTAIM offering access to non atomic reactivity local index, for example on lone pairs. We apply this "subvalence" concept to study the cation selectivity in enzymes involved in blood coagulation (GLA domains of three coagulation factors). We show that the calcium ions are clearly able to form partially covalent charge transfer networks between the subdomain of the metal ion and the carboxylate oxygen lone pairs whereas magnesium does not have such ability. Our analysis also explains the different role of two groups (high affinity and low affinity cation binding sites) present in GLA domains. If the presence of Ca(II) is mandatory in the central "high affinity" region to conserve a proper folding and a charge transfer network, external sites are better stabilised by Mg(II), rather than Ca(II), in agreement with experiment. The central role of discrete water molecules is also discussed in order to understand the stabilities of the observed X-rays structures of the Gla domain. Indeed, the presence of explicit water molecules generating indirect cation-protein interactions through water networks is shown to be able to reverse the observed electronic selectivity occuring when cations directly interact with the Gla domain without the need of water.