Quantum chemistry study of the multiphoton absorption in enhanced green fluorescent protein at the single amino acid residue level

The chromophore (CRO) of fluorescent proteins (FPs) is embedded in a complex environment that is a source of specific interactions with the CRO. Understanding how these interactions influence FPs spectral properties is important for a directed design of novel markers with desired characteristics. In this work, we apply computational chemistry methods to gain insight into one-, two- and three-photon absorption (1PA, 2PA, 3PA) tuning in enhanced green fluorescent protein (EGFP). To achieve this goal, we built EGFP models differing in: i) number and position of hydrogen-bonds (h-bonds) donors to the CRO and ii) the electric field, as approximated by polarizable force  field, acting on the CRO. We  find that h-bonding to the CRO's phenolate oxygen results in stronger one- and multiphoton absorption. The brighter absorption can be also achieved by creating more positive electric field near the CRO's phenolate moiety. Interestingly, while individual CRO-environment h-bonds usually enhance 1PA and 2PA, it takes a few h-bond donors to enhance 3PA. Clearly, response of the absorption intensity to many-body effects depends on the excitation mechanism.  We further employ symmetry-adapted perturbation theory (SAPT) to reveal excellent (2PA) and good (3PA) correlation of multiphoton intensity with electrostatic and induction interaction energies. This points to importance of accounting for mutual CRO-environment polarization in quantitative calculations of absorption spectra in FPs.

Harmonic Infrared and Raman Spectra in Molecular Environments Using the Polarizable Embedding Model

We present a fully analytic approach to calculate infrared (IR) and Raman spectra of molecules embedded in complex molecular environments modeled using the fragment-based polarizable embedding (PE) model. We provide the theory for the calculation of analytic second-order geometric derivatives of molecular energies and first-order geometric derivatives of electric dipole moments and dipole-dipole polarizabilities within the PE model. The derivatives are implemented using a general open-ended response theory framework, thus allowing for an extension to higher-order derivatives. The embedding-potential parameters used to describe the environment in the PE model are derived through first-principles calculations, thus allowing a wide variety of systems to be modeled, including solvents, proteins, and other large and complex molecular environments. Here, we present proof-of-principle calculations of IR and Raman spectra of acetone in different solvents. This work is an important step toward calculating accurate vibrational spectra of molecules embedded in realistic environments.

The Challenge of Accurate Computation of Two-Photon Absorption Properties of Organic Chromophores in the Condensed Phase

Two strategies are applied to evaluate the effect of the environment on the two-photon absorption (TPA) cross sections for two characteristic excited states of C₂H₄ upon complexation with H₂O. The supermolecular strategy provides the reference complexation-induced shifts and uses either the EOM-CCSD or ADC(2) method. The embedding strategy is based on frozen-density-embedding theory (FDET) and uses only fundamental constants. The TPA cross sections from high-level supermolecular calculations are extremely basis-set-sensitive. Literature data and the present study indicate that accuracy of the absolute TPA cross sections below 100 atomic units and their shifts below 10 atomic units remains a challenge. The obtained FDET results show a similar basis-set behavior. For the largest basis set (d-aug-cc-pVQZ), TPA cross sections obtained from these two strategies are in excellent agreement. The complexation-induced shifts have the correct sign of the effect and a small (12-33%) relative error in magnitude. The deviations of the FDET-derived shifts from the reference are of similar magnitude as the reliability threshold of the reference shifts.

What is the Optimal Size of the Quantum Region in Embedding Calculations of Two-Photon Absorption Spectra of Fluorescent Proteins?

We systematically investigate an impact of the size and content of a quantum (QM) region, treated at the density functional theory level, in embedding calculations on one- (OPA) and two-photon absorption (TPA) spectra of the following fluorescent proteins (FPs) models: Aequorea victoria green FP (avGFP) with neutral (avGFP-n) and anionic (avGFP-a) chromophore as well as Citrine FP. We find that amino acid (a.a.) residues as well as water molecules hydrogen-bonded (h-bonded) to the chromophore usually boost both OPA and TPA processes intensity. The presence of hydrophobic a.a. residues in the quantum region also non-negligibly affects both absorption spectra but decreases absorption intensity. We conclude that to reach a quantitative description of OPA and TPA spectra in multiscale modeling of FPs, the quantum region should consist of a chromophore and most of a.a. residues and water molecules in a radius of 0.30–0.35 nm (ca. 200–230 atoms) when the remaining part of the system is approximated by the electrostatic point-charges. The optimal size of the QM region can be reduced to 80–100 atoms by utilizing a more advanced polarizable embedding model. We also find components of the QM region that are specific to a FP under study. We propose that the F165 a.a. residue is important in tuning the TPA spectrum of avGFP-n but not other investigated FPs. In the case of Citrine, Y203 and M69 a.a. residues must definitely be part of the QM subsystem. Furthermore, we find that long-range electrostatic interactions between the QM region and the rest of the protein cannot be neglected even for the most extensive QM regions (ca. 350 atoms).

Two-Photon Absorption Cross-Sections in Fluorescent Proteins Containing Non-canonical Chromophores Using Polarizable QM/MM

Multi-photon absorption properties, particularly two-photon absorption (2PA), of fluorescent proteins (FPs) have made them attractive tools in deep-tissue clinical imaging. Although the diversity of photophysical properties for FPs is wide, there are some caveats predominant among the existing FP variants that need to be overcome, such as low quantum yields and small 2PA cross-sections. From a computational perspective, Salem et al. (2016) suggested the inclusion of non-canonical amino acids in the chromophore of the red fluorescent protein DsRed, through the replacement of the tyrosine amino acid. The 2PA properties of these new non-canonical chromophores (nCCs) were determined in vacuum, i.e., without taking into account the protein environment. However, in the computation of response properties, such as 2PA cross-sections, the environment plays an important role. To account for environment and protein-chromophore coupling effects, quantum mechanical/molecular mechanical (QM/MM) schemes can be useful. In this work, the polarizable embedding (PE) model is employed along with time-dependent density functional theory to describe the 2PA properties of a selected set of chromophores made from non-canonical amino acids as they are embedded in the DsRed protein matrix. The objective is to provide insights to determine whether or not the nCCs could be developed and, thus, generate a new class of FPs. Results from this investigation show that within the DsRed environment, the nCC 2PA cross-sections are diminished relative to their values in vacuum. However, further studies toward understanding the 2PA limit of these nCCs using different protein environments are needed.

CPPE: An Open-Source C++ and Python Library for Polarizable Embedding

We present a modular open-source library for polarizable embedding (PE) named CPPE. The library is implemented in C++, and it additionally provides a Python interface for rapid prototyping and experimentation in a high-level scripting language. Our library integrates seamlessly with existing quantum chemical program packages through an intuitive and minimal interface. Until now, CPPE has been interfaced to three packages, Q-Chem, Psi4, and PySCF. Furthermore, we show CPPE in action using all three program packages for a computational spectroscopy application. With CPPE, host program interfaces only require minor programming effort, paving the way for new combined methodologies and broader availability of the PE model.

Steady-State Linear and Non-linear Optical Spectroscopy of Organic Chromophores and Bio-macromolecules

Bio-macromolecules as DNA, lipid membranes and (poly)peptides are essential compounds at the core of biological systems. The development of techniques and methodologies for their characterization is therefore necessary and of utmost interest, even though difficulties can be experienced due to their intrinsic complex nature. Among these methods, spectroscopies, relying on optical properties are especially important to determine their macromolecular structures and behaviors, as well as the possible interactions and reactivity with external dyes—often drugs or pollutants—that can (photo)sensitize the bio-macromolecule leading to eventual chemical modifications, thus damages. In this review, we will focus on the theoretical simulation of electronic spectroscopies of bio-macromolecules, considering their secondary structure and including their interaction with different kind of (photo)sensitizers. Namely, absorption, emission and electronic circular dichroism (CD) spectra are calculated and compared with the available experimental data. Non-linear properties will be also taken into account by two-photon absorption, a highly promising technique (i) to enhance absorption in the red and infra-red windows and (ii) to enhance spatial resolution. Methodologically, the implications of using implicit and explicit solvent, coupled to quantum and thermal samplings of the phase space, will be addressed. Especially, hybrid quantum mechanics/molecular mechanics (QM/MM) methods are explored for a comparison with solely QM methods, in order to address the necessity to consider an accurate description of environmental effects on spectroscopic properties of biological systems.

Open-ended formulation of self-consistent field response theory with the polarizable continuum model for solvation

The study of high-order absorption properties of molecules is a field of growing importance. Quantum-chemical studies can help design chromophores with desirable characteristics. Given that most experiments are performed in solution, it is important to devise a cost-effective strategy to include solvation effects in quantum-chemical studies of these properties. We here present an open-ended formulation of self-consistent field (SCF) response theory for a molecular solute coupled to a polarizable continuum model (PCM) description of the solvent. Our formulation relies on the open-ended, density matrix-based quasienergy formulation of SCF response theory of Thorvaldsen, et al., [J. Chem. Phys., 2008, 129, 214108] and the variational formulation of the PCM, as presented by Lipparini et al., [J. Chem. Phys., 2010, 133, 014106]. Within the PCM approach to solvation, the mutual solute-solvent polarization is represented by means of an apparent surface charge (ASC) spread over the molecular cavity defining the solute-solvent boundary. In the variational formulation, the ASC is an independent, variational degree of freedom. This allows us to formulate response theory for molecular solutes in the fixed-cavity approximation up to arbitrary order and with arbitrary perturbation operators. For electric dipole perturbations, pole and residue analyses of the response functions naturally lead to the identification of excitation energies and transition moments. We document the implementation of this approach in the Dalton program package using a recently developed open-ended response code and the PCMSolver libraries and present results for one-, two-, three-, four- and five-photon absorption processes of three small molecules in solution.