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

Unravelling the ultrafast dynamics of a N-BODIPY compound

Although the photophysics of BODIPY compounds has been widely investigated in the last few years, their analogues N-BODIPY, with nitrogen substitution at the boron center, did not receive comparable attention. In this work we report the synthesis and photochemical characterization of a substituted N-BODIPY compound, by means of a combined theoretical and spectroscopic approach. Compared to a standard BODIPY, the compound under investigation presents a lower fluorescence quantum yield (QY) in the visible region. The excited state relaxation dynamics of the dye was studied in different solvents, showing further fluorescence quenching in polar solvents, and excited state decay rates strongly dependent on the environment polarity. The role of the pendant moieties and the involvement of charge transfer states in the excited state dynamics was experimentally addressed by transient absorption spectroscopy, and further analyzed with TD-DFT calculations, which allowed precise assignment of the transient signals to the correspondent electronic configuration. The complete picture of the N-BODIPY behavior shows the presence of both charge transfer and localized states, influencing the observed photophysics to different amounts, depending on the excitation conditions and the surrounding environment

Fine control of chlorophyll-carotenoid interactions defines the functionality of light-harvesting proteins in plants

V.B. and C.D.P.D. acknowledge the support from the Leverhulme Trust RPG-2015-337. This research utilized Queen Mary’s MidPlus computational facilities, supported by QMUL Research-IT and funded by EPSRC grant EP/K000128/1. W.P.B acknowledges support from the Photosynthetic Antenna Research Center (PARC), an Energy Frontier Research Center funded by the U.S. Department of Energy, Office of Science, Office of Basic Energy Sciences under Award DE-SC0001035 for initial development of the TDC calculation code, as well as support from Army Research Office (ARO-MURI) Award W911NF1210420 for further development

A polarisable QM/MM description of NMR chemical shifts of a photoreceptor protein

We present a polarisable QM/MM investigation of NMR chemical shifts of a photoreceptor protein belonging to the Blue Light-Using Flavin family. Two different structures have been proposed for this photoreceptor which show a large variability in terms of the position and orientation of the protein residues around the flavin chromophore. Here, the two structures have been investigated with our multiscale approach using both DFT and MP2 level of theory. The picture that comes out comparing the (Formula presented.) H chemical shifts of the flavin and the most strongly interacting protein residues with the available experimental data, indicates a different behaviour of the two structures, with one showing a better correlation with NMR measurements. This shows that hybrid quantum chemical-classical simulations of NMR chemical shifts can indeed become a valuable tool to investigate the structure of complex biosystems

EXAT: EXcitonic analysis tool

We introduce EXcitonic Analysis Tool (EXAT), a program able to compute optical spectra of large excitonic systems directly from the output of quantum mechanical calculations performed with the popular Gaussian 16 package. The software is able to combine in an excitonic scheme the single-chromophore properties and exciton couplings to simulate energies, coefficients, and excitonic spectra (UV-vis, CD, and LD). The effect of the environment can also be included using a Polarizable Continuum Model. EXAT also presents a simple graphical user interface, which shows on-screen both site and exciton properties. To show the potential of the method, we report two applications on a a chiral perturbed BODIPY system and DNA G-quadruplexes, respectively. The program is available online at . (c) 2017 Wiley Periodicals, Inc

Coupling to Charge Transfer States is the Key to Modulate the Optical Bands for Efficient Light Harvesting in Purple Bacteria

The photosynthetic apparatus of purple bacteria uses exciton delocalization and static disorder to modulate the position and broadening of its absorption bands, leading to efficient light harvesting. Its main antenna complex, LH2, contains two rings of identical bacteriochlorophyll pigments, B800 and B850, absorbing at 800 and 850 nm, respectively. It has been an unsolved problem why static disorder of the strongly coupled B850 ring is several times larger than that of the B800 ring. Here we show that mixing between excitons and charge transfer states in the B850 ring is responsible for the effect. The linear absorption spectrum of the LH2 system is simulated by using a multiscale approach with an exciton Hamiltonian generalized to include the charge transfer states that involve adjacent pigment pairs, with static disorder modeled microscopically by molecular dynamics simulations. Our results show that sufficient inhomogeneous broadening of the B850 band, needed for efficient light harvesting, is only obtained by utilizing static disorder in the coupling between local excited and interpigment charge transfer states

The modeling of the absorption lineshape for embedded molecules through a polarizable QM/MM approach

We present a computational strategy to simulate the absorption lineshape of a molecule embedded in a complex environment by using a polarizable QM/MM approach. This strategy is presented in two alternative formulations, one based on a molecular dynamics simulation of the structural fluctuations of the system and the other using normal modes and harmonic frequencies calculated on optimized geometries. The comparison for the case of a chromophore within a strongly inhomogeneous and structured environment, namely the intercalation pocket of DNA, shows that the MD-based approach is able to reproduce the experimental spectral bandshape. In contrast, the static approach overestimates the vibronic coupling, resulting in a much broader band

Single-chain self-folding in an amphiphilic copolymer: An integrated experimental and computational study

An amphiphilic random copolymer of hydrophilic poly(ethylene glycol) methyl ether methacrylate with hydrophobic perfluorohexylethyl acrylate, PEGMA77-co-FA23, was synthesized by ATRP and used to investigate self-assembling into nanostructures in water and chloroform solutions, both experimentally and computationally. The dynamic light scattering measurements on water solutions of the copolymer at room temperature evidenced the presence of nanoassemblies with hydrodynamic diameter D h = 4 ± 1 nm. The behavior of fluorescence emission intensity of water solutions with added ethidium bromide suggested confinement of the molecular rotor within a hydrophobic environment of the copolymer. Moreover, these nanoassemblies were thermoresponsive and reversibly collapsed into much larger, multi-chain aggregates with D h = 390 ± 20 nm at a critical temperature of 55 °C (5 mg mL −1 ). Molecular dynamics simulations revealed the formation of single-chain, prolate globular nanoassemblies with a structural variability in water solution at room temperature. The evolutions of the simulated radius of gyration (R g ), asphericity, prolateness and solvent-accessible surface area were analyzed along the folding trajectories. Thus, self-folding appeared to result from the interplay between hydrophobic interactions and structural constraints which leads to rather complex nanostructures (R g = 20–25 Å 200 ns simulation). By contrast, folding in much more open polymer conformations (R g = 30–40 Å) was predicted for chloroform solutions