Fragment Molecular Orbital Study on Electron Tunneling Mechanisms in Bacterial Photosynthetic Reaction Center

Toward Accurate Prediction of Ion Mobility in Organic Semiconductors by Atomistic Simulation

A multiscale scheme (MLMS: Multi-Level Multi-Scale) to predict the ion mobility (μ) of amorphous organic semiconductors is proposed, which was successfully applied to the hole mobility predictions of 14 organic systems. An inverse relationship between μ and reorganization energy is observed due to local polaronic distortions. Another moderate inverse correlation between μ and distribution of site energy change exists, representing the effects of geometric flexibility. The former and the latter represent the intramolecular and intermolecular geometric effects, respectively. In addition, a linear correlation between transfer coupling and μ is observed, showing the importance of orbital overlaps between monomers. Especially, the highest hole mobility of C6-2TTN is due to its large transfer coupling. On the other hand, another high hole mobility of CBP turned out to come from the high first neighbor density (ρ_FND) of its first self-solvation, emphasizing the proper description of amorphous structural configurations with a sufficiently large number of monomers. In general, systems with either unusually high transfer coupling or high first neighbor density can potentially have high μ regardless of geometric effects. Especially, the newly suggested design parameter, ρ_FND, is pointing to a new direction as opposed to the traditional π-conjugation strategy.

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.

Recognition of repulsive and attractive regions of selected serotonin receptor binding site using FMO-EDA approach

The complexes of selected long-chain arylpiperazines with homology models of 5-HT_1A, 5-HT_2A, and 5-HT₇ receptors were investigated using quantum mechanical methods. The molecular geometries of the ligand-receptor complexes were firstly optimized with the Our own N-layered Integrated molecular Orbital and molecular Mechanics (ONIOM) method. Next, the fragment molecular orbitals method with an energy decomposition analysis scheme (FMO-EDA) was employed to estimate the interaction energies in binding sites. The results clearly showed that orthosteric binding sites of studied serotonin receptors have both attractive and repulsive regions. In the case of 5-HT_1A and 5-HT_2A two repulsive areas, located in the lower part of the binding pocket, and one large area of attraction engaging many residues at the top of all helices were identified. Additionally, for the 5-HT₇ receptor, the third area of destabilization located at the extracellular end of the helix 6 was found.

Electron Transfer Pathways of Cyclobutane Pyrimidine Dimer Photolyase Revisited

The photoinduced electron transfer (ET) reaction of cyclobutane pyrimidine dimer (CPD) photolyase plays an essential role in its DNA repair reaction, and the molecular mechanism of the ET reaction has attracted a large number of experimental and theoretical studies. We investigated the quantum mechanical nature of their ET reactions, characterized by multiple ET pathways of the CPD photolyase derived from Anacystis nidulans. Using the generalized Mulliken-Hush (GMH) method and the bridge green function (GF) methods, we estimated the electronic coupling matrix element, T_DA, to be 36 ± 30 cm^-1 from the donor (FADH^-) to the acceptor (CPD). The estimated ET time was 386 ps, in good agreement with the experimental value (250 ps) in the literature. Furthermore, we performed the molecular dynamics (MD) simulations and ab initio molecular orbital (MO) calculations, and explored the electron tunneling pathway. We examined 20 different structures during the MD trajectory and quantitatively evaluated the electron tunneling currents for each of them. As a result, we demonstrated that the ET route via Asn349 was the dominant pathway among the five major routes via (Adenine/Asn349), (Adenine/Glu283), (Adenine/Glu283/Asn349/Met353), (Met353/Asn349), and (Asn349), indicating that Asn349 is an essential amino acid residue in the ET reaction.

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

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

Photosynthetic diode: electron transport rectification by wetting the quinone cofactor

We report 11 μs of molecular dynamics simulations of the electron-transfer reaction between primary and secondary quinone cofactors in the bacterial reaction center. The main question addressed here is the mechanistic reason for unidirectional electron transfer between chemically identical cofactors. We find that electron is trapped at the secondary quinone by wetting of the protein pocket following electron transfer on the time-scale shorter than the backward transition. This mechanism provides effective rectification of the electron transport, making the reaction center a molecular diode operating by cyclic charge-induced electrowetting.

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.

Electron-correlated fragment-molecular-orbital calculations for biomolecular and nano systems

Recent developments in the fragment molecular orbital (FMO) method for theoretical formulation, implementation, and application to nano and biomolecular systems are reviewed. The FMO method has enabled ab initio quantum-mechanical calculations for large molecular systems such as protein-ligand complexes at a reasonable computational cost in a parallelized way. There have been a wealth of application outcomes from the FMO method in the fields of biochemistry, medicinal chemistry and nanotechnology, in which the electron correlation effects play vital roles. With the aid of the advances in high-performance computing, the FMO method promises larger, faster, and more accurate simulations of biomolecular and related systems, including the descriptions of dynamical behaviors in solvent environments. The current status and future prospects of the FMO scheme are addressed in these contexts.

The nonheme iron in photosystem II

Photosystem II (PSII), the light-driven water:plastoquinone (PQ) oxidoreductase of oxygenic photosynthesis, contains a nonheme iron (NHI) at its electron acceptor side. The NHI is situated between the two PQs QA and QB that serve as one-electron transmitter and substrate of the reductase part of PSII, respectively. Among the ligands of the NHI is a (bi)carbonate originating from CO2, the substrate of the dark reactions of oxygenic photosynthesis. Based on recent advances in the crystallography of PSII, we review the structure of the NHI in PSII and discuss ideas concerning its function and the role of bicarbonate along with a comparison to the reaction center of purple bacteria and other enzymes containing a mononuclear NHI site.