Chromophore-Protein Coupling beyond Nonpolarizable Models: Understanding Absorption in Green Fluorescent Protein

The nature of the coupling of the photoexcited chromophore with the environment in a prototypical system like green fluorescent protein (GFP) is to date not understood, and its description still defies state-of-the-art multiscale approaches. To identify which theoretical framework of the chromophore protein complex can realistically capture its essence, we employ here a variety of electronic-structure methods, namely, time-dependent density functional theory (TD-DFT), multireference perturbation theory (NEVPT2 and CASPT2), and quantum Monte Carlo (QMC) in combination with static point charges (QM/MM), DFT embedding (QM/DFT), and classical polarizable embedding through induced dipoles (QM/MMpol). Since structural modifications can significantly affect the photophysics of GFP, we also account for thermal fluctuations through extensive molecular dynamics simulations. We find that a treatment of the protein through static point charges leads to significantly blue-shifted excitation energies and that including thermal fluctuations does not cure the coarseness of the MM description. While TDDFT calculations on large cluster models indicate the need of a responsive protein, this response is not simply electrostatic: An improved description of the protein in the ground state or in response to the excitation of the chromophore via ground-state or state-specific DFT and MMpol embedding does not significantly modify the results obtained with static point charges. Through the use of QM/MMpol in a linear response formulation, a different picture in fact emerges in which the main environment response to the chromophore excitation is the one coupling the transition density and the corresponding induced dipoles. Such interaction leads to significant red-shifts and a satisfactory agreement with full QM cluster calculations at the same level of theory. Our findings demonstrate that, ultimately, faithfully capturing the effects of the environment in GFP requires a quantum treatment of large photoexcited regions but that a QM/classical model can be a useful approximation when extended beyond the electrostatic-only formulation

Introducing QMC/MMpol: Quantum Monte Carlo in Polarizable Force Fields for Excited States

We present for the first time a quantum mechanics/molecular mechanics scheme which combines quantum Monte Carlo with the reaction field of classical polarizable dipoles (QMC/MMpol). In our approach, the optimal dipoles are self-consistently generated at the variational Monte Carlo level and then used to include environmental effects in diffusion Monte Carlo. We investigate the performance of this hybrid model in describing the vertical excitation energies of prototypical small molecules solvated in water, namely, methylenecyclopropene and s-trans acrolein. Two polarization regimes are explored where either the dipoles are optimized with respect to the ground-state solute density (polGS) or different sets of dipoles are separately brought to equilibrium with the states involved in the electronic transition (polSS). By comparing with reference supermolecular calculations where both solute and solvent are treated quantum mechanically, we find that the inclusion of the response of the environment to the excitation of the solute leads to superior results than the use of a frozen environment (point charges or polGS), in particular, when the solute-solvent coupling is dominated by electrostatic effects which are well recovered in the polSS condition. QMC/MMpol represents therefore a robust scheme to treat important environmental effects beyond static point charges, combining the accuracy of QMC with the simplicity of a classical approach

Electronic charge and spin density distribution in a quantum ring with spin-orbit and Coulomb interactions

Charge and spin density distributions are studied within a nano-ring structure endowed with Rashba and Dresselhaus spin orbit coupling (SOI). For a small number of interacting electrons, in the presence of an external magnetic field, the energy spectrum of the system is calculated through an exact numerical diagonalization procedure. The eigenstates thus determined are used to estimate the charge and spin densities around the ring. We find that when more than two electrons are considered, the charge-density deformations induced by SOI are dramatically flattened by the Coulomb repulsion, while the spin density ones are amplified.Comment: 7 pages, 7 figure

Understanding Conformational Dynamics of Complex Lipid Mixtures Relevant to Biology

This is a perspective article entitled “Frontiers in computational biophysics: understanding conformational dynamics of complex lipid mixtures relevant to biology” which is following a CECAM meeting with the same name.Fil: Friedman, Ran. Linnæus University; ArgentinaFil: Khalid, Syma. University of Southampton; Reino UnidoFil: Aponte Santamaría, Camilo. Ruprecht-Karls-Universität Heidelberg; Alemania. Universidad de los Andes; ColombiaFil: Arutyunova, Elena. University of Alberta; CanadáFil: Becker, Marlon. Westfälische Wilhelms Universität; AlemaniaFil: Boyd, Kevin J.. University of Connecticut; Estados UnidosFil: Christensen, Mikkel. University Aarhus; DinamarcaFil: Coimbra, João T. S.. Universidad de Porto; PortugalFil: Concilio, Simona. Universita di Salerno; ItaliaFil: Daday, Csaba. Heidelberg Institute for Theoretical Studies; AlemaniaFil: Eerden, Floris J. van. University of Groningen; Países BajosFil: Fernandes, Pedro A.. Universidad de Porto; PortugalFil: Gräter, Frauke. Heidelberg University; Alemania. Heidelberg Institute for Theoretical Studies; AlemaniaFil: Hakobyan, Davit. Westfälische Wilhelms Universität; AlemaniaFil: Heuer, Andreas. Westfälische Wilhelms Universität; AlemaniaFil: Karathanou, Konstantina. Freie Universität Berlin; AlemaniaFil: Keller, Fabian. Westfälische Wilhelms Universität; AlemaniaFil: Lemieux, M. Joanne. University of Alberta; CanadáFil: Marrink, Siewert J.. University of Groningen; Países BajosFil: May, Eric R.. University of Connecticut; Estados UnidosFil: Mazumdar, Antara. University of Groningen; Países BajosFil: Naftalin, Richard. Colegio Universitario de Londres; Reino UnidoFil: Pickholz, Mónica Andrea. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Instituto de Física de Buenos Aires. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Instituto de Física de Buenos Aires; ArgentinaFil: Piotto, Stefano. Universita di Salerno; ItaliaFil: Pohl, Peter. Johannes Kepler University; AustriaFil: Quinn, Peter. Colegio Universitario de Londres; Reino UnidoFil: Ramos, Maria J.. Universidad de Porto; PortugalFil: Schiøtt, Birgit. University Aarhus; DinamarcaFil: Sengupta, Durba. National Chemical Laboratory India; IndiaFil: Sessa, Lucia. Universita di Salerno; ItaliaFil: Vanni, Stefano. University Of Fribourg;Fil: Zeppelin, Talia. University Aarhus; DinamarcaFil: Zoni, Valeria. University of Fribourg; SuizaFil: Bondar, Ana-Nicoleta. Freie Universität Berlin; AlemaniaFil: Domene, Carmen. University of Oxford; Reino Unido. University of Bath; Reino Unid

Coulomb and Spin-Orbit Interaction Effects in a Mesoscopic Ring

We study a ring structure on a two-dimensional electron gas with both Rashba and Dresselhaus spin-orbit interaction in an external magnetic field. The competition between the two kinds of SOI leads to a deformation of the charge and spin densities around the ring. After examining the one-electron eigenstates we investigate Coulomb effects. We use two different ring models: a continuous model with analytical single particle states (Tan-Inkson) and a discrete model (tight-binding) which is more convenient for many-body calculations. We show that for more than two electrons the deformation of the charge density is smeared out by the Coulomb interaction, whilst the deformation of the spin density is amplified.Rannsóknasjóðu

Chromophore-protein coupling beyond nonpolarizable models: understanding absorption in green fluorescent protein

The nature of the coupling of the photoexcited chromophore with the environment in a prototypical system like green fluorescent protein (GFP) is to date not understood, and its description still defies state-of-the-art multiscale approaches. To identify which theoretical framework of the chromophore–protein complex can realistically capture its essence, we employ here a variety of electronic-structure methods, namely, time-dependent density functional theory (TD-DFT), multireference perturbation theory (NEVPT2 and CASPT2), and quantum Monte Carlo (QMC) in combination with static point charges (QM/MM), DFT embedding (QM/DFT), and classical polarizable embedding through induced dipoles (QM/MMpol). Since structural modifications can significantly affect the photophysics of GFP, we also account for thermal fluctuations through extensive molecular dynamics simulations. We find that a treatment of the protein through static point charges leads to significantly blue-shifted excitation energies and that including thermal fluctuations does not cure the coarseness of the MM description. While TDDFT calculations on large cluster models indicate the need of a responsive protein, this response is not simply electrostatic: An improved description of the protein in the ground state or in response to the excitation of the chromophore via ground-state or state-specific DFT and MMpol embedding does not significantly modify the results obtained with static point charges. Through the use of QM/MMpol in a linear response formulation, a different picture in fact emerges in which the main environment response to the chromophore excitation is the one coupling the transition density and the corresponding induced dipoles. Such interaction leads to significant red-shifts and a satisfactory agreement with full QM cluster calculations at the same level of theory. Our findings demonstrate that, ultimately, faithfully capturing the effects of the environment in GFP requires a quantum treatment of large photoexcited regions but that a QM/classical model can be a useful approximation when extended beyond the electrostatic-only formulation

State-Specific Embedding Potentials for Excitation-Energy Calculations

Embedding potentials are frequently used to describe the effect of an environment on the electronic structure of molecules in larger systems, including their excited states. If such excitations are accompanied by significant rearrangements in the electron density of the embedded molecule, large differential polarization effects may take place, which in turn can require state-specific embedding potentials for an accurate theoretical description. We outline here how to extend wave function in density functional theory (WF/DFT) methods to compute the excitation energies of a molecule in a responsive environment through the use of state-specific density-based embedding potentials constructed within a modified subsystem DFT approach. We evaluate the general expression of the ground- and excited-state energy difference of the total system both with the use of state-independent and state-dependent embedding potentials and propose some practical recipes to construct the approximate excited-state DFT density of the active part used to polarize the environment. We illustrate these concepts with the state-independent and state-dependent WF/DFT computation of the excitation energies of <i>p</i>-nitroaniline, acrolein, methylenecyclopropene, and <i>p</i>-nitrophenolate in various solvents