Fluorescent membrane probes' behavior in lipid bilayers: insights from molecular dynamics simulations

Interaction of Hoechst 33342 with POPC Membranes at Different pH Values

Hoechst 33342 (H33342) is a fluorescent probe that is commonly used to stain the DNA of living cells. To do so, it needs to interact with and permeate through cell membranes, despite its high overall charge at physiological pH values. In this work, we address the effect of pH in the association of H33342 with lipid bilayers using a combined experimental and computational approach. The partition of H33342 to 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) lipid membranes was experimentally quantified using fluorescence spectroscopy and isothermal titration calorimetry (ITC) measurements. Quantum chemical calculations were performed to select the most stable isomer of H33342 for the overall charges 0, +1, and +2, expected to predominate across the 5 < pH < 10 range. The interaction of these isomers with POPC bilayers was then studied by both unrestrained and umbrella sampling molecular dynamics (MD) simulations. Both experimental results and computational free energy profiles indicate that the partition coefficient of H33342 displays a small variation over a wide pH range, not exceeding one order of magnitude. The enthalpy variation upon partition to the membrane suggests efficient hydrogen bonding between the probe and the lipid, namely, for the protonated +2 form, which was confirmed in the MD simulation studies. The relatively high lipophilicity obtained for the charged species contrasts with the decrease in their general hydrophobicity as estimated from octanol/water partition. This highlights the distinction between lipophilicity and hydrophobicity, as well as the importance of considering the association with lipid bilayers when predicting the affinity for biomembranes.

Fluorescent Probes cis- and trans-Parinaric Acids in Fluid and Gel Lipid Bilayers: A Molecular Dynamics Study

Fluorescence probes are indispensable tools in biochemical and biophysical membrane studies. Most of them possess extrinsic fluorophores, which often constitute a source of uncertainty and potential perturbation to the host system. In this regard, the few available intrinsically fluorescent membrane probes acquire increased importance. Among them, cis- and trans-parinaric acids (c-PnA and t-PnA, respectively) stand out as probes of membrane order and dynamics. These two compounds are long-chained fatty acids, differing solely in the configurations of two double bonds of their conjugated tetraene fluorophore. In this work, we employed all-atom and coarse-grained molecular dynamics simulations to study the behavior of c-PnA and t-PnA in lipid bilayers of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) and 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC), representative of the liquid disordered and solid ordered lipid phases, respectively. All-atom simulations indicate that the two probes show similar location and orientation in the simulated systems, with the carboxylate facing the water/lipid interface and the tail spanning the membrane leaflet. The two probes establish interactions with the solvent and lipids to a similar degree in POPC. However, the almost linear t-PnA molecules have tighter lipid packing around them, especially in DPPC, where they also interact more with positively charged lipid choline groups. Probably for these reasons, while both probes show similar partition (assessed from computed free energy profiles across bilayers) to POPC, t-PnA clearly partitions more extensively than c-PnA to the gel phase. t-PnA also displays more hindered fluorophore rotation, especially in DPPC. Our results agree very well with experimental fluorescence data from the literature and allow deeper understanding of the behavior of these two reporters of membrane organization.

Interactions between Rhodamine Dyes and Model Membrane Systems—Insights from Molecular Dynamics Simulations

Background: rhodamines are dyes widely used as fluorescent tags in cell imaging, probing of mitochondrial membrane potential, and as P-glycoprotein model substrates. In all these applications, detailed understanding of the interaction between rhodamines and biomembranes is fundamental. Methods: we combined atomistic molecular dynamics (MD) simulations and fluorescence spectroscopy to characterize the interaction between rhodamines 123 and B (Rh123 and RhB, respectively) and POPC bilayers. Results: while the xanthene moiety orients roughly parallel to the membrane plane in unrestrained MD simulations, variations on the relative position of the benzoic ring (below the xanthene for Rh123, above it for RhB) were observed, and related to the structure of the two dyes and their interactions with water and lipids. Subtle distinctions were found among different ionization forms of the probes. Experimentally, RhB displayed a lipid/water partition coefficient more than two orders of magnitude higher than Rh123, in agreement with free energy profiles obtained from umbrella sampling MD. Conclusions: this work provided detailed insights on the similarities and differences in the behavior of bilayer-inserted Rh123 and RhB, related to the structure of the probes. The much higher affinity of RhB for the membranes increases the local concentration and explains its higher apparent affinity for P-glycoprotein reconstituted in model membranes.

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.

Rational Design of Nile Red Analogs for Sensing in Membranes

Development of next-generation fluorescent probes is a key element in the quest for a greater understanding of complex biological environments (e.g., membranes) by bioimaging. Such fluorescence-based techniques rely on specialized small molecules that possess excellent fluorescent properties but also do not perturb the native biological environment in which they reside. Herein we present a theoretical/computational strategy for the design of novel optical probes for sensing in membranes based on the parent chromophore Nile Red. Using a combination of time-dependent density functional theory (TD-DFT) and molecular dynamics (MD), we have studied the optical properties and accommodation in a model membrane of Nile Red and eight analogs. Special attention has been given to the design of probes with improved solvatochromism and two-photon absorption (2PA) without altering the membrane properties. Of the eight studied analogs, two probes were found to possess attractive probe features and are hence suggested to be taken forward to chemical synthesis and experimental exploration.

Meningeal Foam Cells and Ependymal Cells in Axolotl Spinal Cord Regeneration

A previously unreported population of foam cells (foamy macrophages) accumulates in the invasive fibrotic meninges during gap regeneration of transected adult Axolotl spinal cord (salamander Ambystoma mexicanum) and may act beneficially. Multinucleated giant cells (MNGCs) also occurred in the fibrotic meninges. Actin-label localization and transmission electron microscopy showed characteristic foam cell and MNGC podosome and ruffled border-containing sealing ring structures involved in substratum attachment, with characteristic intermediate filament accumulations surrounding nuclei. These cells co-localized with regenerating cord ependymal cell (ependymoglial) outgrowth. Phase contrast-bright droplets labeled with Oil Red O, DiI, and DyRect polar lipid live cell label showed accumulated foamy macrophages to be heavily lipid-laden, while reactive ependymoglia contained smaller lipid droplets. Both cell types contained both neutral and polar lipids in lipid droplets. Foamy macrophages and ependymoglia expressed the lipid scavenger receptor CD36 (fatty acid translocase) and the co-transporter toll-like receptor-4 (TLR4). Competitive inhibitor treatment using the modified fatty acid Sulfo-N-succinimidyl Oleate verified the role of the lipid scavenger receptor CD36 in lipid uptake studies in vitro. Fluoromyelin staining showed both cell types took up myelin fragments in situ during the regeneration process. Foam cells took up DiI-Ox-LDL and DiI-myelin fragments in vitro while ependymoglia took up only DiI-myelin in vitro. Both cell types expressed the cysteine proteinase cathepsin K, with foam cells sequestering cathepsin K within the sealing ring adjacent to the culture substratum. The two cell types act as sinks for Ox-LDL and myelin fragments within the lesion site, with foamy macrophages showing more Ox-LDL uptake activity. Cathepsin K activity and cellular localization suggested that foamy macrophages digest ECM within reactive meninges, while ependymal cells act from within the spinal cord tissue during outgrowth into the lesion site, acting in complementary fashion. Small MNGCs also expressed lipid transporters and showed cathepsin K activity. Comparison of ³H-glucosamine uptake in ependymal cells and foam cells showed that only ependymal cells produce glycosaminoglycan and proteoglycan-containing ECM, while the cathepsin studies showed both cell types remove ECM. Interaction of foam cells and ependymoglia in vitro supported the dispersion of ependymal outgrowth associated with tissue reconstruction in Axolotl spinal cord regeneration.

Temperature Dependence of the Structure and Dynamics of a Dye-Labeled Lipid in a Planar Phospholipid Bilayer: A Computational Study

Fluorescent probes are widely employed to label lipids for the investigation of structural and dynamic properties of model and cell membranes through optical microscopy techniques. Although the effect of tagging a lipid with an organic dye is generally assumed to be negligible, optically modified lipids can nonetheless affect the local lipid structure and, in turn, the lipid lateral mobility. To better assess this potential issue, all-atom (MD) molecular dynamics simulations have been performed to study structural and dynamic effects in a model DOPC membrane in the presence of a standard Rhodamine B-labeled DOPE lipid (RHB) as a function of temperature, i.e., 293 K, 303 K, and 320 K. As the temperature is increased, we observe similar changes in the structural properties of both pure DOPC and RHB-DOPC lipid bilayers: an increase of the area per lipid, a reduction of the membrane thickness and a decrease of lipid order parameters. The partial density profile of the RHB headgroups and their orientation within the lipid bilayer confirm the amphiphilic nature of the RHB fluorescent moiety, which mainly partitions in the DOPC glycerol backbone region at each temperature. Moreover, at all temperatures, our results on lipid lateral diffusion support a non-neutral role of the dye with respect to the unlabeled lipid mobility, thus suggesting important implications for optical microscopy studies of lipid membranes.

Diphenylhexatriene membrane probes DPH and TMA-DPH: A comparative molecular dynamics simulation study

Fluorescence spectroscopy and microscopy have been utilized as tools in membrane biophysics for decades now. Because phospholipids are non-fluorescent, the use of extrinsic membrane probes in this context is commonplace. Among the latter, 1,6-diphenylhexatriene (DPH) and its trimethylammonium derivative (TMA-DPH) have been extensively used. It is widely believed that, owing to its additional charged group, TMA-DPH is anchored at the lipid/water interface and reports on a bilayer region that is distinct from that of the hydrophobic DPH. In this study, we employ atomistic MD simulations to characterize the behavior of DPH and TMA-DPH in 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) and POPC/cholesterol (4:1) bilayers. We show that although the dynamics of TMA-DPH in these membranes is noticeably more hindered than that of DPH, the location of the average fluorophore of TMA-DPH is only ~3-4Å more shallow than that of DPH. The hindrance observed in the translational and rotational motions of TMA-DPH compared to DPH is mainly not due to significant differences in depth, but to the favorable electrostatic interactions of the former with electronegative lipid atoms instead. By revealing detailed insights on the behavior of these two probes, our results are useful both in the interpretation of past work and in the planning of future experiments using them as membrane reporters.

Dye lipophilicity and retention in lipid membranes: implications for single-molecule spectroscopy

Fluorescence studies of individual lipid vesicles rely on the proper positioning of probes in the lipid milieu. This is true for both positional tags and chemoselective fluorogenic probes that undergo chemical modification following reaction with an analyte of interest within the lipid environment. The present report describes lipophilicity and localization estimations for a series of BODIPY dyes bearing substituents of varying hydrophobicity. We also studied fluorogenic trap-reporter probes that undergo fluorescence emission enhancement upon trapping of reactive oxygen species (ROS), including lipid peroxyl radicals. We show that caution has to be taken to extrapolate ensemble partition measurements of dyes to the single-molecule regime as a result of the dramatically different lipid concentration prevailing in ensemble versus single-molecule experiments. We show that the mole fraction of dyes that remains embedded in liposomes during a typical single-molecule experiment may be accurately determined from a ratiometric single-particle imaging analysis. We further demonstrate that fluorescence correlation spectroscopy (FCS) provides a very rapid and reliable estimate of the lipophilic nature of a given dye under highly dilute single-molecule-like conditions. Our combined single-particle spectroscopy and FCS experiments suggest that the minimal mole fraction of membrane-associated dyes (x(m)) as determined from FCS experiments is about 0.5 for adequate dye retention during single-molecule imaging in lipid membranes. Our work further highlights the dramatic effect that chemical modifications can have on chemoselective fluorogenic probe localization.

Limited perturbation of a DPPC bilayer by fluorescent lipid probes: a molecular dynamics study

The properties of lipid bilayer nanometer-scale domains could be crucial for understanding cell membranes. Fluorescent probes are often used to study bilayers, yet their effects on host lipids are not well understood. We used molecular dynamics simulations to investigate perturbations in a fluid DPPC bilayer upon incorporation of three indocarbocyanine probes: DiI-C18:0, DiI-C18:2, or DiI-C12:0. We find a 10-12% decrease in chain order for DPPC in the solvation shell nearest the probe but smaller effects in subsequent shells, indicating that the probes significantly alter only their local environment. We also observe order perturbations of lipids directly across from the probe in the opposite leaflet. Additionally, the DPPC headgroup phosphorus-to-nitrogen vector of lipids nearest the probe exhibits preferential orientation pointing away from the DiI. We show that, while DiI probes perturb their local environment, they do not strongly influence the average properties of "nanoscopic" domains containing a few hundred lipids.

Simple estimation of Förster Resonance Energy Transfer (FRET) orientation factor distribution in membranes

Because of its acute sensitivity to distance in the nanometer scale, Förster resonance energy transfer (FRET) has found a large variety of applications in many fields of chemistry, physics, and biology. One important issue regarding the correct usage of FRET is its dependence on the donor-acceptor relative orientation, expressed as the orientation factor κ². Different donor/acceptor conformations can lead to κ² values in the 0 ≤ κ² ≤ 4 range. Because the characteristic distance for FRET, R₀, is proportional to (κ²)^1/6, uncertainties in the orientation factor are reflected in the quality of information that can be retrieved from a FRET experiment. In most cases, the average value of κ² corresponding to the dynamic isotropic limit (<κ²> = 2/3) is used for computation of R₀ and hence donor-acceptor distances and acceptor concentrations. However, this can lead to significant error in unfavorable cases. This issue is more critical in membrane systems, because of their intrinsically anisotropic nature and their reduced fluidity in comparison to most common solvents. Here, a simple numerical simulation method for estimation of the probability density function of κ² for membrane-embedded donor and acceptor fluorophores in the dynamic regime is presented. In the simplest form, the proposed procedure uses as input the most probable orientations of the donor and acceptor transition dipoles, obtained by experimental (including linear dichroism) or theoretical (such as molecular dynamics simulation) techniques. Optionally, information about the widths of the donor and/or acceptor angular distributions may be incorporated. The methodology is illustrated for special limiting cases and common membrane FRET pairs.

Lateral distribution of NBD-PC fluorescent lipid analogs in membranes probed by molecular dynamics-assisted analysis of Förster Resonance Energy Transfer (FRET) and fluorescence quenching

Förster resonance energy transfer (FRET) is a powerful tool used for many problems in membrane biophysics, including characterization of the lateral distribution of lipid components and other species of interest. However, quantitative analysis of FRET data with a topological model requires adequate choices for the values of several input parameters, some of which are difficult to obtain experimentally in an independent manner. For this purpose, atomistic molecular dynamics (MD) simulations can be potentially useful as they provide direct detailed information on transverse probe localization, relative probe orientation, and membrane surface area, all of which are required for analysis of FRET data. This is illustrated here for the FRET pairs involving 1,6-diphenylhexatriene (DPH) as donor and either 1-palmitoyl,2-(6-[N-(7-nitrobenz-2-oxa-1,3-diazol-4-yl)amino] hexanoyl)- sn-glycero-3-phosphocholine (C6-NBD-PC) or 1-palmitoyl,2-(12-[N-(7-nitrobenz-2-oxa-1,3-diazol-4-yl)amino]dodecanoyl)-sn-glycero-3-phosphocholine (C12-NBD-PC) as acceptors, in fluid vesicles of 1,2-dipalmitoyl-sn-3-glycerophosphocholine (DPPC, 50 °C). Incorporation of results from MD simulations improves the statistical quality of model fitting to the experimental FRET data. Furthermore, the decay of DPH in the presence of moderate amounts of C12-NBD-PC (>0.4 mol%) is consistent with non-random lateral distribution of the latter, at variance with C6-NBD-PC, for which aggregation is ruled out up to 2.5 mol% concentration. These conclusions are supported by analysis of NBD-PC fluorescence self-quenching. Implications regarding the relative utility of these probes in membrane studies are discussed.

Recent developments in molecular dynamics simulations of fluorescent membrane probes

Due to their sensitivity and versatility, the use of fluorescence techniques in membrane biophysics is widespread. Because membrane lipids are non-fluorescent, extrinsic membrane probes are widely used. However, the behaviour of these probes when inserted in the bilayer is often poorly understood, and it can be hard to distinguish between legitimate membrane properties and perturbation resulting from probe incorporation. Atomistic molecular dynamics simulations present a convenient way to address these issues and have been increasingly used in recent years in this context. This article reviews the application of molecular dynamics to the study of fluorescent membrane probes, focusing on recent work with complex design fluorophores and ordered bilayer systems.