Long-Range Electron Transfer in Biomolecules. Tunneling or Hopping?

Molecular fragmentation as a way to reveal early electron dynamics induced by attosecond pulses

We present a theoretical study of the electron and nuclear dynamics that would arise in an attosecond two-color XUV-pump/XUV-probe experiment in glycine. In this scheme, the broadband pump pulse suddenly ionizes the molecule and creates an electronic wave packet that subsequently evolves under the influence of the nuclear motion until it is finally probed by the second XUV pulse. To describe the different steps of such an experiment, we have combined a multi-reference static-exchange scattering method with a trajectory surface hopping approach. We show that by changing the central frequency of the pump pulse, i.e., by engineering the initial electronic wave packet with the pump pulse, one can drive the cation dynamics into a specific fragmentation pathway. Reminiscence of this early electron dynamics can be observed in specific fragmentation channels (not all of them) as a function of the pump-probe delay and in time-resolved photoelectron spectra at specific photoelectron energies. The optimum conditions to visualize the initial electronic coherence in photoelectron and photo-ion spectra depend very much on the characteristics of the pump pulse as well as on the electronic structure of the molecule under study.

Enhancement of direct electron transfer in graphene bioelectrodes containing novel cytochrome c553 variants with optimized heme orientation

The highly efficient bioelectrodes based on single layer graphene (SLG) functionalized with pyrene self-assembled monolayer and novel cytochromec₅₅₃(cytc₅₅₃)peptide linker variants were rationally designed to optimize the direct electron transfer (DET) between SLG and the heme group of cyt. Through a combination of photoelectrochemical and quantum mechanical (QM/MM) approaches we show that the specific amino acid sequence of a short peptide genetically inserted between the cytc₅₅₃holoprotein and thesurface anchoring C-terminal His₆-tag plays a crucial role in ensuring the optimal orientation and distance of the heme group with respect to the SLG surface. Consequently, efficient DET occurring between graphene and cyt c₅₅₃ leads to a 20-fold enhancement of the cathodic photocurrent output compared to the previously reported devices of a similar type. The QM/MM modeling implies that a perpendicular or parallel orientation of the heme group with respect to the SLG surface is detrimental to DET, whereas the tilted orientation favors the cathodic photocurrent generation. Our work confirms the possibility of fine-tuning the electronic communication within complex bio-organic nanoarchitectures and interfaces due to optimization of the tilt angle of the heme group, its distance from the SLG surface and optimal HOMO/LUMO levels of the interacting redox centers.

Molecular mechanism of direct electron transfer in the robust cytochrome-functionalised graphene nanosystem

Molecular mechanism of DET between graphene and cytochrome c depends on the metal in the bio-organic interface: Co enhances the cathodic current via electron hopping from graphene to haem, whereas Ni exerts the opposite effect via tunnelling.

Charge-Transfer Knowledge Graph among Amino Acids Derived from High-Throughput Electronic Structure Calculations for Protein Database

The charge-transfer coupling is an important component in tight-binding methods. Because of the highly complex chemical structure of biomolecules, the anisotropic feature of charge-transfer couplings in realistic proteins cannot be ignored. In this work, we have performed the first large-scale quantitative assessment of charge-transfer preference by calculating the charge-transfer couplings in all 20 × 20 possible amino acid side-chain combinations, which are extracted from available high-quality structures of thousands of protein complexes. The charge-transfer database quantitatively shows distinct features of charge-transfer couplings among millions of amino acid side-chain combinations. The overall distribution of charge-transfer couplings reveals that only one average or representative structure cannot be regarded as the typical charge-transfer preference in realistic proteins. This work provides us an alternative route to comprehensively understand the charge-transfer couplings for the overall distribution of realistic proteins in the foreseen big data scenario.

Predicting and Interpreting the Structure of Type IV Pilus of Electricigens by Molecular Dynamics Simulations

Nanowires that transfer electrons to extracellular acceptors are important in organic matter degradation and nutrient cycling in the environment. Geobacter pili of the group of Type IV pilus are regarded as nanowire-like biological structures. However, determination of the structure of pili remains challenging due to the insolubility of monomers, presence of surface appendages, heterogeneity of the assembly, and low-resolution of electron microscopy techniques. Our previous study provided a method to predict structures for Type IV pili. In this work, we improved on our previous method using molecular dynamics simulations to optimize structures of Neisseria gonorrhoeae (GC), Neisseria meningitidis and Geobacter uraniireducens pilus. Comparison between the predicted structures for GC and Neisseria meningitidis pilus and their native structures revealed that proposed method could predict Type IV pilus successfully. According to the predicted structures, the structural basis for conductivity in G.uraniireducens pili was attributed to the three N-terminal aromatic amino acids. The aromatics were interspersed within the regions of charged amino acids, which may influence the configuration of the aromatic contacts and the rate of electron transfer. These results will supplement experimental research into the mechanism of long-rang electron transport along pili of electricigens.

Two Aromatic Rings Coupled a Sulfur-Containing Group to Favor Protein Electron Transfer by Instantaneous Formations of π∴S:π↔π:S∴π or π∴π:S↔π:π∴S Five-Electron Bindings

The cooperative interactions among two aromatic rings with a S-containing group are described, which may participate in electron hole transport in proteins. Ab initio calculations reveal the possibility for the formations of the π∴S:π↔π:S∴π and π∴π:S↔π:π∴S five-electron bindings in the corresponding microsurrounding structures in proteins, both facilitating electron hole transport as efficient relay stations. The relay functionality of these two special structures comes from their low local ionization energies and proper binding energies, which varies with the different aromatic amino acids, S-containing residues, and the arrangements of the same aromatic rings according to the local microsurroundings in proteins.

A New Type of Electron Relay Station in Proteins: Three-Piece S:Π∴S↔S∴Π:S Resonance Structure

A type of relay station for electron transfer in proteins, three-piece five-electron bonding, is introduced in this paper, which is also first proposed here. The ab initio calculations predict the formation of S:Π∴S↔S∴Π:S resonance binding with an aromatic ring located in the middle of two sulfur-containing groups, which may participate in electron-hole transport in proteins. These special structures can lower the local ionization energies to capture electron holes efficiently and may be easily formed and broken because of their proper binding energies. In addition, the UV-vis spectra provide evidence of the formations of the three-piece five-electron binding. The cooperation of three adjacent pieces may be advantage to promote electron transfer a longer distance.

Hole transfer kinetics of DNA

Not long after the discovery of the double-helical structure of DNA in 1952, researchers proposed that charge transfer along a one-dimensional π-array of nucleobases might be possible. At the end of the 1990s researchers discovered that a positive charge (a hole) generated in DNA migrates more than 200 Å along the structure, a discovery that ignited interest in the charge-transfer process in DNA. As a result, DNA became an interesting potential bottom-up material for constructing nanoelectronic sensors and devices because DNA can form various complex two-dimensional and three-dimensional structures, such as smiley faces and cubes. From the fundamental aspects of the hole transfer process, DNA is one of the most well-studied organic molecules with many reports on the synthesis of artificial nucleobase analogues. Thus, DNA offers a unique system to study how factors such as the HOMO energy and molecular flexibility affect hole transfer kinetics. Understanding the hole transfer mechanism requires a discussion of the hole transfer rate constants (kHT). This Account reviews the kHT values determined by our group and by Lewis and Wasielewski's group, obtained by a combination of the synthesis of modified DNA and time-resolved spectroscopy. DNA consists of G/C and A/T base pairs; the HOMO localizes on the purine bases G and A, and G has a lower oxidation potential and a higher energy HOMO. Typically, long-range hole transfer proceeded via sequential hole transfer between G/C's. The kinetics of this process in DNA sequences, including those with mismatches, is reproducible via kinetic modeling using the determined kHT for each hole transfer step between G/C's. We also determined the distance dependence parameter (β), which describes the steepness of the exponential decrease of kHT. Because of this value, >0.6 Å(-1) for hole transfer in DNA, DNA itself does not serve as a molecular wire. Interestingly, hole transfer proceeded exceptionally fast for some sequences in which G/C's are located close to each other, an observation that we cannot explain by a simple sequential hole transfer between G/C's but rather through hole delocalization over the nucleobases. To further investigate and refine the factors that affect kHT, we examined various artificial nucleobases. We clearly demonstrated that kHT depends strongly on the HOMO energy gap between the bases (ΔHOMO), and that kHT can be increased with decreasing ΔHOMO. We reduced ΔHOMO between the two type of base pairs by replacing adenines (A's) with deazaadenines ((z)A's) or diaminopurines (D's) and showed that the hole transfer rate through the G/C and A/T mix sequence increased by more than 3 orders of magnitude. We also investigated how DNA flexibility affects kHT. Locked nucleic acid (LNA) modification, which makes DNA more rigid, lowered kHT by more than 2 orders of magnitude. On the other hand, 5-Me-2'-deoxyzebularine (B) modification, which increases DNA flexibility, increased kHT by more than 1 order of magnitude. These new insights in hole transfer kinetics obtained from modified DNAs may aid in the design of new molecular-scale conducting materials.

Single-molecule studies of disulfide bond reduction pathways used by human thioredoxin

Disulfide bond reduction pathways used by human thioredoxin (hTrx) are studied at the single molecule level using a recombinant protein (I27SS)8. (I27SS)8 contains eight tandem repeats of identical immunoglobulin-like modules with one disulfide bond in each module. Single (I27SS)8 molecules are stretched at constant force applied by a cantilever in a force-clamp mode of atomic force microscopy (FC-AFM). Disulfide reduction events are accurately detected from stepwise increases in the end-to-end length of (I27SS)8. Earlier FC-AFM studies observed one disulfide reduction pathway used by hTrx and suggested an additional electron tunneling mechanism. Here, a very large set of unbiased FC-AFM data is collected in a range of clamping forces. By analyzing the data using exponential fits and dwell time histograms two disulfide reduction pathways used by hTrx are resolved. Based on previous studies one of these pathways is attributed to force-dependent Michaelis-Menten catalysis. The latter reduction pathway is weakly force-inhibited and occurs sporadically. Bimolecular nucleophilic substitutions (SN2) and electron tunneling (ET) mechanisms are discussed to explain the second pathway. Direct SN2 and ET mechanisms cannot be discounted, but a hypothetical E2-SN2 mechanism involving a hydride reducing a disulfide bond provides an interesting alternative, which needs to be verified in future experiments.

HOMO energy gap dependence of hole-transfer kinetics in DNA

DNA consists of two type of base-pairs, G-C and A-T, in which the highest occupied molecular orbital (HOMO) localizes on the purine bases G and A. While the hole transfer through consecutive Gs or As occurs faster than 10(9) s(-1), a significant drop in the hole transfer rate was observed for G-C and A-T mixed random sequences. In this study, by using various natural and artificial nucleobases having different HOMO levels, the effect of the HOMO-energy gap between bases (Δ(HOMO)) on the hole-transfer kinetics in DNA was investigated. The results demonstrated that the hole transfer rate can be increased by decreasing the Δ(HOMO) and can be finely tuned over 3 orders of magnitude by varying the Δ(HOMO).