Adalgisa Sinicropi

I - Computational Photochemistry

This chapter discusses photochemical reaction path concept and its use in mechanistic investigations. The field of computational photochemistry is a relatively young field, especially when applied to the study of ultrafast reactions, but it is now established as a branch of computational chemistry and as a powerful, sometimes unique, way to simulate the molecular mechanism underlying fundamental chemical and biological events such as vision, primitive photosynthesis, phototropism, photochromism, bleaching, fluorescence, phosphorescence. These days, computational strategies are available for locating conical intersection and singlet/triplet crossing points and for constructing inter-state “photochemical” reaction pathways. These tools comprise methodologies for the optimisation of low-lying crossings between pair of potential energy surfaces and the computation of relaxation paths from a photoexcited reactant (For example, from theFranck-Condon (FC) structure) to a deactivation channel.

An artificial molecular switch that mimics the visual pigment and completes its photocycle in picoseconds

Single molecules that act as light-energy transducers (e.g., converting the energy of a photon into atomic-level mechanical motion) are examples of minimal molecular devices. Here, we focus on a molecular switch designed by merging a conformationally locked diarylidene skeleton with a retinal-like Schiff base and capable of mimicking, in solution, different aspects of the transduction of the visual pigment Rhodopsin. Complementary ab initio multiconfigurational quantum chemistry-based computations and time-resolved spectroscopy are used to follow the light-induced isomerization of the switch in methanol. The results show that, similar to rhodopsin, the isomerization occurs on a 0.3-ps time scale and is followed by <10-ps cooling and solvation. The entire (2-photon-powered) switch cycle was traced by following the evolution of its infrared spectrum. These measurements indicate that a full cycle can be completed within 20 ps.

Bathochromic Shift in Green Fluorescent Protein: A Puzzle for QM/MM Approaches

We present an extensive investigation of the vertical excitations of the anionic and neutral forms of wild-type green fluorescent protein using time-dependent density functional theory (TDDFT), multiconfigurational perturbation theory (CASPT2), and quantum Monte Carlo (QMC) methods within a quantum mechanics/molecular mechanics (QM/MM) scheme. The protein models are constructed via room-temperature QM/MM molecular dynamics simulations based on DFT and are representative of an average configuration of the chromophore-protein complex. We thoroughly verify the reliability of our structures through simulations with an extended QM region, different nonpolarizable force fields, as well as partial reoptimization with the CASPT2 approach. When computing the excitations, we find that wave function as well as density functional theory methods with long-range corrected functionals agree in the gas phase with the extrapolation of solution experiments but fail in reproducing the bathochromic shift in the protein, which should be particularly significant in the neutral case. In particular, while all methods correctly predict a shift in the absorption between the anionic and neutral forms of the protein, the location of the theoretical absorption maxima is significantly blue-shifted and too close to the gas-phase values. These results point to either an intrinsic limitation of nonpolarizable force-field embedding in the computation of the excitations or to the need to explore alternative protonation states of amino acids in the close vicinity of the chomophore.

Chromophore–Protein Coupling beyond Nonpolarizable Models: Understanding Absorption in Green Fluorescent Protein

The nature of the coupling of the photoexcited chromophore with the environment in a prototypical system like green fluorescent protein (GFP) is to date not understood, and its description still defies state-of-the-art multiscale approaches. To identify which theoretical framework of the chromophore-protein complex can realistically capture its essence, we employ here a variety of electronic-structure methods, namely, time-dependent density functional theory (TD-DFT), multireference perturbation theory (NEVPT2 and CASPT2), and quantum Monte Carlo (QMC) in combination with static point charges (QM/MM), DFT embedding (QM/DFT), and classical polarizable embedding through induced dipoles (QM/MMpol). Since structural modifications can significantly affect the photophysics of GFP, we also account for thermal fluctuations through extensive molecular dynamics simulations. We find that a treatment of the protein through static point charges leads to significantly blue-shifted excitation energies and that including thermal fluctuations does not cure the coarseness of the MM description. While TDDFT calculations on large cluster models indicate the need of a responsive protein, this response is not simply electrostatic: An improved description of the protein in the ground state or in response to the excitation of the chromophore via ground-state or state-specific DFT and MMpol embedding does not significantly modify the results obtained with static point charges. Through the use of QM/MMpol in a linear response formulation, a different picture in fact emerges in which the main environment response to the chromophore excitation is the one coupling the transition density and the corresponding induced dipoles. Such interaction leads to significant red-shifts and a satisfactory agreement with full QM cluster calculations at the same level of theory. Our findings demonstrate that, ultimately, faithfully capturing the effects of the environment in GFP requires a quantum treatment of large photoexcited regions but that a QM/classical model can be a useful approximation when extended beyond the electrostatic-only formulation.

Catalytic surface radical in dye-decolorizing peroxidase: a computational, spectroscopic and site-directed mutagenesis study

Dye-decolorizing peroxidase (DyP) of Auricularia auricula-judae has been expressed in Escherichia coli as a representative of a new DyP family, and subjected to mutagenic, spectroscopic, crystallographic and computational studies. The crystal structure of DyP shows a buried haem cofactor, and surface tryptophan and tyrosine residues potentially involved in long-range electron transfer from bulky dyes. Simulations using PELE (Protein Energy Landscape Exploration) software provided several binding-energy optima for the anthraquinone-type RB19 (Reactive Blue 19) near the above aromatic residues and the haem access-channel. Subsequent QM/MM (quantum mechanics/molecular mechanics) calculations showed a higher tendency of Trp-377 than other exposed haem-neighbouring residues to harbour a catalytic protein radical, and identified the electron-transfer pathway. The existence of such a radical in H2O2-activated DyP was shown by low-temperature EPR, being identified as a mixed tryptophanyl/tyrosyl radical in multifrequency experiments. The signal was dominated by the Trp-377 neutral radical contribution, which disappeared in the W377S variant, and included a tyrosyl contribution assigned to Tyr-337 after analysing the W377S spectra. Kinetics of substrate oxidation by DyP suggests the existence of high- and low-turnover sites. The high-turnover site for oxidation of RB19 (kcat> 200 s−1) and other DyP substrates was assigned to Trp-377 since it was absent from the W377S variant. The low-turnover site/s (RB19 kcat ~20 s−1) could correspond to the haem access-channel, since activity was decreased when the haem channel was occluded by the G169L mutation. If a tyrosine residue is also involved, it will be different from Tyr-337 since all activities are largely unaffected in the Y337S variant.

Structure of the intersection space associated with Z/E photoisomerization of retinal in rhodopsin proteins

In this paper we employ a CASSCF/AMBER quantum-mechanics/molecular-mechanics tool to map the intersection space (IS) of a protein. In particular, we provide evidence that the S1 excited-state potential-energy surface of the visual photoreceptor rhodopsin is spanned by an IS segment located right at the bottom of the surface. Analysis of the molecular structures of the protein chromophore (a protonated Schiff base of retinal) along IS reveals a type of geometrical deformation not observed in vacuo. Such a structure suggests that conical intersections mediating different photochemical reactions reside along the same intersection space. This conjecture is investigated by mapping the intersection space of the rhodopsin chromophore model 2-Z-hepta-2,4,6-trieniminium cation and of the conjugated hydrocarbon 3-Z-deca-1,3,5,6,7-pentaene.

Modeling, Preparation, and Characterization of a Dipole Moment Switch Driven by Z/E Photoisomerization

We report the results of a multidisciplinary research effort where the methods of computational photochemistry and retrosynthetic analysis/synthesis have contributed to the preparation of a novel N-alkylated indanylidene-pyrroline Schiff base featuring an exocyclic double bond and a permanent zwitterionic head. We show that, due to its large dipole moment and efficient photoisomerization, such a system may constitute the prototype of a novel generation of electrostatic switches achieving a reversible light-induced dipole moment change on the order of 30 D. The modeling of a peptide fragment incorporating the zwitterionic head into a conformationally rigid side chain shows that the switch can effectively modulate the fluorescence of a tryptophan probe.

Excited State Geometries and Vertical Emission Energies of Solvated Dyes for DSSC: A PCM/TD-DFT Benchmark Study

The ability of Time-Dependent Density Functional Theory (TD-DFT) to provide excited state geometries and reproduce emission energies of organic D-π-A dyes designed for DSSC applications is evaluated. The performance of six functionals (CAM-B3LYP, MPW1K, ωB97X-D, LC-BLYP, LC-ωPBE, and M06-HF) in combination with three basis sets (cc-pVDZ, 6-31+G(d,p), and 6-311+G(2d,p)) has been analyzed. Solvent effects have been taken into account by means of a Polarizable Continuum Model in both LR and SS formalisms. Our LR-PCM/TD-DFT results show that accurate emission energies are obtained only when solvent effects are included in the computation of excited state geometries and when a range separated hybrid functional is used. Vertical emission energies are reproduced with a mean absolute error of at most 0.2 eV. The accuracy is further improved using the SS-PCM formalism.

Modeling the Fluorescence of Protein-Embedded Tryptophans with ab Initio Multiconfigurational Quantum Chemistry: The Limiting Cases of Parvalbumin and Monellin

We show that a quantum-mechanics/molecular-mechanics strategy based on ab initio (i.e., first principle) multiconfigurational perturbation theory can reproduce the spectral properties of a tryptophan residue embedded in the contrasting hydrophobic and hydrophilic environments of parvalbumin and monellin, respectively. We show that the observed absorption and emission energies can be reproduced with a less than 3 kcal mol(-1) error. The analysis of the computed emission energies based on a protein disassembly scheme and protein electrostatic potential mapping allows for a detailed understanding of the factors modulating the tryptophan emission. It is shown that for monellin, where the tryptophan is exposed to the solvent, the fluorescence wavelength is controlled not only by the distribution of the point charges of the protein-solvent environment but also by specific hydrogen bonds and, most important, by the environment-induced change in chromophore structure. In contrast, in parvalbumin, where the chromophore is embedded in the protein core, the structure and emission maxima are the same as those of an isolated 3-methylindole fluorophore. Consistently, we find that in parvalbumin the solvation does not change significantly the computed emission energy.

Characterization of the conical intersection of the visual pigment rhodopsin at the CASPT2//CASSCF/AMBER level of theory

The branching plane associated with the conical intersection controlling the photochemical Z → E isomerization reaction of rhodopsin has been mapped using a CASPT2//CASSCF/AMBER quantum mechanics–molecular mechanics method. The nature of the derivative coupling and gradient difference vectors spanning the branching plane has been investigated, showing that the conical intersection is not only associated with the isomerization process but also to a charge transfer along the retinal backbone. Using a simple Landau–Zener model, the paper discusses the possible effects of the documented conical intersection topologies on the efficiency of the reactive process. It is argued that the peculiar shape of the conical intersection favours decay at structures that are geometrically displaced towards the photoproduct bathorhodopsin.