Unraveling the similarity of the photoabsorption of deprotonated p-coumaric acid in the gas phase and within the photoactive yellow protein

Using advanced QM/MM methods, the surprisingly negligible shift of the lowest-lying bright electronic excitation of the deprotonated p-coumaric acid (pCA(-)) within the photoactive yellow protein (PYP) is shown to stem from a subtle balance between hypsochromic and bathochromic effects. More specifically, it is found that the change in the excitation energy as a consequence of the disruption of the planarity of pCA(-) inside PYP is nearly canceled out by the shift induced by the intermolecular interactions of the chromophore and the protein as a whole. These results provide important insights about the primary absorption and the tuning of the chromophore by the protein environment in PYP.     

Quantum chemical comparison of vertical,adiabatic, and 0‐0 excitation energies: The PYP and GFP chromophores

A number of benchmark studies investigating the performance of quantum chemical methods for calculating vertical excitation energies are today available in the literature. However, less established is the variation between methods in their estimates of the differences between vertical, adiabatic, and 0‐0 excitation energies. To this end, such excitation energies are here calculated for the bright S₁ states of the anionic chromophores of the photoactive yellow protein (PYP) and the green fluorescent protein (GFP) in the gas phase using configuration interaction singles, complete active space self‐consistent field, coupled‐cluster singles and doubles, and time‐dependent density functional theory methods. Although the estimates of the excitation energies vary by more than 1 eV between the methods, the differences between the different types of excitation energies are found to be relatively method‐insensitive, varying by ～0.1 eV only for these particular chromophores. Specifically, the adiabatic energies are uniformly 0.10–0.17 (PYP) and 0.06–0.17 eV (GFP) lower than the vertical energies, and the 0‐0 energies are similarly 0.09–0.14 (PYP) and 0.07–0.17 eV (GFP) lower than the adiabatic energies. © 2012 Wiley Periodicals, Inc.     

Controlling Radical Formation in the Photoactive Yellow Protein Chromophore

To understand how photoactive proteins function, it is necessary to understand the photoresponse of the chromophore. Photoactive yellow protein (PYP) is a prototypical signaling protein. Blue light triggers trans–cis isomerization of the chromophore covalently bound within PYP as the first step in a photocycle that results in the host bacterium moving away from potentially harmful light. At higher energies, photoabsorption has the potential to create radicals and free electrons; however, this process is largely unexplored. Here, we use photoelectron spectroscopy and quantum chemistry calculations to show that the molecular structure and conformation of the isolated PYP chromophore can be exploited to control the competition between trans–cis isomerization and radical formation. We also find evidence to suggest that one of the roles of the protein is to impede radical formation in PYP by preventing torsional motion in the electronic ground state of the chromophore.     

Photoisomerization and proton transfer in photoactive yellow protein

The photoactive yellow protein (PYP) is a bacterial photosensor containing a para-coumaryl thioester chromophore that absorbs blue light, initiating a photocycle involving a series of conformational changes. Here, we present computational studies to resolve uncertainties and controversies concerning the correspondence between atomic structures and spectroscopic measurements on early photocycle intermediates. The initial nanoseconds of the PYP photocycle are examined using time-dependent density functional theory (TDDFT) to calculate the energy profiles for chromophore photoisomerization and proton transfer, and to calculate excitation energies to identify photocycle intermediates. The calculated potential energy surface for photoisomerization matches key, experimentally determined, spectral parameters. The calculated excitation energy of the photocycle intermediate cryogenically trapped in a crystal structure by Genick et al. [Genick, U. K.; Soltis, S. M.; Kuhn, P.; Canestrelli, I. L.; Getzoff, E. D. Nature 1998, 392, 206-209] supports its assignment to the PYP(B) (I(0)) intermediate. Differences between the time-resolved room temperature (298 K) spectrum of the PYP(B) intermediate and its low temperature (77 K) absorbance are attributed to a predominantly deprotonated chromophore in the former and protonated chromophore in the latter. This contrasts with the widely held belief that chromophore protonation does not occur until after the PYP(L) (I(1) or pR) intermediate. The structure of the chromophore in the PYP(L) intermediate is determined computationally and shown to be deprotonated, in agreement with experiment. Calculations based on our PYP(B) and PYP(L) models lead to insights concerning the PYP(BL) intermediate, observed only at low temperature. The results suggest that the proton is more mobile between Glu46 and the chromophore than previously realized. The findings presented here provide an example of the insights that theoretical studies can contribute to a unified analysis of experimental structures and spectra.     

Impact of the computational method on the geometric and electronic properties of oligo(phenylene vinylene) radical cations

We report on a quantum-chemical study of the electronic and optical properties of unsubstituted oligo(phenylene vinylene) (OPV) radical cations. Our goal is to distinguish the impact of the choice of molecular geometry from the impact of the choice of quantum-chemical method, on the calculated optical transition energies. The geometry modifications upon ionization of the OPV chains are found to depend critically on the theoretical formalism: Hartree-Fock (HF) geometry optimizations lead to self-localization of the charged defects while pure density functional theory (DFT) results in a complete delocalization of the geometric modifications over the whole conjugated backbone. The electronic structure and vertical transition energy associated with the lowest excited state of the radical cations have been calculated at the post-Hartree-Fock level within a configuration interaction (HF-CI) scheme and using the time-dependent DFT (TD-DFT) formalism for different radical cation geometries. Interestingly, the changes in the calculated optical properties obtained when using different geometric structures are less important within a given method than the differences between methods for a given structure. The optical excitation is localized with HF-CI and delocalized with TD-DFT, almost irrespective of the molecular geometry; as a result, HF-CI excitation energies tend to saturate as the chain length increases, in contrast to the results from TD-DFT.     

Calculating absorption shifts for retinal proteins: computational challenges

Rhodopsins can modulate the optical properties of their chromophores over a wide range of wavelengths. The mechanism for this spectral tuning is based on the response of the retinal chromophore to external stress and the interaction with the charged, polar, and polarizable amino acids of the protein environment and is connected to its large change in dipole moment upon excitation, its large electronic polarizability, and its structural flexibility. In this work, we investigate the accuracy of computational approaches for modeling changes in absorption energies with respect to changes in geometry and applied external electric fields. We illustrate the high sensitivity of absorption energies on the ground-state structure of retinal, which varies significantly with the computational method used for geometry optimization. The response to external fields, in particular to point charges which model the protein environment in combined quantum mechanical/molecular mechanical (QM/MM) applications, is a crucial feature, which is not properly represented by previously used methods, such as time-dependent density functional theory (TDDFT), complete active space self-consistent field (CASSCF), and Hartree-Fock (HF) or semiempirical configuration interaction singles (CIS). This is discussed in detail for bacteriorhodopsin (bR), a protein which blue-shifts retinal gas-phase excitation energy by about 0.5 eV. As a result of this study, we propose a procedure which combines structure optimization or molecular dynamics simulation using DFT methods with a semiempirical or ab initio multireference configuration interaction treatment of the excitation energies. Using a conventional QM/MM point charge representation of the protein environment, we obtain an absorption energy for bR of 2.34 eV. This result is already close to the experimental value of 2.18 eV, even without considering the effects of protein polarization, differential dispersion, and conformational sampling.     

Initial photoinduced dynamics of the photoactive yellow protein.

The photoactive yellow protein (PYP) is the photoreceptor protein responsible for initiating the blue-light repellent response of the Halorhodospira halophila bacterium. Optical excitation of the intrinsic chromophore in PYP, p-coumaric acid, leads to the initiation of a photocycle that comprises several distinct intermediates. The dynamical processes responsible for the initiation of the PYP photocycle have been explored with several time-resolved techniques, which include ultrafast electronic and vibrational spectroscopies. Ultrafast electronic spectroscopies, such as pump-visible probe, pump-dump-visible probe, and fluorescence upconversion, are useful in identifying the timescales and connectivity of the transient intermediates, while ultrafast vibrational spectroscopies link these intermediates to dynamic structures. Herein, we present the use of these techniques for exploring the initial dynamics of PYP, and show how these techniques provide the basis for understanding the complex relationship between protein and chromophore, which ultimately results in biological function.     

Effect of microhydration on the electronic structure of the chromophores of the photoactive yellow and green fluorescent proteins

Electronic structure calculations of microhydrated model chromophores (in their deprotonated anionic forms) of the photoactive yellow and green fluorescent proteins (PYP and GFP) are reported. Electron-detachment and excitation energies as well as binding energies of mono- and dihydrated isomers are computed and analyzed. Microhydration has different effects on the excited and ionized states. In lower-energy planar isomers, the interaction with one water molecule blueshifts the excitation energies by 0.1-0.2 eV, whereas the detachment energies increase by 0.4-0.8 eV. The important consequence is that microhydration by just one water molecule converts the resonance (autoionizing) excited states of the bare chromophores into bound states. In the lower-energy microhydrated clusters, interactions with water have negligible effect on the chromophore geometry; however, we also identified higher-energy dihydrated clusters of PYP in which two water molecules form hydrogen-bonding network connecting the carboxylate and phenolate moieties and the chromophore is strongly distorted resulting in a significant shift of excitation energies (up to 0.6 eV).     

Excited-State Electronic Structure with Configuration Interaction Singles and Tamm-Dancoff Time-Dependent Density Functional Theory on Graphical Processing Units

Excited-state calculations are implemented in a development version of the GPU-based TeraChem software package using the configuration interaction singles (CIS) and adiabatic linear response Tamm-Dancoff time-dependent density functional theory (TDA-TDDFT) methods. The speedup of the CIS and TDDFT methods using GPU-based electron repulsion integrals and density functional quadrature integration allows full ab initio excited-state calculations on molecules of unprecedented size. CIS/6-31G and TD-BLYP/6-31G benchmark timings are presented for a range of systems, including four generations of oligothiophene dendrimers, photoactive yellow protein (PYP), and the PYP chromophore solvated with 900 quantum mechanical water molecules. The effects of double and single precision integration are discussed, and mixed precision GPU integration is shown to give extremely good numerical accuracy for both CIS and TDDFT excitation energies (excitation energies within 0.0005 eV of extended double precision CPU results).     

Electronic excitations of glycine, alanine, and cysteine conformers from first-principles calculations

The electronic and optical properties are studied for three conformers of amino acid molecules using gradient-corrected (spin-) density functional theory within a projector-augmented wave scheme and the supercell method. We investigate single-particle excitations such as ionization energies and electron affinities as well as pair excitations. By comparing eigenvalues resulting from several local and nonlocal energy functionals, the influence of treatment of exchange and correlation is demonstrated. The excitations are described within the Delta-self-consistent field method with an occupation number constraint to obtain excitation energies and Stokes shifts. The results are used to also discuss the optical absorption properties. In contrast to the lowest single- and two-particle excitation energies, remarkable changes are found in absorption spectra in dependence on the conformation of the molecule geometry.     

Optoelectronic Properties of Carbon Nanorings: Excitonic Effects from Time-Dependent Density Functional Theory

The electronic structure and size-scaling of optoelectronic properties in cycloparaphenylene carbon nanorings are investigated using time-dependent density functional theory (TDDFT). The TDDFT calculations on these molecular nanostructures indicate that the lowest excitation energy surprisingly becomes larger as the carbon nanoring size is increased, in contradiction with typical quantum confinement effects. In order to understand their unusual electronic properties, I performed an extensive investigation of excitonic effects by analyzing electron-hole transition density matrices and exciton binding energies as a function of size in these nanoring systems. The transition density matrices allow a global view of electronic coherence during an electronic excitation, and the exciton binding energies give a quantitative measure of electron-hole interaction energies in the nanorings. Based on overall trends in exciton binding energies and their spatial delocalization, I find that excitonic effects play a vital role in understanding the unique photoinduced dynamics in these carbon nanoring systems.     

Tuning the photophysical properties of pyrene-based systems: a theoretical study

Recently new molecular systems based on the pyrene moiety were developed for photovoltaic applications. Here we present the results of a quantum chemical study focused on the effects induced by some different substituents on the electronic properties of pyrene, to obtain general hints for the molecular design of new pyrene-based systems. In particular, a series of electron-donating (hydroxy, amino, acetylamino) and electron-withdrawing (cyano, carbamoyl, formyl, ethynyl, ethenyl) groups were considered. Furthermore, in addition to the single pyrene molecule, two pyrene units linked by ethenylene, ethynylene, 2,5-thienylene, and ethynylene-p-phenylene containing chains of different lengths were taken into account. For all of the model structures presented, the ground state geometries have been optimized using the density functional approach, while the vertical transition energies were calculated using the time-dependent density functional theory. We will show that the tuning of the lowest electronic excitation energy (i.e., the HOMO-LUMO energy gap) as well as the localization of the spatial distributions of the frontier molecular orbitals (i.e., the nature of the electron-hole pair, generated by photon absorption) can be obtained through the analysis of the pyrene frontier molecular orbitals. This approach allows to evaluate the most suitable position of the substituents on the pyrene moiety giving rise to enhanced electronic effects also in function of their electronic nature. In this way, pyrene-structures with tailored electronic properties could be modeled. Our screening shows that promising candidates for photovoltaic applications could be molecular structures formed by two pyrene units joined/linked by a short conjugated bridge containing double or triple bonds (henceforth pyrene-linked dimers). As far as the single pyrene units are considered, the most significant reduction of the transition energy of the lowest optical electronic excitation is obtained with disubstituted pyrenes with push-pull character.     

Structural evolution of the chromophore in the primary stages of trans/cis isomerization in photoactive yellow protein

We have studied the structural changes induced by optical excitation of the chromophore in wild-type photoactive yellow protein (PYP) in liquid solution with a combined approach of polarization-sensitive ultrafast infrared spectroscopy and density functional theory calculations. We identify the nuC8-C9 marker modes for solution phase PYP in the P and I0 states, from which we derive that the first intermediate state I0 that appears with a 3 ps time constant can be characterized to have a cis geometry. This is the first unequivocal demonstration that the formation of I0 correlates with the conversion from the trans to the cis state. For the P and I0 states we compare the experimentally measured vibrational band patterns and anisotropies with calculations and find that for both trans and cis configurations the planarity of the chromophore has a strong influence. The C7=C8-(C9=O)-S moiety of the chromophore in the dark P state has a trans geometry with the C=O group slightly tilted out-of-plane, in accordance with the earlier reported structure obtained in an X-ray diffraction study of PYP crystals. In the case of I0, experiment and theory are only in agreement when the C7=C8-(C9=O)-S moiety has a planar configuration. We find that the carboxylic side group of Glu46 that is hydrogen-bonded to the chromophore phenolate oxygen does not alter its orientation on going from the electronic ground P state, via the electronic excited P state to the intermediate I0 state, providing conclusive experimental evidence that the primary stages of PYP photoisomerization involve flipping of the enone thioester linkage without significant relocation of the phenolate moiety.     

Spectral tuning in visual pigments: an ONIOM(QM:MM) study on bovine rhodopsin and its mutants

We have investigated geometries and excitation energies of bovine rhodopsin and some of its mutants by hybrid quantum mechanical/molecular mechanical (QM/MM) calculations in ONIOM scheme, employing B3LYP and BLYP density functionals as well as DFTB method for the QM part and AMBER force field for the MM part. QM/MM geometries of the protonated Schiff-base 11- cis-retinal with B3LYP and DFTB are very similar to each other. TD-B3LYP/MM excitation energy calculations reproduce the experimental absorption maximum of 500 nm in the presence of native rhodopsin environment and predict spectral shifts due to mutations within 10 nm, whereas TD-BLYP/MM excitation energies have red-shift error of at least 50 nm. In the wild-type rhodopsin, Glu113 shifts the first excitation energy to blue and accounts for most of the shift found. Other amino acids individually contribute to the first excitation energy but their net effect is small. The electronic polarization effect is essential for reproducing experimental bond length alternation along the polyene chain in protonated Schiff-base retinal, which correlates with the computed first excitation energy. It also corrects the excitation energies and spectral shifts in mutants, more effectively for deprotonated Schiff-base retinal than for the protonated form. The protonation state and conformation of mutated residues affect electronic spectrum significantly. The present QM/MM calculations estimate not only the experimental excitation energies but also the source of spectral shifts in mutants.     

Photo absorption of ‐coumaric acid in aqueous solution: RISM‐SCF‐SEDD theory approach

Photo absorption properties of p‐coumaric acid, the chromophore of photoactive yellow protein, in aqueous solution were investigated by means of reference interaction site model self‐consistent field with spatial electron density distribution (RISM‐SCF‐SEDD) method. RISM‐SCF‐SEDD is a combination methodology of electronic structure theory and statistical mechanics for molecular liquids. Here, time‐dependent density functional theory was coupled with RISM equation to study the electronic structure of p‐coumaric acid in aqueous system. Excitation energies of the chromophore in its neutral, two monoanionic and dianionic forms were computed to elucidate the effect of the deprotonation and solvation on the spectroscopic properties. We found that solvation strongly affects the excitation character of the chromophore, especially for phenolate anion and dianion. The free energy difference among the four protonation states is also discussed. © 2017 Wiley Periodicals, Inc.     

Protonation of the chromophore in the photoactive yellow protein

The photoactive yellow protein (PYP) acts as a light sensor to its bacterial host: it responds to light by changing shape. After excitation by blue light, PYP undergoes several transformations, to partially unfold into its signaling state. One of the crucial steps in this photocycle is the protonation of p-coumaric acid after excitation and isomerization of this chromophore. Experimentalists still debate on the nature of the proton donor and on whether it donates the hydrogen directly or indirectly. To obtain better knowledge of the mechanism, we studied this proton transfer using Car-Parrinello molecular dynamics, classical molecular dynamics, and computer simulations combining these two methods (quantum mechanics/molecular mechanics, QMMM). The simulations reproduce the chromophore structure and hydrogen-bond network of the protein measured by X-ray crystallography and NMR. When the chromophore is protonated, it leaves the assumed proton donor, glutamic acid 46, with a negative charge in a hydrophobic environment. We show that the stabilization of this charge is a very important factor in the mechanism of protonation. Protonation frequently occurs in simplified ab initio simulations of the chromophore binding pocket in vacuum, where amino acids can easily hydrogen bond to Glu46. When the complete protein environment is incorporated in a QMMM simulation on the complete protein, no proton transfer is observed within 14 ps. The hydrogen-bond rearrangements in this time span are not sufficient to stabilize the new protonation state. Force field molecular dynamics simulations on a much longer time scale have shown which internal rearrangements of the protein are needed. Combining these simulations with more QMMM calculations enabled us to check the stability of protonation states and clarify the initial requirements for the proton transfer in PYP.     

Numerical simulation of nonadiabatic electron excitation in the strong field regime. 2. Linear polyene cations

Time-dependent Hartree-Fock theory has been used to study the electronic optical response of a series of linear polyene cations (+1 and +2) in strong laser fields. The interaction of ethylene, butadiene, and hexatriene, with pulsed and CW fields corresponding to 8.75 x 10(13) W/cm(2) and 760 nm, have been calculated using the 6-31G(d,p) basis set. Nonadiabatic processes including nonlinear response of the dipole moment to the field and non-resonant energy deposition into excited states were more pronounced for the monocations in comparison with dications. For a given charge state and geometry, the nonadiabatic effects in the charge distribution and instantaneous dipole increased with the length of the polyene. For pulsed fields, the instantaneous dipole continued to oscillate after the field returned to zero and corresponded to a non-resonant electronic excitation involving primarily the lowest electronic transition. For a given molecule and fixed charge state, the degree of nonadiabatic coupling and excitation was greater for geometries with lower excitation energies.     

Excited state properties, fluorescence energies, and lifetime of a poly(fluorene-pyridine) copolymer, based on TD-DFT investigation

The structural and electronic properties of the fluorene-pyridine copolymer (FPy)(n), (n = 1-4) were investigated theoretically by means of quantum mechanical calculations based on density functional theory (DFT) and time-dependent DFT (TD-DFT) using the B3LYP functional. Geometry optimizations of these oligomers were performed for the ground state and the lowest excited state. It was found that (FPy)(n) is nonplanar in its ground state, whereas a more pronounced trend toward planarity is observed in the S(1) state. Absorption and fluorescence energies have been extrapolated to infinite chain length making use of their good linearity with respect to 1/n. An extrapolated value of 2.64 eV is obtained for vertical excitation energy. The S(1)<--S(0) electronic excitation is characterized as a highest occupied molecular orbital to lowest unoccupied molecular orbital transition and is dominating in terms of oscillator strength. Fluorescence energies and radiative lifetime were calculated as well. The obtained results indicate that the fluorescence energy and radiative lifetime of (FPy)(n) are 2.16 eV and 0.38 ns, respectively. The decrease of fluorescence energy and radiative lifetime with the increase in the chain length is discussed.