Investigation into Biological Environments through (Non)linear Optics: A Multiscale Study of Laurdan Derivatives

Generalized polarization and time-resolved fluorescence provide evidence for different populations of Laurdan in lipid vesicles

The solvatochromic dye Laurdan is widely used in sensing the lipid packing of both model and biological membranes. The fluorescence emission maximum shifts from about 440 nm (blue channel) in condensed membranes (So) to about 490 nm (green channel) in the liquid-crystalline phase (Lα). Although the fluorescence intensity based generalized polarization (GP) is widely used to characterize lipid membranes, the fluorescence lifetime of Laurdan, in the blue and the green channel, is less used for that purpose. Here we explore the correlation between GP and fluorescence lifetimes by spectroscopic measurements on the So and Lα phases of large unilamellar vesicles of DMPC and DPPC. A positive correlation between GP and the lifetimes is observed in each of the optical channels for the two lipid phases. Microfluorimetric determinations on giant unilamellar vesicles of DPPC and DOPC at room temperature are performed under linearly polarized two-photon excitation to disentangle possible subpopulations of Laurdan at a scale below the optical resolution. Fluorescence intensities, GP and fluorescence lifetimes depend on the angle between the orientation of the linear polarization of the excitation light and the local normal to the membrane of the optical cross-section. This angular variation depends on the lipid phase and the emission channel. GP and fluorescence intensities in the blue and green channel in So and in the blue channel in Lα exhibit a minimum near 90o. Surprisingly, the intensity in the green channel in Lα reaches a maximum near 90o. The fluorescence lifetimes in the two optical channels also reach a pronounced minimum near 90o in So and Lα, apart from the lifetime in the blue channel in Lα where the lifetime is short with minimal angular variation. To our knowledge, these experimental observations are the first to demonstrate the existence of a bent conformation of Laurdan in lipid membranes, as previously suggested by molecular dynamics calculations.

Multiscale Modeling Unravels the Influence of Biomembranes on the Photochemical Properties of Embedded Anti-Oxidative Polyphenolic and Phenanthroline Chelating Dyes

The embedding of caffeate methyl ester, the flavonoids luteolin and quercetin, and the o-phenanthroline and neocuproine in a liquid disordered lipid bilayer has been studied through extensive atomistic calculations. The location and the orientation of these bio-active antioxidants are explained and analyzed. While the two phenanthrolines strongly associate with the lipid tail region, the other three compounds are rather found among the head groups. The simulations showcase conformational changes of the flavonoids. Through the use of a hybrid quantum mechanics-molecular mechanics scheme and supported by a profound benchmarking of the electronic excited-state method for these compounds, the influence of the anisotropic environment on the compounds' optical properties is analyzed. Influences of surrounding water molecules and of the polar parts of the lipids on the transition dipole moments and excited-state dipole moments are weighted with respect to a change in conformation. The current study highlights the importance of the mapping of molecular interactions in model membranes and pinpoints properties, which can be biomedically used to discriminate and detect different lipid environments.

Fluorescent Multifunctional Organic Nanoparticles for Drug Delivery and Bioimaging: A Tutorial Review

Fluorescent organic nanoparticles (FONs) are a large family of nanostructures constituted by organic components that emit light in different spectral regions upon excitation, due to the presence of organic fluorophores. FONs are of great interest for numerous biological and medical applications, due to their high tunability in terms of composition, morphology, surface functionalization, and optical properties. Multifunctional FONs combine several functionalities in a single nanostructure (emission of light, carriers for drug-delivery, functionalization with targeting ligands, etc.), opening the possibility of using the same nanoparticle for diagnosis and therapy. The preparation, characterization, and application of these multifunctional FONs require a multidisciplinary approach. In this review, we present FONs following a tutorial approach, with the aim of providing a general overview of the different aspects of the design, preparation, and characterization of FONs. The review encompasses the most common FONs developed to date, the description of the most important features of fluorophores that determine the optical properties of FONs, an overview of the preparation methods and of the optical characterization techniques, and the description of the theoretical approaches that are currently adopted for modeling FONs. The last part of the review is devoted to a non-exhaustive selection of some recent biomedical applications of FONs.

Benchmarking Density Functional Approximations for Excited-State Properties of Fluorescent Dyes

This study presents an extensive analysis of the predictive power of time-dependent density functional theory in determining the excited-state properties of two groups of important fluorescent dyes, difluoroboranes and hydroxyphenylimidazo[1,2-a]pyridine derivatives. To ensure statistically meaningful results, the data set is comprised of 85 molecules manifesting diverse photophysical properties. The vertical excitation energies and dipole moments (in the electronic ground and excited states) of the aforementioned dyes were determined using the RI-CC2 method (reference) and with 18 density functional approximations (DFA). The set encompasses DFAs with varying amounts of exact exchange energy (EEX): from 0% (e.g., SVWN, BLYP), through a medium (e.g., TPSSh, B3LYP), up to a major contribution of EEX (e.g., BMK, MN15). It also includes range-separated hybrids (CAM-B3LYP, LC-BLYP). Similar error profiles of vertical energy were obtained for both dye groups, although the errors related to hydroxyphenylimidazopiridines are significantly larger. Overall, functionals including 40-55% of EEX (SOGGA11-X, BMK, M06-2X) ensure satisfactory agreement with the reference vertical excitation energies obtained using the RI-CC2 method; however, MN15 significantly outperforms them, providing a mean absolute error of merely 0.04 eV together with a very high correlation coefficient (R2 = 0.98). Within the investigated set of functionals, there is no single functional that would equally accurately determine ground- and excited-state dipole moments of difluoroboranes and hydroxyphenylimidazopiridine derivatives. Depending on the chosen set of dyes, the most accurate μGS predictions were delivered by MN15 incorporating a major EEX contribution (difluoroboranes) and by PBE0 containing a minor EEX fraction (hydroxyphenylimidazopiridines). Reverse trends are observed for μES, i.e., for difluoroboranes the best results were obtained with functionals including a minor fraction of EEX, specifically PBE0, while in the case of hydroxyphenylimidazopiridines, much more accurate predictions were provided by functionals incorporating a major EEX contribution (BMK, MN15).

First hyperpolarizability of water at the air-vapor interface: a QM/MM study questions standard experimental approximations

Surface Second-Harmonic Generation (S-SHG) experiments provide a unique approach to probe interfaces. One important issue for S-SHG is how to interpret the S-SHG intensities at the molecular level. Established frameworks commonly assume that each molecule emits light according to an average molecular hyperpolarizability tensor β(-2ω,ω,ω). However, for water molecules, this first hyperpolarizability is known to be extremely sensitive to their environment. We have investigated the molecular first hyperpolarizability of water molecules within the liquid-vapor interface, using a quantum description with explicit, inhomogeneous electrostatic embedding. The resulting average molecular first hyperpolarizability tensor depends on the distance relative to the interface, and it practically respects the Kleinman symmetry everywhere in the liquid. Within this numerical approach, based on the dipolar approximation, the water layer contributing to the Surface Second Harmonic Generation (S-SHG) intensity is less than a nanometer. The results reported here question standard interpretations based on a single, averaged hyperpolarizability for all molecules at the interface. Not only the molecular first hyperpolarizability tensor significantly depends on the distance relative to the interface, but it is also correlated to the molecular orientation. Such hyperpolarizability fluctuations may impact the S-SHG intensity emitted by an aqueous interface.

Unravelling the Effects of Cholesterol on the Second-Order Nonlinear Optical Responses of Di-8-ANEPPS Dye Embedded in Phosphatidylcholine Lipid Bilayers

Cholesterol is known for its role in maintaining the correct fluidity and rigidity of the animals cell membranes and thus their functions. Assessing the content and the role of cholesterol in lipid bilayers is therefore of crucial importance for a deeper understanding and control of membrane functioning. In this computational work, we investigate bilayers built from three types of glycerophospholipid phosphatidylcholine (PC) lipids, namely dipalmitoylphosphatidylcholine (DPPC), 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC), and dioleoylphosphatidylcholine (DOPC), and containing different amounts of cholesterol by analyzing the second-harmonic generation (SHG) nonlinear optical (NLO) response of a probe molecule, di-8-ANEPPS, inserted into the membranes. This molecular property presents the advantage to be specific to interfacial regions such as lipid bilayers. To unravel these effects, Molecular Dynamics (MD) simulations have been performed on both DPPC and DOPC lipids by varying the cholesterol mole fraction (from 0 to 0.66), while POPC was only considered as a pure bilayer. In the case of the structural properties of the bilayers, all the analyses converge toward the same conclusion: as the mole fraction of cholesterol increases, the systems become more rigid, confirming the condensing effect of cholesterol. In addition, the chromophore is progressively more aligned with respect to the normal to the bilayer. On the contrary, addition of unsaturation disorders the lipid bilayers, with barely no impact on the alignment of the chromophore. Then, using the frames obtained from the MD simulations, the first hyperpolarizability β of the dye in its environment has been computed at the TDDFT level. On the one hand, the addition of cholesterol induces a progressive increase of the diagonal component the β tensor parallel to the bilayer normal. On the other hand, larger β values have been calculated for the unsaturated than for the saturated lipid systems. In summary, this study illustrates the relationship between the composition and structure of the bilayers and the NLO responses of the embedded dye.

Computational and photophysical characterization of a Laurdan malononitrile derivative

The malononitrile group is considered one of the strongest natural electron-withdrawing groups in a chemist's arsenal. However, surprisingly little is known about how this group functions in push-pull fluorophores. In a recent computational study, we reported that replacing the ketone group of the traditional push-pull dye Laurdan with a malononitrile group significantly improves the optical properties while retaining the membrane behavior of the parent molecule Laurdan. Motivated by these results, we report here the synthesis and photophysical characterization of the said compound, 6-(1-undecyl-2,2-dicyanovinyl)-N,N-dimethyl-2-naphthylamine (CN-Laurdan). To our surprise, this new CN-Laurdan probe is found to be much less bright than the parent Laurdan due to a large drop in the fluorescence quantum yield. Using computational methods, we determine that the origin of this low quantum yield is related to the existence of a non-radiative decay pathway related to a rotation of the malononitrile moiety, suggesting that the molecule could nonetheless function very well as a molecular rotor. We confirm experimentally that CN-Laurdan functions as a molecular rotor by measuring the quantum yield in methanol/glycerol mixtures of increasing viscosity. Specifically, we found a consistent increase in the quantum yield across the entire range of tested viscosities.

Computational Characterization of Novel Malononitrile Variants of Laurdan with Improved Photophysical Properties for Sensing in Membranes

Fluorescent probes are powerful tools for improving our understanding of cellular membranes and other complex biological environments. Using simulations, we gain atomistic and electronic insights into the effectiveness of the probes. In the current work, we have used various computational approaches to comprehensively investigate the properties of the fluorescent probe Laurdan and two Laurdan-like probes: AADAL and ECL. In addition, we propose the development of their corresponding novel malononitrile variants, which are computationally characterized herein. For the candidate probes, electronic structure calculations were used to rationalize their optical properties, including their ability for two-photon activation, and molecular dynamics simulations were used to unravel atomistic details of their functioning within lipid bilayers. While Laurdan, AADAL, and ECL were found to have very similar optical and membrane partitioning profiles, their malononitrile variants were found to show significantly improved optical properties, especially in regard to two-photon cross sections, and they appear to retain the desired membrane characteristics of the parent Laurdan molecule.

The influence of lipid membranes on fluorescent probes' optical properties

BACKGROUND

Organic fluorophores embedded in lipid bilayers can nowadays be described by a multiscale computational approach. Combining different length and time scales, a full characterization of the probe localization and optical properties led to novel insight into the effect of the environments.

SCOPE OF REVIEW

Following an introduction on computational advancements, three relevant probes are reviewed that delineate how a multiscale approach can lead to novel insight into the probes' (non) linear optical properties. Attention is paid to the quality of the theoretical description of the optical techniques.

MAJOR CONCLUSIONS

Computation can assess a priori novel probes' optical properties and guide the analysis and interpretation of experimental data in novel studies. The properties can be used to gain information on the phase and condition of the surrounding biological environment.

GENERAL SIGNIFICANCE

Computation showed that a canonical view on some of the probes should be revisited and adapted.

Influence of Membrane Phase on the Optical Properties of DPH

The fluorescent molecule diphenylhexatriene (DPH) has been often used in combination with fluorescence anisotropy measurements, yet little is known regarding the non-linear optical properties. In the current work, we focus on them and extend the application to fluorescence, while paying attention to the conformational versatility of DPH when it is embedded in different membrane phases. Extensive hybrid quantum mechanics/molecular mechanics calculations were performed to investigate the influence of the phase- and temperature-dependent lipid environment on the probe. Already, the transition dipole moments and one-photon absorption spectra obtained in the liquid ordered mixture of sphingomyelin (SM)-cholesterol (Chol) (2:1) differ largely from the ones calculated in the liquid disordered DOPC and solid gel DPPC membranes. Throughout the work, the molecular conformation in SM:Chol is found to differ from the other environments. The two-photon absorption spectra and the ones obtained by hyper-Rayleigh scattering depend strongly on the environment. Finally, a stringent comparison of the fluorescence anisotropy decay and the fluorescence lifetime confirm the use of DPH to gain information upon the surrounding lipids and lipid phases. DPH might thus open the possibility to detect and analyze different biological environments based on its absorption and emission properties.

The Secret Lives of Fluorescent Membrane Probes as Revealed by Molecular Dynamics Simulations

Fluorescent probes have been employed for more than half a century to study the structure and dynamics of model and biological membranes, using spectroscopic and/or microscopic experimental approaches. While their utilization has led to tremendous progress in our knowledge of membrane biophysics and physiology, in some respects the behavior of bilayer-inserted membrane probes has long remained inscrutable. The location, orientation and interaction of fluorophores with lipid and/or water molecules are often not well known, and they are crucial for understanding what the probe is actually reporting. Moreover, because the probe is an extraneous inclusion, it may perturb the properties of the host membrane system, altering the very properties it is supposed to measure. For these reasons, the need for independent methodologies to assess the behavior of bilayer-inserted fluorescence probes has been recognized for a long time. Because of recent improvements in computational tools, molecular dynamics (MD) simulations have become a popular means of obtaining this important information. The present review addresses MD studies of all major classes of fluorescent membrane probes, focusing in the period between 2011 and 2020, during which such work has undergone a dramatic surge in both the number of studies and the variety of probes and properties accessed.

Cyanine dyes with tail length asymmetry enhance photoselection: A multiscale study on DiD probes in a liquid disordered membrane

Visualization of membrane domains like lipid rafts in natural or artificial membranes is a crucial task for cell biology. For this purpose, fluorescence microscopy is often used. Since fluorescing probes in lipid membranes partition specifically in e.g. local liquid disordered or liquid ordered environments, the consequent changes in their orientation and location are both theoretically and experimentally of interest. Here we focused on a liquid disordered membrane phase and performed molecular dynamics (MD) simulations of the indocarbocyanine DiD probes by varying the length of the attached alkyl tails and also the length of the cyanine backbone. From the probed compounds in a DOPC lipid bilayer at ambient temperature, a varying orientation of the transition dipole moment was observed, which is crucial for fluorescence microscopy and which, through photoselection, was found to be surprisingly more effective for asymmetric probes than for the symmetric ones. Furthermore, we observed that the orientation of the probes was dependent on the tail length; with the methyls or propyls attached, DiD oriented with its tails facing the water, contrary to the ones with longer tails. With advanced hybrid QM/MM calculations we show that the different local environment for differently oriented probes affected the one-photon absorption spectra, that was blue-shifted for the short-tailed DiD with respect to the DiDs with longer tails. We show here that the presented probes can be successfully used for fluorescence microscopy and we believe that the described properties bring further insight for the experimental use of these probes.

Conformational Changes as Driving Force for Phase Recognition: The Case of Laurdan

The development of a universal probe to assess the phase of a lipid membrane is one of the most ambitious goals for fluorescence spectroscopy. The ability of a well-known molecule as Laurdan to reach this aim is here exploited as the behavior of the probe is fully characterized in a dipalmitoylphosphatidylcholine (DPPC) solid gel (So) phase by means of molecular dynamics simulations. Laurdan can take two conformations, depending on whether the carbonyl oxygen points toward the β-position of the naphthalene core (Conf-I) or to the α-position (Conf-II). We observe that Conf-I has an elongated form in this environment, whereas Conf-II takes an L-shape. Interestingly, our theoretical calculations show that these two conformations behave in an opposite way from what is reported in the literature for a DPPC membrane in a liquid disordered (Ld) phase, where Conf-I assumes an L-shape and Conf-II is elongated. Moreover, our results show that in DPPC (So) no intermixing between the conformations is present, whereas it has been seen in a fluid environment such as DOPC (Ld). Through a careful analysis of angle distributions and by means of the rotational autocorrelation function, we predict that the two conformers of Laurdan behave differently in different membrane environments.

Atomistic Picture of Fluorescent Probes with Hydrocarbon Tails in Lipid Bilayer Membranes: An Investigation of Selective Affinities and Fluorescent Anisotropies in Different Environmental Phases

By reverting to spectroscopy, changes in the biological environment of a fluorescent probe can be monitored and the presence of various phases of the surrounding lipid bilayer membranes can be detected. However, it is currently not always clear in which phase the probe resides. The well-known orange 1,1'-dioctadecyl-3,3,3',3'-tetramethylindodicarbo-cyanine perchlorate (DiI-C18(5)) fluorophore, for instance, and the new, blue BODIPY (4,4-difluoro-4-bora-3 a,4 a-diaza- s-indacene) derivative were experimentally seen to target and highlight identical parts of giant unilamellar vesicles of various compositions, comprising mixtures of dipalmitoylphosphatidylcholine (DPPC), dioleoylphosphatidylcholine (DOPC), sphingomyelin (SM), and cholesterol (Chol). However, it was not clear which of the coexisting membrane phases were visualized (Bacalum et al., Langmuir. 2016, 32, 3495). The present study addresses this issue by utilizing large-scale molecular dynamics simulations and the z-constraint method, which allows evaluating Gibbs free-energy profiles. The current calculations give an indication why, at room temperature, both BODIPY and DiI-C18(5) probes prefer the gel (S_o) phase in DOPC/DPPC (2:3 molar ratio) and the liquid-ordered (L_o) phase in DOPC/SM/Chol (1:2:1 molar ratio) mixtures. This study highlights the important differences in orientation and location and therefore in efficiency between the probes when they are used in fluorescence microscopy to screen various lipid bilayer membrane phases. Dependent on the lipid composition, the angle between the transition-state dipole moments of both probes and the normal to the membrane is found to deviate clearly from 90°. It is seen that the DiI-C18(5) probe is located in the headgroup region of the SM/Chol mixture, in close contact with water molecules. A fluorescence anisotropy study also indicates that DiI-C18(5) gives rise to a distinctive behavior in the SM/Chol membrane compared to the other considered membranes. The latter behavior has not been seen for the studied BODIPY probe, which is located deeper in the membrane.

Orientation of Laurdan in Phospholipid Bilayers Influences Its Fluorescence: Quantum Mechanics and Classical Molecular Dynamics Study

Fluidity of lipid membranes is known to play an important role in the functioning of living organisms. The fluorescent probe Laurdan embedded in a lipid membrane is typically used to assess the fluidity state of lipid bilayers by utilizing the sensitivity of Laurdan emission to the properties of its lipid environment. In particular, Laurdan fluorescence is sensitive to gel vs liquid–crystalline phases of lipids, which is demonstrated in different emission of the dye in these two phases. Still, the exact mechanism of the environment effects on Laurdan emission is not understood. Herein, we utilize dipalmitoylphosphatidylcholine (DPPC) and dioleoylphosphatidylcholine (DOPC) lipid bilayers, which at room temperature represent gel and liquid–crystalline phases, respectively. We simulate absorption and emission spectra of Laurdan in both DOPC and DPPC bilayers with quantum chemical and classical molecular dynamics methods. We demonstrate that Laurdan is incorporated in heterogeneous fashion in both DOPC and DPPC bilayers, and that its fluorescence depends on the details of this embedding.

Dynamics of the phospholipid shell of microbubbles: a fluorescence photoselection and spectral phasor approach

The lipid organization of microbubbles is important in many applications.