Spectroscopic characterization and assignment of reduction potentials in the tetraheme cytochrome C554 from Nitrosomonas europaea

The tetraheme cytochrome c(554) (cyt c(554)) from Nitrosomonas europaea is an essential electron transfer component in the biological oxidation of ammonia. The protein contains one 5-coordinate heme and three bis-His coordinated hemes in a 3D arrangement common to a newly characterized class of multiheme proteins. The ligand binding, electrochemical properties, and heme-heme interactions are investigated with M?ssbauer and X- and Q-band (parallel/perpendicular mode) EPR spectroscopy. The results indicate that the 5-coordinate heme will not bind the common heme ligands, CN(-), F(-), CO, and NO in a wide pH range. Thus, cyt c(554) functions only in electron transfer. Analysis of a series of electrochemically poised and chemically reduced samples allows assignment of reduction potentials for heme 1 through 4 of +47, +47, -147, and -276 mV, respectively. The spectroscopic results indicate that the hemes are weakly exchange-coupled (J approximately -0.5 cm(-)(1)) in two separate pairs and in accordance with the structure: hemes 2/4 (high-spin/low-spin), hemes 1/3 (low-spin/low-spin). There is no evidence of exchange coupling between the pairs. A comparison of the reduction potentials between homologous hemes of cyt c(554) and other members of this new class of multiheme proteins is discussed. Heme 1 has a unique axial N(delta)-His coordination which contributes to a higher potential relative to the homologous hemes of hydroxylamine oxidoreductase (HAO) and the split-Soret cytochrome. Heme 2 is 300 mV more positive than heme 4 of HAO, which is attributed to hydroxide coordination to heme 4 of HAO.     

Characterization of electronic structure and properties of a Bis(histidine) heme model complex

Ferric and ferrous hemes, such as those present in electron transfer proteins, often have low-lying spin states that are very close in energy. To explore the relationship between spin state, geometry, and cytochrome electron transfer, we investigate, using density functional theory, the relative energies, electronic structure, and optimized geometries for a high- and low-spin ferric and ferrous heme model complex. Our model consists of an iron-porphyrin axially ligated by two imidazoles, which model the interaction of a heme with histidine residues. Using the B3LYP hybrid functional, we found that, in the ferric model heme complex, the doublet is lower in energy than the sextet by 8.4 kcal/mol and the singlet ferrous heme is 6.7 kcal/mol more stable than the quintet. The difference between the high-spin ferric and ferrous model heme energies yields an adiabatic electron affinity (AEA) of 5.24 eV, and the low-spin AEA is 5.17 eV. Both values are large enough to ensure electron trapping, and electronic structure analysis indicates that the iron d(pi) orbital is involved in the electron transfer between hemes. M?ssbauer parameters calculated to verify the B3LYP electronic structure correlate very well with experimental values. Isotropic hyperfine coupling constants for the ligand nitrogen atoms were also evaluated. The optimized geometries of the ferric and ferrous hemes are consistent with structures from X-ray crystallography and reveal that the iron-imidazole distances are significantly longer in the high-spin hemes, which suggests that the protein environment, modeled here by the imidazoles, plays an important role in regulating the spin state. Iron-imidazole dissociation energies, force constants, and harmonic frequencies were calculated for the ferric and ferrous low-spin and high-spin hemes. In both the ferric and the ferrous cases, a single imidazole ligand is more easily dissociated from the high-spin hemes.     

Redox dependent conformational changes in the mixed valence form of the cytochrome c oxidase from p. The reorganization of glutamic acid 278 is coupled to the electron transfer from/to heme a and the binuclear center. denitrificans

In this work we present the separation of FTIR difference signals induced by electron transfer to/from the redox centers of the cytochrome c oxidase from P. denitrificans and compare electrochemically induced FTIR difference spectra with those induced by CO photolysis. FTIR difference spectra of rebinding of CO to the half reduced (mixed valence) form of the cytochrome c oxidase after photolysis reflect the conformational changes induced by the rebinding of CO and by electron transfer reactions from heme a3 to heme a and further on to CUA. During this process, heme a3 (and CUB) are oxidized, whereas heme a and CuA are reduced. By subtracting these difference spectra from an electrochemically induced FTIR difference spectrum, where all four cofactors are reduced, the contributions for heme a3 (and CuB) could be separated. Correspondingly, the spectral contributions of heme a and CuA have been separated. The comparison of these spectra with the spectra calculated for the hemes on the basis of their redox dependent changes previously published in Hellwig et al., (Biochemistry 38, (1999) 1685-1694) show a high degree of similarity, except for additional signals coupled to the reorganization of the binuclear center upon CO rebinding. The separated spectra clearly show that the signals attributed to Glu278, an amino acid discussed to be crucial for proton pumping, is coupled to electron transfer to/from heme a and the binuclear heme a3-CuB center.     

Dynamic docking of cytochrome b5 with myoglobin and alpha-hemoglobin: heme-neutralization "squares" and the binding of electron-transfer-reactive configurations

Intracomplex electron transfer (ET) occurs most often in intrinsically transient, low affinity complexes. As a result, the means by which adequate specificity and reactivity are obtained to support effective ET is still poorly understood. We report here on two such ET complexes: cytochrome b5 (cyt b5) in reaction with its physiological partners, myoglobin (Mb) and hemoglobin (Hb). These complexes obey the Dynamic Docking (DD) paradigm: a large ensemble of weakly bound protein-protein configurations contribute to binding in the rapid-exchange limit, but only a few are ET-active. We report the ionic-strength dependence of the second-order rate constant, k2, for photoinitiated ET from within all four combinations of heme-neutralized Zn deuteroporphyrin-substituted Mb/alphaHb undergoing ET with cyt b5, the four "corners" of a "heme-neutralization square". These experiments provide insights into the relative importance of both global and local electrostatic contributions to the binding of reactive configurations, which are too few to be observed directly. To interpret the variations of k2 arising from heme neutralization, we have developed a procedure by which comparisons of the ET rate constants for a heme-neutralization square permit us to decompose the free energy of reactive binding into individual local electrostatic contributions associated with interactions between (i) the propionates of the two hemes and (ii) the heme of each protein with the polypeptide of its partner. Most notably, we find the contribution from the repulsion between propionates of partner hemes to the reactive binding free energy to be surprisingly small, DeltaG(Hb) approximately +1 kcal/mol at ambient temperature, 18 mM ionic strength, and we speculate about possible causes of this observation. To confirm the fundamental assumption of these studies, that the structure of a heme-neutralized protein is unaltered either by substitution of Zn or by heme neutralization, we have obtained the X-ray structure of ZnMb prepared with the porphyrin dimethyl ester and find it to be nearly isostructural with the native protein.     

Spatially Resolved Confocal Resonant Raman Microscopic Analysis of Anode‐Grown Geobacter sulfurreducens Biofilms

When grown on the surface of an anode electrode, Geobacter sulfurreducens forms a multi‐cell thick biofilm in which all cells appear to couple the oxidation of acetate with electron transport to the anode, which serves as the terminal metabolic electron acceptor. Just how electrons are transported through such a biofilm from cells to the underlying anode surface over distances that can exceed 20 microns remains unresolved. Current evidence suggests it may occur by electron hopping through a proposed network of redox cofactors composed of immobile outer membrane and/or extracellular multi‐heme c‐type cytochromes. In the present work, we perform a spatially resolved confocal resonant Raman (CRR) microscopic analysis to investigate anode‐grown Geobacter biofilms. The results confirm the presence of an intra‐biofilm redox gradient whereby the probability that a heme is in the reduced state increases with increasing distance from the anode surface. Such a gradient is required to drive electron transport toward the anode surface by electron hopping via cytochromes. The results also indicate that at open circuit, when electrons are expected to accumulate in redox cofactors involved in electron transport due to the inability of the anode to accept electrons, nearly all c‐type cytochrome hemes detected in the biofilm are oxidized. The same outcome occurs when a comparable potential to that measured at open circuit (?0.30 V vs. SHE) is applied to the anode, whereas nearly all hemes are reduced when an exceedingly negative potential (?0.50 V vs. SHE) is applied to the anode. These results suggest that nearly all c‐type cytochrome hemes detected in the biofilm can be electrochemically accessed by the electrode, but most have oxidation potentials too negative to transport electrons originating from acetate metabolism. The results also reveal a lateral heterogeneity (x–y dimensions) in the type of c‐type cytochromes within the biofilm that may affect electron transport to the electrode.     

Correlations of structure and electronic properties from EPR spectroscopy of hydroxylamine oxidoreductase.

Hydroxylamine oxidoreductase (HAO) from the autotrophic nitrifying bacterium Nitrosomonas europaea catalyzes the oxidation of NH(2)OH to HNO(2). The enzyme contains eight hemes per subunit which participate in catalytic function and electron transport. The structure of the enzyme shows a unique spatial arrangement of the eight hemes, subsets of which are now observed in four other proteins. The spatial arrangement displays three types of diheme pairing motifs. At least four of the eight hemes are electronically coupled in two distinguishable pairs and one of these pairs is at the active site of the enzyme. Here, the use of quantitative simulation of the EPR signals allows determination of exchange couplings, and assignments of signals and reduction potentials to hemes of the crystal structure. The absence of any obvious heme-to-heme bonding pathway in the crystal structure suggests that the observed exchange interactions are derived from direct electronic overlap of porphyrin orbitals. This provides evidence for heme pairs which function as biological two-electron redox centers in electron-transfer processes.     

Direct Electrochemistry with Nitrate Reductase in Chitosan Films

Stable films made from chitosan (CS) on pyrolytic graphite (PG) electrode gave direct electrochemistry for incorporated enzyme nitrate reductase (NR). Cyclic voltammetry (CV) of NR‐CS films showed a pair of well‐defined and nearly reversible redox peaks at about ?0.430 V vs. SCE at pH 7.0 phosphate buffers, which was considered as the redox of FAD and heme‐ion. The electron transfer between NR and PG electrode was greatly facilitated in CS films. Integration of reduction peaks at different scan rates from 0.01 to 0.5 V s^?;1 gave nearly constant charge values, which is characteristic of thin‐larger‐electrochemical behavior. The pH of the solution strongly affected the direct electron transfer of NR‐CS films. EDTA accelerated the electron transfer, and it was proposed as a stimulant for the system. Reflectance absorption infrared spectra demonstrated that NR retained a nearly native conformation in CS films.     

Does the electron spin affect the rates of electron tunneling in electrochemical systems?

An effect of spin degeneracy of electron energy levels in the metal electrode on the observable characteristics of electrochemical systems in the absence of magnetic field is discussed. Single-electrode outer-sphere electron transfer reactions are considered as well as redox-mediated electron tunneling in electrochemical contacts. Particular attention is paid to the difference between the spin-less model and the limit of infinitely large Coulomb repulsion of the electrons occupying the same valence orbital in the redox group. Adiabatic and non-adiabatic regimes of the transitions are studied and the expressions for the tunnel current are obtained.     

Protein dynamics control of electron transfer in photosynthetic reaction centers from Rps. sulfoviridis

In the cycle of photosynthetic reaction centers, the initially oxidized special pair of bacteriochlorophyll molecules is subsequently reduced by an electron transferred over a chain of four hemes of the complex. Here, we examine the kinetics of electron transfer between the proximal heme c-559 of the chain and the oxidized special pair in the reaction center from Rps. sulfoviridis in the range of temperatures from 294 to 40 K. The experimental data were obtained for three redox states of the reaction center, in which one, two, or three nearest hemes of the chain are reduced prior to special pair oxidation. The experimental kinetic data are analyzed in terms of a Sumi-Marcus-type model developed in our previous paper,1 in which similar measurements were reported on the reaction centers from Rps. viridis. The model allows us to establish a connection between the observed nonexponential electron-transfer kinetics and the local structural relaxation dynamics of the reaction center protein on the microsecond time scale. The activation energy for relaxation dynamics of the protein medium has been found to be around 0.1 eV for all three redox states, which is in contrast to a value around 0.4-0.6 eV in Rps. viridis.1 The possible nature of the difference between the reaction centers from Rps. viridis and Rps. sulfoviridis, which are believed to be very similar, is discussed. The role of the protein glass transition at low temperatures and that of internal water molecules in the process are analyzed.     

Electronic structures of heme a of cytochrome c oxidase in the redox states—Charge density migration to the propionate groups of heme a

The electronic structures of heme a of cytochrome c oxidase in the redox states were studied, using hybrid density functional theory with a polarizable continuum model and a point charge model. We found that the most stable electronic configurations of the d electrons of the Fe ion are determined by the orbital interactions of the d orbitals of the Fe ion with the π orbitals of the porphyrin ring and the His residues. The redox reaction of the Fe ion influences the charge density on the formyl group through the π conjugation of the porphyrin ring. In addition, we found the charge transfer from the Fe ion to the propionate group of heme a in the redox change despite the lack of the π‐conjugation. We elucidated that the charge propagation originates from the heme a structure itself and that the origin of the charge delocalization to the heme propionate is the orbital interactions between the d orbital of the Fe ion and the p orbitals of the carboxylate part of the heme propionate via the π conjugation of the porphyrin ring and the σ* orbital of the C? C bond of the propionate group. The electrostatic effect by surrounding proteins enhances the charge transfer from the Fe ion to the propionate group. These results indicate that heme propionate groups serve electron mediators in electron transfer as well as electrostatic anchors, and that proteins surrounding the active site reinforce the congenital abilities of the cofactors. © 2009 Wiley Periodicals, Inc. J Comput Chem 2010     

Electrochemical transistor based on bridge tunneling contact containing two redox groups

The method of kinetic equations was used to show that a bridge tunneling contact containing two redox groups in the sequential configuration immersed into a solution of electrolyte at the room temperature features pronounced transistor properties at a given set of physical system parameters. Valence electrons of redox groups interact strongly with the classical phonon subsystem of the liquid medium. Debye screening of the electric field in the tunneling gap and Coulomb repulsion between electrons in different redox groups are taken into account. The case of nonadiabatic electron transfer both between redox groups and between electrodes and redox groups is considered in the limit of infinitely high Coulomb repulsion between electrons in a redox group. For sufficiently high absolute values of difference δ between unperturbed energy levels of redox groups, the system features voltammetric characteristics typical for a transistor. The amplification effect appears due to a strong dependence of tunneling current on overpotential. The emphasis is upon the peculiaritiespeculiarities of voltammetric characteristics in the case of asymmetric tunneling contacts.     

Surface‐Tuned Electron Transfer and Electrocatalysis of Hexameric Tyrosine‐Coordinated Heme Protein

Molecular modeling, electrochemical methods, and quartz crystal microbalance were used to characterize immobilized hexameric tyrosine‐coordinated heme protein (HTHP) on bare carbon or on gold electrodes modified with positively and negatively charged self‐assembled monolayers (SAMs), respectively. HTHP binds to the positively charged surface but no direct electron transfer (DET) is found due to the long distance of the active sites from the electrode surfaces. At carboxyl‐terminated surfaces, the neutrally charged bottom of HTHP can bind to the SAM. For this “disc” orientation all six hemes are close to the electrode and their direct electron transfer should be efficient. HTHP on all negatively charged SAMs showed a quasi‐reversible redox behavior with rate constant k_s values between 0.93 and 2.86 s^?1 and apparent formal potentials ${E{{0{^{\prime }}\hfill \atop {\rm app}\hfill}}}$ between ‐131.1 and ‐249.1 mV. On the MUA/MU‐modified electrode, the maximum surface concentration corresponds to a complete monolayer of the hexameric HTHP in the disc orientation. HTHP electrostatically immobilized on negatively charged SAMs shows electrocatalysis of peroxide reduction and enzymatic oxidation of NADH.     

Structure–Function of the Cytochrome b6f Complex†

The structure and function of the cytochrome b₆ f complex is considered in the context of recent crystal structures of the complex as an eight subunit, 220 kDa symmetric dimeric complex obtained from the thermophilic cyanobacterium, Mastigocladus laminosus, and the green alga, Chlamydomonas reinhardtii. A major problem confronted in crystallization of the cyanobacterial complex, proteolysis of three of the subunits, is discussed along with initial efforts to identify the protease. The evolution of these cytochrome complexes is illustrated by conservation of the hydrophobic heme‐binding transmembrane domain of the cyt b polypeptide between b₆ f and bc₁ complexes, and the rubredoxin‐like membrane proximal domain of the Rieske [2Fe‐2S] protein. Pathways of coupled electron and proton transfer are discussed in the framework of a modified Q cycle, in which the heme c_n, not found in the bc₁ complex, but electronically tightly coupled to the heme b_n of the b₆ f complex, is included. Crystal structures of the cyanobacterial complex with the quinone analogue inhibitors, NQNO or tridecyl‐stigmatellin, show the latter to be ligands of heme c_n, implicating heme c_n as an n‐side plastoquinone reductase. Existing questions include (a) the details of the shuttle of: (i) the [2Fe‐2S] protein between the membrane‐bound PQH₂ electron/H⁺ donor and the cytochrome f acceptor to complete the p‐side electron transfer circuit; (ii) PQ/PQH₂ between n‐ and p‐sides of the complex across the intermonomer quinone exchange cavity, through the narrow portal connecting the cavity with the p‐side [2Fe‐2S] niche; (b) the role of the n‐side of the b₆ f complex and heme c_n in regulation of the relative rates of noncyclic and cyclic electron transfer. The likely presence of cyclic electron transport in the b₆ f complex, and of heme c_n in the firmicute bc complex suggests the concept that hemes b_n‐c_n define a branch point in bc complexes that can support electron transport pathways that differ in detail from the Q cycle supported by the bc₁ complex.     

Redox Multifunctionality in a Series of PtII Dithiolene Complexes of a Tetrathiafulvalene‐Based Diphosphine Ligand

The redox‐active and chelating diphosphine, 3,4‐dimethyl‐3′,4′‐bis(diphenylphosphino)‐tetrathiafulvalene, denoted as P2 , is engaged in a series of platinum complexes, [(P2)Pt(dithiolene)], with different dithiolate ligands, such as 1,2‐benzenedithiolate (bdt), 1,3‐dithiole‐2‐thione‐4,5‐dithiolate (dmit), and 5,6‐dihydro‐1,4‐dithiin‐2,3‐dithiolate (dddt). The complexes are structurally characterized by X‐ray diffraction, together with a model compound derived from bis(diphenylphosphino)ethane, namely, [(dppe)Pt(dddt)] . Four successive reversible electron‐transfer processes are found for the [(P2)Pt(dddt)] complex, associated with the two covalently linked but electronically uncoupled electrophores, that is, the TTF core and the platinum dithiolene moiety. The assignments of the different redox processes to either one or the other electrophore is made thanks to the electrochemical properties of the model compound [(dppe)Pt(dddt)] lacking the TTF redox core, and with the help of theoretical calculations (DFT) to understand the nature and energy of the frontier orbitals of the [(P2)Pt(dithiolene)] complexes in their different oxidation states. The first oxidation of the highly electron‐rich [(P2)Pt(dddt)] complex can be unambiguously assigned to the redox process affecting the Pt(dddt) moiety rather than the TTF core, a rare example in the coordination chemistry of tetrathiafulvalenes acting as ligands.     

Artificial cytochrome b: computer modeling and evaluation of redox potentials.

We generated atomic coordinates of an artificial protein that was recently synthesized to model the central part of the native cytochrome b (Cb) subunit consisting of a four-helix bundle with two hemes. Since no X-ray structure is available, the structural elements of the artificial Cb were assembled from scratch using all known chemical and structural information available and avoiding strain as much as possible. Molecular dynamics (MD) simulations applied to this model protein exhibited root-mean-square deviations as small as those obtained from MD simulations starting with the crystal structure of the native Cb subunit. This demonstrates that the modeled structure of the artificial Cb is relatively rigid and strain-free. The model structure of the artificial Cb was used to determine the redox potentials of the two hemes by calculating the electrostatic energies from the solution of the linearized Poisson-Boltzmann equation (LPBE). The calculated redox potentials agree within 20 meV with the experimentally measured values. The dependence of the redox potentials of the hemes on the protein environment was analyzed. Accordingly, the total shift in the redox potentials is mainly due to the low dielectric medium of the protein, the protein backbone charges, and the salt bridges formed between the arginines and the propionic acid groups of the hemes. The difference in the shift of the redox potentials is due to the interactions with the hydrophilic side chains and the salt bridges formed with the propionic acids of the hemes. For comparison and to test the computational procedure, the redox potentials of the two hemes in the native Cb from the cytochrome bc(1) (Cbc(1)) complex were also calculated. Also in this case the computed redox potentials agree well with experiments.     

Control of the Electronic Ground State on an Electron‐Transfer Copper Site by Second‐Sphere Perturbations

The Cu_A center is a dinuclear copper site that serves as an optimized hub for long‐range electron transfer in heme–copper terminal oxidases. Its electronic structure can be described in terms of a σ_u* ground‐state wavefunction with an alternative, less populated ground state of π_u symmetry, which is thermally accessible. It is now shown that second‐sphere mutations in the Cu_A containing subunit of Thermus thermophilus ba₃ oxidase perturb the electronic structure, which leads to a substantial increase in the population of the π_u state, as shown by different spectroscopic methods. This perturbation does not affect the redox potential of the metal site, and despite an increase in the reorganization energy, it is not detrimental to the electron‐transfer kinetics. The mutations were achieved by replacing the loops that are involved in protein–protein interactions with cytochrome c, suggesting that transient protein binding could also elicit ground‐state switching in the oxidase, which enables alternative electron‐transfer pathways.     

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.     

A Series of Diamagnetic Pyridine Monoimine Rhenium Complexes with Different Degrees of Metal‐to‐Ligand Charge Transfer: Correlating 13C NMR Chemical Shifts with Bond Lengths in Redox‐Active Ligands

A set of pyridine monoimine (PMI) rhenium(I) tricarbonyl chlorido complexes with substituents of different steric and electronic properties was synthesized and fully characterized. Spectroscopic (NMR and IR) and single‐crystal X‐ray diffraction analyses of these complexes showed that the redox‐active PMI ligands are neutral and that the overall electronic structure is little affected by the choices of the substituent at the ligand backbone. One‐ and two‐electron reduction products were prepared from selected starting compounds and could also be characterized by multiple spectroscopic methods and X‐ray diffraction. The final product of a one‐electron reduction in THF is a diamagnetic metal–metal‐bonded dimer after loss of the chlorido ligand. Bond lengths in and NMR chemical shifts of the PMI ligand backbone indicate partial electron transfer to the ligand. Two‐electron reduction in THF also leads to the loss of the chlorido ligand and a pentacoordinate complex is obtained. The comparison with reported bond lengths and ¹³C NMR chemical shifts of doubly reduced free pyridine monoaldimine ligands indicates that both redox equivalents in the doubly reduced rhenium complex investigated here are located in the PMI ligand. With diamagnetic complexes varying over three formal reduction stages at the PMI ligand we were, for the first time, able to establish correlations of the ¹³C NMR chemical shifts with the relevant bond lengths in redox‐active ligands over a full redox series.     

Peculiarity of Hexamethylenetetratellurafulvalene (HMTTeF) Charge Transfer Complexes of Donor-Acceptor (D-A) Type

The complex formation of hexamethylenetetratellurafulvalene (HMTTeF) with 28 kinds of organic electron acceptors yielded 31 charge transfer (CT) complexes. The infrared and ultra-violet-visible-near-infrared spectra of the complexes were examined to study the ionicity of their ground states in solid. A plot of CT transition energies and the difference of redox potentials; ΔE(DA) of donor (D) and acceptor (A) molecules indicated that four complexes have a neutral ground state. Four other complexes exhibit characteristic features of a fully ionic ground state based on the vibrational spectra. Notably, the HCBD, F₄TCNQ and DDQ complexes indicate both a relatively low first CT band and high conductivity in a solid in spite of the fully ionic character being very plausible. Twenty-three complexes having a partially ionic ground state have a CT band below 4×10 cm⁻¹ and are highly conductive. The preparation of good single crystals of the HMTTeF complexes for structural analysis was only successful with Et₂TCNQ and BTDA-TCNQ, which have the structure of DA alternately stacking. These two complexes indicate high conductivities in spite of their disadvantageous packing manner. The intermolecular interactions are found to be strongly enhanced by both the bulky molecular orbital of HMTTeF and the decreased on-site Coulomb repulsion in the HMTTeF complexes. These two factors in particular seem to prevent both the fully ionic and the DA alternating HMTTeF complexes from becoming insulators, even though the redox parameters and the crystal structures predict them to be insulating.     

Thermodynamics of electron flow in the bacterial deca-heme cytochrome MtrF

Electron-transporting multi-heme cytochromes are essential to the metabolism of microbes that inhabit soils and carry out important biogeochemical processes. Recently the first crystal structure of a prototype bacterial deca-heme cytochrome (MtrF) has been resolved and its electrochemistry characterized. However, the molecular details of electron transport along heme chains in the cytochrome are difficult to access via experiment due to the nearly identical chemical nature of the heme cofactors. Here we employ large-scale molecular dynamics simulations to compute the redox potentials of the 10 hemes of MtrF in aqueous solution. We find that as a whole they fall within a range of ~0.3 V, in agreement with experiment. Individual redox potentials give rise to a free energy profile for electron transport that is approximately symmetric with respect to the center of the protein. Our calculations indicate that there is no significant potential bias along the orthogonal octa- and tetra-heme chains, suggesting that under aqueous conditions MtrF is a nearly reversible two-dimensional conductor.