Photophysics in Biomembranes: Computational Insight into the Interaction between Lipid Bilayers and Chromophores

Light is ubiquitously available to probe the structure and dynamics of biomolecules and biological tissues. Generally, this cannot be done directly with visible light, because of the absence of absorption by those biomolecules. This problem can be overcome by incorporating organic molecules (chromophores) that show an optical response in the vicinity of those biomolecules. Since those optical properties are strongly dependent on the chromophore's environment, time-resolved spectroscopic studies can provide a wealth of information on biosystems at the molecular scale in a nondestructive way. In this work, we give an overview on the multiscale computational strategy developed by us in the last eight years and prove that theoretical studies and simulations are needed to explain, guide, and predict observations in fluorescence experiments. As we challenge the accepted views on existing probes, we discover unexplored abilities that can discriminate surrounding lipid bilayers and their temperature-dependent as well as solvent-dependent properties. We focus on three archetypal chromophores: diphenylhexatriene (DPH), Laurdan, and azobenzene. Our method shows that conformational changes should not be neglected for the prototype rod-shaped molecule DPH. They determine its position and orientation in a liquid-ordered (Lo) sphingomyelin/cholesterol (SM/Chol) bilayer and are responsible for a strong differentiation of its absorption spectra and fluorescence decay times in dioleoylphosphatidylcholine (DOPC) and dipalmitoylphosphatidylcholine (DPPC) membranes, which are at room temperature in liquid-disordered (Ld) and solid-gel (So) phases, respectively. Thanks to its pronounced first excited state dipole moment, Laurdan has long been known as a solvatochromic probe. Since this molecule has however two conformers, we prove that they exhibit different properties in different lipid membrane phases. We see that the two conformers are only blocked in one phase but not in another. Supported by fluorescence anisotropy decay simulations, Laurdan can therefore be regarded as a molecular rotor. Finally, the conformational versatility of azobenzene in saturated Ld lipid bilayers is simulated, along with its photoisomerization pathways. By means of nonadiabatic QM/MM surface hopping analyses (QM/MM-SH), a dual mechanism is found with a torsional mechanism and a slow conversion for trans-to-cis. For cis-to-trans, simulations show a much higher quantum yield and a so-called "pedal-like" mechanism. The differences are related to the different potential energy surfaces as well as the interactions with the surrounding alkyl chains. When tails of increased length are attached to this probe, cis is pushed toward the polar surface, while trans is pulled toward the center of the membrane.