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

Observation of Giant Conductance Fluctuations in a Protein

Proteins are insulating molecular solids, yet even those containing easily reduced or oxidized centers can have single-molecule electronic conductances that are too large to account for with conventional transport theories. Here, we report the observation of remarkably high electronic conductance states in an electrochemically-inactive protein, the ~200 kD α_Vβ₃ extracelluar domain of human integrin. Large current pulses (up to nA) were observed for long durations (many ms, corresponding to many pC of charge transfer) at large gap (>5nm) distances in an STM when the protein was bound specifically by a small peptide ligand attached to the electrodes. The effect is greatly reduced when a homologous, weakly-binding protein (α₄β₁) is used as a control. In order to overcome the limitations of the STM, the time- and voltage-dependence of the conductance were further explored using a fixed-gap (5 nm) tunneling junction device that was small enough to trap a single protein molecule at any one time. Transitions to a high conductance (~ nS) state were observed, the protein being "on" for times from ms to tenths of a second. The high-conductance states only occur above ~ 100mV applied bias, and thus are not an equilibrium property of the protein. Nanoamp two-level signals indicate the specific capture of a single molecule in an electrode gap functionalized with the ligand. This offers a new approach to label-free electronic detection of single protein molecules. Electronic structure calculations yield a distribution of energy level spacings that is consistent with a recently proposed quantum-critical state for proteins, in which small fluctuations can drive transitions between localized and band-like electronic states.

Using the fragment molecular orbital method to investigate agonist-orexin-2 receptor interactions

The understanding of binding interactions between any protein and a small molecule plays a key role in the rationalization of affinity and selectivity and is essential for an efficient structure-based drug discovery (SBDD) process. Clearly, to begin SBDD, a structure is needed, and although there has been fantastic progress in solving G-protein-coupled receptor (GPCR) crystal structures, the process remains quite slow and is not currently feasible for every GPCR or GPCR-ligand complex. This situation significantly limits the ability of X-ray crystallography to impact the drug discovery process for GPCR targets in 'real-time' and hence there is still a need for other practical and cost-efficient alternatives. We present here an approach that integrates our previously described hierarchical GPCR modelling protocol (HGMP) and the fragment molecular orbital (FMO) quantum mechanics (QM) method to explore the interactions and selectivity of the human orexin-2 receptor (OX2R) and its recently discovered nonpeptidic agonists. HGMP generates a 3D model of GPCR structures and its complexes with small molecules by applying a set of computational methods. FMO allowsab initioapproaches to be applied to systems that conventional QM methods would find challenging. The key advantage of FMO is that it can reveal information on the individual contribution and chemical nature of each residue and water molecule to the ligand binding that normally would be difficult to detect without QM. We illustrate how the combination of both techniques provides a practical and efficient approach that can be used to analyse the existing structure-function relationships (SAR) and to drive forward SBDD in a real-world example for which there is no crystal structure of the complex available.

Stereoselective Fluorescence Quenching in the Electron Transfer Photooxidation of Nucleobase-Related Azetidines by Cyanoaromatics

Electron transfer involving nucleic acids and their derivatives is an important field in bioorganic chemistry, specifically in connection with its role in the photo-driven DNA damage and repair. Four-membered ring heterocyclic oxetanes and azetidines have been claimed to be the intermediates involved in the repair of DNA (6-4) photoproduct by photolyase. In this context, we examine here the redox properties of the two azetidine isomers obtained from photocycloaddition between 6-aza-1,3-dimethyluracil and cyclohexene. Steady-state and time-resolved fluorescence experiments using a series of photoreductants and photooxidants have been run to evaluate the efficiency of the electron transfer process. Analysis of the obtained quenching kinetics shows that the azetidine compounds can act as electron donors. Additionally, it appears that the cis isomer is more easily oxidized than its trans counterpart. This result is in agreement with electrochemical studies performed on both azetidine derivatives.

Coherent control of long-range photoinduced electron transfer by stimulated X-ray Raman processes

We show that X-ray pulses resonant with selected core transitions can manipulate electron transfer (ET) in molecules with ultrafast and atomic selectivity. We present possible protocols for coherently controlling ET dynamics in donor-bridge-acceptor (DBA) systems by stimulated X-ray resonant Raman processes involving various transitions between the D, B, and A sites. Simulations presented for a Ru(II)-Co(III) model complex demonstrate how the shapes, phases and amplitudes of the X-ray pulses can be optimized to create charge on demand at selected atoms, by opening up otherwise blocked ET pathways.

Mean First-Passage Times in Biology

Many biochemical processes, such as charge hopping or protein folding, can be described by an average timescale to reach a final state, starting from an initial state. Here, we provide a pedagogical treatment of the mean first-passage time (MFPT) of a physical process, which depends on the number of intervening states between the initial state and the target state. Our aim in this tutorial review is to provide a clear development of the mean first-passage time formalism and to show some of its practical utility. The MFPT treatment can provide a useful link between microscopic rates and the average timescales often probed by experiment.

Charge transport in molecular junctions: From tunneling to hopping with the probe technique

We demonstrate that a simple phenomenological approach can be used to simulate electronic conduction in molecular wires under thermal effects induced by the surrounding environment. This "Landauer-Büttiker's probe technique" can properly replicate different transport mechanisms, phase coherent nonresonant tunneling, ballistic behavior, and hopping conduction. Specifically, our simulations with the probe method recover the following central characteristics of charge transfer in molecular wires: (i) the electrical conductance of short wires falls off exponentially with molecular length, a manifestation of the tunneling (superexchange) mechanism. Hopping dynamics overtakes superexchange in long wires demonstrating an ohmic-like behavior. (ii) In off-resonance situations, weak dephasing effects facilitate charge transfer, but under large dephasing, the electrical conductance is suppressed. (iii) At high enough temperatures, kBT/ϵB > 1/25, with ϵB as the molecular-barrier height, the current is enhanced by a thermal activation (Arrhenius) factor. However, this enhancement takes place for both coherent and incoherent electrons and it does not readily indicate on the underlying mechanism. (iv) At finite-bias, dephasing effects may impede conduction in resonant situations. We further show that memory (non-Markovian) effects can be implemented within the Landauer-Büttiker's probe technique to model the interaction of electrons with a structured environment. Finally, we examine experimental results of electron transfer in conjugated molecular wires and show that our computational approach can reasonably reproduce reported values to provide mechanistic information.

Control of electrochemical signals from quantum dots conjugated to organic materials by using DNA structure in an analog logic gate

Various bio-logic gates have been studied intensively to overcome the rigidity of single-function silicon-based logic devices arising from combinations of various gates. Here, a simple control tool using electrochemical signals from quantum dots (QDs) was constructed using DNA and organic materials for multiple logic functions. The electrochemical redox current generated from QDs was controlled by the DNA structure. DNA structure, in turn, was dependent on the components (organic materials) and the input signal (pH). Independent electrochemical signals from two different logic units containing QDs were merged into a single analog-type logic gate, which was controlled by two inputs. We applied this electrochemical biodevice to a simple logic system and achieved various logic functions from the controlled pH input sets. This could be further improved by choosing QDs, ionic conditions, or DNA sequences. This research provides a feasible method for fabricating an artificial intelligence system.

Relation between Dephasing Time and Energy Gap Fluctuations in Biomolecular Systems

Excitation energy and charge transfer are fundamental processes in biological systems. Because of their quantum nature, the effect of dephasing on these processes is of interest especially when trying to understand their efficiency. Moreover, recent experiments have shown quantum coherences in such systems. As a first step toward a better understanding, we studied the relationship between dephasing time and energy gap fluctuations of the individual molecular subunits. A larger set of molecular simulations has been investigated to shed light on this dependence. This set includes bacterio-chlorophylls in Fenna-Matthews-Olson complexes, the PE545 aggregate, the LH2 complexes, DNA, photolyase, and cryptochromes. For the individual molecular subunits of these aggregates it has been confirmed quantitatively that an inverse proportionality exists between dephasing time and average gap energy fluctuation. However, for entire complexes including the respective intermolecular couplings, such a relation still needs to be verified.

Fragment Molecular Orbital Method Applied to Lead Optimization of Novel Interleukin-2 Inducible T-Cell Kinase (ITK) Inhibitors

Inhibition of inducible T-cell kinase (ITK), a nonreceptor tyrosine kinase, may represent a novel treatment for allergic asthma. In our previous reports, we described the discovery of sulfonylpyridine (SAP), benzothiazole (BZT), indazole (IND), and tetrahydroindazole (THI) series as novel ITK inhibitors and how computational tools such as dihedral scans and docking were used to support this process. X-ray crystallography and modeling were applied to provide essential insight into ITK-ligand interactions. However, "visual inspection" traditionally used for the rationalization of protein-ligand affinity cannot always explain the full complexity of the molecular interactions. The fragment molecular orbital (FMO) quantum-mechanical (QM) method provides a complete list of the interactions formed between the ligand and protein that are often omitted from traditional structure-based descriptions. FMO methodology was successfully used as part of a rational structure-based drug design effort to improve the ITK potency of high-throughput screening hits, ultimately delivering ligands with potency in the subnanomolar range.

High-Energy Long-Lived Mixed Frenkel-Charge-Transfer Excitons: From Double Stranded (AT)n to Natural DNA

The electronic excited states populated upon absorption of UV photons by DNA are extensively studied in relation to the UV-induced damage to the genetic code. Here, we report a new unexpected relaxation pathway in adenine-thymine double-stranded structures (AT)n . Fluorescence measurements on (AT)n hairpins (six and ten base pairs) and duplexes (20 and 2000 base pairs) reveal the existence of an emission band peaking at approximately 320 nm and decaying on the nanosecond time scale. Time-dependent (TD)-DFT calculations, performed for two base pairs and exploring various relaxation pathways, allow the assignment of this emission band to excited states resulting from mixing between Frenkel excitons and adenine-to-thymine charge-transfer states. Emission from such high-energy long-lived mixed (HELM) states is in agreement with their fluorescence anisotropy (0.03), which is lower than that expected for π-π* states (≥0.1). An increase in the size of the system quenches π-π* fluorescence while enhancing HELM fluorescence. The latter process varies linearly with the hypochromism of the absorption spectra, both depending on the coupling between π-π* and charge-transfer states. Subsequently, we identify the common features between the HELM states of (AT)n structures with those reported previously for alternating (GC)n : high emission energy, low fluorescence anisotropy, nanosecond lifetimes, and sensitivity to conformational disorder. These features are also detected for calf thymus DNA in which HELM states could evolve toward reactive π-π* states, giving rise to delayed fluorescence.

Radical transfer in E. coli ribonucleotide reductase: a NH2Y731/R411A-α mutant unmasks a new conformation of the pathway residue 731

Ribonucleotide reductases (RNRs) catalyze the conversion of ribonucleotides to deoxyribonucleotides in all living organisms. The catalytic cycle of E. coli RNR involves a long-range proton-coupled electron transfer (PCET) from a tyrosyl radical (Y122˙) in subunit β2 to a cysteine (C439) in the active site of subunit α2, which subsequently initiates nucleotide reduction. This oxidation occurs over 35 Å and involves a specific pathway of redox active amino acids (Y122 ↔ [W48?] ↔ Y356 in β2 to Y731 ↔ Y730 ↔ C439 in α2). The mechanisms of the PCET steps at the interface of the α2β2 complex remain puzzling due to a lack of structural information for this region. Recently, DFT calculations on the 3-aminotyrosyl radical (NH2Y731˙)-α2 trapped by incubation of NH2Y731-α2/β2/CDP(substrate)/ATP(allosteric effector) suggested that R411-α2, a residue close to the α2β2 interface, interacts with NH2Y731˙ and accounts in part for its perturbed EPR parameters. To examine its role, we further modified NH2Y731-α2 with a R411A substitution. NH2Y731˙/R411A generated upon incubation of NH2Y731/R411A-α2/β2/CDP/ATP was investigated using multi-frequency (34, 94 and 263 GHz) EPR, 34 GHz pulsed electron-electron double resonance (PELDOR) and electron-nuclear double resonance (ENDOR) spectroscopies. The data indicate a large conformational change in NH2Y731˙/R411A relative to the NH2Y731˙ single mutant. Particularly, the inter-spin distance from NH2Y731˙/R411A in one αβ pair to Y122˙ in a second αβ pair decreases by 3 Å in the presence of the R411A mutation. This is the first experimental evidence for the flexibility of pathway residue Y731-α2 in an α2β2 complex and suggests a role for R411 in the stacked Y731/Y730 conformation involved in collinear PCET. Furthermore, NH2Y731˙/R411A serves as a probe of the PCET process across the subunit interface.

Quantum effects in ultrafast electron transfers within cryptochromes

Cryptochromes and photolyases are flavoproteins that may undergo ultrafast charge separation upon electronic excitation of their flavin cofactors.

Electron transport via a soluble photochromic photoreceptor

Electron transport properties via a photochromic biological photoreceptor have been studied in junctions of monolayer assemblies in solid-state configurations.

The Fragment Molecular Orbital Method Reveals New Insight into the Chemical Nature of GPCR–Ligand Interactions

Our interpretation of ligand-protein interactions is often informed by high-resolution structures, which represent the cornerstone of structure-based drug design. However, visual inspection and molecular mechanics approaches cannot explain the full complexity of molecular interactions. Quantum Mechanics approaches are often too computationally expensive, but one method, Fragment Molecular Orbital (FMO), offers an excellent compromise and has the potential to reveal key interactions that would otherwise be hard to detect. To illustrate this, we have applied the FMO method to 18 Class A GPCR-ligand crystal structures, representing different branches of the GPCR genome. Our work reveals key interactions that are often omitted from structure-based descriptions, including hydrophobic interactions, nonclassical hydrogen bonds, and the involvement of backbone atoms. This approach provides a more comprehensive picture of receptor-ligand interactions than is currently used and should prove useful for evaluation of the chemical nature of ligand binding and to support structure-based drug design.

Building blocks for bioinspired electrets: molecular-level approach to materials for energy and electronics

In biology, an immense diversity of protein structural and functional motifs originates from only 20 common proteinogenic native amino acids arranged in various sequences. Is it possible to attain the same diversity in electronic materials based on organic macromolecules composed of non-native residues with different characteristics? This publication describes the design, preparation and characterization of non-native aromatic β-amino acid residues, i.e. derivatives of anthranilic acid, for polyamides that can efficiently mediate hole transfer. Chemical derivatization with three types of substituents at two positions of the aromatic ring allows for adjusting the energy levels of the frontier orbitals of the anthranilamide residues over a range of about one electronvolt. Most importantly, the anthranilamide residues possess permanent electric dipoles, adding to the electronic properties of the bioinspired conjugates they compose, making them molecular electrets.

Conformationally Gated Charge Transfer in DNA Three-Way Junctions

Molecular structures that direct charge transport in two or three dimensions possess some of the essential functionality of electrical switches and gates. We use theory, modeling, and simulation to explore the conformational dynamics of DNA three-way junctions (TWJs) that may control the flow of charge through these structures. Molecular dynamics simulations and quantum calculations indicate that DNA TWJs undergo dynamic interconversion among "well stacked" conformations on the time scale of nanoseconds, a feature that makes the junctions very different from linear DNA duplexes. The studies further indicate that this conformational gating would control charge flow through these TWJs, distinguishing them from conventional (larger size scale) gated devices. Simulations also find that structures with polyethylene glycol linking groups ("extenders") lock conformations that favor CT for 25 ns or more. The simulations explain the kinetics observed experimentally in TWJs and rationalize their transport properties compared with double-stranded DNA.

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.

Probing and Exploiting the Interplay between Nuclear and Electronic Motion in Charge Transfer Processes

The Born-Oppenheimer approximation refers to the assumption that the nuclear and electronic wave functions describing a molecular system evolve and can be determined independently. It is now well-known that this approximation often breaks down and that nuclear-electronic (vibronic) coupling contributes greatly to the ultrafast photophysics and photochemistry observed in many systems ranging from simple molecules to biological organisms. In order to probe vibronic coupling in a time-dependent manner, one must use spectroscopic tools capable of correlating the motions of electrons and nuclei on an ultrafast time scale. Recent developments in nonlinear multidimensional electronic and vibrational spectroscopies allow monitoring both electronic and structural factors with unprecedented time and spatial resolution. In this Account, we present recent studies from our group that make use of different variants of frequency-domain transient two-dimensional infrared (T-2DIR) spectroscopy, a pulse sequence combining electronic and vibrational excitations in the form of a UV-visible pump, a narrowband (12 cm(-1)) IR pump, and a broadband (400 cm(-1)) IR probe. In the first example, T-2DIR is used to directly compare vibrational dynamics in the ground and relaxed electronic excited states of Re(Cl)(CO)3(4,4'-diethylester-2,2'-bipyridine) and Ru(4,4'-diethylester-2,2'-bipyridine)2(NCS)2, prototypical charge transfer complexes used in photocatalytic CO2 reduction and electron injection in dye-sensitized solar cells. The experiments show that intramolecular vibrational redistribution (IVR) and vibrational energy transfer (VET) are up to an order of magnitude faster in the triplet charge transfer excited state than in the ground state. These results show the influence of electronic arrangement on vibrational coupling patterns, with direct implications for vibronic coupling mechanisms in charge transfer excited states. In the second example, we show unambiguously that electronic and vibrational movement are coupled in a donor-bridge-acceptor complex based on a Pt(II) trans-acetylide design motif. Time-resolved IR (TRIR) spectroscopy reveals that the rate of electron transfer (ET) is highly dependent on the amount of excess energy localized on the bridge following electronic excitation. Using an adaptation of T-2DIR, we are able to selectively perturb bridge-localized vibrational modes during charge separation, resulting in the donor-acceptor charge separation pathway being completely switched off, with all excess energy redirected toward the formation of a long-lived intraligand triplet state. A series of control experiments reveal that this effect is mode specific: it is only when the high-frequency bridging C≡C stretching mode is pumped that radical changes in photoproduct yields are observed. These experiments therefore suggest that one may perturb electronic movement by stimulating structural motion along the reaction coordinate using IR light. These studies add to a growing body of evidence suggesting that controlling the pathways and efficiency of charge transfer may be achieved through synthetic and perturbative approaches aiming to modulate vibronic coupling. Achieving such control would represent a breakthrough for charge transfer-based applications such as solar energy conversion and molecular electronics.

Proton-coupled electron transfer dynamics in the catalytic mechanism of a [NiFe]-hydrogenase

The movement of protons and electrons is common to the synthesis of all chemical fuels such as H2. Hydrogenases, which catalyze the reversible reduction of protons, necessitate transport and reactivity between protons and electrons, but a detailed mechanism has thus far been elusive. Here, we use a phototriggered chemical potential jump method to rapidly initiate the proton reduction activity of a [NiFe] hydrogenase. Coupling the photochemical initiation approach to nanosecond transient infrared and visible absorbance spectroscopy afforded direct observation of interfacial electron transfer and active site chemistry. Tuning of intramolecular proton transport by pH and isotopic substitution revealed distinct concerted and stepwise proton-coupled electron transfer mechanisms in catalysis. The observed heterogeneity in the two sequential proton-associated reduction processes suggests a highly engineered protein environment modulating catalysis and implicates three new reaction intermediates; Nia-I, Nia-D, and Nia-SR(-). The results establish an elementary mechanistic understanding of catalysis in a [NiFe] hydrogenase with implications in enzymatic proton-coupled electron transfer and biomimetic catalyst design.

Electron-transfer acceleration investigated by time resolved infrared spectroscopy

Ultrafast electron transfer (ET) processes are important primary steps in natural and artificial photosynthesis, as well as in molecular electronic/photonic devices. In biological systems, ET often occurs surprisingly fast over long distances of several tens of angströms. Laser-pulse irradiation is conveniently used to generate strongly oxidizing (or reducing) excited states whose reactions are then studied by time-resolved spectroscopic techniques. While photoluminescence decay and UV-vis absorption supply precise kinetics data, time-resolved infrared absorption (TRIR) and Raman-based spectroscopies have the advantage of providing additional structural information and monitoring vibrational energy flows and dissipation, as well as medium relaxation, that accompany ultrafast ET. We will discuss three cases of photoinduced ET involving the Re(I)(CO)3(N,N) moiety (N,N = polypyridine) that occur much faster than would be expected from ET theories. [Re(4-N-methylpyridinium-pyridine)(CO)3(N,N)](2+) represents a case of excited-state picosecond ET between two different ligands that remains ultrafast even in slow-relaxing solvents, beating the adiabatic limit. This is caused by vibrational/solvational excitation of the precursor state and participation of high-frequency quantum modes in barrier crossing. The case of Re-tryptophan assemblies demonstrates that excited-state Trp → *Re(II) ET is accelerated from nanoseconds to picoseconds when the Re(I)(CO)3(N,N) chromophore is appended to a protein, close to a tryptophan residue. TRIR in combination with DFT calculations and structural studies reveals an interaction between the N,N ligand and the tryptophan indole. It results in partial electronic delocalization in the precursor excited state and likely contributes to the ultrafast ET rate. Long-lived vibrational/solvational excitation of the protein Re(I)(CO)3(N,N)···Trp moiety, documented by dynamic IR band shifts, could be another accelerating factor. The last discussed process, back-ET in a porphyrin-Re(I)(CO)3(N,N) dyad, demonstrates that formation of a hot product accelerates highly exergonic ET in the Marcus inverted region. Overall, it follows that ET can be accelerated by enhancing the electronic interaction and by vibrational excitation of the reacting system and its medium, stressing the importance of quantum nuclear dynamics in ET reactivity. These effects are experimentally accessible by time-resolved vibrational spectroscopies (IR, Raman) in combination with quantum chemical calculations. It is suggested that structural dynamics play different mechanistic roles in light-triggered ET involving electronically excited donors or acceptors than in ground-state processes. While TRIR spectroscopy is well suitable to elucidate ET processes on a molecular-level, transient 2D-IR techniques combining optical and two IR (or terahertz) laser pulses present future opportunities for investigating, driving, and controlling ET.

Bidirectional Photoinduced Electron Transfer in Ruthenium(II)-Tris-bipyridyl-Modified PpcA, a Multi-heme c-Type Cytochrome from Geobacter sulfurreducens

PpcA, a tri-heme cytochrome c7 from Geobacter sulfurreducens, was investigated as a model for photosensitizer-initiated electron transfer within a multi-heme "molecular wire" protein architecture. Escherichia coli expression of PpcA was found to be tolerant of cysteine site-directed mutagenesis, demonstrated by the successful expression of natively folded proteins bearing cysteine mutations at a series of sites selected to vary characteristically with respect to the three -CXXCH- heme binding domains. The introduced cysteines readily reacted with Ru(II)-(2,2'-bpy)2(4-bromomethyl-4'-methyl-2,2'-bipyridine) to form covalently linked constructs that support both photo-oxidative and photo-reductive quenching of the photosensitizer excited state, depending upon the initial heme redox state. Excited-state electron-transfer times were found to vary from 6 × 10(-12) to 4 × 10(-8) s, correlated with the distance and pathways for electron transfer. The fastest rate is more than 10(3)-fold faster than previously reported for photosensitizer-redox protein constructs using amino acid residue linking. Clear evidence for inter-heme electron transfer within the multi-heme protein is not detected within the lifetimes of the charge-separated states. These results demonstrate an opportunity to develop multi-heme c-cytochromes for investigation of electron transfer in protein "molecular wires" and to serve as frameworks for metalloprotein designs that support multiple-electron-transfer redox chemistry.

Vibrational control of electron-transfer reactions: a feasibility study for the fast coherent transfer regime

Electron donors are connectedvialeft and right bridges to electron acceptors. Following electron-transfer initiation, the IR excitation of selected bridge vibrational modes can tune the directionality of electron transfer.