An interesting possibility of forming special hole stepping stones with high-stacking aromatic rings in proteins: three-π five-electron and four-π seven-electron resonance bindings

Long-range hole transfer of proteins plays an important role in many biological processes of living organisms. Therefore, it is highly useful to examine the possible hole stepping stones, which can facilitate hole transfer in proteins. However, the structures of stepping stones are diverse because of the complexity of the protein structures. In the present work, we proposed a series of special stepping stones, which are instantaneously formed by three and four packing aromatic side chains of amino acids to capture a hole, corresponding to three-π five-electron (π:π∴π↔π∴π:π) and four-π seven-electron (π:π∴π:π↔π:π:π∴π) resonance bindings with appropriate binding energies. The aromatic amino acids include histidine (His), phenylalanine (Phe), tyrosine (Tyr) and tryptophan (Trp). The formations of these special stepping stones can effectively reduce the local ionization potential of the high π-stacking region to efficiently capture the migration hole. The quick formations and separations of them promote the efficient hole transfer in proteins. More interestingly, we revealed that a hole cannot delocalize over infinite aromatic rings along the high π-π packing structure at the same time and the micro-surroundings of proteins can modulate the formations of π:π∴π↔π∴π:π and π:π∴π:π↔π:π:π∴π bindings. These results may contribute a new avenue to better understand the potential hole transfer pathway in proteins.

Electron transfer pathways from quantum dynamics simulations

This work explores the possibility of simulating an electron transfer process between a donor and an acceptor in real time using time-dependent density functional theory electron dynamics. To achieve this objective, a central issue to resolve is the definition of the initial state. This must be a non-equilibrium electronic state able to trigger the charge transfer dynamics; here, two schemes are proposed to prepare such states. One is based on the combination of the density matrices of the donor and acceptor converged separately with appropriate charges (for example, -1 for the donor and +1 for the acceptor). The second approach relied on constrained DFT to localize the charge on each fragment. With these schemes, electron transfer processes are simulated in different model systems of increasing complexity: an atomic hydrogen dimer, a polyacetylene chain, and the active site of the T. cruzi hybrid type A heme peroxidase, for which two possible electron transfer paths have been postulated. For the latter system, the present methodology applied in a hybrid Quantum Mechanics - Molecular Mechanics framework allows us to establish the relative probabilities of each path and provides insight into the inhibition of the electron transfer provoked by the substitution of tryptophan by phenylalanine in the W233F mutant.

Charge Transfer and Biocompatibility Aspects in Conducting Polymer-Based Enzymatic Biosensors and Biofuel Cells

Charge transfer (CT) is a very important issue in the design of biosensors and biofuel cells. Some nanomaterials can be applied to facilitate the CT in these bioelectronics-based devices. In this review, we overview some CT mechanisms and/or pathways that are the most frequently established between redox enzymes and electrodes. Facilitation of indirect CT by the application of some nanomaterials is frequently applied in electrochemical enzymatic biosensors and biofuel cells. More sophisticated and still rather rarely observed is direct charge transfer (DCT), which is often addressed as direct electron transfer (DET), therefore, DCT/DET is also targeted and discussed in this review. The application of conducting polymers (CPs) for the immobilization of enzymes and facilitation of charge transfer during the design of biosensors and biofuel cells are overviewed. Significant attention is paid to various ways of synthesis and application of conducting polymers such as polyaniline, polypyrrole, polythiophene poly(3,4-ethylenedioxythiophene). Some DCT/DET mechanisms in CP-based sensors and biosensors are discussed, taking into account that not only charge transfer via electrons, but also charge transfer via holes can play a crucial role in the design of bioelectronics-based devices. Biocompatibility aspects of CPs, which provides important advantages essential for implantable bioelectronics, are discussed.

Tunneling matrix element and tunneling pathways of protein electron transfer calculated with a fragment molecular orbital method

Practical ways to calculate the tunneling matrix elements and analyze the tunneling pathways for protein electron-transfer (ET) reactions with a fragment molecular orbital (FMO) method are presented. The straightforward use of minimal basis sets only for the atoms involved in the covalent bond detachment in FMO can properly describe the ETs through the protein main-chains with the cost-effective two-body corrections (FMO2) without losing the quality of double-zeta basis sets. The current FMO codes have been interfaced with density functional theory, polarizable continuum model, and model core potentials, with which the FMO-based protein ET calculations can consider the effects of electron correlation, solvation, and transition-metal redox centers. The reasonable performance of the FMO-based ET calculations is demonstrated for three different sets of protein-ET model molecules: (1) hole transfer between two tryptophans covalently bridged by a polyalanine linker in the ideal α-helix and β-strand conformations, (2) ET between two plastoquinones covalently bridged by a polyalanine linker in the ideal α-helix and β-strand conformations, and (3) hole transfer between ruthenium (Ru) and copper (Cu) complexes covalently bridged by a stretch of a polyglycine linker as a model for Ru-modified derivatives of azurin.

Artificial, Photoinduced Activation of Nitrogenase Using Directed and Mediated Electron Transfer Processes

Nitrogenase, a bacteria-based enzyme, is the sole enzyme that is able to generate ammonia by atmospheric nitrogen fixation. Thus, improved understanding of its utilization and developing methods to artificially activate it may contribute to basic research, as well as to the design of future artificial systems. Here, we present methods to artificially activate nitrogenase using photoinduced reactions. Two nitrogenase variants originating from Azotobacter vinelandii were examined using photoactivated CdS nanoparticles (NPs) capped with thioglycolic acid (TGA) or 2-mercaptoethanol (ME) ligands. The effect of methyl viologen (MV) as a redox mediator of hydrogen and ammonia generation was tested and analyzed. We further determined the NPs conductive band edges and their effect on the nitrogenase photoactivation. The nano-biohybrid systems comprising CdS NPs and nitrogenase were further imaged by transmission electron microscopy, visualizing their formation for the first time. Our results show that the ME-capped CdS NPs–nitrogenase enzyme biohybrid system with added MV as a redox mediator leads to a five-fold increase in the production of ammonia compared with the non-mediated biohybrid system; nevertheless, it stills lag behind the natural process rate. On the contrary, a maximal hydrogen generation amount was achieved by the αL158C MoFe-P and the ME-capped CdS NPs.

Ergodicity Breaking in Thermal Biological Electron Transfer? Cytochrome C Revisited II

It was recently suggested that cytochrome c operates in an ergodicity-breaking regime characterized by unusually large energy gap thermal fluctuations and associated reorganization free energies for heme oxidation of up to 3.0 eV. The large fluctuations were reported to lower activation free energy for oxidation of the heme cofactor by almost a factor of 2 compared to the case where ergodicity is maintained. Our group has recently investigated this claim computationally at several levels of theory and found no evidence for such large energy gap fluctuations. Here we address the points of our earlier work that have raised criticism and we also extend our previous investigation by considering a simple linear polarizability model for cytochrome c oxidation. We find very consistent results among all our computational approaches, ranging from classical molecular dynamics, to the linear polarizability model to QM(PMM)/MM to full QM(DFT)/MM electrostatic emdedding. None of them support the notion of very large energy gap fluctuations or ergodicity breaking. The deviation between our simulations and the ones reported in [ J. Phys. Chem. B 2017, 121, 4958] is traced back to rather large electric fields at the Fe site of the heme c cofactor in that study, not seen in our simulations, neither with the AMBER nor with the CHARMM force field. While ergodicity breaking effects may well occur in other biological ET, our numerical evidence suggests that this is not the case for cytochrome c.

Printed-Circuit-Board-Based Two-Electrode System for Electronic Characterization of Proteins

Proteins have been increasingly suggested as suitable candidates for the fabrication of biological computers and other biomolecular-based electronic devices mainly due to their interesting structure-related intrinsic electrical properties. These natural biopolymers are environmentally friendly substitutes for conventional inorganic materials and find numerous applications in bioelectronics. Effective manipulation of protein biomolecules allows for accurate fabrication of nanoscaled device dimensions for miniaturized electronics. The prerequisite, however, demands an interrogation of its various electronic properties prior to understanding the complex charge transfer mechanisms in protein molecules, the knowledge of which will be crucial toward development of such nanodevices. One significantly preferred method in recent times involves the utilization of solid-state sensors where interactions of proteins could be investigated upon contact with metals such as gold. Therefore, in this work, proteins (hemoglobin and collagen) were integrated within a two-electrode system, and the resulting electronic profiles were investigated. Interestingly, structure-related electronic profiles representing semiconductive-like behaviors were observed. These characteristic electronic profiles arise from the metal (Au)-semiconductor (protein) junction, clearly demonstrating the formation of a Schottky junction. Further interpretation of the electronic behavior of proteins was done by the calculation of selected solid-state parameters. For example, the turn-on voltage of hemoglobin was measured to occur at a lower turn-on voltage, indicating the possible influence of the hem group present as a cofactor in each subunit of this tetrameric protein.

Competing ultrafast photoinduced electron transfer and intersystem crossing of [Re(CO) 3 (Dmp)(His124)(Trp122)] + in Pseudomonas aeruginosa azurin: a nonadiabatic dynamics study

We present a computational study of sub-picosecond nonadiabatic dynamics in a rhenium complex coupled electronically to a tryptophan (Trp) side chain of Pseudomonas aeruginosa azurin, a prototypical protein used in the study of electron transfer in proteins. To gain a comprehensive understanding of the photoinduced processes in this system, we have carried out vertical excitation calculations at the TDDFT level of theory as well as nonadiabatic dynamics simulations using the surface hopping including arbitrary couplings (SHARC) method coupled to potential energy surfaces represented with a linear vibronic coupling model. The results show that the initial photoexcitation populates both singlet metal-to-ligand charge transfer (MLCT) and singlet charge-separated (CS) states, where in the latter an electron was transferred from the Trp amino acid to the complex. Subsequently, a complex mechanism of simultaneous intersystem crossing and electron transfer leads to the sub-picosecond population of triplet MLCT and triplet CS states. These results confirm the assignment of the sub-ps time constants of previous experimental studies and constitute the first computational evidence for the ultrafast formation of the charge-separated states in Re-sensitized azurin.

Ergodicity-Breaking in Thermal Biological Electron Transfer? Cytochrome C Revisited

It was recently suggested that certain redox proteins operate in an ergodicity-breaking regime to facilitate biological electron transfer (ET). A signature for this is a large variance reorganization free energy (several electronvolts) but a significantly smaller Stokes reorganization free energy due to incomplete protein relaxation on the time scale of the ET event. Here we investigate whether this picture holds for oxidation of cytochrome c in aqueous solution, at various levels of theory including classical molecular dynamics with two additive and one electronically polarizable force field, and QM/MM calculations with the QM region treated by full electrostatic DFT embedding and by the perturbed matrix method. Sampling the protein and energy gap dynamics over more than 250 ns, we find no evidence for ergodicity-breaking effects. In particular, the inclusion of electronic polarizability of the heme group at QM/MM levels did not induce nonergodic effects, contrary to previous reports by Matyushov et al. The well-known problem of overestimation of reorganization free energies with additive force fields is cured when the protein and solvent are treated as electronically polarizable. Ergodicity-breaking effects may occur in other redox proteins, and our results suggest that long simulations, ideally on the ET time scale, with electronically polarizable force fields are required to obtain strong numerical evidence for them.

Quantum biological tunnel junction for electron transfer imaging in live cells

Quantum biological electron transfer (ET) essentially involves in virtually all important biological processes such as photosynthesis, cellular respiration, DNA repair, cellular homeostasis, and cell death. However, there is no real-time imaging method to capture biological electron tunnelling in live cells to date. Here, we report a quantum biological electron tunnelling (QBET) junction and its application in real-time optical detection of QBET and the dynamics of ET in mitochondrial cytochrome c during cell life and death process. QBET junctions permit to see the behaviours of electron tunnelling through barrier molecules with different barrier widths. Using QBET spectroscopy, we optically capture real-time ET in cytochrome c redox dynamics during cellular apoptosis and necrosis in living cells. The non-invasive real-time QBET spectroscopic imaging of ET in live cell open a new era in life sciences and medicine by providing a way to capture spatiotemporal ET dynamics and to reveal the quantum biological mechanisms.

Although quantum biological electron transfer is important in many biological processes, imaging of the events in live cells has remained challenging. Here, the authors demonstrate real-time optical detection of quantum biological electron tunnelling between nanoparticles and cytochrome c inside living cells.

Towards direct enzyme wiring: a theoretical investigation of charge carrier transfer mechanisms between glucose oxidase and organic semiconductors

In this work, a general theoretical and numerical approach based on semiconductor theory, which could be applied to a study of direct enzyme wiring, has been discussed. Marcus-Hush theory was applied to evaluate the potential transfer of charge carriers (holes and electrons) between glucose oxidase (GOx) and organic semiconductors. Two mechanisms of multistep hopping of charge to/from the oxidised/reduced flavin-based moiety through residues of aromatic amino acids located in GOx and long range charge direct tunnelling from the cofactor to the organic semiconductor surface have been proposed and evaluated. It was determined that the hole-hopping mechanism is possible and proceeds at a low ionization potential of the organic semiconductor. The calculations reveal that hopping of electrons is blocked, but direct electron tunnelling between the cofactor and the organic semiconductor is still probable. The most optimal conditions and tunable characteristics of GOx-based biosensors such as the ionization potential, electron affinity of organic semiconductors and distance between the enzyme and surface were estimated for the first time.

Computational Design of Protein Function

The computational design of enzymes is a tremendous challenge for both chemistry and biochemistry. The ability to design stable and functional biocatalysts that could operate under different conditions to perform chemical reactions without precedent in nature, allowing the large-scale production of chemicals à la carte, would revolutionise both synthetic, pharmacologic and materials chemistry. Despite the great advances achieved, this highly multidisciplinary area of research is still in its infancy. This chapter describes the ‘inside-out’ protocol for computational enzyme design and both the achievements and limitations of the current technology are highlighted. Furthermore, molecular dynamics simulations have proved to be invaluable in the enzyme design process, constituting an important tool for discovering elusive catalytically relevant conformations of the engineered or designed enzyme. As a complement to the ‘inside-out’ design protocol, different examples where hybrid QM/MM approaches have been directly applied to discover beneficial mutations in rational computational enzyme design are described.

Genetically encoded conductive protein nanofibers secreted by engineered cells

Bacterial biofilms are promising tools for functional applications as bionanomaterials.

A theoretical multiscale treatment of protein-protein electron transfer: The ferredoxin/ferredoxin-NADP(+) reductase and flavodoxin/ferredoxin-NADP(+) reductase systems

In the photosynthetic electron transfer (ET) chain, two electrons transfer from photosystem I to the flavin-dependent ferredoxin-NADP(+) reductase (FNR) via two sequential independent ferredoxin (Fd) electron carriers. In some algae and cyanobacteria (as Anabaena), under low iron conditions, flavodoxin (Fld) replaces Fd as single electron carrier. Extensive mutational studies have characterized the protein-protein interaction in FNR/Fd and FNR/Fld complexes. Interestingly, even though Fd and Fld share the interaction site on FNR, individual residues on FNR do not participate to the same extent in the interaction with each of the protein partners, pointing to different electron transfer mechanisms. Despite of extensive mutational studies, only FNR/Fd X-ray structures from Anabaena and maize have been solved; structural data for FNR/Fld remains elusive. Here, we present a multiscale modelling approach including coarse-grained and all-atom protein-protein docking, the QM/MM e-Pathway analysis and electronic coupling calculations, allowing for a molecular and electronic comprehensive analysis of the ET process in both complexes. Our results, consistent with experimental mutational data, reveal the ET in FNR/Fd proceeding through a bridge-mediated mechanism in a dominant protein-protein complex, where transfer of the electron is facilitated by Fd loop-residues 40-49. In FNR/Fld, however, we observe a direct transfer between redox cofactors and less complex specificity than in Fd; more than one orientation in the encounter complex can be efficient in ET.

Progress and challenges in simulating and understanding electron transfer in proteins

This Review presents an overview of the most common numerical simulation approaches for the investigation of electron transfer (ET) in proteins. We try to highlight the merits of the different approaches but also the current limitations and challenges. The article is organized into three sections. Section 2 deals with direct simulation algorithms of charge migration in proteins. Section 3 summarizes the methods for testing the applicability of the Marcus theory for ET in proteins and for evaluating key thermodynamic quantities entering the reaction rates (reorganization energies and driving force). Recent studies interrogating the validity of the theory due to the presence of non-ergodic effects or of non-linear responses are also described. Section 4 focuses on the tunneling aspects of electron transfer. How can the electronic coupling between charge transfer states be evaluated by quantum chemistry approaches and rationalized? What interesting physics regarding the impact of protein dynamics on tunneling can be addressed? We will illustrate the different sections with examples taken from the literature to show what types of system are currently manageable with current methodologies. We also take care to recall what has been learned on the biophysics of ET within proteins thanks to the advent of atomistic simulations.