Theoretical Studies on Excited States of a Phenolate Anion in the Environment of Photoactive Yellow Protein

QM/MM Study of Substituent and Solvent Effects on the Excited State Dynamics of the Photoactive Yellow Protein Chromophore

Substituent and solvent effects on the excited state dynamics of the Photoactive Yellow Protein chromophore are studied using the average solvent electrostatic potential from molecular dynamics (ASEP/MD) method. Four molecular models were considered: the ester and thioester derivatives of the p-coumaric acid anion and their methylated derivatives. We found that the solvent produces dramatic modifications on the free energy profile of the S1 state: 1) Two twisted structures that are minima in the gas phase could not be located in aqueous solution. 2) Conical intersections (CIs) associated with the rotation of the single bond adjacent to the phenyl group are found for the four derivatives in water solution but only for thio derivatives in the gas phase. 3) The relative stability of minima and CIs is reverted with respect to the gas phase values, affecting the prevalent de-excitation paths. As a consequence of these changes, three competitive de-excitation channels are open in aqueous solution: the fluorescence emission from a planar minimum on S1, the trans-cis photoisomerization through a CI that involves the rotation of the vinyl double bond, and the nonradiative, nonreactive, de-excitation through the CI associated with the rotation of the single bond adjacent to the phenyl group. In the gas phase, the minima are the structures with the lower energy, while in solution these are the conical intersections. In solution, the de-excitation prevalent path seems to be the photoisomerization for oxo compounds, while thio compounds return to the initial trans ground state without emission.

Solvent effects on de-excitation channels in the p-coumaric acid methyl ester anion, an analogue of the photoactive yellow protein (PYP) chromophore

Environmental effects on the deactivation channels of the PYP chromophore in the gas phase and water solution are compared at the CASPT2//CASSCF/cc-pVDZ level.

Electronic Absorption Spectra from MM and ab initio QM/MM Molecular Dynamics: Environmental Effects on the Absorption Spectrum of Photoactive Yellow Protein

We describe a new interface of the GPU parallelized TeraChem electronic structure package and the Amber molecular dynamics package for quantum mechanical (QM) and mixed QM and molecular mechanical (MM) molecular dynamics simulations. This QM/MM interface is used for computation of the absorption spectra of the photoactive yellow protein (PYP) chromophore in vacuum, aqueous solution, and protein environments. The computed excitation energies of PYP require a very large QM region (hundreds of atoms) covalently bonded to the chromophore in order to achieve agreement with calculations that treat the entire protein quantum mechanically. We also show that 40 or more surrounding water molecules must be included in the QM region in order to obtain converged excitation energies of the solvated PYP chromophore. These results indicate that large QM regions (with hundreds of atoms) are a necessity in QM/MM calculations.

Electronic Structure of the PYP Chromophore in Its Native Protein Environment

We report on supermolecular ab initio calculations which clarify the role of the local amino acid environment in determining the unique electronic structure properties of the photoactive yellow protein (PYP) chromophore. The extensive ab initio calculations, at the level of the CC2 and EOM-CCSD methods, allow us to explicitly address how the interactions between the deprotonated p-coumaric thio-methyl ester (pCTM-) chromophore and the surrounding amino acids act together to create a specifically stabilized pCTM- species. Particularly noteworthy is the role of the Arg52 amino acid in stabilizing the chromophore against autoionization, and the role of the Tyr42 and Glu46 amino acids in determining the hydrogen-bonding properties that carry the dominant energetic effects.

Absorption spectra of photoactive yellow protein chromophores in vacuum

The absorption spectra of two photoactive yellow protein model chromophores have been measured in vacuum using an electrostatic ion storage ring. The absorption spectrum of the isolated chromophore is an important reference for deducing the influence of the protein environment on the electronic energy levels of the chromophore and separating the intrinsic properties of the chromophore from properties induced by the protein environment. In vacuum the deprotonated trans-thiophenyl-p-coumarate model chromophore has an absorption maximum at 460 nm, whereas the photoactive yellow protein absorbs maximally at 446 nm. The protein environment thus only slightly blue-shifts the absorption. In contrast, the absorption of the model chromophore in aqueous solution is significantly blue-shifted (lambda(max) = 395 nm). A deprotonated trans-p-coumaric acid has also been studied to elucidate the effect of thioester formation and phenol deprotonation. The sum of these two changes on the chromophore induces a red shift both in vacuum and in aqueous solution.

Impact of Sulfur vs Oxygen on the Low-Lying Excited States oftrans-p-Coumaric Acid andtrans-p-Coumaric Thio Acid

The low-lying excited singlet states of trans-p-coumaric acid (CA) and trans-p-coumaric thio acid (CTA) are investigated in view of characterizing the chromophore of the photoactive yellow protein (PYP), with particular regard to the impact of sulfur on the chromophore's electronic structure. The comparative ab initio study, performed with the highly accurate EOM-CCSD method, shows that the electronic state ordering upon vertical excitation and following in-plane geometry relaxation indeed depends in a very sensitive fashion on the presence of either sulfur or oxygen. The study identifies three relevant excited singlet states, two of which are of pi-pi type while the third state is of n-pi character. The study highlights the role of the latter n-pi state which is shown to be the lowest-lying excited state of CTA at all in-plane geometries under consideration, whereas this is not the case for CA.

Ultrafast Dynamics of Isolated Model Photoactive Yellow Protein Chromophores: “Chemical Perturbation Theory” in the Laboratory

Pump-probe and pump-dump probe experiments have been performed on several isolated model chromophores of the photoactive yellow protein (PYP). The observed transient absorption spectra are discussed in terms of the spectral signatures ascribed to solvation, excited-state twisting, and vibrational relaxation. It is observed that the protonation state has a profound effect on the excited-state lifetime of p-coumaric acid. Pigments with ester groups on the coumaryl tail end and charged phenolic moieties show dynamics that are significantly different from those of other pigments. Here, an unrelaxed ground-state intermediate could be observed in pump-probe signals. A similar intermediate could be identified in the sinapinic acid and in isomerization-locked chromophores by means of pump-dump probe spectroscopy; however, in these compounds it is less pronounced and could be due to ground-state solvation and/or vibrational relaxation. Because of strong protonation-state dependencies and the effect of electron donor groups, it is argued that charge redistribution upon excitation determines the twisting reaction pathway, possibly through interaction with the environment. It is suggested that the same pathway may be responsible for the initiation of the photocycle in native PYP.

Electronic excitations of biomolecules studied by quantum chemistry

Significant methodological advances have been made over the past ten years in developing reliable quantum chemical methods for the treatment of electronically excited states. These methods can nowadays be used routinely by the experienced researcher to accurately compute excitation spectra of medium-sized organic molecules; results have been reported for several popular photobiological systems, including green fluorescent protein. First steps are currently being taken to account for the solvochromic shifts of chromophore excitations caused by particular protein environments and to dynamically simulate photochemical reactions in the excited states.

On the absorbance changes in the photocycle of the photoactive yellow protein: a quantum-chemical analysis

Spectral changes in the photocycle of the photoactive yellow protein (PYP) are investigated by using ab initio multiconfigurational second-order perturbation theory at the available structures experimentally determined. Using the dark ground-state crystal structure [Genick, U. K., Soltis, S. M., Kuhn, P., Canestrelli, I. L. & Getzoff, E. D. (1998) Nature (London) 392, 206-209], the pipi* transition to the lowest excited state is related to the typical blue-light absorption observed at 446 nm. The different nature of the second excited state (npi*) is consistent with the alternative route detected at 395-nm excitation. The results suggest the low-temperature photoproduct PYP(HL) as the most plausible candidate for the assignment of the cryogenically trapped early intermediate (Genick et al.). We cannot establish, however, a successful correspondence between the theoretical spectrum for the nanosecond time-resolved x-ray structure [Perman, B., Srajer, V., Ren, Z., Teng, T., Pradervand, C., et al. (1998) Science 279, 1946-1950] and any of the spectroscopic photoproducts known up to date. It is fully confirmed that the colorless light-activated intermediate recorded by millisecond time-resolved crystallography [Genick, U. K., Borgstahl, G. E. O., Ng, K., Ren, Z., Pradervand, C., et al. (1997) Science 275, 1471-1475] is protonated, nicely matching the spectroscopic features of the photoproduct PYP(M). The overall contribution demonstrates that a combined analysis of high-level theoretical results and experimental data can be of great value to perform assignments of detected intermediates in a photocycle.