The key to the yellow-to-cyan tuning in the green fluorescent protein family is polarisation

Computational approaches have to date failed to fully capture the large (about 0.4 eV) excitation energy tuning displayed by the nearly identical anionic chromophore in different green fluorescent protein (GFP) variants. Here, we present a thorough comparative study of a set of proteins in this sub-family, including the most red- (phiYFP) and blue-shifted (mTFP0.7) ones. We employ a classical polarisable embedding through induced dipoles and combine it with time-dependent density functional theory and multireference perturbation theory in order to capture both state-specific induction contributions and the coupling of the polarisation of the protein to the chromophore transition density. The obtained results show that only upon inclusion of both these two effects generated by the mutual polarisation between the chromophore and the protein can the full spectral tuning be replicated. We finally discuss how this mutual polarisation affects the correlation between excitation energies, dipole moment variation, and molecular electrostatic field

Exciton properties and optical spectra of light harvesting complex II from a fully atomistic description

We present a fully atomistic simulation of linear optical spectra (absorption, fluorescence and circular dichroism) of the Light Harvesting Complex II (LHCII) trimer using a hybrid approach, which couples a quantum chemical description of the chlorophylls with a classical model for the protein and the external environment (membrane and water). The classical model uses a polarizable Molecular Mechanics force field, thus allowing mutual polarization effects in the calculations of the excitonic properties. The investigation is performed both on the crystal structure and on structures generated by a μs long classical molecular dynamics simulation of the complex within a solvated membrane. The results show that this integrated approach not only provides a good description of the excitonic properties and optical spectra without the need for additional refinements of the excitonic parameters, but it also allows an atomistic investigation of the relative importance of electronic, structural and environment effects in determining the optical spectra

An enhanced sampling QM/AMOEBA approach: The case of the excited state intramolecular proton transfer in solvated 3-hydroxyflavone

We present an extension of the polarizable quantum mechanical (QM)/AMOEBA approach to enhanced sampling techniques. This is achieved by connecting the enhanced sampling PLUMED library to the machinery based on the interface of Gaussian and Tinker to perform QM/AMOEBA molecular dynamics. As an application, we study the excited state intramolecular proton transfer of 3-hydroxyflavone in two solvents: methanol and methylcyclohexane. By using a combination of molecular dynamics and umbrella sampling, we find an ultrafast component of the transfer, which is common to the two solvents, and a much slower component, which is active in the protic solvent only. The mechanisms of the two components are explained in terms of intramolecular vibrational redistribution and intermolecular hydrogen-bonding, respectively. Ground and excited state free energies along an effective reaction coordinate are finally obtained allowing for a detailed analysis of the solvent mediated mechanism

Energy, Structures, and Response Properties with a Fully Coupled QM/AMOEBA/ddCOSMO Implementation

We present the implementation of a fully coupled polarizable QM/MM/continuum model based on the AMOEBA polarizable force field and the domain decomposition implementation of the conductor-like screening model. Energies, response properties, and analytical gradients with respect to both QM and MM nuclear positions are available, and a generic, atomistic cavity can be employed. The model is linear scaling in memory requirements and computational cost with respect to the number of classical atoms and is therefore suited to model large, complex systems. Using three variants of the green-fluorescent protein, we investigate the overall computational cost of such calculations and the effect of the continuum model on the convergence of the computed properties with respect to the size of the embedding. We also demonstrate the fundamental role of polarization effects by comparing polarizable and nonpolarizable embeddings to fully QM ones

Structure of the stress-related LHCSR1 complex determined by an integrated computational strategy

Light-harvesting complexes (LHCs) are pigment-protein complexes whose main function is to capture sunlight and transfer the energy to reaction centers of photosystems. In response to varying light conditions, LH complexes also play photoregulation and photoprotection roles. In algae and mosses, a sub-family of LHCs, light-harvesting complex stress-related (LHCSR), is responsible for photoprotective quenching. Despite their functional and evolutionary importance, no direct structural information on LHCSRs is available that can explain their unique properties. In this work, we propose a structural model of LHCSR1 from the moss P. patens, obtained through an integrated computational strategy that combines homology modeling, molecular dynamics, and multiscale quantum chemical calculations. The model is validated by reproducing the spectral properties of LHCSR1. Our model reveals the structural specificity of LHCSR1, as compared with the CP29 LH complex, and poses the basis for understanding photoprotective quenching in mosses

Towards accurate solvation free energies of large biological systems

Continuum solvation models like PCM or COSMO are the standard tool to calculate solvation free energies in a quantum level, but have been typically limited to small biological molecules due to its large computational cost. Recently, a new implementation of COSMO based on a domain decomposition strategy (ddCOSMO) [1] has been presented, which speeds up calculations by several orders of magnitude, thus paving the way for its application to very large systems. Here, we report the parameterization of ddCOSMO to the prediction of hydration free energies based on the MST solvation model developed in Barcelona, [2][3]. The parameterization is based on the PM6 semi-empirical Hamiltonian, on a set of over 200 experimental hydration free energies. The new model opens the way to the accurate prediction of hydration free energies of very large biomolecules, thus going beyond the usual classical MM-PBSA or MM-GBSA approaches

Ultrafast Transient Infrared Spectroscopy of Photoreceptors with Polarizable QM/MM Dynamics

Ultrafast transient infrared (TRIR) spectroscopy is widely used to measure the excitation-induced structural changes of protein-bound chromophores. Here, we design a novel and general strategy to compute TRIR spectra of photoreceptors by combining μs-long MM molecular dynamics with ps-long QM/AMOEBA Born-Oppenheimer molecular dynamics (BOMD) trajectories for both ground and excited electronic states. As a proof of concept, the strategy is here applied to AppA, a blue-light-utilizing flavin (BLUF) protein, found in bacteria. We first analyzed the short-time evolution of the embedded flavin upon excitation revealing that its dynamic Stokes shift is ultrafast and mainly driven by the internal reorganization of the chromophore. A different normal-mode representation was needed to describe ground- and excited-state IR spectra. In this way, we could assign all of the bands observed in the measured transient spectrum. In particular, we could characterize the flavin isoalloxazine-ring region of the spectrum, for which a full and clear description was missing

A different perspective for nonphotochemical quenching in plant antenna complexes

Light-harvesting complexes of plants exert a dual function of light-harvesting (LH) and photoprotection through processes collectively called nonphotochemical quenching (NPQ). While LH processes are relatively well characterized, those involved in NPQ are less understood. Here, we characterize the quenching mechanisms of CP29, a minor LHC of plants, through the integration of two complementary enhanced-sampling techniques, dimensionality reduction schemes, electronic calculations and the analysis of cryo-EM data in the light of the predicted conformational ensemble. Our study reveals that the switch between LH and quenching state is more complex than previously thought. Several conformations of the lumenal side of the protein occur and differently affect the pigments’ relative geometries and interactions. Moreover, we show that a quenching mechanism localized on a single chlorophyll-carotenoid pair is not sufficient but many chlorophylls are simultaneously involved. In such a diffuse mechanism, short-range interactions between each carotenoid and different chlorophylls combined with a protein-mediated tuning of the carotenoid excitation energies have to be considered in addition to the commonly suggested Coulomb interactions

A different perspective for nonphotochemical quenching in plant antenna complexes

Light-harvesting complexes of plants exert a dual function of light-harvesting (LH) and photoprotection through processes collectively called nonphotochemical quenching (NPQ). While LH processes are relatively well characterized, those involved in NPQ are less understood. Here, we characterize the quenching mechanisms of CP29, a minor LHC of plants, through the integration of two complementary enhanced-sampling techniques, dimensionality reduction schemes, electronic calculations and the analysis of cryo-EM data in the light of the predicted conformational ensemble. Our study reveals that the switch between LH and quenching state is more complex than previously thought. Several conformations of the lumenal side of the protein occur and differently affect the pigments’ relative geometries and interactions. Moreover, we show that a quenching mechanism localized on a single chlorophyll-carotenoid pair is not sufficient but many chlorophylls are simultaneously involved. In such a diffuse mechanism, short-range interactions between each carotenoid and different chlorophylls combined with a protein-mediated tuning of the carotenoid excitation energies have to be considered in addition to the commonly suggested Coulomb interactions