Partitioning and localization of environment-sensitive 2-(2'-pyridyl)- and 2-(2'-pyrimidyl)-indoles in lipid membranes: a joint refinement using fluorescence measurements and molecular dynamics simulations

Fluorescent Probes and Quenchers in Studies of Protein Folding and Protein-Lipid Interactions

Fluorescence spectroscopy provides numerous methodological tools for structural and functional studies of biological macromolecules and their complexes. All fluorescence-based approaches require either existence of an intrinsic probe or an introduction of an extrinsic one. Moreover, studies of complex systems often require an additional introduction of a specific quencher molecule acting in combination with a fluorophore to provide structural or thermodynamic information. Here, we review the fundamentals and summarize the latest progress in applications of different classes of fluorescent probes and their specific quenchers, aimed at studies of protein folding and protein-membrane interactions. Specifically, we discuss various environment-sensitive dyes, FRET probes, probes for short-distance measurements, and several probe-quencher pairs for studies of membrane penetration of proteins and peptides. The goals of this review are: (a) to familiarize the readership with the general concept that complex biological systems often require both a probe and a quencher to decipher mechanistic details of functioning and (b) to provide example of the immediate applications of the described methods.

Temperature-sensitive emission of dialkylaminostyrylhetarene dyes and their incorporation into phospholipid aggregates: Applicability for thermal sensing and cellular uptake behavior

A series of dialkylaminostyrylhetarene dyes constructed from electron-rich and electron-deficient moieties of various structures connected via vinylene π-bridges are introduced as temperature-sensitive luminophores. The temperature dependent emission of the dyes in the acidified dichloromethane solutions derives from temperature-induced shift of the equilibrium between neutral and protonated forms of the dyes. The heating-induced blue shift and intensification of emission of neutral form of the dyes make them a promising basis for development of nanoparticles exhibiting temperature-sensitivity in aqueous solutions at pH typical of biological liquids. Hydrophobicity-driven incorporation of the water insoluble dyes into L-α-phosphatidylcholine(PC)-based bilayers allows to obtain water dispersible dye-PC aggregates, and to follow their emission in the aqueous solutions. Structure of the dyes has strong impact on the efficacy of the dyes incorporation into the PC-based bilayers, temperature sensitivity of emission of the dye-PC aggregates and its reversibility under the heating/cooling cycles. This enables structural optimization of the dyes in order to obtain the dye-PC species demonstrating maximal temperature dependence and reversibility of their luminescence in aqueous solutions. The selected leader exhibits low cytotoxicity exemplified for M-HeLa and Chang Liver cell lines, while the efficient cell internalization of the dye, manifested in the staining of the cell cytoplasm, opens further opportunities for biosensing applications.

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.

Atomistic molecular dynamics simulations of the LCST conformational transition in poly(N-vinylcaprolactam) in water

Thermoresponsive poly(N-vinylcaprolactam) (PVCL) has received growing interest due to a temperature-induced phase transition, which switches its solubility in aqueous solutions. However, the lower critical solution temperature (LCST) of PVCL is greatly influenced by the molecular weight, morphology and the environment. Therefore, despite of numerous experimental studies of the thermal response of PVCL, a driving force and a molecular origin of conformation transitions in solution remain far less studied. To get a better understanding of the coil-to-globule conformation transition of PVCL in aqueous solution, we examined the structure and conformation dynamics of a single-chain PVCL₃₀ in a temperature range of 280-360 K by using atomistic molecular dynamics (MD) simulations. The united-atom GROMOS G53a6 force field was re-parameterized and fine-tuned by DFT calculations to reproduce the experimental LCST transition of PVCL. Our MD model reproduces the LCST transition of PVCL₃₀ to occur within a temperature range of 34.6-38.5°. MD simulation results suggest a significant difference between the hydration state of the carbonyl group of PVCL below and above the LCST threshold. The analysis of the number of hydrogen bonds of PVCL with water molecules demonstrates that dehydration of the polymer plays an important role and drives the temperature-induced polymer collapse. Finally, the developed MD model and FF parameters were successfully tested for large-scale systems, such as mixture PVCL₃₀ oligomer and single-chain PVCL₈₁₆ polymer, respectively.

Association and sequestered dissociation of an anticancer drug from liposome membrane: Role of hydrophobic hydration

Herein, the interaction of a potent anticancer drug (Sanguinarine, SG) with dimyristoyl-l-α-phosphatidylglycerol (DMPG) liposome membrane has been investigated at physiological pH. The spectroscopic fluorescence decay results demonstrate a modification of the photophysics of SG within DMPG-encapsulated state leading to preferential stabilization of the iminium ion over the alkanolamine form. This suggests a key role of electrostatic force underlying the interaction. The complex dependence of the thermodynamic parameters on temperature yields a unique finding of a positive heat capacity change (ΔC_p) indicating the signature of hydrophobic hydration. The study also demonstrates the application of β-cyclodextrin (βCD) as a prospective host system resulting in release of the DMPG-bound drug. A calorimetric exploration of the DMPG-βCD interaction reveals an intrinsically complex thermodynamics of the process leading to ΔC_p > 0 and thus marking the instrumental role of hydrophobic hydration which follows that the DMPG-βCD interaction is accompanied with burial of polar molecular surfaces. A systematic investigation of the diffusion of the drug within various microheterogeneous environments by Fluorescence Correlation Spectroscopy (FCS) categorically reinforces our arguments.

Poly(vinyl alcohol) as a water protecting agent for silver nanoparticles: the role of polymer size and structure

Chemical modification of silver nanoparticles (AgNPs) with a stabilizing agent, such as poly(vinyl alcohol) (PVA), plays an important role in shape-controlled seeded-growth and colloidal stability. However, theoretical aspects of the stabilizing mechanism of PVA are still poorly understood. To gain a better understanding of the role of PVA in water protecting effects for silver nanoparticles, we developed an atomistic model of a AgNP grafted with single-chain PVA of various lengths. Our model, designed for classical molecular dynamics (MD) simulations, approximates the AgNP as a quasi-spherical silver nanocrystal with 3.9 nm diameter and uses a united-atom representation for PVA with its polymer chain length varying from 220 up to 1540 repeating units. We found that PVA adsorbs onto the AgNP surface through multiple non-covalent interactions, among which non-covalent bonding of the hydroxyl groups plays a key role. The analysis of adsorption isotherms by using the Hill, Scatchard, and McGhee & von Hippel models exhibits evidence for positive binding cooperativity with the cooperativity parameter varying from 1.55 to 2.12. Our results indicate that the size of the PVA polymer rather than its structure plays a crucial role in providing water protecting effects for the AgNP core, varying from 40% up to 91%. The water-protecting efficiency was well approximated by the Langmuir-Freundlich equation, allowing us to predict that the saturated coverage of the nanoparticle of a given diameter of 3.9 nm should occur when the PVA molecular weight approaches 115 kDa, which corresponds to the number of vinyl alcohol monomers being equal to 3100 units.

Computational refinement of spectroscopic FRET measurements

This article supplies raw data related to a research article entitled “Joint refinement of FRET measurements using spectroscopic and computational tools” (Kyrychenko et al., 2017) [1], in which we demonstrate the use of molecular dynamics simulations to estimate FRET orientational factors in a benchmark donor-linker-acceptor system of enhanced cyan (ECFP) and enhanced yellow (EYFP) fluorescent proteins. This can improve the recalculation of donor-acceptor distance information from single-molecule FRET measurements.

Non-typical fluorescence studies of excited and ground state proton and hydrogen transfer

Fluorescence studies of tautomerization have been carried out for various systems that exhibit single and double proton or hydrogen translocation in various environments, such as liquid and solid condensed phases, ultracold supersonic jets, and finally, polymer matrices with single emitters. We focus on less explored areas of application of fluorescence for tautomerization studies, using porphycene, a porphyrin isomer, as an example. Fluorescence anisotropy techniques allow investigations of self-exchange reactions, where the reactant and product are formally identical. Excitation with polarized light makes it possible to monitor tautomerization in single molecules and to detect their three-dimensional orientation. Analysis of fluorescence from single vibronic levels of jet-isolated porphycene not only demonstrates coherent tunneling of two internal protons, but also indicates that the process is vibrational mode-specific. Next, we present bifunctional proton donor-acceptor systems, molecules that are able, depending on the environment, to undergo excited state single intramolecular or double intermolecular proton transfer. For molecules that have donor and acceptor groups located in separate moieties linked by a single bond, excited state tautomerization can be coupled to mutual twisting of the two subunits.

Molecular modeling of lipid probes and their influence on the membrane

In this review a number of Molecular Dynamics simulation studies are discussed which focus on the understanding of the behavior of lipid probes in biomembranes. Experiments often use specialized probe moieties or molecules to report on the behavior of a membrane and try to gain information on the membrane as a whole from the probe lipids as these probes are the only things an experiment sees. Probes can be used to make NMR, EPR and fluorescence accessible to the membrane and use fluorescent or spin-active moieties for this purpose. Clearly membranes with and without probes are not identical which makes it worthwhile to elucidate the differences between them with detailed atomistic simulations. In almost all cases these differences are confined to the local neighborhood of the probe molecules which are sparsely used and generally present as single molecules. In general, the behavior of the bulk membrane lipids can be qualitatively understood from the probes but in most cases their properties cannot be directly quantitatively deduced from the probe behavior. This article is part of a Special Issue entitled: Biosimulations edited by Ilpo Vattulainen and Tomasz Róg.

Octadecyl ferulate behavior in 1,2-Dioleoylphosphocholine liposomes

Octadecyl ferulate was prepared using solid acid catalyst, monitored using Supercritical Fluid Chromatography and purified to a 42% yield. Differential scanning calorimetry measurements determined octadecyl ferulate to have melting/solidification phase transitions at 67 and 39°C, respectively. AFM imaging shows that 5-mol% present in a lipid bilayer induced domains to form. Phase behavior measurements confirmed that octadecyl ferulate increased transition temperature of phospholipids. Fluorescence measurements demonstrated that octadecyl ferulate stabilized liposomes against leakage, maintained antioxidant capacity within liposomes, and oriented such that the feruloyl moiety remained in the hydrophilic region of the bilayer. Molecular modeling calculation indicated that antioxidant activity was mostly influenced by interactions within the bilayer.

Hydroxytyrosol and tyrosol esters partitioning into, location within, and effect on DOPC liposome bilayer behavior

The phenols hydroxytyrosol and tyrosol made abundantly available through olive oil processing were enzymatically transesterified into effective lipophilic antioxidants with cuphea oil. The hydroxytyrosyl and tyrosyl esters made from cuphea oil were assessed for their ability to partition into, locate within and effect the bilayer behavior of 1,2-dioloeoylphosphatidylcholine liposomes and compared to their counterparts made from decanoic acid. Partitioning into liposomes was on the same scale for both hydroxytyrosyl derivatives and both tyrosyl derivatives. All were found to locate nearly at the same depth within the bilayer. Each was found to affect bilayer behavior in a distinct manner.

Calibration of Distribution Analysis of the Depth of Membrane Penetration Using Simulations and Depth-Dependent Fluorescence Quenching

Determination of the depth of membrane penetration provides important information for studies of membrane protein folding and protein-lipid interactions. Here, we use a combination of molecular dynamics (MD) simulations and depth-dependent fluorescence quenching to calibrate the methodology for extracting quantitative information on membrane penetration. In order to investigate the immersion depth of the fluorescent label in lipid bilayer, we studied 7-nitrobenz-2-oxa-1,3-diazole (NBD) attached to the lipid headgroup in NBD-PE incorporated into POPC bilayer. The immersion depth of NBD was estimated by measuring steady-state and time-resolved fluorescence quenching with spin-labeled lipids co-incorporated into lipid vesicles. Six different spin-labeled lipids were utilized: one with headgroup-attached Tempo probe (Tempo-PC) and five with acyl chain-labeled n-Doxyl moieties (n-Doxyl-PC where n is a chain labeling position equal to 5, 7, 10, 12, and 14, respectively). The Stern-Volmer analysis revealed that NBD quenching in membranes occurs by both static and dynamic collisional quenching processes. Using the methodology of Distribution Analysis, the immersion depth and the apparent half-width of the transversal distributions of the NBD moiety were estimated to be 14.7 and 6.7 Å, respectively, from the bilayer center. This position is independently validated by atomistic MD simulations of NBD-PE lipids in a POPC bilayer (14.4 Å). In addition, we demonstrate that MD simulations of the transverse overlap integrals between dye and quencher distributions can be used for proper analysis of the depth-dependent quenching profile. Finally, we illustrate the application of this methodology by determining membrane penetration of site selectively labeled mutants of diphtheria toxin T-domain.

Molecular dynamics simulations of depth distribution of spin-labeled phospholipids within lipid bilayer

Spin-labeled lipids are commonly used as fluorescence quenchers in studies of membrane penetration of dye-labeled proteins and peptides using depth-dependent quenching. Accurate calculations of depth of the fluorophore rely on the use of several spin labels placed in the membrane at various positions. The depth of the quenchers (spin probes) has to be determined independently; however, experimental determination of transverse distributions of spin probe depths is difficult. In this Article, we use molecular dynamics (MD) simulations to study the membrane behavior and depth distributions of spin-labeled phospholipids in a 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) bilayer. To probe different depths within the bilayer, a series containing five Doxyl-labeled lipids (n-Doxyl PC) has been studied, in which a spin moiety was covalently attached to nth carbon atoms (where n = 5, 7, 10, 12, and 14) of the sn-2 stearoyl chain of the host phospholipid. Our results demonstrate that the chain-attached spin labels are broadly distributed across the model membrane and their environment is characterized by a high degree of mobility and structural heterogeneity. Despite the high thermal disorder, the depth distributions of the Doxyl labels were found to correlate well with their attachment positions, indicating that the distribution of the spin label within the model membrane is dictated by the depth of the nth lipid carbon atom and not by intrinsic properties of the label. In contrast, a much broader and heterogeneous distribution was observed for a headgroup-attached Tempo spin label of Tempo-PC lipids. MD simulations reveal that, due to the hydrophobic nature, a Tempo moiety favors partitioning from the headgroup region deeper into the membrane. Depending on the concentration of Tempo-PC lipids, the probable depth of the Tempo moiety could span a range from 14.4 to 18.2 Å from the membrane center. Comparison of the MD-estimated immersion depths of Tempo and n-Doxyl labels with their suggested experimental depth positions allows us to review critically the possible sources of error in depth-dependent fluorescence quenching studies.

Validation of depth-dependent fluorescence quenching in membranes by molecular dynamics simulation of tryptophan octyl ester in POPC bilayer

Depth-dependent fluorescence quenching is an important tool for studying the penetration of proteins and peptides into lipid bilayers. Extracting quantitative information from quenching data is, however, complicated by (1) a limited number of experimentally available quenchers and (2) thermal disorder resulting in broad distributions of the transverse positions of both quenchers and fluorophores. Here we validate and refine a general approach to determining the location of a fluorescent probe along the bilayer normal from quenching data, based on a molecular dynamics (MD) simulation of a model compound, tryptophan octyl ester (TOE), in a 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) bilayer. The TOE ring was found to lie deeply within the bilayer (most probable position of 13.3 Å and center-of-weight of the distribution of 14.8 Å from the bilayer center), and it was very broadly distributed (with 9 Å depth distribution width), which is consistent with previous experimental observations. The depth-dependent quenching profiles were simulated by treating carbon atoms of the lipid acyl chain of POPC as "pseudo-quenchers" and calculating appropriate transverse overlaps and collision rates with indole atoms of TOE. These simulated quenching profiles were well fitted by a Gaussian function of depth, as is routinely done with experimental data subjected to the distribution analysis procedure [Methods Enzymol. 1997, 278, 462-473]. Comparison of the collisional pseudoquenching profiles with the actual profiles of the indole moiety of TOE allows for testing of the validity of the data analysis and identification of the possible sources of error in calculating depths of membrane penetration from quenching data.

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.

Partitioning of 2,6-Bis(1H-Benzimidazol-2-yl)pyridine fluorophore into a phospholipid bilayer: complementary use of fluorescence quenching studies and molecular dynamics simulations

Successful use of fluorescence sensing in elucidating the biophysical properties of lipid membranes requires knowledge of the distribution and location of an emitting molecule in the bilayer. We report here that 2,6-bis(1H-benzimidazol-2-yl)pyridine (BBP), which is almost non-fluorescent in aqueous solutions, reveals a strong emission enhancement in a hydrophobic environment of a phospholipid bilayer, making it interesting for fluorescence probing of water content in a lipid membrane. Comparing the fluorescence behavior of BBP in a wide variety of solvents with those in phospholipid vesicles, we suggest that the hydrogen bonding interactions between a BBP fluorophore and water molecules play a crucial role in the observed "light switch effect". Therefore, the loss of water-induced fluorescence quenching inside a membrane are thought to be due to deep penetration of BBP into the hydrophobic, water-free region of a bilayer. Characterized by strong quenching by transition metal ions in solution, BBP also demonstrated significant shielding from the action of the quencher in the presence of phospholipid vesicles. We used the increase in fluorescence intensity, measured upon titration of probe molecules with lipid vesicles, to estimate the partition constant and the Gibbs free energy (ΔG) of transfer of BBP from aqueous buffer into a membrane. Partitioning BBP revealed strongly favorable ΔG, which depends only slightly on the lipid composition of a bilayer, varying in a range from -6.5 to -7.0kcal/mol. To elucidate the binding interactions of the probe with a membrane on the molecular level, a distribution and favorable location of BBP in a POPC bilayer were modeled via atomistic molecular dynamics (MD) simulations using two different approaches: (i) free, diffusion-driven partitioning of the probe molecules into a bilayer and (ii) constrained umbrella sampling of a penetration profile of the dye molecule across a bilayer. Both of these MD approaches agreed with regard to the preferred location of a BBP fluorophore within the interfacial region of a bilayer, located between the hydrocarbon acyl tails and the initial portion of the lipid headgroups. MD simulations also revealed restricted permeability of water molecules into this region of a POPC bilayer, determining the strong fluorescence enhancement observed experimentally for the membrane-partitioned form of BBP.