Dataset Information

Simulation of X-ray absorption spectra with orthogonality constrained density functional theory.

ABSTRACT: Orthogonality constrained density functional theory (OCDFT) [F. A. Evangelista, P. Shushkov and J. C. Tully, J. Phys. Chem. A, 2013, 117, 7378] is a variational time-independent approach for the computation of electronic excited states. In this work we extend OCDFT to compute core-excited states and generalize the original formalism to determine multiple excited states. Benchmark computations on a set of 13 small molecules and 40 excited states show that unshifted OCDFT/B3LYP excitation energies have a mean absolute error of 1.0 eV. Contrary to time-dependent DFT, OCDFT excitation energies for first- and second-row elements are computed with near-uniform accuracy. OCDFT core excitation energies are insensitive to the choice of the functional and the amount of Hartree-Fock exchange. We show that OCDFT is a powerful tool for the assignment of X-ray absorption spectra of large molecules by simulating the gas-phase near-edge spectrum of adenine and thymine.

SUBMITTER: Derricotte WD

PROVIDER: S-EPMC4449318 | biostudies-literature | 2015 Jun

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Simulation of X-ray absorption spectra with orthogonality constrained density functional theory.

Derricotte Wallace D WD Evangelista Francesco A FA

Physical chemistry chemical physics : PCCP 20150601 22

Orthogonality constrained density functional theory (OCDFT) [F. A. Evangelista, P. Shushkov and J. C. Tully, J. Phys. Chem. A, 2013, 117, 7378] is a variational time-independent approach for the computation of electronic excited states. In this work we extend OCDFT to compute core-excited states and generalize the original formalism to determine multiple excited states. Benchmark computations on a set of 13 small molecules and 40 excited states show that unshifted OCDFT/B3LYP excitation energies ...[more]

PMID: 25690350

Similar Datasets

Project description:A detailed understanding and interpretation of absorption spectra of molecular systems, especially in condensed phases, requires computational models that allow their structural and electronic features to be connected to the observed macroscopic spectra. This work is focused on modeling the electronic absorption spectrum of a fluorescent probe, namely, the 9-(4-((bis(2-((2-(ethylthio)ethyl)thio)ethyl)amino)methyl)phenyl)-6-(pyrrolidin-1-yl)-3H-xanthen-3-one molecule, depicted by a combined classical-quantum chemical approach. Particularly, first classical molecular dynamics (MD) has been used to explore the configurational space, and next, the absorption spectrum has been reconstructed by averaging the results of time-dependent density functional theory (TD-DFT) calculations performed on equispaced molecular conformations extracted from MD to properly sample the configurational space explored at finite temperature. To verify the effect of molecular conformation on the spectral profile, the generated electronic absorption spectra were compared with those obtained considering a single structure corresponding to the optimized one, an approach also referred to as static. This comparison allows one to highlight a sizable though small shift between the maxima of the corresponding reconstructed absorption spectra, highlighting the importance of conformational sampling in the case of this rather flexible molecule. Four different exchange and correlation functionals (PBE, BLYP, PBE0, B3LYP) were considered to compute vertical transition via TD-DFT calculations. From the results obtained in gas and in condensed, here solution, phases, it appears that the magnitude of the shift is actually more affected by the phase in which the system is found than by the functional used. This fact underlines the central importance of conformational mobility, that is flexibility, of this molecule. From a more quantitative point of view, a comparison with available experimental data shows that hybrid functionals, such as PBE0 and B3LYP, enable one to faithfully reproduce the observed absorption maxima.

Project description:Superoxide reductase (SOR) is a non-heme iron enzyme that reduces superoxide to peroxide at a diffusion-controlled rate. Sulfur K-edge X-ray absorption spectroscopy (XAS) is used to investigate the ground-state electronic structure of the resting high-spin and CN- bound low-spin FeIII forms of the 1Fe SOR from Pyrococcus furiosus. A computational model with constrained imidazole rings (necessary for reproducing spin states), H-bonding interaction to the thiolate (necessary for reproducing Fe-S bond covalency of the high-spin and low-spin forms), and H-bonding to the exchangeable axial ligand (necessary to reproduce the ground state of the low-spin form) was developed and then used to investigate the enzymatic reaction mechanism. Reaction of the resting ferrous site with superoxide and protonation leading to a high-spin FeIII-OOH species and its subsequent protonation resulting in H2O2 release is calculated to be the most energetically favorable reaction pathway. Our results suggest that the thiolate acts as a covalent anionic ligand. Replacing the thiolate with a neutral noncovalent ligand makes protonation very endothermic and greatly raises the reduction potential. The covalent nature of the thiolate weakens the FeIII bond to the proximal oxygen of this hydroperoxo species, which raises its pKa by an additional 5 log units relative to the pKa of a primarily anionic ligand, facilitating its protonation. A comparison with cytochrome P450 indicates that the stronger equatorial ligand field from the porphyrin results in a low-spin FeIII-OOH species that would not be capable of efficient H2O2 release due to a spin-crossing barrier associated with formation of a high-spin 5C FeIII product. Additionally, the presence of the dianionic porphyrin pi ring in cytochrome P450 allows O-O heterolysis, forming an FeIV-oxo porphyrin radical species, which is calculated to be extremely unfavorable for the non-heme SOR ligand environment. Finally, the 5C FeIII site that results from the product release at the end of the O2- reduction cycle is calculated to be capable of reacting with a second O2-, resulting in superoxide dismutase (SOD) activity. However, in contrast to FeSOD, the 5C FeIII site of SOR, which is more positively charged, is calculated to have a high affinity for binding a sixth anionic ligand, which would inhibit its SOD activity.

Dataset Information

Simulation of X-ray absorption spectra with orthogonality constrained density functional theory.

Publications

Simulation of X-ray absorption spectra with orthogonality constrained density functional theory.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets