A 300-fold conductivity increase in microbial cytochrome nanowires due to temperature-induced restructuring of hydrogen bonding networks

On the fundamentals of quantum rate theory and the long-range electron transport in respiratory chains

It has been shown that both the electron-transfer rate constant of an electrochemical reaction and the conductance quantum are correlated with the concept of quantum capacitance. This simple association between the two separate concepts has an entirely quantum rate basis that encompasses the electron-transfer rate theory as originally proposed by Rudolph A. Marcus whether statistical mechanics is appropriately taken into account. I have prepared a concise review of the quantum mechanical rate theory principles focused on its quantum electrodynamics character to demonstrate that it can reconcile the conflicting views established on attempting to use the super-exchange (supported on electron transfer) or 'metallic-like' (supported on conductance quantum) mechanisms separately to explain the highly efficient long-range electron transport observed in the respiratory processes of living cells. The unresolved issues related to long-range electron transport are clarified in light of the quantum rate theory with a discussion focused on Geobacter sulfurreducens films as a reference standard of the respiration chain. Theoretical analyses supported by experimental data suggest that the efficiency of respiration within a long-range electron transport path is intrinsically a quantum mechanical event that follows relativistic quantum electrodynamics principles as addressed by quantum rate theory.

To be or not to be a cytochrome: electrical characterizations are inconsistent with Geobacter cytochrome 'nanowires'

Geobacter sulfurreducens profoundly shapes Earth's biogeochemistry by discharging respiratory electrons to minerals and other microbes through filaments of a two-decades-long debated identity. Cryogenic electron microscopy has revealed filaments of redox-active cytochromes, but the same filaments have exhibited hallmarks of organic metal-like conductivity under cytochrome denaturing/inhibiting conditions. Prior structure-based calculations and kinetic analyses on multi-heme proteins are synthesized herein to propose that a minimum of ~7 cytochrome 'nanowires' can carry the respiratory flux of a Geobacter cell, which is known to express somewhat more (≥20) filaments to increase the likelihood of productive contacts. By contrast, prior electrical and spectroscopic structural characterizations are argued to be physiologically irrelevant or physically implausible for the known cytochrome filaments because of experimental artifacts and sample impurities. This perspective clarifies our mechanistic understanding of physiological metal-microbe interactions and advances synthetic biology efforts to optimize those interactions for bioremediation and energy or chemical production.

Widespread extracellular electron transfer pathways for charging microbial cytochrome OmcS nanowires via periplasmic cytochromes PpcABCDE

Extracellular electron transfer (EET) via microbial nanowires drives globally-important environmental processes and biotechnological applications for bioenergy, bioremediation, and bioelectronics. Due to highly-redundant and complex EET pathways, it is unclear how microbes wire electrons rapidly (>106 s-1) from the inner-membrane through outer-surface nanowires directly to an external environment despite a crowded periplasm and slow (<105 s-1) electron diffusion among periplasmic cytochromes. Here, we show that Geobacter sulfurreducens periplasmic cytochromes PpcABCDE inject electrons directly into OmcS nanowires by binding transiently with differing efficiencies, with the least-abundant cytochrome (PpcC) showing the highest efficiency. Remarkably, this defined nanowire-charging pathway is evolutionarily conserved in phylogenetically-diverse bacteria capable of EET. OmcS heme reduction potentials are within 200 mV of each other, with a midpoint 82 mV-higher than reported previously. This could explain efficient EET over micrometres at ultrafast (<200 fs) rates with negligible energy loss. Engineering this minimal nanowire-charging pathway may yield microbial chassis with improved performance.

Optimizing electrochemically active microorganisms as a key player in the bioelectrochemical system: Identification methods and pathways to large-scale implementation

The rapid global economic growth driven by industrialization and population expansion has resulted in significant issues, including reliance on fossil fuels, energy scarcity, water crises, and environmental emissions. To address these issues, bioelectrochemical systems (BES) have emerged as a dual-purpose solution, harnessing electrochemical processes and the capabilities of electrochemically active microorganisms (EAM) to simultaneously recover energy and treat wastewater. This review examines critical performance factors in BES, including inoculum selection, pretreatment methods, electrodes, and operational conditions. Further, authors explore innovative approaches to suppress methanogens and simultaneously enhance the EAM in mixed cultures. Additionally, advanced techniques for detecting EAM are discussed. The rapid detection of EAM facilitates the selection of suitable inoculum sources and optimization of enrichment strategies in BESs. This optimization is essential for facilitating the successful scaling up of BES applications, contributing substantially to the realization of clean energy and sustainable wastewater treatment. This analysis introduces a novel viewpoint by amalgamating contemporary research on the selective enrichment of EAM in mixed cultures. It encompasses identification and detection techniques, along with methodologies tailored for the selective enrichment of EAM, geared explicitly toward upscaling applications in BES.

Electron transfer in a crystalline cytochrome with four hemes

Diffusion of electrons over distances on the order of 100 μm has been observed in crystals of a small tetraheme cytochrome (STC) from Shewanella oneidensis [J. Huang et al. J. Am. Chem. Soc. 142, 10459-10467 (2020)]. Electron transfer between hemes in adjacent subunits of the crystal is slower and more strongly dependent on temperature than had been expected based on semiclassical electron-transfer theory. We here explore explanations for these findings by molecular-dynamics simulations of crystalline and monomeric STC. New procedures are developed for including time-dependent quantum mechanical energy differences in the gap between the energies of the reactant and product states and for evaluating fluctuations of the electronic-interaction matrix element that couples the two hemes. Rate constants for electron transfer are calculated from the time- and temperature-dependent energy gaps, coupling factors, and Franck-Condon-weighted densities of states using an expression with no freely adjustable parameters. Back reactions are considered, as are the effects of various protonation states of the carboxyl groups on the heme side chains. Interactions with water are found to dominate the fluctuations of the energy gap between the reactant and product states. The calculated rate constant for electron transfer from heme IV to heme Ib in a neighboring subunit at 300 K agrees well with the measured value. However, the calculated activation energy of the reaction in the crystal is considerably smaller than observed. We suggest two possible explanations for this discrepancy. The calculated rate constant for transfer from heme I to II within the same subunit of the crystal is about one-third that for monomeric STC in solution.

Long-Range Electron Transport Rates Depend on Wire Dimensions in Cytochrome Nanowires

The ability to redirect electron transport to new reactions in living systems opens possibilities to store energy, generate new products, or probe physiological processes. Recent work by Huang et al. showed that 3D crystals of small tetraheme cytochromes (STC) can transport electrons over nanoscopic to mesoscopic distances by an electron hopping mechanism, making them promising materials for nanowires. However, fluctuations at room temperature may distort the nanostructure, hindering efficient electron transport. Classical molecular dynamics simulations of these fluctuations at the nano- and mesoscopic scales allowed us to develop a graph network representation to estimate maximum electron flow that can be driven through STC wires. In longer nanowires, transient structural fluctuations at protein-protein interfaces tended to obstruct efficient electron transfer, but these blockages are ameliorated in thicker crystals where alternative electron transfer pathways become more efficient. The model implies that more flexible proteinprotein interfaces limit the required minimum diameter to carry currents commensurate with conventional electronics.

Heavy Water Reduces the Electronic Conductance of Protein Wires via Deuteron Interactions with Aromatic Residues

Proteins are versatile, self-assembling nanoelectronic components, but their hopping conductivity is expected to be influenced by solvent fluctuations. The role of the solvent was investigated by measuring the single molecule conductance of several proteins in both H2O and D2O. The conductance of a homologous series of protein wires decreases more rapidly with length in D2O, indicating a 6-fold decrease in carrier diffusion constant relative to the same protein in H2O. The effect was found to depend on the specific aromatic amino acid composition. A tryptophan zipper protein showed a decrease in conductance similar to that of the protein wires, whereas a phenylalanine zipper protein was insensitive to solvent changes. Tryptophan contains an indole amine, whereas the phenylalanine aromatic ring has no exchangeable protons, so the effect of heavy water on conductance is a consequence of specific D- or H-interactions with the aromatic residues.

Long-Range Conductivity in Proteins Mediated by Aromatic Residues

Single-molecule measurements show that many proteins, lacking any redox cofactors, nonetheless exhibit electrical conductance on the order of a nanosiemen over 10 nm distances, implying that electrons can transit an entire protein in less than a nanosecond when subject to a potential difference of less than 1 V. This is puzzling because, for fast transport (i.e., a free energy barrier of zero), the hopping rate is determined by the reorganization energy of approximately 0.8 eV, and this sets the time scale of a single hop to at least 1 μs. Furthermore, the Fermi energies of typical metal electrodes are far removed from the energies required for sequential oxidation and reduction of the aromatic residues of the protein, which should further reduce the hopping current. Here, we combine all-atom molecular dynamics (MD) simulations of non-redox-active proteins (consensus tetratricopeptide repeats) with an electron transfer theory to demonstrate a molecular mechanism that can account for the unexpectedly fast electron transport. According to our MD simulations, the reorganization energy produced by the energy shift on charging (the Stokes shift) is close to the conventional value of 0.8 eV. However, the non-ergodic sampling of molecular configurations by the protein results in reaction-reorganization energies, extracted directly from the distribution of the electrostatic energy fluctuations, that are only ∼0.2 eV, which is small enough to enable long-range conductivity, without invoking quantum coherent transport. Using the MD values of the reorganization energies, we calculate a current decay with distance that is in agreement with experiment.

Structural Determinants of Redox Conduction Favor Robustness over Tunability in Microbial Cytochrome Nanowires

Structural determinants of a 103-fold variation in electrical conductivity for helical homopolymers of tetra-, hexa-, and octa-heme cytochromes (named Omc- E, S, and Z, respectively) from Geobacter sulfurreducens are investigated with the Pathways model for electron tunneling, classical molecular dynamics, and hybrid quantum/classical molecular mechanics. Thermally averaged electronic couplings for through-space heme-to-heme electron transfer in the "nanowires" computed with density functional theory are ≤0.015 eV. Pathways analyses also indicate that couplings match within a factor of 5 for all "nanowires", but some alternative tunneling routes are found involving covalent protein backbone bonds (Omc- S and Z) or propionic acid-ligating His H-bonds on adjacent hemes (OmcZ). Reorganization energies computed from electrostatic vertical energy gaps or a version of the Marcus continuum expression parameterized on the total (donor + acceptor) solvent-accessible surface area typically agree within 20% and fall within the range 0.48-0.98 eV. Reaction free energies in all three "nanowires" are ≤|0.28| eV, even though Coulombic interactions primarily tune the site redox energies by 0.7-1.2 eV. Given the conserved energetic parameters, redox conductivity differs by < 103-fold among the cytochrome "nanowires". Redox currents do not exceed 3.0 × 10-3 pA at a physiologically relevant 0.1 V bias, with the slowest electron transfers being on a (μs) timescale much faster than typical (ms) enzymatic turnovers. Thus, the "nanowires" are proposed to be functionally robust to variations in structure that provide a habitat-customized protein interface. The 30 pA to 30 nA variation in conductivity previously reported from atomic force microscopy experiments is not intrinsic to the structures and/or does not result from the physiologically relevant redox conduction mechanism.

Electron Tunneling in Biology: When Does it Matter?

Electrons can tunnel between cofactor molecules positioned along biological electron transport chains up to a distance of ≃ 20 Å on the millisecond time scale of enzymatic turnover. This tunneling range determines the design of biological energy chains facilitating the cross-membrane transport of electrons. Tunneling distance and cofactors' redox potentials become the main physical parameters affecting the rate of electron transport. In addition, universal charge-transport properties are assigned to all proteins, making protein identity, flexibility, and dynamics insignificant. This paradigm is challenged by dynamical models of electron transfer, showing that the electron hopping rate is constant within the crossover distance R* ≃ 12 Å, followed with an exponential falloff at longer distances. If this hypothesis is fully confirmed, natural and man-made energy chains for electron transport should be best designed by placing redox cofactors near the crossover distance R*. Protein flexibility and dynamics affect the magnitude of the maximum hopping rate within the crossover distance. Changes in protein flexibility between forward and backward transitions contribute to vectorial charge transport. For biological energy chains, charge transport through proteins is not defined by universal parameters, and protein identity matters.

Extracellular cytochrome nanowires appear to be ubiquitous in prokaryotes

Electrically conductive appendages from the anaerobic bacterium Geobacter sulfurreducens, recently identified as extracellular cytochrome nanowires (ECNs), have received wide attention due to numerous potential applications. However, whether other organisms employ similar ECNs for electron transfer remains unknown. Here, using cryoelectron microscopy, we describe the atomic structures of two ECNs from two major orders of hyperthermophilic archaea present in deep-sea hydrothermal vents and terrestrial hot springs. Homologs of Archaeoglobus veneficus ECN are widespread among mesophilic methane-oxidizing Methanoperedenaceae, alkane-degrading Syntrophoarchaeales archaea, and in the recently described megaplasmids called Borgs. The ECN protein subunits lack similarities in their folds; however, they share a common heme arrangement, suggesting an evolutionarily optimized heme packing for efficient electron transfer. The detection of ECNs in archaea suggests that filaments containing closely stacked hemes may be a common and widespread mechanism for long-range electron transfer in both prokaryotic domains of life.

Hybrid Quantum Mechanical/Molecular Mechanical Methods For Studying Energy Transduction in Biomolecular Machines

Hybrid quantum mechanical/molecular mechanical (QM/MM) methods have become indispensable tools for the study of biomolecules. In this article, we briefly review the basic methodological details of QM/MM approaches and discuss their applications to various energy transduction problems in biomolecular machines, such as long-range proton transports, fast electron transfers, and mechanochemical coupling. We highlight the particular importance for these applications of balancing computational efficiency and accuracy. Using several recent examples, we illustrate the value and limitations of QM/MM methodologies for both ground and excited states, as well as strategies for calibrating them in specific applications. We conclude with brief comments on several areas that can benefit from further efforts to make QM/MM analyses more quantitative and applicable to increasingly complex biological problems.

Determinants of synergistic cell-cell interactions in bacteria

Bacteria are ubiquitous and colonize virtually every conceivable habitat on earth. To achieve this, bacteria require different metabolites and biochemical capabilities. Rather than trying to produce all of the needed materials by themselves, bacteria have evolved a range of synergistic interactions, in which they exchange different commodities with other members of their local community. While it is widely acknowledged that synergistic interactions are key to the ecology of both individual bacteria and entire microbial communities, the factors determining their establishment remain poorly understood. Here we provide a comprehensive overview over our current knowledge on the determinants of positive cell-cell interactions among bacteria. Taking a holistic approach, we review the literature on the molecular mechanisms bacteria use to transfer commodities between bacterial cells and discuss to which extent these mechanisms favour or constrain the successful establishment of synergistic cell-cell interactions. In addition, we analyse how these different processes affect the specificity among interaction partners. By drawing together evidence from different disciplines that study the focal question on different levels of organisation, this work not only summarizes the state of the art in this exciting field of research, but also identifies new avenues for future research.

Living electronics: A catalogue of engineered living electronic components

Biology leverages a range of electrical phenomena to extract and store energy, control molecular reactions and enable multicellular communication. Microbes, in particular, have evolved genetically encoded machinery enabling them to utilize the abundant redox-active molecules and minerals available on Earth, which in turn drive global-scale biogeochemical cycles. Recently, the microbial machinery enabling these redox reactions have been leveraged for interfacing cells and biomolecules with electrical circuits for biotechnological applications. Synthetic biology is allowing for the use of these machinery as components of engineered living materials with tuneable electrical properties. Herein, we review the state of such living electronic components including wires, capacitors, transistors, diodes, optoelectronic components, spin filters, sensors, logic processors, bioactuators, information storage media and methods for assembling these components into living electronic circuits.

Structure of Geobacter cytochrome OmcZ identifies mechanism of nanowire assembly and conductivity

OmcZ nanowires produced by Geobacter species have high electron conductivity (>30 S cm-1). Of 111 cytochromes present in G. sulfurreducens, OmcZ is the only known nanowire-forming cytochrome essential for the formation of high-current-density biofilms that require long-distance (>10 µm) extracellular electron transport. However, the mechanisms underlying OmcZ nanowire assembly and high conductivity are unknown. Here we report a 3.5-Å-resolution cryogenic electron microscopy structure for OmcZ nanowires. Our structure reveals linear and closely stacked haems that may account for conductivity. Surface-exposed haems and charge interactions explain how OmcZ nanowires bind to diverse extracellular electron acceptors and how organization of nanowire network re-arranges in different biochemical environments. In vitro studies explain how G. sulfurreducens employ a serine protease to control the assembly of OmcZ monomers into nanowires. We find that both OmcZ and serine protease are widespread in environmentally important bacteria and archaea, thus establishing a prevalence of nanowire biogenesis across diverse species and environments.

Effect of Water Deuteration on Protein Electron Transfer

Traditional theories of long-range protein electron transfer describe the reaction rate in terms of the tunneling distance and the reaction free energy. They do not recognize two physical effects: (i) local wetting of the active site by hydration water and (ii) protein identity affecting the rate through dynamics and flexibility. We find, by molecular dynamics simulations, a significant, ∼25 times, slowing down of the rate of protein electron transfer upon deuteration. H/D substitution changes the rate constant pre-exponential factor in the regime of electron transfer controlled by medium dynamics. Switching from light to heavy water increases the effective medium relaxation time. The effect is caused by both a global change in the flexibility of the protein backbone and locally stronger hydrogen bonds to charged residues.

Structural Determinants of Redox Conduction Favor Robustness over Tunability in Microbial Cytochrome Nanowires

Helical homopolymers of multiheme cytochromes catalyze biogeochemically significant electron transfers with a reported 10 3 -fold variation in conductivity. Herein, classical molecular dynamics and hybrid quantum/classical molecular mechanics are used to elucidate the structural determinants of the redox potentials and conductivities of the tetra-, hexa-, and octaheme outer-membrane cytochromes E, S, and Z, respectively, from Geobacter sulfurreducens . Second-sphere electrostatic interactions acting on minimally polarized heme centers are found to regulate redox potentials over a computed 0.5-V range. However, the energetics of redox conduction are largely robust to the structural diversity: Single-step electronic couplings (⟨H mn ⟩), reaction free energies , and reorganization energies (λ mn ) are always respectively <|0.026|, <|0.26|, and between 0.5 - 1.0 eV. With these conserved parameter ranges, redox conductivity differed by less than a factor of 10 among the 'nanowires' and is sufficient to meet the demands of cellular respiration if 10 2 - 10 3 'nanowires' are expressed. The 'nanowires' are proposed to be differentiated by the protein packaging to interface with a great variety of environments, and not by conductivity, because the rate-limiting electron transfers are elsewhere in the respiratory process. Conducting-probe atomic force microscopy measurements that find conductivities 10 3 -10 6 -fold more than cellular demands are suggested to report on functionality that is either not used or not accessible under physiological conditions. The experimentally measured difference in conductivity between Omc- S and Z is suggested to not be an intrinsic feature of the CryoEM-resolved structures.

Theory of Protein Charge Transfer: Electron Transfer between Tryptophan Residue and Active Site of Azurin

One reaction step in the conductivity relay of azurin, electron transfer between the Cu-based active site and the tryptophan residue, is studied theoretically and by classical molecular dynamics simulations. Oxidation of tryptophan results in electrowetting of this residue. This structural change makes the free energy surfaces of electron transfer nonparabolic as described by the Q-model of electron transfer. We analyze the medium dynamical effect on protein electron transfer produced by coupled Stokes-shift dynamics and the dynamics of the donor-acceptor distance modulating electron tunneling. The equilibrium donor-acceptor distance falls in the plateau region of the rate constant, where it is determined by the protein-water dynamics, and the probability of electron tunneling does not affect the rate. The crossover distance found here puts most intraprotein electron-transfer reactions under the umbrella of dynamical control. The crossover between the medium-controlled and tunneling-controlled kinetics is combined with the effect of the protein-water medium on the activation barrier to formulate principles of tunability of protein-based charge-transfer chains. The main principle in optimizing the activation barrier is the departure from the Gaussian-Gibbsian statistics of fluctuations promoting activated transitions. This is achieved either by incomplete (nonergodic) sampling, breaking the link between the Stokes-shift and variance reorganization energies, or through wetting-induced structural changes of the enzyme's active site.

Assessing Thermal Response of Redox Conduction for Anti-Arrhenius Kinetics in a Microbial Cytochrome Nanowire

A micrometers-long helical homopolymer of the outer-membrane cytochrome type S (OmcS) from Geobacter sulfurreducens is proposed to transport electrons to extracellular acceptors in an ancient respiratory strategy of biogeochemical and technological significance. OmcS surprisingly exhibits higher conductivity upon cooling (anti-Arrhenius kinetics), an effect previously attributed to H-bond restructuring and heme redox potential shifts. Herein, the temperature sensitivity of redox conductivity is more thoroughly examined with conventional and constant-redox and -pH molecular dynamics and quantum mechanics/molecular mechanics. A 30 K drop in temperature constituted a weak perturbation to electron transfer energetics, changing electronic couplings (⟨H_mn⟩), reaction free energies (ΔG_mn), reorganization energies (λ_mn), and activation energies (E_a) by at most |0.002|, |0.050|, |0.120|, and |0.045| eV, respectively. Changes in ΔG_mn reflected -0.07 ± 0.03 V shifts in redox potentials that were caused in roughly equal measure by altered electrostatic interactions with the solvent and protein. Changes in intraprotein H-bonding reproduced the earlier observations. Single-particle diffusion and multiparticle steady-state flux models, parametrized with Marcus theory rates, showed that biologically relevant incoherent hopping cannot qualitatively or quantitatively describe electrical conductivity measured by atomic force microscopy in filamentous OmcS. The discrepancy is attributed to differences between solution-phase simulations and solid-state measurements and the need to model intra- and intermolecular vibrations explicitly.

Microbial biofilms as living photoconductors due to ultrafast electron transfer in cytochrome OmcS nanowires

Light-induced microbial electron transfer has potential for efficient production of value-added chemicals, biofuels and biodegradable materials owing to diversified metabolic pathways. However, most microbes lack photoactive proteins and require synthetic photosensitizers that suffer from photocorrosion, photodegradation, cytotoxicity, and generation of photoexcited radicals that are harmful to cells, thus severely limiting the catalytic performance. Therefore, there is a pressing need for biocompatible photoconductive materials for efficient electronic interface between microbes and electrodes. Here we show that living biofilms of Geobacter sulfurreducens use nanowires of cytochrome OmcS as intrinsic photoconductors. Photoconductive atomic force microscopy shows up to 100-fold increase in photocurrent in purified individual nanowires. Photocurrents respond rapidly (<100 ms) to the excitation and persist reversibly for hours. Femtosecond transient absorption spectroscopy and quantum dynamics simulations reveal ultrafast (~200 fs) electron transfer between nanowire hemes upon photoexcitation, enhancing carrier density and mobility. Our work reveals a new class of natural photoconductors for whole-cell catalysis.

Despite enormous potential of solar-driven biocatalysis, most living systems lack photoactive proteins and require toxic and expensive synthetic materials limiting the performance. Here, a class of natural photoconductors is demonstrated through sub-picosecond heme-to-heme electron transfer in bacteria-produced protein nanowires.

Structure of Geobacter OmcZ filaments suggests extracellular cytochrome polymers evolved independently multiple times

While early genetic and low-resolution structural observations suggested that extracellular conductive filaments on metal-reducing organisms such as Geobacter were composed of type IV pili, it has now been established that bacterial c-type cytochromes can polymerize to form extracellular filaments capable of long-range electron transport. Atomic structures exist for two such cytochrome filaments, formed from the hexaheme cytochrome OmcS and the tetraheme cytochrome OmcE. Due to the highly conserved heme packing within the central OmcS and OmcE cores, and shared pattern of heme coordination between subunits, it has been suggested that these polymers have a common origin. We have now used cryo-electron microscopy (cryo-EM) to determine the structure of a third extracellular filament, formed from the Geobacter sulfurreducens octaheme cytochrome, OmcZ. In contrast to the linear heme chains in OmcS and OmcE from the same organism, the packing of hemes, heme:heme angles, and between-subunit heme coordination is quite different in OmcZ. A branched heme arrangement within OmcZ leads to a highly surface exposed heme in every subunit, which may account for the formation of conductive biofilm networks, and explain the higher measured conductivity of OmcZ filaments. This new structural evidence suggests that conductive cytochrome polymers arose independently on more than one occasion from different ancestral multiheme proteins.