Chlorophyll and pheophytin protonated and deprotonated ions: Observation and theory.

Pheophytin a and chlorophyll a have been investigated by electrospray mass spectrometry in the positive and negative modes, in view of the importance of the knowledge of their properties in photosynthesis. Pheophytin and chlorophyll are both observed intensely in the protonated mode, and their main fragmentation route is the loss of their phytyl chain. Pheophytin is observed intact in the negative mode, while under collisions, it is primarily cleaved beyond the phytyl chain and loses the attaching propionate group. Chlorophyll is not detected in normal conditions in the negative mode, but addition of methanol solvent molecule is detected. Fragmentation of this adduct primarily forms a product (-30 amu) that dissociates into dephytyllated deprotonated chlorophyll. Semi-empirical molecular dynamics calculations show that the phytyl chain is unfolded from the chlorin cycle in pheophytin a and folded in chlorophyll a. Density functional theory calculations have been conducted to locate the charges on protonated and deprotonated pheophytin a and chlorophyll a and have found the major location sites that are notably more stable in energy by more than 0.5 eV than the others. The deprotonation site is found identical for pheophytin a and the chlorophyll a-methanol adduct. This is in line with experiment and calculation locating the addition of methanol on a double bond of deprotonated chlorophyll a.

QM/MM study of l-lactate oxidation by flavocytochrome b2.

In this work, we have performed molecular dynamics simulations using a hybrid Quantum Mechanics/Molecular Mechanics (QM/MM) scheme to study the mechanism of l-lactate oxidation by flavocytochrome b2 (Fcb2). Our results obtained at the QM(AM1)/MM level have been improved by single-point corrections using density functional theory (DFT) methods. Free energy surfaces have been calculated in the framework of the hydride transfer hypothesis. This mechanism involves the transfer of the lactate hydroxyl proton to H373 while the substrate αH atom is transferred as a hydride to the flavin mononucleotide (FMN) prosthetic group anchored in the active site. Four different systems have been modeled: wild-type enzyme considering R289 in a distal or a proximal conformation observed in crystal structures and the D282N and Y254L variants (with R289 in a distal position). Simulation results highlight the influence of the environment on the catalytic mechanism by describing a step-wise process in the WT enzyme with R289 in a distal position and a concerted mechanism for the other systems. In the step-wise mechanism, the hydride transfer to flavin can occur only after a proton transfer from substrate to H373. Modifications of the electrostatic field around l-lactate or H373 disfavor the highly charged complex resulting from this proton transfer. Simulations of the Y254L variant also reveal some effect of steric changes.

Electron transfer within a reaction path model calibrated by constrained DFT calculations: application to mixed-valence organic compounds.

The quantum dynamics of electron transfer in mixed-valence organic compounds is investigated using a reaction path model calibrated by constrained density functional theory (cDFT). Constrained DFT is used to define diabatic states relevant for describing the electron transfer, to obtain equilibrium structures for each of these states and to estimate the electronic coupling between them. The harmonic analysis at the diabatic minima yields normal modes forming the dissipative bath coupled to the electronic states. In order to decrease the system-bath coupling, an effective one dimensional vibronic Hamiltonian is constructed by partitioning the modes into a linear reaction path which connects both equilibrium positions and a set of secondary vibrational modes, coupled to this reaction coordinate. Using this vibronic model Hamiltonian, dissipative quantum dynamics is carried out using Redfield theory, based on a spectral density which is determined from the cDFT results. In a first benchmark case, the model is applied to a series of mixed-valence organic compounds formed by two 1,4-dimethoxy-3-methylphenylene fragments linked by an increasing number of phenylene bridges. This allows us to examine the coherent electron transfer in extreme situations leading to a ground adiabatic state with or without a barrier and therefore to the trapping of the charge or to an easy delocalization.

Study of the docking of competitive inhibitors at a model of tyrosinase active site: insights from joint broken-symmetry/Spin-Flip DFT computations and ELF topological analysis.

Following our previous study (Piquemal et al., New J. Chem., 2003, 27, 909), we present here a DFT study of the inhibition of the Tyrosinase enzyme. Broken-symmetry DFT computations are supplemented with Spin-Flip TD-DFT calculations, which, for the first time, are applied to such a dicopper enzyme. The chosen biomimetic model encompasses a dioxygen molecule, two Cu(II) cations, and six imidazole rings. The docking energy of a natural substrate, namely phenolate, together with those of several inhibitor and non-inhibitor compounds, are reported and show the ability of the model to rank the most potent inhibitors in agreement with experimental data. With respect to broken-symmetry calculations, the Spin-Flip TD-DFT approach reinforces the possibility for theory to point out potent inhibitors: the need for the deprotonation of the substrates, natural or inhibitors, is now clearly established. Moreover, Electron Localization Function (ELF) topological analysis computations are used to deeply track the particular electronic distribution of the Cu-O-Cu three-center bonds involved in the enzymatic Cu(2)O(2) metallic core (Piquemal and Pilmé, J. Mol. Struct.: Theochem, 2006, 77, 764). It is shown that such bonds exhibit very resilient out-of-plane density expansions that play a key role in docking interactions: their 3D-orientation could be the topological electronic signature of oxygen activation within such systems.