Influence of DPH on the structure and dynamics of a DPPC bilayer

We have conducted extensive molecular dynamics (MD) simulations together with differential scanning calorimetry (DSC) and nuclear magnetic resonance (NMR) experiments to quantify the influence of free 1,6-diphenyl-1,3,5-hexatriene (DPH) fluorescent probes on the structure and dynamics of a dipalmitoylphosphatidylcholine bilayer. Atomistic MD simulations show that in the membrane-water interface the influence of DPH is minor, whereas in the acyl-chain region DPH gives rise to major perturbations. In the latter case, DPH is found to influence a wide range of membrane properties, such as the packing and ordering of hydrocarbon tails and the lateral diffusion of lipid molecules. The effects are prominent but of local nature, i.e., the changes observed in the properties of lipid molecules are significant in the vicinity of DPH, but reduce rapidly as the distance from the probe increases. Long-range perturbations due to DPH are hence not expected. Detailed DSC and (2)H NMR measurements support this view. DSC shows only subtle perturbation to the cooperative behavior of the membrane system in the presence of DPH, and (2)H NMR shows that DPH gives rise to a slight increase in the lipid chain order, in agreement with MD simulations. Potential effects of other probes such as pyrene are briefly discussed.     

Extracting Curvature Preferences of Lipids Assembled in Flat Bilayers Shows Possible Kinetic Windows for Genesis of Bilayer Asymmetry and Domain Formation in Biological Membranes

Studies on the assembly of pure lipid components allow mechanistic insights toward understanding the structural and functional aspects of biological membranes. Molecular dynamic (MD) simulations on membrane systems provide molecular details on membrane dynamics that are difficult to obtain experimentally. A large number of MD studies have remained somewhat disconnected from a key concept of amphipathic assembly resulting in membrane structures—shape parameters of lipid molecules in those structures in aqueous environments. This is because most of the MD studies have been done on flat lipid membranes. With the above in view, we analyzed MD simulations of 26 pure lipid patches as a function of (1) lipid type(s) and (2) time of MD simulations along with 35–40 ns trajectories of five pure lipids. We report, for the first time, extraction of curvature preferences of lipids from MD simulations done on flat bilayers. Our results may lead to mechanistic insights into the possible origins of bilayer asymmetries and domain formation in biological membranes.     

Setting up and running molecular dynamics simulations of membrane proteins

Molecular dynamics simulations have become a popular and powerful technique to study lipids and membrane proteins. We present some general questions and issues that should be considered prior to embarking on molecular dynamics simulation studies of membrane proteins and review common simulation methods. We suggest a practical approach to setting up and running simulations of membrane proteins, and introduce two new (related) methods to embed a protein in a lipid bilayer. Both methods rely on placing lipids and the protein(s) on a widely spaced grid and then 'shrinking' the grid until the bilayer with the protein has the desired density, with lipids neatly packed around the protein. When starting from a grid based on a single lipid structure, or several potentially different lipid structures (method 1), the bilayer will start well-packed but requires more equilibration. When starting from a pre-equilibrated bilayer, either pure or mixed, most of the structure of the bilayer stays intact, reducing equilibration time (method 2). The main advantages of these methods are that they minimize equilibration time and can be almost completely automated, nearly eliminating one time consuming step in MD simulations of membrane proteins.     

All-atom molecular dynamics simulations using orientational constraints from anisotropic NMR samples

Orientational constraints obtained from solid state NMR experiments on anisotropic samples are used here in molecular dynamics (MD) simulations for determining the structure and dynamics of several different membrane-bound molecules. The new MD technique is based on the inclusion of orientation dependent pseudo-forces in the COSMOS-NMR force field. These forces drive molecular rotations and re-orientations in the simulation, such that the motional time-averages of the tensorial NMR properties approach the experimentally measured parameters. The orientational-constraint-driven MD simulations are universally applicable to all NMR interaction tensors, such as chemical shifts, dipolar couplings and quadrupolar interactions. The strategy does not depend on the initial choice of coordinates, and is in principle suitable for any flexible molecule. To test the method on three systems of increasing complexity, we used as constraints some deuterium quadrupolar couplings from the literature on pyrene, cholesterol and an antimicrobial peptide embedded in oriented lipid bilayers. The MD simulations were able to reproduce the NMR parameters within experimental error. The alignment of the three membrane-bound molecules and some aspects of their conformation were thus derived from the NMR data, in good agreement with previous analyses. Furthermore, the new approach yielded for the first time the distribution of segmental orientations with respect to the membrane and the order parameter tensors of all three systems.     

NMR-Based Simulation Studies of Pf1 Coat Protein in Explicit Membranes

As time- and ensemble-averaged measures, NMR observables contain information about both protein structure and dynamics. This work represents a computational study to extract such information for membrane proteins from orientation-dependent NMR observables: solid-state NMR chemical shift anisotropy and dipolar coupling, and solution NMR residual dipolar coupling. We have performed NMR-restrained molecular dynamics simulations to refine the structure of the membrane-bound form of Pf1 coat protein in explicit lipid bilayers using the recently measured chemical shift anisotropy, dipolar coupling, and residual dipolar coupling data. From the simulations, we have characterized detailed protein-lipid interactions and explored the dynamics. All simulations are stable and the NMR restraints are well satisfied. The C-terminal transmembrane (TM) domain of Pf1 finds its optimal position in the membrane quickly (within 6 ns), illustrating efficient solvation of TM domains in explicit bilayer environments. Such rapid convergence also leads to well-converged interaction patterns between the TM helix and the membrane, which clearly show the interactions of interfacial membrane-anchoring residues with the lipids. For the N-terminal periplasmic helix of Pf1, we identify a stable, albeit dynamic, helix orientation parallel to the membrane surface that satisfies the amphiphatic nature of the helix in an explicit lipid bilayer. Such detailed information cannot be obtained solely from NMR observables. Therefore, the present simulations illustrate the usefulness of NMR-restrained MD refinement of membrane protein structure in explicit membranes.     

NMR-Based Simulation Studies of Pf1 Coat Protein in Explicit Membranes

As time- and ensemble-averaged measures, NMR observables contain information about both protein structure and dynamics. This work represents a computational study to extract such information for membrane proteins from orientation-dependent NMR observables: solid-state NMR chemical shift anisotropy and dipolar coupling, and solution NMR residual dipolar coupling. We have performed NMR-restrained molecular dynamics simulations to refine the structure of the membrane-bound form of Pf1 coat protein in explicit lipid bilayers using the recently measured chemical shift anisotropy, dipolar coupling, and residual dipolar coupling data. From the simulations, we have characterized detailed protein-lipid interactions and explored the dynamics. All simulations are stable and the NMR restraints are well satisfied. The C-terminal transmembrane (TM) domain of Pf1 finds its optimal position in the membrane quickly (within 6 ns), illustrating efficient solvation of TM domains in explicit bilayer environments. Such rapid convergence also leads to well-converged interaction patterns between the TM helix and the membrane, which clearly show the interactions of interfacial membrane-anchoring residues with the lipids. For the N-terminal periplasmic helix of Pf1, we identify a stable, albeit dynamic, helix orientation parallel to the membrane surface that satisfies the amphiphatic nature of the helix in an explicit lipid bilayer. Such detailed information cannot be obtained solely from NMR observables. Therefore, the present simulations illustrate the usefulness of NMR-restrained MD refinement of membrane protein structure in explicit membranes.     

Biophysical and computational studies of membrane penetration by the GRP1 pleckstrin homology domain

The pleckstrin homology (PH) domain of the general receptor for phosphoinositides 1 (GRP1) exhibits specific, high-affinity, reversible binding to phosphatidylinositol (3,4,5)-trisphosphate (PI(3,4,5)P(3)) at?the plasma membrane, but the nature and extent of the interaction between this bound complex and the surrounding membrane environment remains unclear. Combining equilibrium and nonequilibrium molecular dynamics (MD) simulations, NMR spectroscopy, and monolayer penetration experiments, we characterize the membrane-associated state of?GRP1-PH. MD simulations show loops flanking the binding site supplement the interaction with PI(3,4,5)P(3) through multiple contacts with the lipid bilayer. NMR data show large perturbations in chemical shift for these loop regions on binding to PI(3,4,5)P(3)-containing DPC micelles. Monolayer penetration experiments and further MD simulations demonstrate that mutating hydrophobic residues to polar residues in the flanking loops reduces membrane penetration. This supports a "dual-recognition" model of binding, with specific GRP1-PH-PI(3,4,5)P(3) interactions supplemented by interactions of loop regions with the lipid bilayer.     

Multiscale modeling of lipids and lipid bilayers

A multiscale modeling approach is applied for simulations of lipids and lipid assemblies on mesoscale. First, molecular dynamics simulation of initially disordered system of lipid molecules in water within all-atomic model was carried out. On the next stage, structural data obtained from the molecular dynamics (MD) simulation were used to build a coarse-grained (ten sites) lipid model, with effective interaction potentials computed by the inverse Monte Carlo method. Finally, several simulations of the coarse-grained model on longer length- and time-scale were performed, both within Monte Carlo and molecular dynamics simulations: a periodical sample of lipid molecules ordered in bilayer, a free sheet of such bilayer without periodic boundary conditions, formation of vesicle from a plain membrane, process of self-assembly of lipids randomly dispersed in volume. It was shown that the coarse-grained model, developed exclusively from all-atomic simulation data, reproduces well all the basic features of lipids in water solution.     

Molecular dynamics simulation of neuropeptide B and neuropeptide W in the dipalmitoylphosphatidylcholine membrane bilayer

GPR7 and GPR8 are recently deorphanized G-protein-coupled receptors that are implicated in the regulation of neuroendocrine function, feeding behavior, and energy homeostasis. Neuropeptide B (NPB) and neuropeptide W (NPW) are two membrane-bound hypothalamic peptides, which specifically antagonize GPR7 and GPR8. Despite years of research, an accurate estimation of structure and molecular recognition of these neuropeptide systems still remains elusive. Herein, we investigated the structure, orientation, and interaction of NPB and NPW in a dipalmitoylphosphatidylcholine bilayer using long-range molecular dynamics (MD) simulation. During 30-ns simulation, membrane-embedded helical axes of NPB and NPW tilted 30 and 15°, respectively, from the membrane normal in order to overcome possible hydrophobic mismatch with the lipid bilayer. The calculation of various structural parameters indicated that NPW is more rigid and compact as compared to NPB. Qualitatively, the peptides exhibited flexible N-terminal (residues 1–12) and rigid C-terminal α-helical parts (residues 13–21), confirming previous NMR data. A strong electrostatic attraction between C-termini and headgroup atoms caused translocation of the peptides towards lower leaflet of the bilayer. The stabilizing hydrogen bonds (H-bonds) between phosphate groups and Trp1, Lys3, and Arg15 of the peptides played important roles for membrane anchoring. MD simulations of Alanine (Ala) mutants revealed that WYK->Ala variant of NPB/NPW lacked crucial H-bond interactions with phospholipid headgroups and also caused severe misfolding in NPB. Altogether, the knowledge of preferred structural fold and interaction of neuropeptides within the membrane bilayer will be useful to develop synthetic agonist or antagonist peptides for GPR7 and GPR8.     

Transmembrane helix-helix interactions: comparative simulations of the glycophorin a dimer

The glycophorin helix dimer is a paradigm for the exploration of helix-helix interactions in integral membrane proteins. Two NMR structures of the dimer are known, one in a detergent micelle and one in a lipid bilayer. Multiple (4 x 50 ns) molecular dynamics simulations starting from each of the two NMR structures, with each structure in either a dodecyl phosphocholine (DPC) micelle or a dimyristoyl phosphatidylcholine (DMPC) bilayer, have been used to explore the conformational dynamics of the helix dimer. Analysis of the helix-helix interaction, mediated by the GxxxG sequence motif, suggests convergence of the simulations to a common model. This is closer to the NMR structure determined in a bilayer than to micelle structure. The stable dimer interface in the final simulation model is characterized by (i) Gly/Gly packing and (ii) Thr/Thr interhelix H-bonds. These results demonstrate the ability of extended molecular dynamics simulations in a lipid bilayer environment to refine membrane protein structures or models derived from experimental data obtained in protein/detergent micelles.     

A molecular dynamics study of the bee venom melittin in aqueous solution, in methanol, and inserted in a phospholipid bilayer

The structural properties of melittin, a small amphipathic peptide found in the bee venom, are investigated in three different environments by molecular dynamics simulation. Long simulations have been performed for monomeric melittin solvated in water, in methanol, and shorter ones for melittin inserted in a dimyristoylphosphatidylcholine bilayer. The resulting trajectories were analysed in terms of structural properties of the peptide and compared to the available NMR data. While in water and methanol solution melittin is observed to partly unfold, the peptide retains its structure when embedded in a lipid bilayer. The latter simulation shows good agreement with the experimentally derived ³J-coupling constants. Generally, it appears that higher the stability of the helical conformation of melittin, lower is the dielectric permittivity of the environment. In addition, peptide-lipid interactions were investigated showing that the C-terminus of the peptide provides an anchor to the lipid bilayer by forming hydrogen bonds with the lipid head groups.     

Arginine in Membranes: The Connection Between Molecular Dynamics Simulations and Translocon-Mediated Insertion Experiments

Several laboratories have carried out molecular dynamics (MD) simulations of arginine interactions with lipid bilayers and found that the energetic cost of placing arginine in lipid bilayers is an order of magnitude greater than observed in molecular biology experiments in which Arg-containing transmembrane helices are inserted across the endoplasmic reticulum membrane by the Sec61 translocon. We attempt here to reconcile the results of the two approaches. We first present MD simulations of guanidinium groups alone in lipid bilayers, and then, to mimic the molecular biology experiments, we present simulations of hydrophobic helices containing single Arg residues at different positions along the helix. We discuss the simulation results in the context of molecular biology results and show that the energetic discrepancy is reduced, but not eliminated, by considering free energy differences between Arg at the interface and at the center of the model helices. The reduction occurs because Arg snorkeling to the interface prevents Arg from residing in the bilayer center where the energetic cost of desolvation is highest. We then show that the problem with MD simulations is that they measure water-to-bilayer free energies, whereas the molecular biology experiments measure the energetics of partitioning from translocon to bilayer, which raises the fundamental question of the relationship between water-to-bilayer and water-to-translocon partitioning. We present two thermodynamic scenarios as a foundation for reconciliation of the simulation and molecular biology results. The simplest scenario is that translocon-to-bilayer partitioning is independent of water-to-bilayer partitioning; there is no thermodynamic cycle connecting the two paths.     

Molecular ordering of interfacially localized tryptophan analogs in ester- and ether-lipid bilayers studied by 2H-NMR.

Perdeuterated indole-d6 and N-methylated indole-d6 were solubilized in lamellar liquid crystalline phases composed of either 1,2-diacyl-glycero-3-phosphocholine (14:0)/water or 1,2-dialkyl-glycero-3-phosphocholine(14:0/water. The molecular ordering of the tryptophan analogs was determined from deuteron quadrupole splittings observed in 2H-NMR spectra on macroscopically aligned lipid bilayers. NMR spectra were recorded with the bilayers oriented perpendicular to or parallel with the external magnetic field, and the values of the splittings differed by a factor of 2 between these distinct orientations, indicating fast rotational motion of the molecules about an axis parallel to the bilayer normal. In all cases the splittings were found to decrease with increasing temperature. Relatively large splittings were observed in all systems, demonstrating that the tryptophans partition into a highly anisotropic environment. Solubilization most likely occurs at the lipid/water interface, as indicated by 1H-NMR chemical shift studies. The 2H-NMR spectra obtained for each analog were found to be rather similar in ester and ether lipids, but with smaller splittings in the ether lipid under similar conditions. The difference was slightly less for the indole molecule. Furthermore, in both lipid systems the positions of the splittings from indole were different from those of N-methyl indole. The results suggest that 1) the tryptophan analogs are solubilized in the interfacial region of the lipid bilayer, 2) the behavior may be modulated by hydrogen bonding in the case of indole, and 3) hydrogen bonding with the lipid carbonyl groups is not likely to play a major role in the solubilization of single indole molecules in the ester lipid bilayer interface.     

Constructing the suitable initial configuration of the membrane-protein system in molecular dynamics simulations

A method for constructing the suitable initial configuration of the membrane-protein system for molecular dynamics (MD) simulations is presented. This method could provide some hydrated initial configurations and help us to determine the best surface area of the system by contracting the surface area and comparing the optimized lowest energy of the system by energy minimization. The gramicidin A (GA) channel in;the fully hydrated dimyristoylphosphatidylcholine (DMPC) bilayer was used as our model. Three configurations with different surface areas were selected and applied for one 400 ps and two 300 ps MD simulations at constant pressure and temperature. All simulations were fairly stable without any constraints. Through analysis of the MD trajectories we found that the system with the best surface area was more stable than the other two systems, whose sizes were changed in the simulations. Further analysis of the bilayer normal length and the order parameters of the lipid alkyl tails indicates that the system with the best surface area shows some characteristics of the L_α phase, while both the smaller and the larger size systems have distinct deviations from the L_α phase that we expect. This illustrates that the correct surface area and the suitable initial configuration have an important influence on the phase of the membrane in the MD simulation. In addition, by comparing the root mean square differences of GA relative to the initial structure and interaction energy between different components of the system for all three systems, we find that the state of the DMPC bilayer has exerted a significant influence on the structure of GA. All these results demonstrate the validity of our method for constructing the initial configuration of the membrane-protein system for MD simulations. Received: 10 September 1998 / Revised version: 19 March 1999 / Accepted: 19 March 1999     

Atomic view of calcium-induced clustering of phosphatidylserine in mixed lipid bilayers

Membranes play key regulatory roles in biological processes, with bilayer composition exerting marked effects on binding affinities and catalytic activities of a number of membrane-associated proteins. In particular, proteins involved in diverse processes such as vesicle fusion, intracellular signaling cascades, and blood coagulation interact specifically with anionic lipids such as phosphatidylserine (PS) in the presence of Ca(2+) ions. While Ca(2+) is suspected to induce PS clustering in mixed phospholipid bilayers, the detailed structural effects of this ion on anionic lipids are not established. In this study, combining magic angle spinning (MAS) solid-state NMR (SSNMR) measurements of isotopically labeled serine headgroups in mixed lipid bilayers with molecular dynamics (MD) simulations of PS lipid bilayers in the presence of different counterions, we provide site-resolved insights into the effects of Ca(2+) on the structure and dynamics of lipid bilayers. Ca(2+)-induced conformational changes of PS in mixed bilayers are observed in both liposomes and Nanodiscs, a nanoscale membrane mimetic of bilayer patches. Site-resolved multidimensional correlation SSNMR spectra of bilayers containing (13)C,(15)N-labeled PS demonstrate that Ca(2+) ions promote two major PS headgroup conformations, which are well resolved in two-dimensional (13)C-(13)C, (15)N-(13)C, and (31)P-(13)C spectra. The results of MD simulations performed on PS lipid bilayers in the presence or absence of Ca(2+) provide an atomic view of the conformational effects underlying the observed spectra.     

Membrane perturbation by the antimicrobial peptide PMAP-23: A fluorescence and molecular dynamics study

Several bioactive peptides exert their biological function by interacting with cellular membranes. Structural data on their location inside lipid bilayers are thus essential for a detailed understanding of their mechanism of action. We propose here a combined approach in which fluorescence spectroscopy and molecular dynamics (MD) simulations were applied to investigate the mechanism of membrane perturbation by the antimicrobial peptide PMAP-23. Fluorescence spectra, depth-dependent quenching experiments, and peptide-translocation assays were employed to determine the location of the peptide inside the membrane. MD simulations were performed starting from a random mixture of water, lipids and peptide, and following the spontaneous self-assembly of the bilayer. Both experimental and theoretical data indicated a peptide location just below the polar headgroups of the membrane, with an orientation essentially parallel to the bilayer plane. These findings, together with experimental results on peptide-induced leakage from large and giant vesicles, lipid flip-flop and peptide exchange between vesicles, support a mechanism of action consistent with the “carpet” model. Furthermore, the atomic detail provided by the simulations suggested the occurrence of an additional, more specific and novel mechanism of bilayer destabilization by PMAP-23, involving the unusual insertion of charged side chains into the hydrophobic core of the membrane.     

Molecular Dynamics and NMR Spectroscopy Studies of E.?coli Lipopolysaccharide Structure and Dynamics

Lipopolysaccharide (LPS), a component of Gram-negative bacterial outer membranes, comprises three regions: lipid A, core oligosaccharide, and O-antigen polysaccharide. Using the CHARMM36 lipid and carbohydrate force fields, we have constructed a model of an Escherichia coli R1 (core) O6 (antigen) LPS molecule. Several all-atom bilayers are built and simulated with lipid A only (LIPA) and varying lengths of 0 (LPS0), 5 (LPS5), and 10 (LPS10) O6 antigen repeating units; a single unit of O6 antigen contains five sugar residues. From ¹H,¹H-NOESY experiments, cross-relaxation rates are obtained from an O-antigen polysaccharide sample. Although some experimental deviations are due to spin-diffusion, the remaining effective proton-proton distances show generally very good agreement between NMR experiments and molecular dynamics simulations. The simulation results show that increasing the LPS molecular length has an impact on LPS structure and dynamics and also on LPS bilayer properties. Terminal residues in a LPS bilayer are more flexible and extended along the membrane normal. As the core and O-antigen are added, per-lipid area increases and lipid bilayer order decreases. In addition, results from mixed LPS0/5 and LPS0/10 bilayer simulations show that the LPS O-antigen conformations at a higher concentration of LPS5 and LPS10 are more orthogonal to the membrane and less flexible. The O-antigen concentration of mixed LPS bilayers does not have a significant effect on per-lipid area and hydrophobic thickness. Analysis of ion and water penetration shows that water molecules can penetrate inside the inner core region, and hydration is critical to maintain the integrity of the bilayer structure.     

Molecular Dynamics and NMR Spectroscopy Studies of E. coli Lipopolysaccharide Structure and Dynamics

Lipopolysaccharide (LPS), a component of Gram-negative bacterial outer membranes, comprises three regions: lipid A, core oligosaccharide, and O-antigen polysaccharide. Using the CHARMM36 lipid and carbohydrate force fields, we have constructed a model of an Escherichia coli R1 (core) O6 (antigen) LPS molecule. Several all-atom bilayers are built and simulated with lipid A only (LIPA) and varying lengths of 0 (LPS0), 5 (LPS5), and 10 (LPS10) O6 antigen repeating units; a single unit of O6 antigen contains five sugar residues. From ¹H,¹H-NOESY experiments, cross-relaxation rates are obtained from an O-antigen polysaccharide sample. Although some experimental deviations are due to spin-diffusion, the remaining effective proton-proton distances show generally very good agreement between NMR experiments and molecular dynamics simulations. The simulation results show that increasing the LPS molecular length has an impact on LPS structure and dynamics and also on LPS bilayer properties. Terminal residues in a LPS bilayer are more flexible and extended along the membrane normal. As the core and O-antigen are added, per-lipid area increases and lipid bilayer order decreases. In addition, results from mixed LPS0/5 and LPS0/10 bilayer simulations show that the LPS O-antigen conformations at a higher concentration of LPS5 and LPS10 are more orthogonal to the membrane and less flexible. The O-antigen concentration of mixed LPS bilayers does not have a significant effect on per-lipid area and hydrophobic thickness. Analysis of ion and water penetration shows that water molecules can penetrate inside the inner core region, and hydration is critical to maintain the integrity of the bilayer structure.     

Effect of membrane tension on the electric field and dipole potential of lipid bilayer membrane

The dipole potential of lipid bilayer membrane controls the difference in permeability of the membrane to oppositely charged ions. We have combined molecular dynamics (MD) simulations and experimental studies to determine changes in electric field and electrostatic potential of 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC) lipid bilayer in response to applied membrane tension. MD simulations based on CHARMM36 force field showed that electrostatic potential of DOPC bilayer decreases by ~45mV in the physiologically relevant range of membrane tension values (0 to 15dyn/cm). The electrostatic field exhibits a peak (~0.8×10(9)V/m) near the water/lipid interface which shifts by 0.9? towards the bilayer center at 15dyn/cm. Maximum membrane tension of 15dyn/cm caused 6.4% increase in area per lipid, 4.7% decrease in bilayer thickness and 1.4% increase in the volume of the bilayer. Dipole-potential sensitive fluorescent probes were used to detect membrane tension induced changes in DOPC vesicles exposed to osmotic stress. Experiments confirmed that dipole potential of DOPC bilayer decreases at higher membrane tensions. These results are suggestive of a potentially new mechanosensing mechanism by which mechanically induced structural changes in the lipid bilayer membrane could modulate the function of membrane proteins by altering electrostatic interactions and energetics of protein conformational states.     

Numerical studies of the membrane fluorescent dyes dynamics in ground and excited states

Fluorescence methods are widely used in studies of biological and model membranes. The dynamics of membrane fluorescent markers in their ground and excited electronic states and correlations with their molecular surrounding within the fully hydrated phospholipid bilayer are still not well understood. In the present work, Quantum Mechanical (QM) calculations and Molecular Dynamics (MD) simulations are used to characterize location and interactions of two membrane polarity probes (Prodan; 6-propionyl-2-dimethylaminonaphthalene and its derivative Laurdan; 2-dimethylamino-6-lauroylnaphthalene) with the dioleoylphosphatidylcholine (DOPC) lipid bilayer model. MD simulations with fluorophores in ground and excited states are found to be a useful tool to analyze the fluorescent dye dynamics and their immediate vicinity. The results of QM calculations and MD simulations are in excellent agreement with available experimental data. The calculation shows that the two amphiphilic dyes initially placed in bulk water diffuse within 10 ns towards their final location in the lipid bilayer. Analysis of solvent relaxation process in the aqueous phase occurs on the picoseconds timescale whereas it takes nanoseconds at the lipid/water interface. Four different relaxation time constants, corresponding to different relaxation processes, where observed when the dyes were embedded into the membrane.