Delineating redox cooperativity in water-soluble and membrane multiheme cytochromes through protein design.

Nature has evolved diverse electron transport proteins and multiprotein assemblies essential to the generation and transduction of biological energy. However, substantially modifying or adapting these proteins for user-defined applications or to gain fundamental mechanistic insight can be hindered by their inherent complexity. De novo protein design offers an attractive route to stripping away this confounding complexity, enabling us to probe the fundamental workings of these bioenergetic proteins and systems, while providing robust, modular platforms for constructing completely artificial electron-conducting circuitry. Here, we use a set of de novo designed mono-heme and di-heme soluble and membrane proteins to delineate the contributions of electrostatic micro-environments and dielectric properties of the surrounding protein medium on the inter-heme redox cooperativity that we have previously reported. Experimentally, we find that the two heme sites in both the water-soluble and membrane constructs have broadly equivalent redox potentials in isolation, in agreement with Poisson-Boltzmann Continuum Electrostatics calculations. BioDC, a Python program for the estimation of electron transfer energetics and kinetics within multiheme cytochromes, also predicts equivalent heme sites, and reports that burial within the low dielectric environment of the membrane strengthens heme-heme electrostatic coupling. We conclude that redox cooperativity in our diheme cytochromes is largely driven by heme electrostatic coupling and confirm that this effect is greatly strengthened by burial in the membrane. These results demonstrate that while our de novo proteins present minimalist, new-to-nature constructs, they enable the dissection and microscopic examination of processes fundamental to the function of vital, yet complex, bioenergetic assemblies.

Oxochlorin frameworks confining a ß-hydroxyketone moiety.

Nominally, meso-hydroxyoxochlorins, like known 5-hydroxy-7-oxo-octaethylchlorin (9), its nickel complex [5-hydroxy-7-oxo-octaethylchlorinato]nickel(II) (9Ni), or the novel 5-hydroxy-7,17-dioxo-octaethylbacteriochlorin (10), incorporate an acetylacetonate (acac)-moiety in the enol form in their chromophore structures. X-Ray diffraction studies of the compounds show the presence of a strong H-bond between the enol and flanking ß-ketone. Like acac, the functionality can be deprotonated. However, unlike regular acac-like moieties, we did not find any indication that this functionality is competent in chelating any of the 3d or 4d transition metal ions tested. Evidently, the conjugation that contributes to the stability of acac as a ligand cannot be expressed in the meso-hydroxyoxochlorins since it would perturb the aromaticity of the porphyrinic chromophores; in other words, the metal binding energies do not offset the loss in aromaticity. The halochromic properties of the molecules provide some more insight into the location of the protonation/deprotonation sites. The interpretation of the findings is supported by computations.

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.

Universal pictures: A lithophane codex helps teenagers with blindness visualize nanoscopic systems.

People with blindness have limited access to the high-resolution graphical data and imagery of science. Here, a lithophane codex is reported. Its pages display tactile and optical readouts for universal visualization of data by persons with or without eyesight. Prototype codices illustrated microscopy of butterfly chitin-from N-acetylglucosamine monomer to fibril, scale, and whole insect-and were given to high schoolers from the Texas School for the Blind and Visually Impaired. Lithophane graphics of Fischer-Spier esterification reactions and electron micrographs of biological cells were also 3D-printed, along with x-ray structures of proteins (as millimeter-scale 3D models). Students with blindness could visualize (describe, recall, distinguish) these systems-for the first time-at the same resolution as sighted peers (average accuracy = 88%). Tactile visualization occurred alongside laboratory training, synthesis, and mentoring by chemists with blindness, resulting in increased student interest and sense of belonging in science.

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.

Oxidation of Dueling Cysteine Promotes Subunit Exchange in SOD1.

The heterodimerization of wild-type (WT) Cu, Zn superoxide dismutase-1 (SOD1) and mutant SOD1 might be a critical step in the pathogenesis of SOD1-linked amyotrophic lateral sclerosis (ALS). Post-translational modifications that accelerate SOD1 heterodimerization remain unidentified. Here, we used capillary electrophoresis to quantify the effect of cysteine-111 oxidation on the rate and free energy of ALS mutant/WT SOD1 heterodimerization. The oxidation of Cys111-ß-SH to sulfinic and sulfonic acid (by hydrogen peroxide) increased rates of heterodimerization (with unoxidized protein) by ∼3-fold. Cysteine oxidation drove the equilibrium free energy of SOD1 heterodimerization by up to ΔΔG = -5.11 ± 0.36 kJ mol-1. Molecular dynamics simulations suggested that this enhanced heterodimerization, between oxidized homodimers and unoxidized homodimers, was promoted by electrostatic repulsion between the two "dueling" Cys111-SO2-/SO3-, which point toward one another in the homodimeric state. Together, these results suggest that oxidation of Cys-111 promotes subunit exchange between oxidized homodimers and unoxidized homodimers, regardless of whether they are mutant or WT dimers.

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.

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.

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 (⟨Hmn⟩), reaction free energies (ΔGmn), reorganization energies (λmn), and activation energies (Ea) by at most |0.002|, |0.050|, |0.120|, and |0.045| eV, respectively. Changes in ΔGmn 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.

[3 + 2]-Cycloadditions with Porphyrin ß,ß'-Bonds: Theoretical Basis of the Counterintuitive meso-Aryl Group Influence on the Rates of Reaction.

Removal of a ß,ß'-bond from meso-tetraarylporphyrin using [3 + 2]-cycloadditions generates meso-tetraarylhydroporphyrins. Literature evidence indicates that meso-tetraphenylporphyrins react more sluggishly with 1,3-dipoles such as ylides and OsO4 (in the presence of pyridine) than meso-tetrakis(pentafluorophenyl)porphyrin. The trend is counterintuitive for the reaction with OsO4, as this formal oxidation reaction is expected to proceed more readily with more electron-rich substrates. This work presents a density functional theory-based computational study of the frontier molecular orbital (FMO) interactions and reaction profile thermodynamics involved in the reaction of archetypical cycloaddition reactions (a simple ylide, OsO4, OsO4·py, OsO4·(py)2, and ozone) with the ß,ß'-double bonds of variously fluorinated meso-arylporphyrins. The trend observed for the Type I cycloaddition of an ylide is straightforward, as lowering the LUMO of the porphyrin with increasing meso-phenyl-fluorination also lowers the reaction barrier. The corresponding simple FMO analyses of Type III cycloadditions do not correctly model the reaction energetics. This is because increasing fluorination leads to lowering of the porphyrin HOMO-2, thus increasing the reaction barrier. However, coordination of pyridine to OsO4 preorganizes the transition state complex; lowering of the energy barrier by the preorganization exceeds the increase in repulsive orbital interactions, overall accelerating the cycloaddition and rationalizing the counterintuitive experimental findings.

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.

Syntheses and Aromaticity Parameters of Hexahydroxypyrrocorphin, Porphotrilactones, and Their Oxidation State Intermediates.

Triple OsO4-mediated dihydroxylation of meso-tetrakis(pentafluorophenyl)porphyrin formed a non-aromatic hexahydroxypyrrocorphin as a single stereo-isomer. A one-step oxidative conversion of all three diol functionalities to lactone moieties generated three out of the four possible porphotrilactone regioisomers that were spectroscopically and structurally characterized. This conversion recovered most of the porphyrinic macrocycle aromatic ring current, as seen in their 1H NMR spectra and modeled using DFT computations. Stepwise OsO4-mediated dihydroxylations of porpho-mono- and -di-lactones generated intermediate oxidation state compounds between the pyrrole-three pyrroline macrocycle of the pyrrocorphin and the pyrrole-three oxazolone chromophore of the trilactones. The aromaticity of these chromophores was reduced with increasing number of oxazolone to pyrroline replacements, showing the importance for the presence of three lactone moieties for the retention of the macrocycle aromaticity in the tris-ß,ß'-modified macrocycles. This work first describes hexahydoxypyrrocorphins, porphotrislactones, and the oxidation state intermediates between them; furthers the understanding of the roles of ß-lactone moieties in the expression of porphyrinic macrocycle aromaticity; and generally broadens access to chemically stable pyrrocorphins and pyrrocorphin analogues.

Data for all: Tactile graphics that light up with picture-perfect resolution.

People who are blind do not have access to graphical data and imagery produced by science. This exclusion complicates learning and data sharing between sighted and blind persons. Because blind people use tactile senses to visualize data (and sighted people use eyesight), a single data format that can be easily visualized by both is needed. Here, we report that graphical data can be three-dimensionally printed into tactile graphics that glow with video-like resolution via the lithophane effect. Lithophane forms of gel electropherograms, micrographs, electronic and mass spectra, and textbook illustrations could be interpreted by touch or eyesight at ≥79% accuracy (n = 360). The lithophane data format enables universal visualization of data by people regardless of their level of eyesight.

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

Although proteins are considered as nonconductors that transfer electrons only up to 1 to 2 nanometers via tunneling, Geobacter sulfurreducens transports respiratory electrons over micrometers, to insoluble acceptors or syntrophic partner cells, via nanowires composed of polymerized cytochrome OmcS. However, the mechanism enabling this long-range conduction is unclear. Here, we demonstrate that individual nanowires exhibit theoretically predicted hopping conductance, at rate (>1010 s-1) comparable to synthetic molecular wires, with negligible carrier loss over micrometers. Unexpectedly, nanowires show a 300-fold increase in their intrinsic conductance upon cooling, which vanishes upon deuteration. Computations show that cooling causes a massive rearrangement of hydrogen bonding networks in nanowires. Cooling makes hemes more planar, as revealed by Raman spectroscopy and simulations, and lowers their reduction potential. We find that the protein surrounding the hemes acts as a temperature-sensitive switch that controls charge transport by sensing environmental perturbations. Rational engineering of heme environments could enable systematic tuning of extracellular respiration.

Stepwise Reduction of ß-Trioxopyrrocorphins: Collapse of the Oxo-Induced Macrocycle Aromaticity.

The diatropic ring current that characterizes the unexpectedly aromatic octaethyltrioxopyrrocorphins gets drastically reduced upon chemical reduction of one and particularly two ketone moieties. With increasing reduction, the chromophores containing one pyrrole, one/two pyrrolinone, and one/two pyrrolines become more similar to regular, nonmacrocycle-aromatic pyrrocorphins (hexahydroporphyrins). Single-crystal diffraction analysis shows the reduction products to be idealized planar. With increasing reduction, their UV-vis spectroscopic signatures are those of conjugated but nonaromatic oligopyrroles. Their diatropic ring currents, as assessed by 1H NMR spectroscopy, showed them to possess largely nonaromatic π-systems. Dihydroxylation of select ß,ß'-dioxobacteriochlorin and ß,ß'-dioxoisobacteriochlorins also resulted in the formation of equivalent mixed pyrrole/two pyrrolinone/pyrroline chromophores. Computations were able to reproduce the experimental trends of the diatropic ring currents and filled in the data for the regioisomers that could not be experimentally accessed. The work further highlights the electronic influence of the ß-oxo-substituents and, more specifically, the origin of the aromaticity of the trioxopyrrocorphins. It also presents a series of chemically robust pyrrocorphins, a chromophore class for which many chemically very sensitive members have been reported.

Electronic energy levels of porphyrins are influenced by the local chemical environment.

Self-assembled islands of 5,10,15,20-tetrakis(pentafluoro-phenyl)porphyrin (2HTFPP) on Au(111) contain two bistable molecular species that differ by shifted electronic energy levels. Interactions with the underlying gold herringbone reconstruction and neighboring 2HTFPP molecules cause approximately 60% of molecules to have shifted electronic energy levels. We observed the packing density decrease from 0.64 ± 0.04 molecules per nm2 to 0.38 ± 0.03 molecules per nm2 after annealing to 200 °C. The molecules with shifted electronic energy levels show longer-range hexagonal packing or are adjacent to molecular vacancies, indicating that molecule-molecule and molecule-substrate interactions contribute to the shifted energies. Multilayers of porphyrins do not exhibit the same shifting of electronic energy levels which strongly suggests that molecule-substrate interactions play a critical role in stabilization of two electronic species of 2HTFPP on Au(111).

Making protons tag along with electrons.

Every living cell needs to get rid of leftover electrons when metabolism extracts energy through the oxidation of nutrients. Common soil microbes such as Geobacter sulfurreducens live in harsh environments that do not afford the luxury of soluble, ingestible electron acceptors like oxygen. Instead of resorting to fermentation, which requires the export of reduced compounds (e.g. ethanol or lactate derived from pyruvate) from the cell, these organisms have evolved a means to anaerobically respire by using nanowires to export electrons to extracellular acceptors in a process called extracellular electron transfer (EET) [ 1]. Since 2005, these nanowires were thought to be pili filaments [ 2]. But recent studies have revealed that nanowires are composed of multiheme cytochromes OmcS [ 3, 4] and OmcZ [ 5] whereas pili remain inside the cell during EET and are required for the secretion of nanowires [ 6]. However, how electrons are passed to these nanowires remains a mystery ( Figure 1A). Periplasmic cytochromes (Ppc) called PpcA-E could be doing the job, but only two of them (PpcA and PpcD) can couple electron/proton transfer - a necessary condition for energy generation. In a recent study, Salgueiro and co-workers selectively replaced an aromatic with an aliphatic residue to couple electron/proton transfer in PpcB and PpcE (Biochem. J. 2021, 478 (14): 2871-2887). This significant in vitro success of their protein engineering strategy may enable the optimization of bioenergetic machinery for bioenergy, biofuels, and bioelectronics applications.

ß-Trioxopyrrocorphins: pyrrocorphins of graded aromaticity.

Octaethyltrioxopyrrocorphins unexpectedly show macrocycle-aromatic properties, even though they contain the macrocyclic π-system of the non-aromatic pyrrocorphins (hexahydroporphyrins). Two of the four possible triketone regioisomers were first reported in 1969 by one-pot oxidation of octaethylporphyrin but remained essentially unexplored since. We detail here the targeted preparation of the remaining two triketone isomers and the optical and NMR spectroscopic properties of all isomers. All four regioisomers possess unique electronic properties, including broadly varying degrees of diatropicity that were experimentally determined using 1H NMR spectroscopy and computationally verified. Structural patterns modulating the aromaticity were recognized. These differences highlight the regioisomerically differentiated influences of the three ß-oxo-functionalities. We also present the solid state structure of the two most common isomers (in their free base form or as zinc complexes), allowing further conclusions to be made about the resonance structures present in these triketones. Remarkably, also, the halochromic properties of the triketones differ sharply from those of regular (hydro)porphyrins, providing further support for the proposed 16-membered, 18 π-electron aromatic ring-current. The work conceptually expands the understanding of tris-modified hydroporphyrinoid analogues and the factors that enable and control porphyrinoid aromaticity.

Computational insights into the membrane fusion mechanism of SARS-CoV-2 at the cellular level.

The membrane fusion mechanism of SARS-CoV-2 offers an attractive target for the development of small molecule antiviral inhibitors. Fusion involves an initial binding of the crown-like trimeric spike glycoproteins of SARS-CoV-2 to the receptor angiotensin II-converting enzyme 2 (ACE2) on the permissive host cellular membrane and a prefusion to post-fusion conversion of the spike trimer. During this conversion, the fusion peptides of the spike trimer are inserted into the host membrane to bring together the host and viral membranes for membrane fusion in highly choreographic events. However, geometric constraints due to interactions with the membranes remain poorly understood. In this study, we build structural models of super-complexes of spike trimer/ACE2 dimers based on the molecular structures of the ACE2/neutral amino acid transporter B(0)AT heterodimer. We determine the conformational constraints due to the membrane geometry on the enzymatic activity of ACE2 and on the viral fusion process. Furthermore, we find that binding three ACE2 dimers per spike trimer is essential to open the central pore as necessary for triggering productive membrane fusion through an elongation of the central stalk. The reported findings thus provide valuable insights for targeting the membrane fusion mechanism for drug design at the molecular level.