Effect of Solvation on Electron Detachment and Excitation Energies of a Green Fluorescent Protein Chromophore Variant

Understanding the Photoprocesses in Biological Systems: Need for Accurate Multireference Treatment

Light-matter interaction is crucial to life itself and revolves around many of the central processes in biology. The need for understanding these photochemical and photophysical processes cannot be overemphasized. Interaction of light with biological systems starts with the absorption of light and subsequent phenomena that occur in the excited states of the system. However, excited states are typically difficult to understand within the mean field approximation of quantum chemical methods. Therefore, suitable multireference methods and methodologies have been developed to understand these phenomena. In this Perspective, we will describe a few methods and methodologies suitable for these descriptions and discuss some persisting difficulties.

Alkylated green fluorescent protein chromophores: dynamics in the gas phase and in aqueous solution

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.

Spatial Decay and Limits of Quantum Solute-Solvent Interactions

Molecular excitations in the liquid-phase environment are renormalized by the surrounding solvent molecules. Herein, we employ the GW approximation to investigate the solvation effects on the ionization energy of phenol in various solvent environments. The electronic effects differ by up to 0.4 eV among the five investigated solvents. This difference depends on both the macroscopic solvent polarizability and the spatial decay of the solvation effects. The latter is probed by separating the electronic subspace and the GW correlation self-energy into fragments. The fragment correlation energy decays with increasing intermolecular distance and vanishes at ∼9 Å, and this pattern is independent of the type of solvent environment. The 9 Å cutoff defines an effective interacting volume within which the ionization energy shift per solvent molecule is proportional to the macroscopic solvent polarizability. Finally, we propose a simple model for computing the ionization energies of molecules in an arbitrary solvent environment.

Data-driven modeling of S0 → S1 excitation energy in the BODIPY chemical space: High-throughput computation, quantum machine learning, and inverse design

Derivatives of BODIPY are popular fluorophores due to their synthetic feasibility, structural rigidity, high quantum yield, and tunable spectroscopic properties. While the characteristic absorption maximum of BODIPY is at 2.5 eV, combinations of functional groups and substitution sites can shift the peak position by ±1 eV. Time-dependent long-range corrected hybrid density functional methods can model the lowest excitation energies offering a semi-quantitative precision of ±0.3 eV. Alas, the chemical space of BODIPYs stemming from combinatorial introduction of-even a few dozen-substituents is too large for brute-force high-throughput modeling. To navigate this vast space, we select 77 412 molecules and train a kernel-based quantum machine learning model providing <2% hold-out error. Further reuse of the results presented here to navigate the entire BODIPY universe comprising over 253 giga (253 × 10⁹) molecules is demonstrated by inverse-designing candidates with desired target excitation energies.

Software for the frontiers of quantum chemistry: An overview of developments in the Q-Chem 5 package

This article summarizes technical advances contained in the fifth major release of the Q-Chem quantum chemistry program package, covering developments since 2015. A comprehensive library of exchange-correlation functionals, along with a suite of correlated many-body methods, continues to be a hallmark of the Q-Chem software. The many-body methods include novel variants of both coupled-cluster and configuration-interaction approaches along with methods based on the algebraic diagrammatic construction and variational reduced density-matrix methods. Methods highlighted in Q-Chem 5 include a suite of tools for modeling core-level spectroscopy, methods for describing metastable resonances, methods for computing vibronic spectra, the nuclear-electronic orbital method, and several different energy decomposition analysis techniques. High-performance capabilities including multithreaded parallelism and support for calculations on graphics processing units are described. Q-Chem boasts a community of well over 100 active academic developers, and the continuing evolution of the software is supported by an "open teamware" model and an increasingly modular design.

Efficient treatment of molecular excitations in the liquid phase environment via stochastic many-body theory

Accurate predictions of charge excitation energies of molecules in the disordered condensed phase are central to the chemical reactivity, stability, and optoelectronic properties of molecules and critically depend on the specific environment. Herein, we develop a stochastic GW method for calculating these charge excitation energies. The approach employs maximally localized electronic states to define the electronic subspace of a molecule and the rest of the system, both of which are randomly sampled. We test the method on three solute-solvent systems: phenol, thymine, and phenylalanine in water. The results are in excellent agreement with the previous high-level calculations and available experimental data. The stochastic calculations for supercells containing up to 1000 electrons representing the solvated systems are inexpensive and require ≤1000 central processing unit hrs. We find that the coupling with the environment accounts for ∼40% of the total correlation energy. The solvent-to-solute feedback mechanism incorporated in the molecular correlation term causes up to 0.6 eV destabilization of the quasiparticle energy. Simulated photo-emission spectra exhibit red shifts, state-degeneracy lifting, and lifetime shortening. Our method provides an efficient approach for an accurate study of excitations of large molecules in realistic condensed phase environments.

An experimental and computational study of the effect of aqueous solution on the multiphoton ionisation photoelectron spectrum of phenol

We revisit the photoelectron spectroscopy of aqueous phenol in an effort to improve our understanding of the impact of inhomogeneous broadening and inelastic scattering on solution-phase photoelectron spectra.

Insights into the halochromism of a bithiophene π* probe

The positive halochromism of the solvatochromic 4-nitro-4'-(N,N-dimethylamino)-2,2'-bithiophene in acetone and 1-butanol was investigated by molecular dynamics simulations and quantum mechanics calculations. Interactions of the dye with the solvated sodium cation were found to have only a minor effect on the observed halochromism, which was ultimately ascribed to a predominant interaction between iodide anions and the N,N-dimethylamino group of the dye.

Discrete Solvent Reaction Field Calculations for One- and Two-Photon Absorptions of Solution-Phase Dimethylaminonitrostilbene Molecule

Based on the configurations generated by molecular dynamics (MD) simulations using the on-the-fly density-functional tight-bonding (DFTB) force field, we investigated performance of the discrete solvent reaction field (DRF) model coupled to time-dependent density functional theory (TD-DFT) for solvatochromic effect of one- and two-photon absorption phenomena. Dimethylaminonitrostilbene (DANS) molecule solvated in chloroform, dichloromethane, and dimethyl sulfoxide solvents was selected as a model system for our research purpose. For every selected MD/DFTB configuration, within the context of the DRF, solute molecule is represented by TD-DFT and solvent molecules are described by atomic charges and polarizabilities. The calculated one-photon absorption energies reproduce well the positive solvatochromic behavior of solvated DANS and are in good agreement with available experimental data. For the two-photon absorption cross section, even though our approach overshot the experimental data by about 20% in absolute magnitude, experimentally observed solvatochromic change was captured qualitatively in this work. At last, we examined the contributions of atomic charges and polarizabilities of solvent molecules to the solvatochromic shifts of properties of interest.

Characterization of the Photophysical Behavior of DFHBI Derivatives: Fluorogenic Molecules that Illuminate the Spinach RNA Aptamer

( Z)-5-(3,5-Difluoro-4-hydroxybenzylidene)-2,3-dimethyl-3,5-dihydro-4 H-imidazol-4-one (DFHBI) and its analogues are fluorogenic molecules that bind the Spinach aptamer (a small RNA molecule), which was selected for imaging RNA. They are extremely weakly fluorescent in liquid solvents. It had been hypothesized that photoisomerization is a very efficient nonradiative process of deactivation. We show, consistent with the results of other studies, that if the isomerization is impeded, the fluorescence signal is enhanced significantly. In addition, we provide a thorough characterization of the photophysical behavior of DFHBI and its derivatives, notably that of ( Z)-5-(3,5-difluoro-4-hydroxybenzylidene)-2-methyl-3-((perfluorophenyl)methyl)-3,5-dihydro-4 H-imidazol-4-one (PFP-DFHBI) in various solvent environments. Solvent-dependent studies were performed with various mixtures of solvents. The results suggest that hydrogen bonding or strong interactions of the solvents with the phenolic-OH group change the absorption band near 420-460 nm and the nature of emission near 430 and 500 nm through various degrees of stabilization and the transformation between the neutral and the anionic species at both ground and excited states. These observations are confirmed by using a methoxy-substituted molecule (( Z)-5-(4-methoxybenzylidene)-2,3-dimethyl-3,5-dihydro-4 H-imidazol-4-one), where the 420-460 nm band is absent in the presence of methanol and the spectra are similar to those of PFP-DFHBI in noninteracting solvents, such as acetonitrile and dichloromethane. Thus, in addition to the major role of photoisomerization as a nonradiative process of deactivation of the excited state, the fluorescence of DFHBI-type molecules is very sensitively dependent upon the pH of the medium as well as upon solvent-specific interactions, such as hydrogen-bonding ability and polarity.

Effect of solvation on the ionization of guanine nucleotide: A hybrid QM/EFP study

Ionization of nucleobases is affected by their biological environment, which includes both the effect of adjacent nucleotides as well as the presence of water around it. Guanine and its nucleotide have the lowest ionization potentials among the various DNA bases. Therefore, the threshold of ionization is dependent on that of guanine and its characterization is crucial to the prediction of interaction of light with DNA. We investigate the effect of solvation on the vertical ionization energies (VIEs) of guanine and its nucleotide. In this work, we have used hybrid quantum mechanics/molecular mechanics (QM/MM) approach with effective fragment potential as the MM method of choice and equation-of-motion coupled-cluster for ionization potential with singles and doubles (EOM-IP-CCSD) as the QM method. The performance of the hybrid scheme with respect to the full QM method shows an accuracy of ≤ 0.02-0.04 eV. The lowest few ionizations of the nucleotide are found to be from different parts of the moiety, that is, the nucleic acid base, phosphate, or sugar, and these ionization energies are very closely spaced giving rise to a very complicated spectrum. Furthermore, microsolvation has large effects on these ionizations and can lead to red or blue shift depending on the position of the water molecule. Even a single water molecule can change the order of ionized states in the nucleotide. The VIEs of the bulk solvated chromophores are predicted and compared to existing experimental spectra. The predominant role of polarization in the solvatochromic shift is noticed. © 2017 Wiley Periodicals, Inc.

Combined quantum-mechanical molecular mechanics calculations with NWChem and AMBER: Excited state properties of green fluorescent protein chromophore analogue in aqueous solution

Combined quantum mechanical molecular mechanics (QM/MM) calculations have become a popular methodology for efficient and accurate description of large molecular systems. In this work we introduce our development of a QM/MM framework based on two well-known codes-NWChem and AMBER. As an initial application area we are focused on excited state properties of small molecules in an aqueous phase using an analogue of the green fluorescent protein (GFP) chromophore as a particular test case. Our approach incorporates high level coupled cluster theory for the analysis of excited states providing a reliable theoretical analysis of effects of an aqueous solvation environment on the photochemical properties of the GFP chromophore. Using a systematic approach, which involves comparison of gas phase and aqueous phase results for different protonation states and conformations, we resolve existing uncertainties regarding the theoretical interpretation of experimental data. We observe that the impact of aqueous environment on charged states generally results in blue shifts of the absorption spectra, but the magnitude of the effect is sensitive to both protonation state and conformation and can be rationalized based on charge movement into the area of higher/lower external electrostatic potentials. At neutral pH levels the experimentally observed absorption signal is most likely coming from the phenol protonated form. Our results also show that the high level electron correlated method is essential for a proper description of excited states of GFP. © 2017 Wiley Periodicals, Inc.

High-level Ab Initio Absorption Spectra Simulations of Neutral, Anionic and Neutral+ Chromophore of Green Fluorescence Protein Chromophore Models in Gas Phase and Solution

Semiclassical ab initio simulations of the absorption spectra of neutral and anionic p-hydroxybenzylidene-2,3-dimethylimidazolinone (p-HBDI), a model chromophore of green fluorescent protein (GFP) and of a positively charged neutral (N+)-HBDI chromophore model, were performed in gas phase with the resolution-of-identity algebraic diagrammatic construction through second-order (RI-ADC(2)) method. The calculated absorption spectra in gas phase are composed of one band centered at 3.51 eV (HBDI), 2.50 eV (HBDI^- ) and 3.02 eV ((N+)-HBDI) owing to the absorption of the first ¹ ππ* transition. Band maxima are redshifted by ~0.1 eV with respect to the corresponding vertical energies. The COSMO-RI-ADC(2) calculations of the first vertical excitation energy of HBDI, HBDI^- and (N+)-HBDI forms in polar solution including microsolvation simulate the observed solvent redshift for neutral HBDI and the solvent blueshift of the HBDI^- and (N+)-HBDI forms. The state-specific solvation approach applied to TDDFT calculations reproduced the experimental solvent shifts for the three HBDI forms, demonstrating a more accurate theoretical description as compared to the linear-response TDDFT approach.

Modifications of the chromophore of Spinach aptamer based on QM:MM calculations

Spinach aptamer was developed as an RNA analog of the green fluorescent protein. The aptamer interacts with its ligand and modifies its electronic spectrum so that it fluoresces brightly at the wavelength of 501 nm. Song et al. investigated modifications of the ligand in their experimental study and found a molecule emitting at 523 nm upon creating a complex with the Spinach aptamer. The crystal structure of the aptamer in complex with its original ligand has been published, which enabled us to study the system computationally. In this article, we suggest several new modifications of the ligand that shift the emission maximum of the complex to even longer wavelengths. Our results are based on combined quantum mechanical/molecular mechanical calculations with DFT method used for geometry optimization and TD-DFT for calculations of absorption and emission energies.