1
|
Wiethorn ZR, Hunter KE, Zuehlsdorff TJ, Montoya-Castillo A. Beyond the Condon limit: Condensed phase optical spectra from atomistic simulations. J Chem Phys 2023; 159:244114. [PMID: 38153146 DOI: 10.1063/5.0180405] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/09/2023] [Accepted: 12/06/2023] [Indexed: 12/29/2023] Open
Abstract
While dark transitions made bright by molecular motions determine the optoelectronic properties of many materials, simulating such non-Condon effects in condensed phase spectroscopy remains a fundamental challenge. We derive a Gaussian theory to predict and analyze condensed phase optical spectra beyond the Condon limit. Our theory introduces novel quantities that encode how nuclear motions modulate the energy gap and transition dipole of electronic transitions in the form of spectral densities. By formulating the theory through a statistical framework of thermal averages and fluctuations, we circumvent the limitations of widely used microscopically harmonic theories, allowing us to tackle systems with generally anharmonic atomistic interactions and non-Condon fluctuations of arbitrary strength. We show how to calculate these spectral densities using first-principles simulations, capturing realistic molecular interactions and incorporating finite-temperature, disorder, and dynamical effects. Our theory accurately predicts the spectra of systems known to exhibit strong non-Condon effects (phenolate in various solvents) and reveals distinct mechanisms for electronic peak splitting: timescale separation of modes that tune non-Condon effects and spectral interference from correlated energy gap and transition dipole fluctuations. We further introduce analysis tools to identify how intramolecular vibrations, solute-solvent interactions, and environmental polarization effects impact dark transitions. Moreover, we prove an upper bound on the strength of cross correlated energy gap and transition dipole fluctuations, thereby elucidating a simple condition that a system must follow for our theory to accurately predict its spectrum.
Collapse
Affiliation(s)
- Zachary R Wiethorn
- Department of Chemistry, University of Colorado Boulder, Boulder, Colorado 80309, USA
| | - Kye E Hunter
- Department of Chemistry, Oregon State University, Corvallis, Oregon 97331, USA
| | - Tim J Zuehlsdorff
- Department of Chemistry, Oregon State University, Corvallis, Oregon 97331, USA
| | | |
Collapse
|
2
|
Ashworth EK, Kao MH, Anstöter CS, Riesco-Llach G, Blancafort L, Solntsev KM, Meech SR, Verlet JRR, Bull JN. Alkylated green fluorescent protein chromophores: dynamics in the gas phase and in aqueous solution. Phys Chem Chem Phys 2023; 25:23626-23636. [PMID: 37649445 DOI: 10.1039/d3cp03250g] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/01/2023]
Abstract
Fluorescent labelling of macromolecular samples, including using the green fluorescent protein (GFP), has revolutionised the field of bioimaging. The ongoing development of fluorescent proteins require a detailed understanding of the photophysics of the biochromophore, and how chemical derivatisation influences the excited state dynamics. Here, we investigate the photophysical properties associated with the S1 state of three alkylated derivatives of the chromophore in GFP, in the gas phase using time-resolved photoelectron imaging, and in water using femtosecond fluorescence upconversion. The gas-phase lifetimes (1.6-10 ps), which are associated with the intrinsic (environment independent) dynamics, are substantially longer than the lifetimes in water (0.06-3 ps), attributed to stabilisation of both twisted intermediate structures and conical intersection seams in the condensed phase. In the gas phase, alkylation on the 3 and 5 positions of the phenyl ring slows the dynamics due to inertial effects, while a 'pre-twist' of the methine bridge through alkylation on the 2 and 6 positions significantly shortens the excited state lifetimes. Formation of a minor, long-lived (≫ 40 ps) excited state population in the gas phase is attributed to intersystem crossing to a triplet state, accessed because of a T1/S1 degeneracy in the so-called P-trap potential energy minimum associated with torsion of the single-bond in the bridging unit connecting to the phenoxide ring. A small amount of intersystem crossing is supported through TD-DFT molecular dynamics trajectories and MS-CASPT2 calculations. No such intersystem crossing occurs in water at T = 300 K or in ethanol at T ≈ 77 K, due to a significantly altered potential energy surface and P-trap geometry.
Collapse
Affiliation(s)
- Eleanor K Ashworth
- School of Chemistry, Norwich Research Park, University of East Anglia, Norwich, NR4 7TJ, UK.
| | - Min-Hsien Kao
- School of Chemistry, Norwich Research Park, University of East Anglia, Norwich, NR4 7TJ, UK.
| | - Cate S Anstöter
- Department of Chemistry, Durham University, Durham, DH1 3LE, UK
| | - Gerard Riesco-Llach
- Institut de Química Computacional i Catàlisi and Departament de Química, Universitat de Girona, C/M.A. Capmany 69, 17003 Girona, Spain
| | - Lluís Blancafort
- Institut de Química Computacional i Catàlisi and Departament de Química, Universitat de Girona, C/M.A. Capmany 69, 17003 Girona, Spain
| | - Kyril M Solntsev
- School of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30332, USA
| | - Stephen R Meech
- School of Chemistry, Norwich Research Park, University of East Anglia, Norwich, NR4 7TJ, UK.
| | - Jan R R Verlet
- Department of Chemistry, Durham University, Durham, DH1 3LE, UK
| | - James N Bull
- School of Chemistry, Norwich Research Park, University of East Anglia, Norwich, NR4 7TJ, UK.
| |
Collapse
|
3
|
Chen MS, Mao Y, Snider A, Gupta P, Montoya-Castillo A, Zuehlsdorff TJ, Isborn CM, Markland TE. Elucidating the Role of Hydrogen Bonding in the Optical Spectroscopy of the Solvated Green Fluorescent Protein Chromophore: Using Machine Learning to Establish the Importance of High-Level Electronic Structure. J Phys Chem Lett 2023; 14:6610-6619. [PMID: 37459252 DOI: 10.1021/acs.jpclett.3c01444] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/28/2023]
Abstract
Hydrogen bonding interactions with chromophores in chemical and biological environments play a key role in determining their electronic absorption and relaxation processes, which are manifested in their linear and multidimensional optical spectra. For chromophores in the condensed phase, the large number of atoms needed to simulate the environment has traditionally prohibited the use of high-level excited-state electronic structure methods. By leveraging transfer learning, we show how to construct machine-learned models to accurately predict the high-level excitation energies of a chromophore in solution from only 400 high-level calculations. We show that when the electronic excitations of the green fluorescent protein chromophore in water are treated using EOM-CCSD embedded in a DFT description of the solvent the optical spectrum is correctly captured and that this improvement arises from correctly treating the coupling of the electronic transition to electric fields, which leads to a larger response upon hydrogen bonding between the chromophore and water.
Collapse
Affiliation(s)
- Michael S Chen
- Department of Chemistry, Stanford University, Stanford, California 94305, United States
| | - Yuezhi Mao
- Department of Chemistry, Stanford University, Stanford, California 94305, United States
| | - Andrew Snider
- Chemistry and Biochemistry, University of California Merced, Merced, California 95343, United States
| | - Prachi Gupta
- Chemistry and Biochemistry, University of California Merced, Merced, California 95343, United States
| | - Andrés Montoya-Castillo
- Department of Chemistry, University of Colorado, Boulder, Boulder, Colorado 80309, United States
| | - Tim J Zuehlsdorff
- Department of Chemistry, Oregon State University, Corvallis, Oregon 97331, United States
| | - Christine M Isborn
- Chemistry and Biochemistry, University of California Merced, Merced, California 95343, United States
| | - Thomas E Markland
- Department of Chemistry, Stanford University, Stanford, California 94305, United States
| |
Collapse
|
4
|
Zhu C. Franck‐Condon
factors within damped displacement harmonic oscillators: Solvent‐enhanced absorption and fluorescence spectra. J CHIN CHEM SOC-TAIP 2022. [DOI: 10.1002/jccs.202200341] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/17/2022]
Affiliation(s)
- Chaoyuan Zhu
- Department of Applied Chemistry and Institute of Molecular Science National Yang Ming Chiao Tung University Hsinchu Taiwan
- Key Laboratory of Theoretical Chemistry of Environment, Ministry of Education School of Environment of South China Normal University Guangzhou PR China
| |
Collapse
|
5
|
Zutterman F, Liégeois V, Champagne B. TDDFT Investigation of the Raman and Resonant Raman Spectra of Fluorescent Protein Chromophore Models. J Phys Chem B 2022; 126:3414-3424. [PMID: 35499480 DOI: 10.1021/acs.jpcb.2c01202] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The off-resonance and resonant Raman spectra have been simulated for models of fluorescent protein chromophores, those of the green fluorescent protein (GFP, called FP1) and of DsRed (called FP2), which presents a longer π-conjugated path, with the aim of providing a systematic investigation of structural but also computational aspects. These were performed at the (time-dependent) density functional theory [(TD)DFT] level. The off-resonance intensities have been calculated from the derivatives of the frequency-dependent polarizability with respect to the normal coordinates while the resonant ones have been evaluated using Huang-Rhys factors determined from the gradients of the excitation energies with respect to the normal coordinates. When applied with the M05 meta-GGA exchange-correlation functional, this simple computational scheme can reproduce most of the experimental Raman signatures of FP1 in its protonated and deprotonated forms, the differences of vibrational signatures of the cis (Z) and trans (E) isomers, as well as their changes as a function of the excitation wavelength. On the other hand, testing the predictions made for FP2 would require new experimental work. It was also observed that simulations with methods that inadequately predict the resonant Raman spectra could nevertheless produce a UV-vis absorption spectrum that is quite similar to the one obtained with better methods, once realistic peak broadening has been applied.
Collapse
Affiliation(s)
- Freddy Zutterman
- Laboratoire de Chimie Théorique (LCT), Unité de Chimie-Physique Théorique et Structurale (UCPTS), NISM (Namur Institute of Structured Matter), Université de Namur, rue de Bruxelles, 61, B-5000 Namur, Belgium
| | - Vincent Liégeois
- Laboratoire de Chimie Théorique (LCT), Unité de Chimie-Physique Théorique et Structurale (UCPTS), NISM (Namur Institute of Structured Matter), Université de Namur, rue de Bruxelles, 61, B-5000 Namur, Belgium
| | - Benoît Champagne
- Laboratoire de Chimie Théorique (LCT), Unité de Chimie-Physique Théorique et Structurale (UCPTS), NISM (Namur Institute of Structured Matter), Université de Namur, rue de Bruxelles, 61, B-5000 Namur, Belgium
| |
Collapse
|
6
|
Lu SY, Zuehlsdorff TJ, Hong H, Aguirre VP, Isborn CM, Shi L. The Influence of Electronic Polarization on Nonlinear Optical Spectroscopy. J Phys Chem B 2021; 125:12214-12227. [PMID: 34726915 DOI: 10.1021/acs.jpcb.1c05914] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The environment surrounding a chromophore can dramatically affect the energy absorption and relaxation process, as manifested in optical spectra. Simulations of nonlinear optical spectroscopy, such as two-dimensional electronic spectroscopy (2DES) and transient absorption (TA), will be influenced by the computational model of the environment. We here compare a fixed point charge molecular mechanics model and a quantum mechanical (QM) model of the environment in computed 2DES and TA spectra of Nile red in water and the chromophore of photoactive yellow protein (PYP) in water and protein environments. In addition to simulating these nonlinear optical spectra, we directly juxtapose the computed excitation energy correlation function to the dynamic Stokes shift function often used to analyze environment dynamics. Overall, we find that for the three systems studied here the mutual electronic polarization provided by the QM environment manifests in broader 2DES signals, as well as a larger reorganization energy and a larger static Stokes shift due to stronger coupling between the chromophore and the environment.
Collapse
Affiliation(s)
- Shao-Yu Lu
- Department of Chemistry and Biochemistry, University of California Merced, Merced, California 95343, United States
| | - Tim J Zuehlsdorff
- Department of Chemistry, Oregon State University, Corvallis, Oregon 97331, United States
| | - Hanbo Hong
- Department of Chemistry and Biochemistry, University of California Merced, Merced, California 95343, United States
| | - Vincent P Aguirre
- Department of Chemistry and Biochemistry, University of California Merced, Merced, California 95343, United States
| | - Christine M Isborn
- Department of Chemistry and Biochemistry, University of California Merced, Merced, California 95343, United States
| | - Liang Shi
- Department of Chemistry and Biochemistry, University of California Merced, Merced, California 95343, United States
| |
Collapse
|
7
|
Borrego-Sánchez A, Zemmouche M, Carmona-García J, Francés-Monerris A, Mulet P, Navizet I, Roca-Sanjuán D. Multiconfigurational Quantum Chemistry Determinations of Absorption Cross Sections (σ) in the Gas Phase and Molar Extinction Coefficients (ε) in Aqueous Solution and Air-Water Interface. J Chem Theory Comput 2021; 17:3571-3582. [PMID: 33974417 PMCID: PMC8444339 DOI: 10.1021/acs.jctc.0c01083] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/14/2020] [Indexed: 11/29/2022]
Abstract
Theoretical determinations of absorption cross sections (σ) in the gas phase and molar extinction coefficients (ε) in condensed phases (water solution, interfaces or surfaces, protein or nucleic acids embeddings, etc.) are of interest when rates of photochemical processes, J = ∫ ϕ(λ) σ(λ) I(λ) dλ, are needed, where ϕ(λ) and I(λ) are the quantum yield of the process and the irradiance of the light source, respectively, as functions of the wavelength λ. Efficient computational strategies based on single-reference quantum-chemistry methods have been developed enabling determinations of line shapes or, in some cases, achieving rovibrational resolution. Developments are however lacking for strongly correlated problems, with many excited states, high-order excitations, and/or near degeneracies between states of the same and different spin multiplicities. In this work, we define and compare the performance of distinct computational strategies using multiconfigurational quantum chemistry, nuclear sampling of the chromophore (by means of molecular dynamics, ab initio molecular dynamics, or Wigner sampling), and conformational and statistical sampling of the environment (by means of molecular dynamics). A new mathematical approach revisiting previous absolute orientation algorithms is also developed to improve alignments of geometries. These approaches are benchmarked through the nπ* band of acrolein not only in the gas phase and water solution but also in a gas-phase/water interface, a common situation for instance in atmospheric chemistry. Subsequently, the best strategy is used to compute the absorption band for the adduct formed upon addition of an OH radical to the C6 position of uracil and compared with the available experimental data. Overall, quantum Wigner sampling of the chromophore with molecular dynamics sampling of the environment with CASPT2 electronic-structure determinations arise as a powerful methodology to predict meaningful σ(λ) and ε(λ) band line shapes with accurate absolute intensities.
Collapse
Affiliation(s)
- Ana Borrego-Sánchez
- Instituto
Andaluz de Ciencias de la Tierra, CSIC-University
of Granada, Av. de las
Palmeras 4, 18100 Armilla, Granada, Spain
| | - Madjid Zemmouche
- MSME,
Univ Gustave Eiffel, CNRS UMR 8208, Univ Paris-Est Créteil 8208, F-77454 Marne-la-Vallée, France
| | - Javier Carmona-García
- Instituto
de Ciencia Molecular, Universitat de València, P.O. Box 22085, València, Spain
| | - Antonio Francés-Monerris
- Université
de Lorraine and CNRS, LPCT UMR 7019, F-54000 Nancy, France
- Departamento
de Química Física, Universitat
de València, C/Dr.
Moliner 50, 46100 Burjassot, Spain
| | - Pep Mulet
- Departamento
de Matemáticas Área de Matemática Aplicada Facultad
de Matemáticas C/Dr. Moliner, 50 46100 Burjassot, Spain
| | - Isabelle Navizet
- MSME,
Univ Gustave Eiffel, CNRS UMR 8208, Univ Paris-Est Créteil 8208, F-77454 Marne-la-Vallée, France
| | - Daniel Roca-Sanjuán
- Instituto
de Ciencia Molecular, Universitat de València, P.O. Box 22085, València, Spain
| |
Collapse
|
8
|
Zuehlsdorff TJ, Shedge SV, Lu SY, Hong H, Aguirre VP, Shi L, Isborn CM. Vibronic and Environmental Effects in Simulations of Optical Spectroscopy. Annu Rev Phys Chem 2021; 72:165-188. [DOI: 10.1146/annurev-physchem-090419-051350] [Citation(s) in RCA: 14] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
Including both environmental and vibronic effects is important for accurate simulation of optical spectra, but combining these effects remains computationally challenging. We outline two approaches that consider both the explicit atomistic environment and the vibronic transitions. Both phenomena are responsible for spectral shapes in linear spectroscopy and the electronic evolution measured in nonlinear spectroscopy. The first approach utilizes snapshots of chromophore-environment configurations for which chromophore normal modes are determined. We outline various approximations for this static approach that assumes harmonic potentials and ignores dynamic system-environment coupling. The second approach obtains excitation energies for a series of time-correlated snapshots. This dynamic approach relies on the accurate truncation of the cumulant expansion but treats the dynamics of the chromophore and the environment on equal footing. Both approaches show significant potential for making strides toward more accurate optical spectroscopy simulations of complex condensed phase systems.
Collapse
Affiliation(s)
- Tim J. Zuehlsdorff
- Department of Chemistry, Oregon State University, Corvallis, Oregon 97331, USA
| | - Sapana V. Shedge
- Department of Chemistry and Chemical Biology, University of California, Merced, California 95343, USA
| | - Shao-Yu Lu
- Department of Chemistry and Chemical Biology, University of California, Merced, California 95343, USA
| | - Hanbo Hong
- Department of Chemistry and Chemical Biology, University of California, Merced, California 95343, USA
| | - Vincent P. Aguirre
- Department of Chemistry and Chemical Biology, University of California, Merced, California 95343, USA
| | - Liang Shi
- Department of Chemistry and Chemical Biology, University of California, Merced, California 95343, USA
| | - Christine M. Isborn
- Department of Chemistry and Chemical Biology, University of California, Merced, California 95343, USA
| |
Collapse
|
9
|
Liu Y, Xu Q, Liu L, Wang L, He D, Zhuang X, Wang M. The spectral-shapes of absorption, emission, ECD and CPL of a fluorene-fused [7]helicene: Vibronic effect and solvent inhomogenous broadening. SPECTROCHIMICA ACTA. PART A, MOLECULAR AND BIOMOLECULAR SPECTROSCOPY 2021; 249:119293. [PMID: 33326916 DOI: 10.1016/j.saa.2020.119293] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/11/2020] [Revised: 11/27/2020] [Accepted: 11/29/2020] [Indexed: 06/12/2023]
Abstract
We compute the line shapes of absorption (ABS), emission (EMI), electronic circular dichroism (ECD), and circularly polarized luminescence (CPL) of a fluorene-fused [7]helicene. To elucidate different physico-chemical factors responsible for the spectral shapes, we account for Franck-Condon (FC) contribution, Herzberg-Teller (HT) effect and Duschinsky rotation, as well as the solvent inhomogenous broadening. The latter is connected to the solvent reorganization energy which is calculated with the state-specific (SS) implementation of Polarizable Continuum Model (PCM). To understand the relative importance of different effects of equilibrium geometry displacements, frequencies changes, and Duschinsky mixing, we carefully compare and discuss the results of Adiabatic Hessian (AH), Adiabatic Shift and Frequencies (ASF), and Adiabatic Shift (AS) models in details. Our results agree with the experimental spectra very well. They show that FC contribution dominates the spectral shapes and HT effect leads to the break of mirror images between absorptive (ABS + ECD) and emissive (EMI + CPL) spectra.
Collapse
Affiliation(s)
- Yanli Liu
- School of Physics and Optoelectronics Engineering, Ludong University, 264025 Yantai, Shandong, PR China
| | - Qiushuang Xu
- School of Physics and Optoelectronics Engineering, Ludong University, 264025 Yantai, Shandong, PR China; School of Physics Engineering, Qufu Normal University, 2673100 Qufu, Shandong, PR China
| | - Lu Liu
- School of Physics and Optoelectronics Engineering, Ludong University, 264025 Yantai, Shandong, PR China
| | - Li Wang
- School of Physics and Optoelectronics Engineering, Ludong University, 264025 Yantai, Shandong, PR China.
| | - Di He
- School of Physics and Optoelectronics Engineering, Ludong University, 264025 Yantai, Shandong, PR China
| | - Xingang Zhuang
- The 41st Research Institute of CETC, 266555 Qingdao, Shandong, PR China
| | - Meishan Wang
- School of Physics and Optoelectronics Engineering, Ludong University, 264025 Yantai, Shandong, PR China.
| |
Collapse
|
10
|
Shedge SV, Zuehlsdorff TJ, Khanna A, Conley S, Isborn CM. Explicit environmental and vibronic effects in simulations of linear and nonlinear optical spectroscopy. J Chem Phys 2021; 154:084116. [PMID: 33639769 DOI: 10.1063/5.0038196] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Accurately simulating the linear and nonlinear electronic spectra of condensed phase systems and accounting for all physical phenomena contributing to spectral line shapes presents a significant challenge. Vibronic transitions can be captured through a harmonic model generated from the normal modes of a chromophore, but it is challenging to also include the effects of specific chromophore-environment interactions within such a model. We work to overcome this limitation by combining approaches to account for both explicit environment interactions and vibronic couplings for simulating both linear and nonlinear optical spectra. We present and show results for three approaches of varying computational cost for combining ensemble sampling of chromophore-environment configurations with Franck-Condon line shapes for simulating linear spectra. We present two analogous approaches for nonlinear spectra. Simulated absorption spectra and two-dimensional electronic spectra (2DES) are presented for the Nile red chromophore in different solvent environments. Employing an average Franck-Condon or 2DES line shape appears to be a promising method for simulating linear and nonlinear spectroscopy for a chromophore in the condensed phase.
Collapse
Affiliation(s)
- Sapana V Shedge
- Chemistry and Chemical Biology, University of California Merced, Merced, California 95343, USA
| | - Tim J Zuehlsdorff
- Department of Chemistry, Oregon State University, Corvallis, Oregon 97331, USA
| | - Ajay Khanna
- Chemistry and Chemical Biology, University of California Merced, Merced, California 95343, USA
| | - Stacey Conley
- Chemistry and Chemical Biology, University of California Merced, Merced, California 95343, USA
| | - Christine M Isborn
- Chemistry and Chemical Biology, University of California Merced, Merced, California 95343, USA
| |
Collapse
|
11
|
Raucci U, Perrella F, Donati G, Zoppi M, Petrone A, Rega N. Ab-initio molecular dynamics and hybrid explicit-implicit solvation model for aqueous and nonaqueous solvents: GFP chromophore in water and methanol solution as case study. J Comput Chem 2020; 41:2228-2239. [PMID: 32770577 DOI: 10.1002/jcc.26384] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/16/2020] [Revised: 06/21/2020] [Accepted: 06/27/2020] [Indexed: 12/15/2022]
Abstract
Solute-solvent interactions are proxies for understanding how the electronic density of a chromophore interacts with the environment in a more exhaustive way. The subtle balance between polarization, electrostatic, and non-bonded interactions need to be accurately described to obtain good agreement between simulations and experiments. First principles approaches providing accurate configurational sampling through molecular dynamics may be a suitable choice to describe solvent effects on solute chemical-physical properties and spectroscopic features, such as optical absorption of dyes. In this context, accurate energy potentials, obtained by hybrid implicit/explicit solvation methods along with employing nonperiodic boundary conditions, are required to represent bulk solvent around a large solute-solvent cluster. In this work, a novel strategy to simulate methanol solutions is proposed combining ab initio molecular dynamics, a hybrid implicit/explicit flexible solvent model, nonperiodic boundary conditions, and time dependent density functional theory. As case study, the robustness of the proposed protocol has been gauged by investigating the microsolvation and electronic absorption of the anionic green fluorescent protein chromophore in methanol and aqueous solution. Satisfactory results are obtained, reproducing the microsolvation layout of the chromophore and, as a consequence, the experimental trends shown by the optical absorption in different solvents.
Collapse
Affiliation(s)
- Umberto Raucci
- Dipartimento di Scienze Chimiche, Università di Napoli Federico II, Complesso Universitario di M.S. Angelo, Naples, Italy
| | - Fulvio Perrella
- Dipartimento di Scienze Chimiche, Università di Napoli Federico II, Complesso Universitario di M.S. Angelo, Naples, Italy
| | - Greta Donati
- Dipartimento di Scienze Chimiche, Università di Napoli Federico II, Complesso Universitario di M.S. Angelo, Naples, Italy.,Dipartimento di Chimica e Biologia "Adolfo Zambelli", Università di Salerno, Fisciano, Italy
| | - Maria Zoppi
- Dipartimento di Scienze Chimiche, Università di Napoli Federico II, Complesso Universitario di M.S. Angelo, Naples, Italy
| | - Alessio Petrone
- Dipartimento di Scienze Chimiche, Università di Napoli Federico II, Complesso Universitario di M.S. Angelo, Naples, Italy
| | - Nadia Rega
- Dipartimento di Scienze Chimiche, Università di Napoli Federico II, Complesso Universitario di M.S. Angelo, Naples, Italy.,Center for Advanced Biomaterials for Healthcare@CRIB, Naples, Italy
| |
Collapse
|
12
|
Zuehlsdorff TJ, Hong H, Shi L, Isborn CM. Nonlinear spectroscopy in the condensed phase: The role of Duschinsky rotations and third order cumulant contributions. J Chem Phys 2020; 153:044127. [PMID: 32752702 DOI: 10.1063/5.0013739] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022] Open
Abstract
First-principles modeling of nonlinear optical spectra in the condensed phase is highly challenging because both environment and vibronic interactions can play a large role in determining spectral shapes and excited state dynamics. Here, we compute two dimensional electronic spectroscopy (2DES) signals based on a cumulant expansion of the energy gap fluctuation operator, with specific focus on analyzing mode mixing effects introduced by the Duschinsky rotation and the role of the third order term in the cumulant expansion for both model and realistic condensed phase systems. We show that for a harmonic model system, the third order cumulant correction captures effects introduced by a mismatch in curvatures of ground and excited state potential energy surfaces, as well as effects of mode mixing. We also demonstrate that 2DES signals can be accurately reconstructed from purely classical correlation functions using quantum correction factors. We then compute nonlinear optical spectra for the Nile red and methylene blue chromophores in solution, assessing the third order cumulant contribution for realistic systems. We show that the third order cumulant correction is strongly dependent on the treatment of the solvent environment, revealing the interplay between environmental polarization and the electronic-vibrational coupling.
Collapse
Affiliation(s)
- Tim J Zuehlsdorff
- Chemistry and Chemical Biology, University of California Merced, Merced, California 95343, USA
| | - Hanbo Hong
- Chemistry and Chemical Biology, University of California Merced, Merced, California 95343, USA
| | - Liang Shi
- Chemistry and Chemical Biology, University of California Merced, Merced, California 95343, USA
| | - Christine M Isborn
- Chemistry and Chemical Biology, University of California Merced, Merced, California 95343, USA
| |
Collapse
|
13
|
Wang CW, Zhu C, Lin SH. Extremely solvent-enhanced absorbance and fluorescence of carbazole interpreted using a damped Franck-Condon simulation. J Chem Phys 2020; 152:104106. [PMID: 32171227 DOI: 10.1063/1.5141037] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022] Open
Abstract
Extremely solvent-enhanced absorption and fluorescence spectra of carbazole were investigated by performing a generalized multi-set damped Franck-Condon spectral simulation. Experimental absorption and fluorescence spectra of carbazole in the gas phase were first well reproduced by performing an un-damped Franck-Condon simulation, but a one-set scaling damped Franck-Condon simulation severely underestimated the intensities of the peaks of experimental absorption and fluorescence spectra of carbazole in n-hexane. Then, a multi-set scaling damped Franck-Condon simulation was proposed and carried out for simulating the extremely solvent-enhanced absorbance and fluorescence, and here, the simulated spectra agreed well with the experimental ones. Five (four) representative solvent-enhanced normal modes corresponding to the combination of ring stretching and ring breathing vibrational motions were determined to be responsible for enhanced absorbance (fluorescence) in n-hexane solution. Furthermore, different scalings were applied to the ground and first-excited states, resulting in different enhancement of absorbance and fluorescence, and this analysis revealed atoms in the carbazole interacting with n-hexane solvent molecules and, hence, leading to different normal-mode vibrational vector patterns in the ground and first-excited states, respectively. Basically, the same conclusion was drawn from a simulation with HF-CIS and the three functionals (TD)B3LYP, (TD)B3LYP-35, and (TD)BHandHLYP. The present multi-set scaling damped Franck-Condon simulation scheme was demonstrated to successfully interpret extremely solvent-enhanced absorbance and fluorescence of carbazole in n-hexane-solvent.
Collapse
Affiliation(s)
- Chen-Wen Wang
- Institute of Molecular Science and Department of Applied Chemistry, National Chiao Tung University, Hsinchu 30010, Taiwan
| | - Chaoyuan Zhu
- Institute of Molecular Science and Department of Applied Chemistry, National Chiao Tung University, Hsinchu 30010, Taiwan
| | - Sheng Hsien Lin
- Institute of Molecular Science and Department of Applied Chemistry, National Chiao Tung University, Hsinchu 30010, Taiwan
| |
Collapse
|
14
|
Zuehlsdorff TJ, Hong H, Shi L, Isborn CM. Influence of Electronic Polarization on the Spectral Density. J Phys Chem B 2019; 124:531-543. [DOI: 10.1021/acs.jpcb.9b10250] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/05/2023]
Affiliation(s)
- Tim J. Zuehlsdorff
- Chemistry and Chemical Biology, University of California Merced, Merced, California 95343, United States
| | - Hanbo Hong
- Chemistry and Chemical Biology, University of California Merced, Merced, California 95343, United States
| | - Liang Shi
- Chemistry and Chemical Biology, University of California Merced, Merced, California 95343, United States
| | - Christine M. Isborn
- Chemistry and Chemical Biology, University of California Merced, Merced, California 95343, United States
| |
Collapse
|
15
|
Lin CY, Romei MG, Oltrogge LM, Mathews II, Boxer SG. Unified Model for Photophysical and Electro-Optical Properties of Green Fluorescent Proteins. J Am Chem Soc 2019; 141:15250-15265. [PMID: 31450887 DOI: 10.1021/jacs.9b07152] [Citation(s) in RCA: 50] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Abstract
Green fluorescent proteins (GFPs) have become indispensable imaging and optogenetic tools. Their absorption and emission properties can be optimized for specific applications. Currently, no unified framework exists to comprehensively describe these photophysical properties, namely the absorption maxima, emission maxima, Stokes shifts, vibronic progressions, extinction coefficients, Stark tuning rates, and spontaneous emission rates, especially one that includes the effects of the protein environment. In this work, we study the correlations among these properties from systematically tuned GFP environmental mutants and chromophore variants. Correlation plots reveal monotonic trends, suggesting that all these properties are governed by one underlying factor dependent on the chromophore's environment. By treating the anionic GFP chromophore as a mixed-valence compound existing as a superposition of two resonance forms, we argue that this underlying factor is defined as the difference in energy between the two forms, or the driving force, which is tuned by the environment. We then introduce a Marcus-Hush model with the bond length alternation vibrational mode, treating the GFP absorption band as an intervalence charge transfer band. This model explains all of the observed strong correlations among photophysical properties; related subtopics are extensively discussed in the Supporting Information. Finally, we demonstrate the model's predictive power by utilizing the additivity of the driving force. The model described here elucidates the role of the protein environment in modulating the photophysical properties of the chromophore, providing insights and limitations for designing new GFPs with desired phenotypes. We argue that this model should also be generally applicable to both biological and nonbiological polymethine dyes.
Collapse
Affiliation(s)
- Chi-Yun Lin
- Department of Chemistry , Stanford University , Stanford , California 94305 , United States
| | - Matthew G Romei
- Department of Chemistry , Stanford University , Stanford , California 94305 , United States
| | - Luke M Oltrogge
- Department of Chemistry , Stanford University , Stanford , California 94305 , United States
| | - Irimpan I Mathews
- Stanford Synchrotron Radiation Lightsource , 2575 Sand Hill Road , Menlo Park , California 94025 , United States
| | - Steven G Boxer
- Department of Chemistry , Stanford University , Stanford , California 94305 , United States
| |
Collapse
|
16
|
Zuehlsdorff TJ, Montoya-Castillo A, Napoli JA, Markland TE, Isborn CM. Optical spectra in the condensed phase: Capturing anharmonic and vibronic features using dynamic and static approaches. J Chem Phys 2019; 151:074111. [PMID: 31438704 DOI: 10.1063/1.5114818] [Citation(s) in RCA: 44] [Impact Index Per Article: 8.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/22/2023] Open
Abstract
Simulating optical spectra in the condensed phase remains a challenge for theory due to the need to capture spectral signatures arising from anharmonicity and dynamical effects, such as vibronic progressions and asymmetry. As such, numerous simulation methods have been developed that invoke different approximations and vary in their ability to capture different physical regimes. Here, we use several models of chromophores in the condensed phase and ab initio molecular dynamics simulations to rigorously assess the applicability of methods to simulate optical absorption spectra. Specifically, we focus on the ensemble scheme, which can address anharmonic potential energy surfaces but relies on the applicability of extreme nuclear-electronic time scale separation; the Franck-Condon method, which includes dynamical effects but generally only at the harmonic level; and the recently introduced ensemble zero-temperature Franck-Condon approach, which straddles these limits. We also devote particular attention to the performance of methods derived from a cumulant expansion of the energy gap fluctuations and test the ability to approximate the requisite time correlation functions using classical dynamics with quantum correction factors. These results provide insights as to when these methods are applicable and able to capture the features of condensed phase spectra qualitatively and, in some cases, quantitatively across a range of regimes.
Collapse
Affiliation(s)
- Tim J Zuehlsdorff
- Chemistry and Chemical Biology, University of California Merced, Merced, California 95343, USA
| | | | - Joseph A Napoli
- Department of Chemistry, Stanford University, Stanford, California 94305, USA
| | - Thomas E Markland
- Department of Chemistry, Stanford University, Stanford, California 94305, USA
| | - Christine M Isborn
- Chemistry and Chemical Biology, University of California Merced, Merced, California 95343, USA
| |
Collapse
|
17
|
Shedge SV, Zuehlsdorff TJ, Servis MJ, Clark AE, Isborn CM. Effect of Ions on the Optical Absorption Spectra of Aqueously Solvated Chromophores. J Phys Chem A 2019; 123:6175-6184. [DOI: 10.1021/acs.jpca.9b03163] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/26/2023]
Affiliation(s)
- Sapana V. Shedge
- Chemistry and Chemical Biology, University of California Merced, Merced, California 95343, United States
| | - Tim J. Zuehlsdorff
- Chemistry and Chemical Biology, University of California Merced, Merced, California 95343, United States
| | - Michael J. Servis
- Department of Chemistry and the Material Science and Engineering Program, Washington State University, Pullman, Washington 99164, United States
| | - Aurora E. Clark
- Department of Chemistry and the Material Science and Engineering Program, Washington State University, Pullman, Washington 99164, United States
- Voiland School of Chemical Engineering and Bioengineering, Washington State University, Pullman, Washington 99164, United States
| | - Christine M. Isborn
- Chemistry and Chemical Biology, University of California Merced, Merced, California 95343, United States
| |
Collapse
|
18
|
Zuehlsdorff TJ, Napoli JA, Milanese JM, Markland TE, Isborn CM. Unraveling electronic absorption spectra using nuclear quantum effects: Photoactive yellow protein and green fluorescent protein chromophores in water. J Chem Phys 2018; 149:024107. [PMID: 30007372 DOI: 10.1063/1.5025517] [Citation(s) in RCA: 27] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
Many physical phenomena must be accounted for to accurately model solution-phase optical spectral line shapes, from the sampling of chromophore-solvent configurations to the electronic-vibrational transitions leading to vibronic fine structure. Here we thoroughly explore the role of nuclear quantum effects, direct and indirect solvent effects, and vibronic effects in the computation of the optical spectrum of the aqueously solvated anionic chromophores of green fluorescent protein and photoactive yellow protein. By analyzing the chromophore and solvent configurations, the distributions of vertical excitation energies, the absorption spectra computed within the ensemble approach, and the absorption spectra computed within the ensemble plus zero-temperature Franck-Condon approach, we show how solvent, nuclear quantum effects, and vibronic transitions alter the optical absorption spectra. We find that including nuclear quantum effects in the sampling of chromophore-solvent configurations using ab initio path integral molecular dynamics simulations leads to improved spectral shapes through three mechanisms. The three mechanisms that lead to line shape broadening and a better description of the high-energy tail are softening of heavy atom bonds in the chromophore that couple to the optically bright state, widening the distribution of vertical excitation energies from more diverse solvation environments, and redistributing spectral weight from the 0-0 vibronic transition to higher energy vibronic transitions when computing the Franck-Condon spectrum in a frozen solvent pocket. The absorption spectra computed using the combined ensemble plus zero-temperature Franck-Condon approach yield significant improvements in spectral shape and width compared to the spectra computed with the ensemble approach. Using the combined approach with configurations sampled from path integral molecular dynamics trajectories presents a significant step forward in accurately modeling the absorption spectra of aqueously solvated chromophores.
Collapse
Affiliation(s)
- Tim J Zuehlsdorff
- Chemistry and Chemical Biology, University of California Merced, Merced, California 95343, USA
| | - Joseph A Napoli
- Department of Chemistry, Stanford University, Stanford, California 94305, USA
| | - Joel M Milanese
- Chemistry and Chemical Biology, University of California Merced, Merced, California 95343, USA
| | - Thomas E Markland
- Department of Chemistry, Stanford University, Stanford, California 94305, USA
| | - Christine M Isborn
- Chemistry and Chemical Biology, University of California Merced, Merced, California 95343, USA
| |
Collapse
|
19
|
Zuehlsdorff TJ, Isborn CM. Combining the ensemble and Franck-Condon approaches for calculating spectral shapes of molecules in solution. J Chem Phys 2018; 148:024110. [PMID: 29331131 DOI: 10.1063/1.5006043] [Citation(s) in RCA: 52] [Impact Index Per Article: 8.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The correct treatment of vibronic effects is vital for the modeling of absorption spectra of many solvated dyes. Vibronic spectra for small dyes in solution can be easily computed within the Franck-Condon approximation using an implicit solvent model. However, implicit solvent models neglect specific solute-solvent interactions on the electronic excited state. On the other hand, a straightforward way to account for solute-solvent interactions and temperature-dependent broadening is by computing vertical excitation energies obtained from an ensemble of solute-solvent conformations. Ensemble approaches usually do not account for vibronic transitions and thus often produce spectral shapes in poor agreement with experiment. We address these shortcomings by combining zero-temperature vibronic fine structure with vertical excitations computed for a room-temperature ensemble of solute-solvent configurations. In this combined approach, all temperature-dependent broadening is treated classically through the sampling of configurations and quantum mechanical vibronic contributions are included as a zero-temperature correction to each vertical transition. In our calculation of the vertical excitations, significant regions of the solvent environment are treated fully quantum mechanically to account for solute-solvent polarization and charge-transfer. For the Franck-Condon calculations, a small amount of frozen explicit solvent is considered in order to capture solvent effects on the vibronic shape function. We test the proposed method by comparing calculated and experimental absorption spectra of Nile red and the green fluorescent protein chromophore in polar and non-polar solvents. For systems with strong solute-solvent interactions, the combined approach yields significant improvements over the ensemble approach. For systems with weak to moderate solute-solvent interactions, both the high-energy vibronic tail and the width of the spectra are in excellent agreement with experiments.
Collapse
Affiliation(s)
- T J Zuehlsdorff
- School of Natural Sciences, University of California Merced, N. Lake Road, Merced, California 95344, USA
| | - C M Isborn
- School of Natural Sciences, University of California Merced, N. Lake Road, Merced, California 95344, USA
| |
Collapse
|
20
|
Zuehlsdorff TJ, Isborn CM. Erratum: “Combining the ensemble and Franck-Condon approaches for calculating spectral shapes of molecules in solution” [J. Chem. Phys. 148, 024110 (2018)]. J Chem Phys 2018; 148:079901. [DOI: 10.1063/1.5024413] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- T. J. Zuehlsdorff
- School of Natural Sciences, University of California Merced, N. Lake Road, Merced, California 95344, USA
| | - C. M. Isborn
- School of Natural Sciences, University of California Merced, N. Lake Road, Merced, California 95344, USA
| |
Collapse
|
21
|
Zutterman F, Liégeois V, Champagne B. Simulation of the UV/Visible Absorption Spectra of Fluorescent Protein Chromophore Models. CHEMPHOTOCHEM 2017. [DOI: 10.1002/cptc.201700002] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/22/2023]
Affiliation(s)
- Freddy Zutterman
- Laboratoire de Chimie Théorique, Unité de Chimie-Physique Théorique et Structurale; Université de Namur; rue de Bruxelles, 61 B-5000 Namur Belgium
| | - Vincent Liégeois
- Laboratoire de Chimie Théorique, Unité de Chimie-Physique Théorique et Structurale; Université de Namur; rue de Bruxelles, 61 B-5000 Namur Belgium
| | - Benoît Champagne
- Laboratoire de Chimie Théorique, Unité de Chimie-Physique Théorique et Structurale; Université de Namur; rue de Bruxelles, 61 B-5000 Namur Belgium
| |
Collapse
|
22
|
Cerezo J, Avila Ferrer FJ, Prampolini G, Santoro F. Modeling Solvent Broadening on the Vibronic Spectra of a Series of Coumarin Dyes. From Implicit to Explicit Solvent Models. J Chem Theory Comput 2015; 11:5810-25. [DOI: 10.1021/acs.jctc.5b00870] [Citation(s) in RCA: 56] [Impact Index Per Article: 6.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023]
Affiliation(s)
- Javier Cerezo
- Istituto
di Chimica dei Composti OrganoMetallici, Consiglio Nazionale delle Richerche (ICCOM-CNR), Area della Ricerca, via G. Moruzzi 1, I-56124 Pisa, Italy
| | - Francisco J. Avila Ferrer
- Istituto
di Chimica dei Composti OrganoMetallici, Consiglio Nazionale delle Richerche (ICCOM-CNR), Area della Ricerca, via G. Moruzzi 1, I-56124 Pisa, Italy
- Physical
Chemistry, Faculty of Science, University of Málaga, Málaga 29071, Spain
| | - Giacomo Prampolini
- Istituto
di Chimica dei Composti OrganoMetallici, Consiglio Nazionale delle Richerche (ICCOM-CNR), Area della Ricerca, via G. Moruzzi 1, I-56124 Pisa, Italy
| | - Fabrizio Santoro
- Istituto
di Chimica dei Composti OrganoMetallici, Consiglio Nazionale delle Richerche (ICCOM-CNR), Area della Ricerca, via G. Moruzzi 1, I-56124 Pisa, Italy
| |
Collapse
|
23
|
Cerezo J, Avila Ferrer FJ, Santoro F. Disentangling vibronic and solvent broadening effects in the absorption spectra of coumarin derivatives for dye sensitized solar cells. Phys Chem Chem Phys 2015; 17:11401-11. [DOI: 10.1039/c5cp00370a] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/14/2023]
Abstract
Individuation of vibronic and solvent contributions to the spectra of a family of coumarin dyes helps to understand the main differences in their lineshapes.
Collapse
Affiliation(s)
- Javier Cerezo
- CNR–Consiglio Nazionale delle Ricerche
- Istituto di Chimica dei Composti Organo Metallici (ICCOM-CNR)
- UOS di Pisa
- Area della Ricerca
- I-56124 Pisa
| | - Francisco J. Avila Ferrer
- CNR–Consiglio Nazionale delle Ricerche
- Istituto di Chimica dei Composti Organo Metallici (ICCOM-CNR)
- UOS di Pisa
- Area della Ricerca
- I-56124 Pisa
| | - Fabrizio Santoro
- CNR–Consiglio Nazionale delle Ricerche
- Istituto di Chimica dei Composti Organo Metallici (ICCOM-CNR)
- UOS di Pisa
- Area della Ricerca
- I-56124 Pisa
| |
Collapse
|