A QM/MM study on the spinach plastocyanin: Redox properties and absorption spectra

Are plastocyanin and ferredoxin specific electron carriers or generic redox capacitors? Classical and murburn perspectives on two photosynthetic proteins

In the light reaction of oxygenic photosynthesis, plastocyanin (PC) and ferredoxins (Fd) are small/diffusible redox-active proteins playing key roles in electron transfer/transport phenomena. In the Z-scheme mechanistic purview, they are considered as specific affinity binding-based electron-relay agents, linking the functions of Cytochrome b₆f (Cyt. b₆f), Photosystem I (PS I) and Fd:NADPH oxidoreductase (FNR). The murburn explanation for photolytic photophosphorylation deems PC/Fd as generic 'redox capacitors', temporally accepting and releasing one-electron equivalents in reaction milieu. Herein, we explore the two theories with respect to structural, distributional and functional aspects of PC/Fd. Amino acid residues located on the surface loci of key patches of PC/Fd vary in electrostatic/contour (topography) signatures. Crystal structures of four different complexes each of Cyt.f-PC and Fd-FNR show little conservation in the contact-surfaces, thereby discrediting 'affinity binding-based electron transfers (ET)' as an evolutionary logic. Further, thermodynamic and kinetic data of wildtype and mutant proteins interactions do not align with Z-scheme. Furthermore, micromolar physiological concentrations of PC and the non-conducive architecture of chloroplasts render the classical model untenable. In the murburn model, as PC is optional, the observation that plants lacking PC survive and grow is justified. Further, the low physiological concentration/distribution of PC in chloroplast lumen/stroma is supported by murburn equilibriums, as higher concentrations would limit electron transfers. Thus, structural evidence, interactive dynamics with redox partners and physiological distribution/role of PC/Fd support the murburn perspective that these proteins serve as generic redox-capacitors in chloroplasts.Communicated by Ramaswamy H. Sarma.

2,3-(Dibenzimidazol-2-yl)quinoxalines: Unexpected Dynamical Effect on Steady-State Electronic Absorption Spectra

We report on the electronic absorption spectra, conformational behavior, and intra- and intermolecular hydrogen bonds of 2,3-(dibenzimidazol-2-yl)-quinoxaline (DBIQ). The experimentally found strong solvent dependence of the absorption spectra of DBIQ solutions cannot be assigned to electronic excitations of the equilibrium ground-state DBIQ structure. Extended consideration including the nonequilibrium structures within the framework of ab initio molecular dynamics (MD) revealed the importance of torsion molecular motions not covered by the static case. The strong impact of solute-solvent hydrogen bonding on stabilization of these nonequilibrium structures and on conformational composition of DBIQ was demonstrated. A presence of twisted nonplanar geometries along the whole MD trajectory was shown to drastically influence not only energies but also characters of electronic excitations, resulting in a change of local π-π* character in a solution of 1,2-dichloroethane to charge-transfer character in polar dimethylsulfoxide.

Charge-Transfer versus Charge-Separated Triplet Excited States of [ReI (dmp)(CO)3 (His124)(Trp122)]+ in Water and in Modified Pseudomonas aeruginosa Azurin Protein

A computational investigation of the triplet excited states of a rhenium complex electronically coupled with a tryptophan side chain and bound to an azurin protein is presented. In particular, by using high-level molecular modeling, evidence is provided for how the electronic properties of the excited-state manifolds strongly depend on coupling with the environment. Indeed, only upon explicitly taking into account the protein environment can two stable triplet states of metal-to-ligand charge transfer or charge-separated nature be recovered. In addition, it is also demonstrated how the rhenium complex plus tryptophan system in an aqueous environment experiences too much flexibility, which prevents the two chromophores from being electronically coupled. This occurrence disables the formation of a charge-separated state. The successful strategy requires a multiscale approach of combining molecular dynamics and quantum chemistry. In this context, the strategy used to parameterize the force fields for the electronic triplet states of the metal complex is also presented.

Steady-State Linear and Non-linear Optical Spectroscopy of Organic Chromophores and Bio-macromolecules

Bio-macromolecules as DNA, lipid membranes and (poly)peptides are essential compounds at the core of biological systems. The development of techniques and methodologies for their characterization is therefore necessary and of utmost interest, even though difficulties can be experienced due to their intrinsic complex nature. Among these methods, spectroscopies, relying on optical properties are especially important to determine their macromolecular structures and behaviors, as well as the possible interactions and reactivity with external dyes—often drugs or pollutants—that can (photo)sensitize the bio-macromolecule leading to eventual chemical modifications, thus damages. In this review, we will focus on the theoretical simulation of electronic spectroscopies of bio-macromolecules, considering their secondary structure and including their interaction with different kind of (photo)sensitizers. Namely, absorption, emission and electronic circular dichroism (CD) spectra are calculated and compared with the available experimental data. Non-linear properties will be also taken into account by two-photon absorption, a highly promising technique (i) to enhance absorption in the red and infra-red windows and (ii) to enhance spatial resolution. Methodologically, the implications of using implicit and explicit solvent, coupled to quantum and thermal samplings of the phase space, will be addressed. Especially, hybrid quantum mechanics/molecular mechanics (QM/MM) methods are explored for a comparison with solely QM methods, in order to address the necessity to consider an accurate description of environmental effects on spectroscopic properties of biological systems.

QM/MM methods for free energies and photochemistry

Hybrid computational methods describing a small region of a biomolecular system with quantum mechanics and the bulk with molecular mechanics, referred to as QM/MM methods, are now a central part of computational biochemistry. This review considers developments in the QM/MM approach that make it easier to calculate free energies using accurate QM-based potential energy expressions. We also describe techniques to treat electronic coupling between the core region and the MM environment. Polarizability of the protein matrix is important but so is electronic coupling. Applications of these new methods, especially to photochemistry, are discussed.

A General Route to Include Pauli Repulsion and Quantum Dispersion Effects in QM/MM Approaches

A methodology to account for nonelectrostatic interactions in Quantum Mechanical (QM)/Molecular Mechanics (MM) approaches is developed. Formulations for Pauli repulsion and dispersion energy, explicitly depending on the QM density, are derived. Such expressions are based on the definition of an auxiliary density on the MM portion and the Tkatchenko-Scheffler (TS) approach, respectively. The developed method is general enough to be applied to any QM/MM method and partition, provided an accurate tuning of a small number of parameters is obtained. The coupling of the method with both nonpolarizable and fully polarizable QM/fluctuating charge (FQ) approaches is reported and applied. A suitable parametrization for the aqueous solution, so that its most representative features are well reproduced, is outlined. Then, the obtained parametrization and method are applied to calculate the nonelectrostatic (repulsion and dispersion) interaction energy of nicotine in aqueous solution.

Modifications of the chromophore of Spinach aptamer based on QM:MM calculations

Spinach aptamer was developed as an RNA analog of the green fluorescent protein. The aptamer interacts with its ligand and modifies its electronic spectrum so that it fluoresces brightly at the wavelength of 501 nm. Song et al. investigated modifications of the ligand in their experimental study and found a molecule emitting at 523 nm upon creating a complex with the Spinach aptamer. The crystal structure of the aptamer in complex with its original ligand has been published, which enabled us to study the system computationally. In this article, we suggest several new modifications of the ligand that shift the emission maximum of the complex to even longer wavelengths. Our results are based on combined quantum mechanical/molecular mechanical calculations with DFT method used for geometry optimization and TD-DFT for calculations of absorption and emission energies.

Controlling charge-transfer properties through a microwave-assisted mono- or bis-annulation of dialkynyl-N-(het)arylpyrroles

Charge-transfer properties dependence on mono or bisannulation of dialkynyl-N-(het)arylpyrroles is reported.

A QM/MM study of the nature of the entatic state in plastocyanin

Plastocyanin is a copper containing protein that is involved in the electron transfer process in photosynthetic organisms. The active site of plastocyanin is described as an entatic state whereby its structure represents a compromise between the structures favored by the oxidized and reduced forms. In this study, the nature of the entatic state is investigated through density functional theory‐based hybrid quantum mechanics/molecular mechanics (QM/MM) molecular dynamics simulations. The strain energy is computed to be 12.8 kcal/mol and 14.5 kcal/mol for the oxidized and reduced forms of the protein, indicating that the active site has an intermediate structure. It is shown that the energy gap between the oxidized and reduced forms varies significantly with the fluctuations in the structure of the active site at room temperature. An accurate determination of the reorganization energy requires averaging over conformation and a large region of the protein around the active site to be treated at the quantum mechanical level. © 2016 The Authors. Journal of Computational Chemistry Published by Wiley Periodicals, Inc.

Kβ Valence to Core X-ray Emission Studies of Cu(I) Binding Proteins with Mixed Methionine - Histidine Coordination. Relevance to the Reactivity of the M- and H-sites of Peptidylglycine Monooxygenase

Biological systems use copper as a redox center in many metalloproteins, where the role of the metal is to cycle between its +1 and +2 oxidation states. This chemistry requires the redox potential to be in a range that can stabilize both Cu(I) and Cu(II) states and often involves protein-derived ligand sets involving mixed histidine-methionine coordination that balance the preferences of both oxidation states. Transport proteins, on the other hand, utilize copper in the Cu(I) state and often contain sites comprised predominately of the cuprophilic residue methionine. The electronic factors that allow enzymes and transporters to balance their redox requirements are complex and are often elusive due to the dearth of spectroscopic probes of the Cu(I) state. Here we present the novel application of X-ray emission spectroscopy to copper proteins via a study of a series of mixed His-Met copper sites where the ligand set varies in a systematic way between the His3 and Met3 limits. The sites are derived from the wild-type peptidylglycine monooxygenase (PHM), two single-site variants which replicate each of its two copper sites (CuM-site and CuH-site), and the transporters CusF and CusB. Clear differences are observed in the Kβ2,5 region at the Met3 and His3 limits. CusB (Met3) has a distinct peak at 8978.4 eV with a broad shoulder at 8975.6 eV, whereas CuH (His3) has two well-resolved features: a more intense feature at 8974.8 eV and a second at 8977.2 eV. The mixed coordination sphere CusF (Met2His) and the PHM CuM variant (Met1His2) have very similar spectra consisting of two features at 8975.2 and 8977.8 eV. An analysis of DFT calculated spectra indicate that the intensity of the higher energy peak near 8978 eV is mediated by mixing of ligand-based orbitals into the Cu d(10) manifold, with S from Met providing more intensity by facilitating increased Cu p-d mixing. Furthermore, reaction of WT PHM with CO (an oxygen analogue) produced the M site CO complex, which showed a unique XES spectrum that could be computationally reproduced by including interactions between Cu(I) and the CO ligand. The study suggests that the valence-to-core (VtC) region can not only serve as a probe of ligand speciation but also offer insight into the coordination geometry, in a fashion similar to XAS pre-edges, and may be sufficiently sensitive to the coordination of exogenous ligands to be useful in the study of reaction mechanisms.

A new record excited state 3MLCT lifetime for metalorganic iron(ii) complexes

The ³MLCT lifetime of iron complexes was pushed to a record lifetime of 26 ps using a new benzimidazolylidene-based C^∧N^∧C ligand opening up new perspectives on the way of applying such complexes to photoelectrochemical processes.

Transition matrices and orbitals from reduced density matrix theory

In this contribution, we report two different methodologies for characterizing the electronic structure reorganization occurring when a chromophore undergoes an electronic transition. For the first method, we start by setting the theoretical background necessary to the reinterpretation through simple tensor analysis of (i) the transition density matrix and (ii) the natural transition orbitals in the scope of reduced density matrix theory. This novel interpretation is made more clear thanks to a short compendium of the one-particle reduced density matrix theory in a Fock space. The formalism is further applied to two different classes of excited states calculation methods, both requiring a single-determinant reference, that express an excited state as a hole-particle mono-excited configurations expansion, to which particle-hole correlation is coupled (time-dependent Hartree-Fock/time-dependent density functional theory) or not (configuration interaction single/Tamm-Dancoff approximation). For the second methodology presented in this paper, we introduce a novel and complementary concept related to electronic transitions with the canonical transition density matrix and the canonical transition orbitals. Their expression actually reflects the electronic cloud polarisation in the orbital space with a decomposition based on the actual contribution of one-particle excitations from occupied canonical orbitals to virtual ones. This approach validates our novel interpretation of the transition density matrix elements in terms of the Euclidean norm of elementary transition vectors in a linear tensor space. A proper use of these new concepts leads to the conclusion that despite the different principles underlying their construction, they provide two equivalent excited states topological analyses. This connexion is evidenced through simple illustrations of (in)organic dyes electronic transitions analysis.

Dye chemistry with time-dependent density functional theory

In this perspective, we present an overview of the determination of excited-state properties of "real-life" dyes, and notably of their optical absorption and emission spectra, performed during the last decade with time-dependent density functional theory (TD-DFT). We discuss the results obtained with both vertical and adiabatic (vibronic) approximations, choosing relevant examples for several series of dyes. These examples include reproducing absorption wavelengths of numerous families of coloured molecules, understanding the specific band shape of amino-anthraquinones, optimising the properties of dyes used in solar cells, mimicking the fluorescence wavelengths of fluorescent brighteners and BODIPY dyes, studying optically active biomolecules and photo-induced proton transfer, as well as improving the properties of photochromes.

2,5-Dithienylpyrrole (DTP) as a donor component in DTP–π–A organic sensitizers: photophysical and photovoltaic properties

DTP–π–A organic dyes efficiently sensitized TiO₂ in DSSCs exhibiting nearly quantitative charge injection yields and regeneration efficiency.