Charge transfer in dynamical biosystems, or the treachery of (static) images

Bioimpedance formalism: A new approach for accessing the health status of cell and tissues

This manuscript describes a novel methodology for studying relaxation dynamics in tissues and cells using characteristic frequency of bioimpedance spectroscopy measurements. The Bioimpedance Formalism allows for the simultaneous study of bioelectrical parameters in the frequency and time domains, providing insight into possible relaxation processes occurring in the tissue or cell of interest. Results from the Cole-Cole analysis showed no multiple relaxation processes associated with heterogeneity, with a visible age group separation in males compared with females. The study of the relaxation dynamic in the time domain revealed that the β parameter can be used to analyse the charge carriers in tissues, cells, or cancer cells, potentially leading to new diagnostic and therapeutic approaches for cancer and other diseases. Overall, this approach presents a promising area of research for gaining insights into the electrical properties of tissues and cells using bioimpedance methods.

How Rigid Are Anthranilamide Molecular Electrets?

As important as molecular electrets are for electronic materials and devices, conformational fluctuations strongly impact their macrodipoles and intrinsic properties. Herein, we employ molecular dynamics (MD) simulations with the polarizable charge equilibrium (PQEq) method to investigate the persistence length (LP) of molecular electrets composed of anthranilamide (Aa) residues. The PQEq-MD dissipates the accepted static notions about Aa macromolecules, and LP represents the shortest Aa rigid segments. The classical model with a single LP value does not describe these oligomers. Introducing multiple LP values for the same macromolecule follows the observed trends and discerns the enhanced rigidity in their middle sections from the reduced stiffness at their terminal regions. Furthermore, LP distinctly depends on solvent polarity. The Aa oligomers maintain extended conformations in nonpolar solvents with LP exceeding 4 nm, while in polar media, increased conformational fluctuations reduce LP to about 2 nm. These characteristics set key guidelines about the utility of Aa conjugates for charge-transfer systems within organic electronics and energy engineering.

Coherence in Chemistry: Foundations and Frontiers

Coherence refers to correlations in waves. Because matter has a wave-particle nature, it is unsurprising that coherence has deep connections with the most contemporary issues in chemistry research (e.g., energy harvesting, femtosecond spectroscopy, molecular qubits and more). But what does the word "coherence" really mean in the context of molecules and other quantum systems? We provide a review of key concepts, definitions, and methodologies, surrounding coherence phenomena in chemistry, and we describe how the terms "coherence" and "quantum coherence" refer to many different phenomena in chemistry. Moreover, we show how these notions are related to the concept of an interference pattern. Coherence phenomena are indeed complex, and ambiguous definitions may spawn confusion. By describing the many definitions and contexts for coherence in the molecular sciences, we aim to enhance understanding and communication in this broad and active area of chemistry.

Efficient simulation of open quantum systems coupled to a reservoir through multiple channels

It is challenging to simulate open quantum systems that are connected to a reservoir through multiple channels. For example, vibrations may induce fluctuations in both energy gaps and electronic couplings, which represent two independent channels of system-bath couplings. Systems of this kind are ubiquitous in the processes of excited state radiationless decay. Combined with density matrix renormalization group (DMRG) and matrix product states (MPS) methods, we develop an interaction-picture chain mapping strategy for vibrational reservoirs to simulate the dynamics of these open systems, resulting in time-dependent spatially local system-bath couplings in the chain-mapped Hamiltonian. This transformation causes the entanglement generated by the system-bath interactions to be restricted within a narrow frequency window of vibrational modes, enabling efficient DMRG/MPS dynamical simulations. We demonstrate the utility of this approach by simulating singlet fission dynamics using a generalized spin-boson Hamiltonian with both diagonal and off-diagonal system-bath couplings. This approach generalizes an earlier interaction-picture chain mapping scheme, allowing for efficient and exact simulation of systems with multi-channel system-bath couplings using matrix product states, which may further our understanding of nonlocal exciton-phonon couplings in exciton transport and the non-Condon effect in energy and electron transfer.

Undulating Free Energy Landscapes Buffer Redox Chains from Environmental Fluctuations

Roller-coaster or undulating free energy landscapes, with alternating high and low potential cofactors, occur frequently in biological redox chains. Yet, there is little understanding of the possible advantages created by these landscapes. We examined the tetraheme subunit associated with Blastochloris viridis reaction centers, comparing the dynamics of the native protein and of hypothetical (in silico) mutants. We computed the variation in the total number of electrons in wild type (WT) and mutant tetrahemes connected to an electron reservoir in the presence of a time-varying potential, as a model for a fluctuating redox environment. We found that roller-coaster free energy landscapes buffer the redox cofactor populations from these fluctuations. The WT roller-coaster landscape slows forward and backward electron transfer in the face of fluctuations, and may offer the advantage of sustaining the reduction of essential cofactors, such as the chlorophyll special pair in photosynthesis, even though an undulating landscape introduces thermodynamically uphill steps in multistep redox chains.

Long-distance electron transport in multicellular freshwater cable bacteria

Filamentous multicellular cable bacteria perform centimeter-scale electron transport in a process that couples oxidation of an electron donor (sulfide) in deeper sediment to the reduction of an electron acceptor (oxygen or nitrate) near the surface. While this electric metabolism is prevalent in both marine and freshwater sediments, detailed electronic measurements of the conductivity previously focused on the marine cable bacteria (Candidatus Electrothrix), rather than freshwater cable bacteria, which form a separate genus (Candidatus Electronema) and contribute essential geochemical roles in freshwater sediments. Here, we characterize the electron transport characteristics of Ca. Electronema cable bacteria from Southern California freshwater sediments. Current-voltage measurements of intact cable filaments bridging interdigitated electrodes confirmed their persistent conductivity under a controlled atmosphere and the variable sensitivity of this conduction to air exposure. Electrostatic and conductive atomic force microscopies mapped out the characteristics of the cell envelope's nanofiber network, implicating it as the conductive pathway in a manner consistent with previous findings in marine cable bacteria. Four-probe measurements of microelectrodes addressing intact cables demonstrated nanoampere currents up to 200 μm lengths at modest driving voltages, allowing us to quantify the nanofiber conductivity at 0.1 S/cm for freshwater cable bacteria filaments under our measurement conditions. Such a high conductivity can support the remarkable sulfide-to-oxygen electrical currents mediated by cable bacteria in sediments. These measurements expand the knowledgebase of long-distance electron transport to the freshwater niche while shedding light on the underlying conductive network of cable bacteria.

Intra-strand phosphate-mediated pathways in microsolvated double-stranded DNA

We argue that dry DNA charge transport in molecular junctions, over distances of tens of nanometers, can take place via independent intra-strand pathways involving the phosphate groups. Such pathways explain recent single-molecule experiments that compare currents in intact and nicked 100 base-pair double-stranded DNA. We explore the conditions that favor independent intra-strand transport channels with the participation of the phosphate groups, as opposed to purely base-mediated transport involving the pi-stacked bases and inter-strand transitions. Our computations demonstrate how long-distance transport pathways in DNA are tuned by the degree of solvation, which affects the level of dynamic disorder in the pi-stacking, and the energies of phosphate-group molecular orbitals.

Excited-State Charge Transfer Coupling from Quasiparticle Energy Density Functional Theory

The recently developed Quasiparticle Energy (QE) scheme, based on a DFT calculation with one more (or less) electron, offers a good description of excitation energies, even with charge transfer characters. In this work, QE is further extended to calculate electron transfer (ET) couplings involving two excited states. We tested it with a donor-acceptor complex, consisting of a furan and a 1,1-dicyanoethylene (DCNE), in which two low lying charge transfer and local excitation states are involved. With generalized Mülliken-Hush and fragment charge-difference schemes, couplings from the QE approach generally agree well with those obtained from TDDFT, except that QE couplings exhibit better exponential distance dependence. Couplings from half-energy gaps with an external field are also calculated and reported. Our results show that the QE scheme is robust in calculating ET couplings with greatly reduced computational time.

Unified Entropy-Ruled Einstein's Relation for Bulk and Low-Dimensional Molecular-Material Systems: A Hopping-to-Band Shift Paradigm

We present a unified paradigm on entropy-ruled Einstein's diffusion-mobility relation (μ/D ratio) for 1D, 2D, and 3D free-electron solid state systems. The localization transport in the extended molecules is well approximated by the continuum time-delayed hopping factor within our unified entropy-ruled transport method of noninteracting quantum systems. Moreover, we generalize an entropy-dependent diffusion relation for 1D, 2D, and 3D systems as defined by D d , h e f f = D d , h e f f = 0 ⁡ exp ( ( d - 1 ) h e f f d + 2 ) , where heff and d are the effective entropy and dimension (d = 1, 2, 3), respectively. This generalized relation is valid for both equilibrium and nonequilibrium transport systems since the parameter heff is closely connected with the nonequilibrium fluctuation theorem-based entropy production rule. Importantly, we herein revisit the Boltzmann approach using an entropy-ruled method for mobility calculation for the universal quantum materials that is expressed as μ d = [ ( d d + 2 ) q d h e f f d η ] v F 2 τ 2 , where v F 2 τ 2 is the diffusion constant for band transport systems and η is the chemical potential. According to our entropy-ruled μ/D relation, the Navamani-Shockley diode equation is transformed.

An Investigation of the Influence of Tyrosine Local Interactions on Electron Hopping in a Model Protein

Multi-step electron transfer reactions are important to the function of many cellular systems. The ways in which such systems have evolved to direct electrons along specific pathways are largely understood, but less so are the ways in which the reduction-oxidation potentials of individual redox sites are controlled. We prepared a series of three new artificial variants of Pseudomonas aeruginosa azurin where a tyrosine (Tyr109) is situated between the native Cu ion and a Ru(II) photosensitizer tethered to a histidine (His107). Arginine, glutamine, or methionine were introduced as position 122, which is near to Tyr109. We investigated the rate of CuI oxidation by a flash-quench generated Ru(III) oxidant over pH values from 5 to 9. While the identity of the residue at position 122 affects some of the physical properties of Tyr109, the rates of CuI oxidation are only weakly dependent on the identity of the residue at 122. The results highlight that more work is still needed to understand how non-covalent interactions of redox active groups are affected in redox proteins.

The energetics and evolution of oxidoreductases in deep time

The core metabolic reactions of life drive electrons through a class of redox protein enzymes, the oxidoreductases. The energetics of electron flow is determined by the redox potentials of organic and inorganic cofactors as tuned by the protein environment. Understanding how protein structure affects oxidation-reduction energetics is crucial for studying metabolism, creating bioelectronic systems, and tracing the history of biological energy utilization on Earth. We constructed ProtReDox (https://protein-redox-potential.web.app), a manually curated database of experimentally determined redox potentials. With over 500 measurements, we can begin to identify how proteins modulate oxidation-reduction energetics across the tree of life. By mapping redox potentials onto networks of oxidoreductase fold evolution, we can infer the evolution of electron transfer energetics over deep time. ProtReDox is designed to include user-contributed submissions with the intention of making it a valuable resource for researchers in this field.

A stronger acceptor decreases the rates of charge transfer: ultrafast dynamics and on/off switching of charge separation in organometallic donor-bridge-acceptor systems

To unravel the role of driving force and structural changes in directing the photoinduced pathways in donor-bridge-acceptor (DBA) systems, we compared the ultrafast dynamics in novel DBAs which share a phenothiazine (PTZ) electron donor and a Pt(ii) trans-acetylide bridge (-C[triple bond, length as m-dash]C-Pt-C[triple bond, length as m-dash]C-), but bear different acceptors conjugated into the bridge (naphthalene-diimide, NDI; or naphthalene-monoimide, NAP). The excited state dynamics were elucidated by transient absorption, time-resolved infrared (TRIR, directly following electron density changes on the bridge/acceptor), and broadband fluorescence-upconversion (FLUP, directly following sub-picosecond intersystem crossing) spectroscopies, supported by TDDFT calculations. Direct conjugation of a strong acceptor into the bridge leads to switching of the lowest excited state from the intraligand 3IL state to the desired charge-separated 3CSS state. We observe two surprising effects of an increased strength of the acceptor in NDI vs. NAP: a ca. 70-fold slow-down of the 3CSS formation-(971 ps)-1vs. (14 ps)-1, and a longer lifetime of the 3CSS (5.9 vs. 1 ns); these are attributed to differences in the driving force ΔGet, and to distance dependence. The 100-fold increase in the rate of intersystem crossing-to sub-500 fs-by the stronger acceptor highlights the role of delocalisation across the heavy-atom containing bridge in this process. The close proximity of several excited states allows one to control the yield of 3CSS from ∼100% to 0% by solvent polarity. The new DBAs offer a versatile platform for investigating the role of bridge vibrations as a tool to control excited state dynamics.

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.

Control of charge transport in electronically active systems towards integrated biomolecular circuits (IbC)

The miniaturization of traditional silicon-based electronics will soon reach its limitation as quantum tunneling and heat become serious problems at the several-nanometer scale. Crafting integrated circuits via self-assembly of electronically active molecules using a "bottom-up" paradigm provides a potential solution to these technological challenges. In particular, integrated biomolecular circuits (IbC) offer promising advantages to achieve this goal, as nature offers countless examples of functionalities entailed by self-assembly and examples of controlling charge transport at the molecular level within the self-assembled structures. To this end, the review summarizes the progress in understanding how charge transport is regulated in biosystems and the key redox-active amino acids that enable the charge transport. In addition, charge transport mechanisms at different length scales are also reviewed, offering key insights for controlling charge transport in IbC in the future.

Single-molecular protein-based bioelectronics via electronic transport: fundamentals, devices and applications

Biomolecular electronics is a rapidly growing multidisciplinary field that combines biology, nanoscience, and engineering to bridge the two important fields of life sciences and molecular electronics. Proteins are remarkable for their ability to recognize molecules and transport electrons, making the integration of proteins into electronic devices a long sought-after goal and leading to the emergence of the field of protein-based bioelectronics, also known as proteotronics. This field seeks to design and create new biomolecular electronic platforms that allow for the understanding and manipulation of protein-mediated electronic charge transport and related functional applications. In recent decades, there have been numerous reports on protein-based bioelectronics using a variety of nano-gapped electrical devices and techniques at the single molecular level, which are not achievable with conventional ensemble approaches. This review focuses on recent advances in physical electron transport mechanisms, device fabrication methodologies, and various applications in protein-based bioelectronics. We discuss the most recent progress of the single or few protein-bridged electrical junction fabrication strategies, summarise the work on fundamental and functional applications of protein bioelectronics that enable high and dynamic electron transport, and highlight future perspectives and challenges that still need to be addressed. We believe that this specific review will stimulate the interdisciplinary research of topics related to protein-related bioelectronics, and open up new possibilities for single-molecule biophysics and biomedicine.

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.

Modified Fermi's golden rule rate expressions

Fermi's golden rule (FGR) serves as the basis for many expressions of spectroscopic observables and quantum transition rates. The utility of FGR has been demonstrated through decades of experimental confirmation. However, there still remain important cases where the evaluation of a FGR rate is ambiguous or ill-defined. Examples are cases where the rate has divergent terms due to the sparsity in the density of final states or time dependent fluctuations of system Hamiltonians. Strictly speaking, assumptions of FGR are no longer valid for such cases. However, it is still possible to define modified FGR rate expressions that are useful as effective rates. The resulting modified FGR rate expressions resolve a long standing ambiguity often encountered in using FGR and offer more reliable ways to model general rate processes. Simple model calculations illustrate the utility and implications of new rate expressions.

Target State Optimized Density Functional Theory for Electronic Excited and Diabatic States

A flexible self-consistent field method, called target state optimization (TSO), is presented for exploring electronic excited configurations and localized diabatic states. The key idea is to partition molecular orbitals into different subspaces according to the excitation or localization pattern for a target state. Because of the orbital-subspace constraint, orbitals belonging to different subspaces do not mix. Furthermore, the determinant wave function for such excited or diabatic configurations can be variationally optimized as a ground state procedure, unlike conventional ΔSCF methods, without the possibility of collapsing back to the ground state or other lower-energy configurations. The TSO method can be applied both in Hartree-Fock theory and in Kohn-Sham density functional theory (DFT). The density projection procedure and the working equations for implementing the TSO method are described along with several illustrative applications. For valence excited states of organic compounds, it was found that the computed excitation energies from TSO-DFT and time-dependent density functional theory (TD-DFT) are of similar quality with average errors of 0.5 and 0.4 eV, respectively. For core excitation, doubly excited states and charge-transfer states, the performance of TSO-DFT is clearly superior to that from conventional TD-DFT calculations. It is shown that variationally optimized charge-localized diabatic states can be defined using TSO-DFT in energy decomposition analysis to gain both qualitative and quantitative insights on intermolecular interactions. Alternatively, the variational diabatic states may be used in molecular dynamics simulation of charge transfer processes. The TSO method can also be used to define basis states in multistate density functional theory for excited states through nonorthogonal state interaction calculations. The software implementing TSO-DFT can be accessed from the authors.

Supramolecular Enhancement of Charge Transport through Pillar[5]arene-Based Self-Assembled Monolayers

Intermolecular charge transport is one of the essential modes for modulating charge transport in molecular electronic devices. Supermolecules are highly promising candidates for molecular devices because of their abundant structures and easy functionalization. Herein, we report an efficient strategy to enhance charge transport through pillar[5]arene self-assembled monolayers (SAMs) by introducing cationic guests. The current density of pillar[5]arene SAMs can be raised up to about 2.1 orders of magnitude by inserting cationic molecules into the cavity of pillar[5]arenes in SAMs. Importantly, we have also observed a positive correlation between the charge transport of pillar[5]arene-based complex SAMs and the binding affinities of the pillar[5]arene-based complexation. Such an enhancement of charge transport is attributed to the efficient host-guest interactions that stabilize the supramolecular complexes and lower the energy gaps for charge transport. This work provides a predictive pattern for the regulation of intermolecular charge transport in guiding the design of next generation switches and functional sensors in supramolecular electronics.

A Method to Investigate the Mechanism of Charge Transport Across Bio-Molecular Junctions with Ferritin

To investigate the mechanisms of charge transport (CT) across biomolecular tunnel junctions, it is required to make electrical contacts by a non-invasive method that leaves the biomolecules unaltered. Although different methods to form biomolecular junctions are available, here we describe the EGaIn-method because it allows us to readily form electrical contacts to monolayers of biomolecules in ordinary laboratory settings and to probe CT as a function of voltage, temperature, or magnetic field. This method relies on a non-Newtonian liquid-metal ally of Ga and In with a few nm thin layer of GaOx floating on its surface giving this material non-Newtonian properties allowing it to be shaped in to cone-shaped tips or stabilized in microchannels. These EGaIn structures form stable contacts to monolayers making it possible to investigate CT mechanisms across biomolecules in great detail.

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.

Generalization on Entropy-Ruled Charge and Energy Transport for Organic Solids and Biomolecular Aggregates

Herein, a generalized version of the entropy-ruled charge and energy transport mechanism for organic solids and biomolecular aggregates is presented. The effects of thermal disorder and electric field on electronic transport in molecular solids have been quantified by entropy, which eventually varies with respect to the typical disorder (static or dynamic). Based on our previous differential entropy (h s )-driven charge transport method, we explore the nonsteady carrier energy flux principle for soft matter systems from small organic solids to macrobiomolecular aggregates. Through this principle, the synergic nature of charge and energy transport in different organic systems is addressed. In this work, entropy is the key parameter to classify whether the carrier dynamics is in a nonsteady or steady state. Besides that, we also propose the formulation for unifying the hopping and band transport, which provides the relaxation time-hopping rate relation and the relaxation time-effective mass ratio. The calculated disorder drift time (or entropy-weighted carrier drift time) for hole transport in an alkyl-substituted triphenylamine (TPA) molecular device is 9.3 × 10-7 s, which illustrates nuclear dynamics-coupled charge transfer kinetics. The existence of nonequilibrium transport is anticipated while the carrier dynamics is in the nonsteady state, which is further examined from the rate of traversing potential in octupolar molecules. Our entropy-ruled Einstein model connects the adiabatic band and nonadiabatic hopping transport mechanisms. The logarithmic current density at different electric field-assisted site energy differences provides information about the typical transport (whether trap-free diffusion or trap-assisted recombination) in molecular devices, which reflects in the Navamani-Shockley diode equation.

Cofactor Dynamics Couples the Protein Surface to the Heme in Cytochrome c, Facilitating Electron Transfer

Electron transport through biomolecules and in biological transport networks is of great importance to bioenergetics and biocatalysis. More generally, it is of crucial importance to understand how the pathways that connect buried metallocofactors to other cofactors, and to protein surfaces, affect the biological chemistry of metalloproteins. In terms of electron transfer (ET), the strongest coupling pathways usually comprise covalent and hydrogen bonded networks, with a limited number of through-space contacts. Herein, we set out to determine the relative roles of hydrogen bonds involved in ET via an established heme-to-surface tunneling pathway in cytochrome (cyt) c (i.e., heme-W59-D60-E61-N62). A series of cyt c variants were produced where a ruthenium tris(diimine) photooxidant was placed at position 62 via covalent modification of the N62C residue. Surprisingly, variants where the H-bonding residues W59 and D60 were replaced (i.e., W59F and D60A) showed no change in ET rate from the ferrous heme to Ru(III). In contrast, changing the composition of an alternative tunneling pathway (i.e., heme-M64-N63-C62) with the M64L substitution shows a factor of 2 decrease in the rate of heme-to-Ru ET. This pathway involves a through-space tunneling step between the heme and M64 residue, and such steps are usually disfavored. To rationalize why the heme-M64-N63-C62 is preferred, molecular dynamics (MD) simulations and Pathways analysis were employed. These simulations show that the change in heme-Ru ET rates is attributed to different conformations with compressed donor-acceptor distances, by ∼2 Å in pathway distance, in the M64-containing protein as compared to the M64L protein. The change in distance is correlated with changes in the electronic coupling that are in accord with the experimentally observed heme-Ru ET rates. Remarkably, the M64L variation at the core of the protein translates to changes in cofactor dynamics at the protein surface. The surface changes identified by MD simulations include dynamic anion-π and dipole-dipole interactions. These interactions influence the strength of tunneling pathways and ET rates by facilitating decreases in through-space tunneling distances in key coupling pathways.

Torsionally broken symmetry assists infrared excitation of biomimetic charge-coupled nuclear motions in the electronic ground state

The concerted interplay between reactive nuclear and electronic motions in molecules actuates chemistry. Here, we demonstrate that out-of-plane torsional deformation and vibrational excitation of stretching motions in the electronic ground...

Chiral Induced Spin Selectivity and Its Implications for Biological Functions

Chirality in life has been preserved throughout evolution. It has been assumed that the main function of chirality is its contribution to structural properties. In the past two decades, however, it has been established that chiral molecules possess unique electronic properties. Electrons that pass through chiral molecules, or even charge displacements within a chiral molecule, do so in a manner that depends on the electron's spin and the molecule's enantiomeric form. This effect, referred to as chiral induced spin selectivity (CISS), has several important implications for the properties of biosystems. Among these implications, CISS facilitates long-range electron transfer, enhances bio-affinities and enantioselectivity, and enables efficient and selective multi-electron redox processes. In this article, we review the CISS effect and some of its manifestations in biological systems. We argue that chirality is preserved so persistently in biology not only because of its structural effect, but also because of its important function in spin polarizing electrons. Expected final online publication date for the Annual Review of Biophysics, Volume 51 is May 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Intramolecular Charge Transfer of 1-Aminoanthraquinone and Ultrafast Solvation Dynamics of Dimethylsulfoxide

The intramolecular charge transfer (ICT) of 1-aminoanthraquinone (AAQ) in the excited state strongly depends on its solvent properties, and the twisted geometry of its amino group has been recommended for the twisted ICT (TICT) state by recent theoretical works. We report the transient Raman spectra of AAQ in a dimethylsulfoxide (DMSO) solution by femtosecond stimulated Raman spectroscopy to provide clear experimental evidence for the TICT state of AAQ. The ultrafast (~110 fs) TICT dynamics of AAQ were observed from the major vibrational modes of AAQ including the ν_C-N + δ_CH and ν_C=O modes. The coherent oscillations in the vibrational bands of AAQ strongly coupled to the nuclear coordinate for the TICT process have been observed, which showed its anharmonic coupling to the low frequency out of the plane deformation modes. The vibrational mode of solvent DMSO, ν_S=O showed a decrease in intensity, especially in the hydrogen-bonded species of DMSO, which clearly shows that the solvation dynamics of DMSO, including hydrogen bonding, are crucial to understanding the reaction dynamics of AAQ with the ultrafast structural changes accompanying the TICT.

Photoinduced hole hopping through tryptophans in proteins

Hole hopping through tryptophan/tyrosine chains enables rapid unidirectional charge transport over long distances. We have elucidated structural and dynamical factors controlling hopping speed and efficiency in two modified azurin constructs that include a rhenium(I) sensitizer, Re(His)(CO)₃(dmp)⁺, and one or two tryptophans (W₁, W₂). Experimental kinetics investigations showed that the two closely spaced (3 to 4 Å) intervening tryptophans dramatically accelerated long-range electron transfer (ET) from Cu^I to the photoexcited sensitizer. In our theoretical work, we found that time-dependent density-functional theory (TDDFT) quantum mechanics/molecular mechanics/molecular dynamics (QM/MM/MD) trajectories of low-lying triplet excited states of Re^I(His)(CO)₃(dmp)⁺-W₁(-W₂) exhibited crossings between sensitizer-localized (*Re) and charge-separated [Re^I(His)(CO)₃(dmp^•-)/(W₁ ^•+ or W₂ ^•+)] (CS1 or CS2) states. Our analysis revealed that the distances, angles, and mutual orientations of ET-active cofactors fluctuate in a relatively narrow range in which the cofactors are strongly coupled, enabling adiabatic ET. Water-dominated electrostatic field fluctuations bring *Re and CS1 states to a crossing where *Re(CO)₃(dmp)⁺←W₁ ET occurs, and CS1 becomes the lowest triplet state. ET is promoted by solvation dynamics around *Re(CO)₃(dmp)⁺(W₁); and CS1 is stabilized by Re(dmp^•-)/W₁ ^•+ electron/hole interaction and enhanced W₁ ^•+ solvation. The second hop, W₁ ^•+←W₂, is facilitated by water fluctuations near the W₁/W₂ unit, taking place when the electrostatic potential at W₂ drops well below that at W₁ ^•+ Insufficient solvation and reorganization around W₂ make W₁ ^•+←W₂ ET endergonic, shifting the equilibrium toward W₁ ^•+ and decreasing the charge-separation yield. We suggest that multiscale TDDFT/MM/MD is a suitable technique to model the simultaneous evolution of photogenerated excited-state manifolds.

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.

TDRD3 promotes DHX9 chromatin recruitment and R-loop resolution

R-loops, which consist of a DNA/RNA hybrid and a displaced single-stranded DNA (ssDNA), are increasingly recognized as critical regulators of chromatin biology. R-loops are particularly enriched at gene promoters, where they play important roles in regulating gene expression. However, the molecular mechanisms that control promoter-associated R-loops remain unclear. The epigenetic ‘reader’ Tudor domain-containing protein 3 (TDRD3), which recognizes methylarginine marks on histones and on the C-terminal domain of RNA polymerase II, was previously shown to recruit DNA topoisomerase 3B (TOP3B) to relax negatively supercoiled DNA and prevent R-loop formation. Here, we further characterize the function of TDRD3 in R-loop metabolism and introduce the DExH-box helicase 9 (DHX9) as a novel interaction partner of the TDRD3/TOP3B complex. TDRD3 directly interacts with DHX9 via its Tudor domain. This interaction is important for recruiting DHX9 to target gene promoters, where it resolves R-loops in a helicase activity-dependent manner to facilitate gene expression. Additionally, TDRD3 also stimulates the helicase activity of DHX9. This stimulation relies on the OB-fold of TDRD3, which likely binds the ssDNA in the R-loop structure. Thus, DHX9 functions together with TOP3B to suppress promoter-associated R-loops. Collectively, these findings reveal new functions of TDRD3 and provide important mechanistic insights into the regulation of R-loop metabolism.

Why Do Most Aromatics Fail to Support Hole Hopping in the Cytochrome c Peroxidase-Cytochrome c Complex?

Electron transport through aromatic species (especially tryptophan and tyrosine) plays a central role in water splitting, redox signaling, oxidative damage protection, and bioenergetics. The cytochrome c peroxidase (CcP)-cytochrome c (Cc) complex (CcP:Cc) is used widely to study interprotein electron transfer (ET) mechanisms. Tryptophan 191 (Trp191) of CcP supports hole hopping charge recombination in the CcP:Cc complex. Experimental studies find that when Trp191 is substituted by tyrosine, phenylalanine, or redox-active aniline derivatives bound in the W191G cavity, enzymatic activity and charge recombination rates both decrease. Theoretical analysis of these CcP:Cc complexes finds that the ET kinetics depend strongly on the chemistry of the modified Trp site. The computed electronic couplings in the W191F and W191G species are orders of magnitude smaller than in the native protein, due largely to the absence of a hopping intermediate and the large tunneling distance. Small molecules bound in the W191G cavity are weakly coupled electronically to the Cc heme, and the structural disorder of the guest molecule in the binding pocket may contribute further to the lack of enzymatic activity. The couplings in W191Y are not substantially weakened compared to the native species, but the redox potential difference for tyrosine vs tryptophan oxidation accounts for the slower rate in the Tyr mutant. Thus, theoretical analysis explains why only the native Trp supports rapid hole hopping in the CcP:Cc complex. Favorable free energies and electronic couplings are essential for establishing an efficient hole hopping relay in this protein-protein complex.

Oxalate decarboxylase uses electron hole hopping for catalysis

The hexameric low-pH stress response enzyme oxalate decarboxylase catalyzes the decarboxylation of the oxalate mono-anion in the soil bacterium Bacillus subtilis. A single protein subunit contains two Mn-binding cupin domains, and catalysis depends on Mn(III) at the N-terminal site. The present study suggests a mechanistic function for the C-terminal Mn as an electron hole donor for the N-terminal Mn. The resulting spatial separation of the radical intermediates directs the chemistry toward decarboxylation of the substrate. A π-stacked tryptophan pair (W96/W274) links two neighboring protein subunits together, thus reducing the Mn-to-Mn distance from 25.9 Å (intrasubunit) to 21.5 Å (intersubunit). Here, we used theoretical analysis of electron hole-hopping paths through redox-active sites in the enzyme combined with site-directed mutagenesis and X-ray crystallography to demonstrate that this tryptophan pair supports effective electron hole hopping between the C-terminal Mn of one subunit and the N-terminal Mn of the other subunit through two short hops of ∼8.5 Å. Replacement of W96, W274, or both with phenylalanine led to a large reduction in catalytic efficiency, whereas replacement with tyrosine led to recovery of most of this activity. W96F and W96Y mutants share the wildtype tertiary structure. Two additional hole-hopping networks were identified leading from the Mn ions to the protein surface, potentially protecting the enzyme from high Mn oxidation states during turnover. Our findings strongly suggest that multistep hole-hopping transport between the two Mn ions is required for enzymatic function, adding to the growing examples of proteins that employ aromatic residues as hopping stations.

A single atom change turns insulating saturated wires into molecular conductors

We present an efficient strategy to modulate tunnelling in molecular junctions by changing the tunnelling decay coefficient, β, by terminal-atom substitution which avoids altering the molecular backbone. By varying X = H, F, Cl, Br, I in junctions with S(CH₂)_(10-18)X, current densities (J) increase >4 orders of magnitude, creating molecular conductors via reduction of β from 0.75 to 0.25 Å⁻¹. Impedance measurements show tripled dielectric constants (ε_r) with X = I, reduced HOMO-LUMO gaps and tunnelling-barrier heights, and 5-times reduced contact resistance. These effects alone cannot explain the large change in β. Density-functional theory shows highly localized, X-dependent potential drops at the S(CH₂)_nX//electrode interface that modifies the tunnelling barrier shape. Commonly-used tunnelling models neglect localized potential drops and changes in ε_r. Here, we demonstrate experimentally that \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\beta \propto 1/\sqrt{{\varepsilon }_{r}}$$\end{document}β∝1/εr, suggesting highly-polarizable terminal-atoms act as charge traps and highlighting the need for new charge transport models that account for dielectric effects in molecular tunnelling junctions.

In molecular junctions, where a molecule is placed between two electrodes, the current passed decays exponentially as a function of length. Here, Chen et al. show that this exponentially attenuation can be controlled by changing a single atom at the end of the molecular wire.

Mechanistic Insights on Heme-to-Heme Transmembrane Electron Transfer Within NADPH Oxydases From Atomistic Simulations

NOX5 is a member of the NADPH oxidase family which is dedicated to the production of reactive oxygen species. The molecular mechanisms governing transmembrane electron transfer (ET) that permits to shuttle electrons over the biological membrane have remained elusive for a long time. Using computer simulations, we report conformational dynamics of NOX5 embedded within a realistic membrane environment. We assess the stability of the protein within the membrane and monitor the existence of cavities that could accommodate dioxygen molecules. We investigate the heme-to-heme electron transfer. We find a reaction free energy of a few tenths of eV (ca. −0.3 eV) and a reorganization free energy of around 1.1 eV (0.8 eV after including electrostatic induction corrections). The former indicates thermodynamically favorable ET, while the latter falls in the expected values for transmembrane inter-heme ET. We estimate the electronic coupling to fall in the range of the μeV. We identify electron tunneling pathways showing that not only the W378 residue is playing a central role, but also F348. Finally, we reveal the existence of two connected O₂₋binding pockets near the outer heme with fast exchange between the two sites on the nanosecond timescale. We show that when the terminal heme is reduced, O₂ binds closer to it, affording a more efficient tunneling pathway than when the terminal heme is oxidized, thereby providing an efficient mechanism to catalyze superoxide production in the final step. Overall, our study reveals some key molecular mechanisms permitting reactive oxygen species production by NOX5 and paves the road for further investigation of ET processes in the wide family of NADPH oxidases by computer simulations.

Delocalization-Assisted Transport through Nucleic Acids in Molecular Junctions

The flow of charge through molecules is central to the function of supramolecular machines, and charge transport in nucleic acids is implicated in molecular signaling and DNA repair. We examine the transport of electrons through nucleic acids to understand the interplay of resonant and nonresonant charge carrier transport mechanisms. This study reports STM break junction measurements of peptide nucleic acids (PNAs) with a G-block structure and contrasts the findings with previous results for DNA duplexes. The conductance of G-block PNA duplexes is much higher than that of the corresponding DNA duplexes of the same sequence; however, they do not display the strong even-odd dependence conductance oscillations found in G-block DNA. Theoretical analysis finds that the conductance oscillation magnitude in PNA is suppressed because of the increased level of electronic coupling interaction between G-blocks in PNA and the stronger PNA-electrode interaction compared to that in DNA duplexes. The strong interactions in the G-block PNA duplexes produce molecular conductances as high as 3% G0, where G0 is the quantum of conductance, for 5 nm duplexes.

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.

On the Search of a Silver Bullet for the Preparation of Bioinspired Molecular Electrets with Propensity to Transfer Holes at High Potentials

Biological structure-function relationships offer incomparable paradigms for charge-transfer (CT) science and its implementation in solar-energy engineering, organic electronics, and photonics. Electrets are systems with co-directionally oriented electric dopes with immense importance for CT science, and bioinspired molecular electrets are polyamides of anthranilic-acid derivatives with designs originating from natural biomolecular motifs. This publication focuses on the synthesis of molecular electrets with ether substituents. As important as ether electret residues are for transferring holes under relatively high potentials, the synthesis of their precursors presents formidable challenges. Each residue in the molecular electrets is introduced as its 2-nitrobenzoic acid (NBA) derivative. Hence, robust and scalable synthesis of ether derivatives of NBA is essential for making such hole-transfer molecular electrets. Purdie-Irvine alkylation, using silver oxide, produces with 90% yield the esters of the NBA building block for iso-butyl ether electrets. It warrants additional ester hydrolysis for obtaining the desired NBA precursor. Conversely, Williamson etherification selectively produces the same free-acid ether derivative in one-pot reaction, but a 40% yield. The high yields of Purdie-Irvine alkylation and the selectivity of the Williamson etherification provide important guidelines for synthesizing building blocks for bioinspired molecular electrets and a wide range of other complex ether conjugates.

In Situ Spectroelectrochemical Investigations of Electrode-Confined Electron-Transferring Proteins and Redox Enzymes

This perspective analyzes recent advances in the spectroelectrochemical investigation of redox proteins and enzymes immobilized on biocompatible or biomimetic electrode surfaces. Specifically, the article highlights new insights obtained by surface-enhanced resonance Raman (SERR), surface-enhanced infrared absorption (SEIRA), protein film infrared electrochemistry (PFIRE), polarization modulation infrared reflection-absorption spectroscopy (PMIRRAS), Förster resonance energy transfer (FRET), X-ray absorption spectroscopy (XAS), electron paramagnetic resonance (EPR), and differential electrochemical mass spectrometry (DMES)-based spectroelectrochemical methods on the structure, orientation, dynamics, and reaction mechanisms for a variety of immobilized species. This includes small heme and copper electron shuttling proteins, large respiratory complexes, hydrogenases, multicopper oxidases, alcohol dehydrogenases, endonucleases, NO-reductases, and dye decolorizing peroxidases, among other enzymes. Finally, I discuss the challenges and foreseeable future developments toward a better understanding of the functioning of these complex macromolecules and their exploitation in technological devices.

Guiding Medicinal Chemistry with Fragment Molecular Orbital (FMO) Method

The understanding of binding interactions between a protein and a small molecule plays a key role in the rationalization of potency and selectivity and in design of new ideas. However, even when a target of interest is structurally enabled, visual inspection and force field-based molecular mechanics calculations cannot always explain the full complexity of the molecular interactions that are critical in drug design. Quantum mechanical methods have the potential to address this shortcoming, but traditionally, computational expense has made the application of these calculations impractical. The fragment molecular orbital (FMO) method offers a solution that combines accuracy, speed, and the ability to characterize important interactions (i.e. its strength in kcal/mol and chemical nature: hydrophobic, electrostatic, etc) that would otherwise be hard to detect. In this chapter, we describe the FMO method and illustrate its application in the discovery of the benzothiazole (BZT) series as novel tyrosine kinase ITK inhibitors for treatment of allergic asthma.

Analyzing GPCR-Ligand Interactions with the Fragment Molecular Orbital (FMO) Method

G-protein-coupled receptors (GPCRs) have enormous physiological and biomedical importance, and therefore it is not surprising that they are the targets of many prescribed drugs. Further progress in GPCR drug discovery is highly dependent on the availability of protein structural information. However, the ability of X-ray crystallography to guide the drug discovery process for GPCR targets is limited by the availability of accurate tools to explore receptor-ligand interactions. Visual inspection and molecular mechanics approaches cannot explain the full complexity of molecular interactions. Quantum mechanics (QM) approaches are often too computationally expensive to be of practical use in time-sensitive situations, but the fragment molecular orbital (FMO) method offers an excellent solution that combines accuracy, speed, and the ability to reveal key interactions that would otherwise be hard to detect. Integration of GPCR crystallography or homology modelling with FMO reveals atomistic details of the individual contributions of each residue and water molecule toward ligand binding, including an analysis of their chemical nature. Such information is essential for an efficient structure-based drug design (SBDD) process. In this chapter, we describe how to use FMO in the characterization of GPCR-ligand interactions.

Correlation between Charge Transport and Base Excision Repair in the MutY-DNA Glycosylase

Experimental evidence suggests that DNA-mediated redox signaling between high-potential [Fe₄S₄] proteins is relevant to DNA replication and repair processes, and protein-mediated charge transfer (CT) between [Fe₄S₄] clusters and nucleic acids is a fundamental process of the signaling and repair mechanisms. We analyzed the dominant CT pathways in the base excision repair glycosylase MutY using molecular dynamics simulations and hole hopping pathway analysis. We find that the adenine nucleobase of the mismatched A·oxoG DNA base pair facilitates [Fe₄S₄]-DNA CT prior to adenine excision by MutY. We also find that the R153L mutation in MutY (linked to colorectal adenomatous polyposis) influences the dominant [Fe4S4]-DNA CT pathways and appreciably decreases their effective CT rates.

Mimicking Natural Photosynthesis: Designing Ultrafast Photosensitized Electron Transfer into Multiheme Cytochrome Protein Nanowires

Efficient nanomaterials for artificial photosynthesis require fast and robust unidirectional electron transfer (ET) from photosensitizers through charge-separation and accumulation units to redox-active catalytic sites. We explored the ultrafast time-scale limits of photo-induced charge transfer between a Ru(II)tris(bipyridine) derivative photosensitizer and PpcA, a 3-heme c-type cytochrome serving as a nanoscale biological wire. Four covalent attachment sites (K28C, K29C, K52C, and G53C) were engineered in PpcA enabling site-specific covalent labeling with expected donor-acceptor (DA) distances of 4-8 Å. X-ray scattering results demonstrated that mutations and chemical labeling did not disrupt the structure of the proteins. Time-resolved spectroscopy revealed three orders of magnitude difference in charge transfer rates for the systems with otherwise similar DA distances and the same number of covalent bonds separating donors and acceptors. All-atom molecular dynamics simulations provided additional insight into the structure-function requirements for ultrafast charge transfer and the requirement of van der Waals contact between aromatic atoms of photosensitizers and hemes in order to observe sub-nanosecond ET. This work demonstrates opportunities to utilize multi-heme c-cytochromes as frameworks for designing ultrafast light-driven ET into charge-accumulating biohybrid model systems, and ultimately for mimicking the photosynthetic paradigm of efficiently coupling ultrafast, light-driven electron transfer chemistry to multi-step catalysis within small, experimentally versatile photosynthetic biohybrid assemblies.

Multifaceted aspects of charge transfer

Charge transfer and charge transport are by far among the most important processes for sustaining life on Earth and for making our modern ways of living possible. Involving multiple electron-transfer steps, photosynthesis and cellular respiration have been principally responsible for managing the energy flow in the biosphere of our planet since the Great Oxygen Event. It is impossible to imagine living organisms without charge transport mediated by ion channels, or electron and proton transfer mediated by redox enzymes. Concurrently, transfer and transport of electrons and holes drive the functionalities of electronic and photonic devices that are intricate for our lives. While fueling advances in engineering, charge-transfer science has established itself as an important independent field, originating from physical chemistry and chemical physics, focusing on paradigms from biology, and gaining momentum from solar-energy research. Here, we review the fundamental concepts of charge transfer, and outline its core role in a broad range of unrelated fields, such as medicine, environmental science, catalysis, electronics and photonics. The ubiquitous nature of dipoles, for example, sets demands on deepening the understanding of how localized electric fields affect charge transfer. Charge-transfer electrets, thus, prove important for advancing the field and for interfacing fundamental science with engineering. Synergy between the vastly different aspects of charge-transfer science sets the stage for the broad global impacts that the advances in this field have.

Backbone charge transport in double-stranded DNA

Understanding charge transport in DNA molecules is a long-standing problem of fundamental importance across disciplines^1,2. It is also of great technological interest due to DNA's ability to form versatile and complex programmable structures. Charge transport in DNA-based junctions has been reported using a wide variety of set-ups^2-4, but experiments so far have yielded seemingly contradictory results that range from insulating^5-8 or semiconducting^9,10 to metallic-like behaviour¹¹. As a result, the intrinsic charge transport mechanism in molecular junction set-ups is not well understood, which is mainly due to the lack of techniques to form reproducible and stable contacts with individual long DNA molecules. Here we report charge-transport measurements through single 30-nm-long double-stranded DNA (dsDNA) molecules with an experimental set-up that enables us to address individual molecules repeatedly and to measure the current-voltage characteristics from 5 K up to room temperature. Strikingly, we observed very high currents of tens of nanoamperes, which flowed through both homogeneous and non-homogeneous base-pair sequences. The currents are fairly temperature independent in the range 5-60 K and show a power-law decrease with temperature above 60 K, which is reminiscent of charge transport in organic crystals. Moreover, we show that the presence of even a single discontinuity ('nick') in both strands that compose the dsDNA leads to complete suppression of the current, which suggests that the backbones mediate the long-distance conduction in dsDNA, contrary to the common wisdom in DNA electronics^2-4.

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.

An Elegant Four-Helical Fold in NOX and STEAP Enzymes Facilitates Electron Transport across Biomembranes-Similar Vehicle, Different Destination

The ferric reductase superfamily comprises several oxidoreductases that use an intracellular electron source to reduce an extracellular acceptor substrate. NADPH oxidases (NOXs) and six-transmembrane epithelial antigen of the prostate enzymes (STEAPs) are iconic members of the superfamily. NOXs produce extracellular reactive oxygen species that exert potent bactericidal activities and trigger redox-signaling cascades that regulate cell division and differentiation. STEAPs catalyze the reduction of extracellular iron and copper which is necessary for the bioavailability of these essential elements. Both NOXs and STEAPs are present as multiple isozymes with distinct regulatory properties and physiological roles. Despite the important roles of NOXs and STEAPs in human physiology and despite their wide involvement in diseases like cancer, their mode of action at the molecular level remained incompletely understood for a long time, in part due to the absence of high-resolution models of the complete enzymes. Our two laboratories have elucidated the three-dimensional structures of NOXs and STEAPs, providing key insight into their mechanisms and evolution. The enzymes share a conserved transmembrane helical domain with an eye-catching hourglass shape. On the extracellular side, a heme prosthetic group is at the bottom of a pocket where the substrate (O₂ in NOX, chelated iron or copper in STEAP) is reduced. On the intracellular side, the inner heme of NOX and the FAD of STEAP are bound to topological equivalent sites. This is a rare case where critical amino acid substitutions and local conformational changes enable a cofactor (heme vs FAD) swap between two structurally and functionally conserved scaffolds. The catalytic core of these enzymes is completed by distinct cytosolic NADPH-binding domains that are topologically unrelated (a ferredoxin reductase-like flavoprotein domain in NOX and a F₄₂₀H₂:NADP⁺-like domain in STEAP), feature different quaternary structures, and underlie specific regulatory mechanisms. Despite their differences, these domains all establish electron-transfer chains that direct the electrons from NADPH to the transmembrane domain. The multistep nature of the process and the chemical nature of the products pose considerable problems in the enzymatic assays. We learned that great care must be exerted in the validation of a candidate inhibitor. Multiple orthogonal assays are required to rule out off-target effects such as ROS-scavenging activities or nonspecific interference with the enzyme redox chain. The structural analysis of STEAP/NOX enzymes led us to further notice that their transmembrane heme-binding topology is shared by other enzymes. We found that the core domain of the cytochrome b subunits of the mitochondrial complex III and photosynthetic cytochrome b6f are closely related to NOXs and STEAPs and likely arose from the same ancestor protein. This observation expands the substrate portfolio of the superfamily since cytochromes b act on ubiquinone. The rigidly packed helices of the NOX/STEAP/cytochrome b domain contrast with the more malleable membrane proteins like ion channels or amino-acid transporters, which undergo large conformational changes to allow passage of relatively large metabolites. This notion of a rigid hourglass scaffold found an unexpected confirmation in the observation, revealed by structural comparisons, that an helical bundle identical to the NOX/STEAP/cytochrome b enzymes is featured by a de novo designed heme-binding protein, PS1. Apparently, nature and protein designers have independently converged to this fold as a versatile scaffold for heme-mediated reactions. The challenge is now to uncover the molecular mechanisms that implement the isozyme-specific regulation of the enzyme functions and develop much needed inhibitors and modulators for chemical biology and drug design studies.

Nanofabrication Techniques in Large-Area Molecular Electronic Devices

The societal impact of the electronics industry is enormous—not to mention how this industry impinges on the global economy. The foreseen limits of the current technology—technical, economic, and sustainability issues—open the door to the search for successor technologies. In this context, molecular electronics has emerged as a promising candidate that, at least in the short-term, will not likely replace our silicon-based electronics, but improve its performance through a nascent hybrid technology. Such technology will take advantage of both the small dimensions of the molecules and new functionalities resulting from the quantum effects that govern the properties at the molecular scale. An optimization of interface engineering and integration of molecules to form densely integrated individually addressable arrays of molecules are two crucial aspects in the molecular electronics field. These challenges should be met to establish the bridge between organic functional materials and hard electronics required for the incorporation of such hybrid technology in the market. In this review, the most advanced methods for fabricating large-area molecular electronic devices are presented, highlighting their advantages and limitations. Special emphasis is focused on bottom-up methodologies for the fabrication of well-ordered and tightly-packed monolayers onto the bottom electrode, followed by a description of the top-contact deposition methods so far used.