Search

When the enhanced electromagnetic field of a confined light mode interacts with photoactive molecules, the system can be driven into the regime of strong coupling, where new hybrid light–matter states, polaritons, are formed. Polaritons, manifested by the Rabi split in the dispersion, have shown potential for controlling the chemistry of the coupled molecules. Here, we show by angle-resolved steady-state experiments accompanied by multi-scale molecular dynamics simulations that the molecular Stokes shift plays a significant role in the relaxation of polaritons formed by organic molecules embedded in a polymer matrix within metallic Fabry–Pérot cavities. Our results suggest that in the case of Rhodamine 6G, a dye with a significant Stokes shift, excitation of the upper polariton leads to a rapid localization of the energy into the fluorescing state of one of the molecules, from where the energy scatters into the lower polariton (radiative pumping), which then emits. In contrast, for excitonic J-aggregates with a negligible Stokes shift, the fluorescing state does not provide an efficient relaxation gateway. Instead, the relaxation is mediated by exchanging energy quanta matching the energy gap between the dark states and lower polariton into vibrational modes (vibrationally assisted scattering). To understand better how the fluorescing state of a molecule that is not strongly coupled to the cavity can transfer its excitation energy to the lower polariton in the radiative pumping mechanism, we performed multi-scale molecular dynamics simulations. The results of these simulations suggest that non-adiabatic couplings between uncoupled molecules and the polaritons are the driving force for this energy transfer process. [ABSTRACT FROM AUTHOR]

The high intensity of free electron lasers opens up the possibility to perform single-shot molecule scattering experiments. However, even for small molecules, radiation damage induced by absorption of high intense x-ray radiation is not yet fully understood. One of the striking effects which occurs under intense x-ray illumination is the creation of double core ionized molecules in considerable quantity. To provide insight into this process, we have studied the dynamics of water molecules in single and double core ionized states by means of electronic transition rate calculations and ab initio molecular dynamics (MD) simulations. From the MD trajectories, photoionization and Auger transition rates were computed based on electronic continuum wavefunctions obtained by explicit integration of the coupled radial Schrödinger equations. These rates served to solve the master equations for the populations of the relevant electronic states. To account for the nuclear dynamics during the core hole lifetime, the calculated electron emission spectra for different molecular geometries were incoherently accumulated according to the obtained time-dependent populations, thus neglecting possible interference effects between different decay pathways. We find that, in contrast to the single core ionized water molecule, the nuclear dynamics for the double core ionized water molecule during the core hole lifetime leaves a clear fingerprint in the resulting electron emission spectra. The lifetime of the double core ionized water was found to be significantly shorter than half of the single core hole lifetime. [ABSTRACT FROM AUTHOR]

X-ray crystallography is an important technique to determine the positions of atoms in a protein crystal. However, because the native environment in which proteins function, is not a crystal, but a solution, it is not a priori clear if the crystal structure represents the functional form of the protein. Because the protein structure and function often depend critically on the pH, the question arises whether proton affinities are affected by crystallization. X-ray diffraction usually does not reveal protons, which makes it difficult to experimentally measure p K a shifts in crystals. Here, we investigate whether this challenge can be addressed by performing in silico titration with constant pH molecular dynamics (MD) simulations. We compare the computed p K a values of proteins between solution and crystal environment and analyze these differences in the context of molecular interactions. For the proteins considered in this work, p K a shifts were mostly found for residues at the crystal interfaces, where the environment is more apolar in the crystal than in water. Although convergence was an obstacle, our simulations suggest that in principle it is possible to apply constant pH MD to protein crystals routinely and assess the effect of crystallization on protein function more systematically than with standard MD simulations. We also highlight technical challenges that need to be addressed to make MD simulations of crystals more reliable.

Collective strong coupling of many molecules to the confined light modes of an optical resonator can influence the photochemistry of these molecules, but the origin of this effect is not yet fully understood. To provide atomistic insights, several approaches have been developed based on quantum chemistry or molecular dynamics methods. However, most of these methods rely on coupling a few molecules (or sometimes only one) to a single cavity mode. To reach the strong coupling regime with such a small number of molecules, much larger vacuum field strengths are employed than in experiments. To keep the vacuum field realistic and avoid potential artefacts, the number of coupled molecules should be significantly increased instead, but that is not always possible due to restrictions on computational hardware and software. To overcome this barrier and model the dynamics of an arbitrarily large ensemble of molecules coupled to realistic cavity fields in atomistic molecular dynamics simulations, we propose to coarse-grain subsets of molecules into one or more effective supermolecules with an enhanced dipole moment and concerted dynamics. To verify the validity of the proposed multiscale model, we performed simulations in which we investigated how the number of molecules that are coupled to the cavity affects excited-state intra-molecular proton transfer, polariton relaxation, and exciton transport., (© 2024 Author(s). Published under an exclusive license by AIP Publishing.)

Molecular dynamics (MD) simulations are widely applied to estimate absolute binding free energies of protein-ligand and protein-protein complexes. A routinely used method for binding free energy calculations with MD is umbrella sampling (US), which calculates the potential of mean force (PMF) along a single reaction coordinate. Surprisingly, in spite of its widespread use, few validation studies have focused on the convergence of the free energy computed along a single path for specific cases, not addressing the reproducibility of such calculations in general. In this work, we therefore investigate the reproducibility and convergence of US along a standard distance-based reaction coordinate for various protein-protein and protein-ligand complexes, following commonly used guidelines for the setup. We show that repeating the complete US workflow can lead to differences of 2-20 kcal/mol in computed binding free energies. We attribute those discrepancies to small differences in the binding pathways. While these differences are unavoidable in the established US protocol, the popularity of the latter could hint at a lack of awareness of such reproducibility problems. To test if the convergence of PMF profiles can be improved if multiple pathways are sampled simultaneously, we performed additional simulations with an adaptive-biasing method, here the accelerated weight histogram (AWH) approach. Indeed, the PMFs obtained from AHW simulations are consistent and reproducible for the systems tested. To the best of our knowledge, our work is the first to attempt a systematic assessment of the pitfalls in one the most widely used protocols for computing binding affinities. We anticipate therefore that our results will provide an incentive for a critical reassessment of the validity of PMFs computed with US, and make a strong case to further benchmark the performance of adaptive-biasing methods for computing binding affinities.

Strong coupling between molecules and confined light modes of optical cavities to form polaritons can alter photochemistry, but the origin of this effect remains largely unknown. While theoretical models suggest a suppression of photochemistry due to the formation of new polaritonic potential energy surfaces, many of these models do not account for the energetic disorder among the molecules, which is unavoidable at ambient conditions. Here, we combine simulations and experiments to show that for an ultra-fast photochemical reaction such thermal disorder prevents the modification of the potential energy surface and that suppression is due to radiative decay of the lossy cavity modes. We also show that the excitation spectrum under strong coupling is a product of the excitation spectrum of the bare molecules and the absorption spectrum of the molecule-cavity system, suggesting that polaritons can act as gateways for channeling an excitation into a molecule, which then reacts normally. Our results therefore imply that strong coupling provides a means to tune the action spectrum of a molecule, rather than to change the reaction., (© 2024. The Author(s).)

High-spin molecules allow for bottom-up qubit design and are promising platforms for magnetic sensing and quantum information science. Optical addressability of molecular electron spins has also been proposed in first-row transition-metal complexes via optically detected magnetic resonance (ODMR) mechanisms analogous to the diamond-nitrogen-vacancy color center. However, significantly less progress has been made on the front of metal-free molecules, which can deliver lower costs and milder environmental impacts. At present, most luminescent open-shell organic molecules are π-diradicals, but such systems often suffer from poor ground-state open-shell characters necessary to realize a stable ground-state molecular qubit. In this work, we use alternancy symmetry to selectively minimize radical-radical interactions in the ground state, generating π-systems with high diradical characters. We call them m -dimers, referencing the need to covalently link two benzylic radicals at their meta carbon atoms for the desired symmetry. Through a detailed electronic structure analysis, we find that the excited states of alternant hydrocarbon m -diradicals contain important symmetries that can be used to construct ODMR mechanisms leading to ground-state spin polarization. The molecular parameters are set in the context of a tris(2,4,6-trichlorophenyl)methyl (TTM) radical dimer covalently tethered at the meta position, demonstrating the feasibility of alternant m -diradicals as molecular color centers.

Molecular dynamics (MD) simulations are routinely performed of biomolecules in solution, because this is their native environment. However, the structures used in such simulations are often obtained with X-ray crystallography, which provides the atomic coordinates of the biomolecule in a crystal environment. With the advent of free electron lasers and time-resolved techniques, X-ray crystallography can now also access metastable states that are intermediates in a biochemical process. Such experiments provide additional data, which can be used, for example, to optimize MD force fields. Doing so requires that the simulation of the biomolecule is also performed in the crystal environment. However, in contrast to simulations of biomolecules in solution, setting up a crystal is challenging. In particular, because not all solvent molecules are resolved in X-ray crystallography, adding a suitable number of solvent molecules, such that the properties of the crystallographic unit cell are preserved in the simulation, can be difficult and typically is a trial-and-error based procedure requiring manual interventions. Such interventions preclude high throughput applications. To overcome this bottleneck, we introduce gmXtal, a tool for setting up crystal simulations for MD simulations with GROMACS. With the information from the protein data bank (rcsb.org) gmXtal automatically (i) builds the crystallographic unit cell; (ii) sets the protonation of titratable residues; (iii) builds missing residues that were not resolved experimentally; and (iv) adds an appropriate number of solvent molecules to the system. gmXtal is available as a standalone tool https://gitlab.com/pbuslaev/gmxtal ., (© 2023. The Author(s).)

The observation that materials can change their properties when placed inside or near an optical resonator has sparked a fervid interest in understanding the effects of strong light-matter coupling on molecular dynamics, and several approaches have been proposed to extend the methods of computational chemistry into this regime. Whereas the majority of these approaches have focused on modeling a single molecule coupled to a single cavity mode, changes to chemistry have so far only been observed experimentally when very many molecules are coupled collectively to multiple modes with short lifetimes. While atomistic simulations of many molecules coupled to multiple cavity modes have been performed with semi-classical molecular dynamics, an explicit description of cavity losses has so far been restricted to simulations in which only a very few molecular degrees of freedom were considered. Here, we have implemented an effective non-Hermitian Hamiltonian to explicitly treat cavity losses in large-scale semi-classical molecular dynamics simulations of organic polaritons and used it to perform both mean-field and surface hopping simulations of polariton relaxation, propagation, and energy transfer., (© 2024 Author(s). Published under an exclusive license by AIP Publishing.)

Constant pH molecular dynamics (MD) is a powerful technique that allows the protonation state of residues to change dynamically, thereby enabling the study of pH dependence in a manner that has not been possible before. Recently, a constant pH implementation was incorporated into the GROMACS MD package. Although this implementation provides good accuracy and performance, manual modification and the preparation of simulation input files are required, which can be complicated, tedious, and prone to errors. To simplify and automate the setup process, we present phbuilder, a tool that automatically prepares constant pH MD simulations for GROMACS by modifying the input structure and topology as well as generating the necessary parameter files. phbuilder can prepare constant pH simulations from both initial structures and existing simulation systems, and it also provides functionality for performing titrations and single-site parametrizations of new titratable group types. The tool is freely available at www.gitlab.com/gromacs-constantph. We anticipate that phbuilder will make constant pH simulations easier to set up, thereby making them more accessible to the GROMACS user community.

The photoisomerization reaction of a fluorescent protein chromophore occurs on the ultrafast timescale. The structural dynamics that result from femtosecond optical excitation have contributions from vibrational and electronic processes and from reaction dynamics that involve the crossing through a conical intersection. The creation and progression of the ultrafast structural dynamics strongly depends on optical and molecular parameters. When using X-ray crystallography as a probe of ultrafast dynamics, the origin of the observed nuclear motions is not known. Now, high-resolution pump-probe X-ray crystallography reveals complex sub-ångström, ultrafast motions and hydrogen-bonding rearrangements in the active site of a fluorescent protein. However, we demonstrate that the measured motions are not part of the photoisomerization reaction but instead arise from impulsively driven coherent vibrational processes in the electronic ground state. A coherent-control experiment using a two-colour and two-pulse optical excitation strongly amplifies the X-ray crystallographic difference density, while it fully depletes the photoisomerization process. A coherent control mechanism was tested and confirmed the wave packets assignment., (© 2023. The Author(s).)

Transport of excitons in organic materials can be enhanced through polariton formation when the interaction strength between these excitons and the confined light modes of an optical resonator exceeds their decay rates. While the polariton lifetime is determined by the Q(uality)-factor of the optical resonator, the polariton group velocity is not. Instead, the latter is solely determined by the polariton dispersion. Yet, experiments suggest that the Q-factor also controls the polariton propagation velocity. To understand this observation, the authors perform molecular dynamics simulations of Rhodamine chromophores strongly coupled to Fabry-Pérot cavities with various Q-factors. The results suggest that propagation in the aforementioned experiments is initially dominated by ballistic motion of upper polariton states at their group velocities, which leads to a rapid expansion of the wavepacket. Cavity decay in combination with non-adiabatic population transfer into dark states, rapidly depletes these bright states, causing the wavepacket to contract. However, because population transfer is reversible, propagation continues, but as a diffusion process, at lower velocity. By controlling the lifetime of bright states, the Q-factor determines the duration of the ballistic phase and the diffusion coefficient in the diffusive regime. Thus, polariton propagation in organic microcavities can be effectively tuned through the Q-factor., (© 2023 The Authors. Advanced Science published by Wiley-VCH GmbH.)

Exciton transport can be enhanced in the strong coupling regime where excitons hybridize with confined light modes to form polaritons. Because polaritons have group velocity, their propagation should be ballistic and long-ranged. However, experiments indicate that organic polaritons propagate in a diffusive manner and more slowly than their group velocity. Here, we resolve this controversy by means of molecular dynamics simulations of Rhodamine molecules in a Fabry-Pérot cavity. Our results suggest that polariton propagation is limited by the cavity lifetime and appears diffusive due to reversible population transfers between polaritonic states that propagate ballistically at their group velocity, and dark states that are stationary. Furthermore, because long-lived dark states transiently trap the excitation, propagation is observed on timescales beyond the intrinsic polariton lifetime. These insights not only help to better understand and interpret experimental observations, but also pave the way towards rational design of molecule-cavity systems for coherent exciton transport., (© 2023. Springer Nature Limited.)

Chromophore cis/trans photoisomerization is a fundamental process in chemistry and in the activation of many photosensitive proteins. A major task is understanding the effect of the protein environment on the efficiency and direction of this reaction compared to what is observed in the gas and solution phases. In this study, we set out to visualize the hula twist (HT) mechanism in a fluorescent protein, which is hypothesized to be the preferred mechanism in a spatially constrained binding pocket. We use a chlorine substituent to break the twofold symmetry of the embedded phenolic group of the chromophore and unambiguously identify the HT primary photoproduct. Through serial femtosecond crystallography, we then track the photoreaction from femtoseconds to the microsecond regime. We observe signals for the photoisomerization of the chromophore as early as 300 fs, obtaining the first experimental structural evidence of the HT mechanism in a protein on its femtosecond-to-picosecond timescale. We are then able to follow how chromophore isomerization and twisting lead to secondary structure rearrangements of the protein β-barrel across the time window of our measurements.

Various approaches have been proposed to include the effect of pH in molecular dynamics (MD) simulations. Among these, the λ-dynamics approach proposed by Brooks and co-workers [Kong, X.; Brooks III, C. L. J. Chem. Phys. 1996 , 105 , 2414-2423] can be performed with little computational overhead and hfor each typeence be used to routinely perform MD simulations at microsecond time scales, as shown in the accompanying paper [Aho, N. et al. J. Chem. Theory Comput. 2022 , DOI: 10.1021/acs.jctc.2c00516]. At such time scales, however, the accuracy of the molecular mechanics force field and the parametrization becomes critical. Here, we address these issues and provide the community with guidelines on how to set up and perform long time scale constant pH MD simulations. We found that barriers associated with the torsions of side chains in the CHARMM36m force field are too high for reaching convergence in constant pH MD simulations on microsecond time scales. To avoid the high computational cost of extending the sampling, we propose small modifications to the force field to selectively reduce the torsional barriers. We demonstrate that with such modifications we obtain converged distributions of both protonation and torsional degrees of freedom and hence consistent p K a estimates, while the sampling of the overall configurational space accessible to proteins is unaffected as compared to normal MD simulations. We also show that the results of constant pH MD depend on the accuracy of the correction potentials. While these potentials are typically obtained by fitting a low-order polynomial to calculated free energy profiles, we find that higher order fits are essential to provide accurate and consistent results. By resolving problems in accuracy and sampling, the work described in this and the accompanying paper paves the way to the widespread application of constant pH MD beyond p K a prediction.

Molecular dynamics (MD) computer simulations are used routinely to compute atomistic trajectories of complex systems. Systems are simulated in various ensembles, depending on the experimental conditions one aims to mimic. While constant energy, temperature, volume, and pressure are rather straightforward to model, pH, which is an equally important parameter in experiments, is more difficult to account for in simulations. Although a constant pH algorithm based on the λ-dynamics approach by Brooks and co-workers [Kong, X.; Brooks III, C. L. J. Chem. Phys. 1996 , 105 , 2414-2423] was implemented in a fork of the GROMACS molecular dynamics program, uptake has been rather limited, presumably due to the poor scaling of that code with respect to the number of titratable sites. To overcome this limitation, we implemented an alternative scheme for interpolating the Hamiltonians of the protonation states that makes the constant pH molecular dynamics simulations almost as fast as a normal MD simulation with GROMACS. In addition, we implemented a simpler scheme, called multisite representation, for modeling side chains with multiple titratable sites, such as imidazole rings. This scheme, which is based on constraining the sum of the λ-coordinates, not only reduces the complexity associated with parametrizing the intramolecular interactions between the sites but also is easily extendable to other molecules with multiple titratable sites. With the combination of a more efficient interpolation scheme and multisite representation of titratable groups, we anticipate a rapid uptake of constant pH molecular dynamics simulations within the GROMACS user community.

Exciton transport in most organic materials is based on an incoherent hopping process between neighboring molecules. This process is very slow, setting a limit to the performance of organic optoelectronic devices. In this Article, we overcome the incoherent exciton transport by strongly coupling localized singlet excitations in a tetracene crystal to confined light modes in an array of plasmonic nanoparticles. We image the transport of the resulting exciton-polaritons in Fourier space at various distances from the excitation to directly probe their propagation length as a function of the exciton to photon fraction. Exciton-polaritons with an exciton fraction of 50% show a propagation length of 4.4 μm, which is an increase by 2 orders of magnitude compared to the singlet exciton diffusion length. This remarkable increase has been qualitatively confirmed with both finite-difference time-domain simulations and atomistic multiscale molecular dynamics simulations. Furthermore, we observe that the propagation length is modified when the dipole moment of the exciton transition is either parallel or perpendicular to the cavity field, which opens a new avenue for controlling the anisotropy of the exciton flow in organic crystals. The enhanced exciton-polariton transport reported here may contribute to the development of organic devices with lower recombination losses and improved performance., Competing Interests: The authors declare no competing financial interest., (© 2022 The Authors. Published by American Chemical Society.)

The strong light-matter coupling regime, in which excitations of materials hybridize with excitations of confined light modes into polaritons, holds great promise in various areas of science and technology. A key aspect for all applications of polaritonic chemistry is the relaxation into the lower polaritonic states. Polariton relaxation is speculated to involve two separate processes: vibrationally assisted scattering (VAS) and radiative pumping (RP), but the driving forces underlying these two mechanisms are not fully understood. To provide mechanistic insights, we performed multiscale molecular dynamics simulations of tetracene molecules strongly coupled to the confined light modes of an optical cavity. The results suggest that both mechanisms are driven by the same molecular vibrations that induce relaxation through nonadiabatic coupling between dark states and polaritonic states. Identifying these vibrational modes provides a rationale for enhanced relaxation into the lower polariton when the cavity detuning is resonant with specific vibrational transitions.

Photoactivation of bacteriophytochrome involves a cis-trans photoisomerization of a biliverdin chromophore, but neither the precise sequence of events nor the direction of the isomerization is known. Here, we used nonadiabatic molecular dynamics simulations on the photosensory protein dimer to resolve the isomerization mechanism in atomic detail. In our simulations the photoisomerization of the D ring occurs in the counterclockwise direction. On a subpicosecond time scale, the photoexcited chromophore adopts a short-lived intermediate with a highly twisted configuration stabilized by an extended hydrogen-bonding network. Within tens of picoseconds, these hydrogen bonds break, allowing the chromophore to adopt a more planar configuration, which we assign to the early Lumi-R state. The isomerization process is completed via helix inversion of the biliverdin chromophore to form the late Lumi-R state. The mechanistic insights into the photoisomerization process are essential to understand how bacteriophytochrome has evolved to mediate photoactivation and to engineer this protein for new applications.

Aptamer selection against novel infections is a complicated and time-consuming approach. Synergy can be achieved by using computational methods together with experimental procedures. This study aims to develop a reliable methodology for a rational aptamer in silico et vitro design. The new approach combines multiple steps: (1) Molecular design, based on screening in a DNA aptamer library and directed mutagenesis to fit the protein tertiary structure; (2) 3D molecular modeling of the target; (3) Molecular docking of an aptamer with the protein; (4) Molecular dynamics (MD) simulations of the complexes; (5) Quantum-mechanical (QM) evaluation of the interactions between aptamer and target with further analysis; (6) Experimental verification at each cycle for structure and binding affinity by using small-angle X-ray scattering, cytometry, and fluorescence polarization. By using a new iterative design procedure, structure- and interaction-based drug design (SIBDD), a highly specific aptamer to the receptor-binding domain of the SARS-CoV-2 spike protein, was developed and validated. The SIBDD approach enhances speed of the high-affinity aptamers development from scratch, using a target protein structure. The method could be used to improve existing aptamers for stronger binding. This approach brings to an advanced level the development of novel affinity probes, functional nucleic acids. It offers a blueprint for the straightforward design of targeting molecules for new pathogen agents and emerging variants., (© 2022 Wiley-VCH GmbH.)

Interactions between an atomically precise gold nanocluster Au 102 ( p -MBA) 44 ( p -MBA = para mercaptobenzoic acid) and a fluorescent organic dye molecule (KU, azadioxatriangulenium) are studied. In solution, the constituents form spontaneously a weakly bound complex leading to quenching of fluorescence of the KU dye via energy transfer. The KU can be separated from the complex by lowering pH, leading to recovery of fluorescence, which forms a basis for an optical reversible pH sensor. However, the sensor is not a stable entity, which could be delivered inside cells. For this purpose, a covalently bound hybrid is synthesized by linking the KU dye to the ligand layer of the cluster via an ester bond. Covalent linking facilitates entry of the cluster-dye hybrids into cells via endocytosis. Inside cells, the hybrids accumulate in endosomes where Au 102 releases its cargo via hydrolysis of the ester bond. Changes of the local pH inside endosomes regulate reversible fluorescence due to variations in the interactions between the Au 102 cluster and the dye. This work presents a concept for delivering reporter molecules into cells by using atomically precise gold nanoclusters as carriers and paves the way for future developments of cluster-reporter sensors for in vivo measurements of e.g. absolute pH values or ion concentrations., Competing Interests: There are no conflicts to declare., (This journal is © The Royal Society of Chemistry.)

Coupling molecules to the confined light modes of an optical cavity is showing great promise for manipulating chemical reactions. However, to fully exploit this principle and use cavities as a new tool for controlling chemistry, a complete understanding of the effects of strong light-matter coupling on molecular dynamics and reactivity is required. While quantum chemistry can provide atomistic insight into the reactivity of uncoupled molecules, the possibilities to also explore strongly coupled systems are still rather limited due to the challenges associated with an accurate description of the cavity in such calculations. Despite recent progress in introducing strong coupling effects into quantum chemistry calculations, applications are mostly restricted to single or simplified molecules in ideal lossless cavities that support a single light mode only. However, even if commonly used planar mirror micro-cavities are characterized by a fundamental mode with a frequency determined by the distance between the mirrors, the cavity energy also depends on the wave vector of the incident light rays. To account for this dependency, called cavity dispersion, in atomistic simulations of molecules in optical cavities, we have extended our multi-scale molecular dynamics model for strongly coupled molecular ensembles to include multiple confined light modes. To validate the new model, we have performed simulations of up to 512 Rhodamine molecules in red-detuned Fabry-Pérot cavities. The results of our simulations suggest that after resonant excitation into the upper polariton at a fixed wave vector, or incidence angle, the coupled cavity-molecule system rapidly decays into dark states that lack dispersion. Slower relaxation from the dark state manifold into both the upper and lower bright polaritons causes observable photo-luminescence from the molecule-cavity system along the two polariton dispersion branches that ultimately evolves toward the bottom of the lower polariton branch, in line with experimental observations. We anticipate that the more realistic cavity description in our approach will help to better understand and predict how cavities can modify molecular properties.

Photosynthetic reaction centres harvest the energy content of sunlight by transporting electrons across an energy-transducing biological membrane. Here we use time-resolved serial femtosecond crystallography 1 using an X-ray free-electron laser 2 to observe light-induced structural changes in the photosynthetic reaction centre of Blastochloris viridis on a timescale of picoseconds. Structural perturbations first occur at the special pair of chlorophyll molecules of the photosynthetic reaction centre that are photo-oxidized by light. Electron transfer to the menaquinone acceptor on the opposite side of the membrane induces a movement of this cofactor together with lower amplitude protein rearrangements. These observations reveal how proteins use conformational dynamics to stabilize the charge-separation steps of electron-transfer reactions.

Organisms use photo-receptors to react to light. The first step is usually the absorption of a photon by a prosthetic group embedded inside the photo-receptor, often a conjugated chromophore. The electronic changes in the chromophore induced by photo-absorption can trigger a cascade of structural or chemical transformations that culminate into a response to light. Understanding how these proteins have evolved to mediate their activation process has remained challenging because the required time and spacial resolutions are notoriously difficult to achieve experimentally. Therefore, mechanistic insights into photoreceptor activation have been predominantly obtained with computer simulations. Here we briefly outline the challenges associated with such computations and review the progress made in this field., (Copyright © 2020 The Authors. Published by Elsevier Ltd.. All rights reserved.)

Phytochrome proteins control the growth, reproduction, and photosynthesis of plants, fungi, and bacteria. Light is detected by a bilin cofactor, but it remains elusive how this leads to activation of the protein through structural changes. We present serial femtosecond X-ray crystallographic data of the chromophore-binding domains of a bacterial phytochrome at delay times of 1 ps and 10 ps after photoexcitation. The data reveal a twist of the D-ring, which leads to partial detachment of the chromophore from the protein. Unexpectedly, the conserved so-called pyrrole water is photodissociated from the chromophore, concomitant with movement of the A-ring and a key signaling aspartate. The changes are wired together by ultrafast backbone and water movements around the chromophore, channeling them into signal transduction towards the output domains. We suggest that the observed collective changes are important for the phytochrome photoresponse, explaining the earliest steps of how plants, fungi and bacteria sense red light., Competing Interests: EC, WW, HT, SP, LC, VK, LH, MP, MC, JK, RN, LI, AN, AC, DM, MM, MK, RB, EN, RT, TT, LF, SI, SO, KM, GG, ES, JI, MS, SW No competing interests declared, (© 2020, Claesson et al.)

When photoactive molecules interact strongly with confined light modes in optical cavities, new hybrid light-matter states form. They are known as polaritons and correspond to coherent superpositions of excitations of the molecules and of the cavity photon. The polariton energies and thus potential energy surfaces are changed with respect to the bare molecules, such that polariton formation is considered a promising paradigm for controlling photochemical reactions. To effectively manipulate photochemistry with confined light, the molecules need to remain in the polaritonic state long enough for the reaction on the modified potential energy surface to take place. To understand what determines this lifetime, we have performed atomistic molecular dynamics simulations of room-temperature ensembles of rhodamine chromophores strongly coupled to a single confined light mode with a 15 fs lifetime. We investigated three popular experimental scenarios and followed the relaxation after optically pumping (i) the lower polariton, (ii) the upper polariton, or (iii) uncoupled molecular states. The results of the simulations suggest that the lifetimes of the optically accessible lower and upper polaritons are limited by (i) ultrafast photoemission due to the low cavity lifetime and (ii) reversible population transfer into the "dark" state manifold. Dark states are superpositions of molecular excitations but with much smaller contributions from the cavity photon, decreasing their emission rates and hence increasing their lifetimes. We find that population transfer between polaritonic modes and dark states is determined by the overlap between the polaritonic and molecular absorption spectra. Importantly, excitation can also be transferred "upward" from the lower polariton into the dark-state reservoir due to the broad absorption spectra of the chromophores, contrary to the common conception of these processes as a "one-way" relaxation from the dark states down to the lower polariton. Our results thus suggest that polaritonic chemistry relying on modified dynamics taking place within the lower polariton manifold requires cavities with sufficiently long lifetimes and, at the same time, strong light-matter coupling strengths to prevent the back-transfer of excitation into the dark states.

The tetrapyrrole chromophore biliverdin IXα (BV) in the bacteriophytochrome from Deinococcus radiodurans (DrBphP) is usually assumed to be fully protonated, but this assumption has not been systematically validated by experiments or extensive computations. Here, we use force field molecular dynamics simulations and quantum mechanics/molecular mechanics calculations with density functional theory and XMCQDPT2 methods to investigate the effect of the five most probable protonation forms of BV on structural stability, binding pocket interactions, and absorption spectra in the two photochromic states of DrBphP. While agreement with X-ray structural data and measured UV/vis spectra suggest that in both states the protonated form of the chromophore dominates, we also find that a minor population with a deprotonated D-ring could contribute to the red-shifted tail in the absorption spectra.

Proton diffusion (PD) across biological membranes is a fundamental process in many biological systems, and much experimental and theoretical effort has been employed for deciphering it. Here, we report on a spectroscopic probe, which can be tightly tethered to the membrane, for following fast (nanosecond) proton transfer events on the surface of membranes. Our probe is composed of a photoacid that serves as our light-induced proton source for the initiation of the PD process. We use our probe to follow PD, and its pH dependence, on the surface of lipid vesicles composed of a zwitterionic headgroup, a negative headgroup, a headgroup that is composed only from the negative phosphate group, or a positive headgroup without the phosphate group. We reveal that the PD kinetic parameters are highly sensitive to the nature of the lipid headgroup, ranging from a fast lateral diffusion at some membranes to the escape of protons from surface to bulk (and vice versa) at others. By referring to existing theoretical models for membrane PD, we found that while some of our results confirm the quasi-equilibrium model, other results are in line with the nonequilibrium model., Competing Interests: The authors declare no conflict of interest.