Atomistic simulation of lipid and DiI dynamics in membrane bilayers under tension

Neutron spin echo shows pHLIP is capable of retarding membrane thickness fluctuations

Cell membranes are responsible for a range of biological processes that require interactions between lipids and proteins. While the effects of lipids on proteins are becoming better understood, our knowledge of how protein conformational changes influence membrane dynamics remains rudimentary. Here, we performed experiments and computer simulations to study the dynamic response of a lipid membrane to changes in the conformational state of pH-low insertion peptide (pHLIP), which transitions from a surface-associated (SA) state at neutral or basic pH to a transmembrane (TM) α-helix under acidic conditions. Our results show that TM-pHLIP significantly slows down membrane thickness fluctuations due to an increase in effective membrane viscosity. Our findings suggest a possible membrane regulatory mechanism, where the TM helix affects lipid chain conformations, and subsequently alters membrane fluctuations and viscosity.

Exploration of lipid bilayer mechanical properties using molecular dynamics simulation

Important biological structures known for their exceptional mechanical qualities, lipid bilayers are essential to many cellular functions. Fluidity, elasticity, permeability, stiffness, tensile strength, compressibility, shear viscosity, line tension, and curvature elasticity are some of the fundamental characteristics affecting their behavior. The purpose of this review is to examine these characteristics in more detail by molecular dynamics simulation, elucidating their importance and the elements that lead to their appearance in lipid bilayers. Comprehending the mechanical characteristics of lipid bilayers is critical for creating medications, drug delivery systems, and biomaterials that interact with biological membranes because it allows one to understand how these materials respond to different stresses and deformations. The influence of mechanical characteristics on important lipid bilayer properties is examined in this review. The mechanical properties of lipid bilayers were clarified through the use of molecular dynamics simulation analysis techniques, including bilayer thickness, stress-strain analysis, lipid bilayer area compressibility, membrane bending rigidity, and time- or ensemble-averaged the area per lipid evaluation. We explain the significance of molecular dynamics simulation analysis methods, providing important new information about the stability and dynamic behavior of the bilayer. In the end, we hope to use molecular dynamics simulation to provide a comprehensive understanding of the mechanical properties and behavior of lipid bilayers, laying the groundwork for further studies and applications. Taken together, careful investigation of these mechanical aspects deepens our understanding of the adaptive capacities and functional roles of lipid bilayers in biological environments.

Bilayer tension-induced clustering of the UPR sensor IRE1

The endoplasmic reticulum acts as a protein quality control center where a range of chaperones and foldases facilitates protein folding. IRE1 is a sensory transmembrane protein that transduces signals of proteotoxic stress by forming clusters and activating a cellular program called the unfolded protein response (UPR). Recently, membrane thickness variation due to membrane compositional changes have been shown to drive IRE1 cluster formation, activating the UPR even in the absence of proteotoxic stress. Here, we demonstrate a direct relationship between bilayer tension and UPR activation based on IRE1 dimer stability. The stability of the IRE1 dimer in a (50%DOPC-50%POPC) membrane at different applied bilayer tensions was analyzed via molecular dynamics simulations. The potential of mean force for IRE1 dimerization predicts a higher concentration of IRE1 dimers for both tensed and compressed ER membranes. This study shows that IRE1 may be a mechanosensitive membrane protein and establishes a direct biophysical relationship between bilayer tension and UPR activation.

Estimation of negative membrane tension in lipid bilayers and its effect on antimicrobial peptide magainin 2-induced pore formation

Positive membrane tension in the stretched plasma membrane of cells and in the stretched lipid bilayer of vesicles has been well analyzed quantitatively, whereas there is limited quantitative information on negative membrane tension in compressed plasma membranes and lipid bilayers. Here, we examined negative membrane tension quantitatively. First, we developed a theory to describe negative membrane tension by analyzing the free energy of lipid bilayers to obtain a theoretical equation for negative membrane tension. This allowed us to obtain an equation describing the negative membrane tension (σosm) for giant unilamellar vesicles (GUVs) in hypertonic solutions due to negative osmotic pressure (Π). Then, we experimentally estimated the negative membrane tension for GUVs in hypertonic solutions by measuring the rate constant (kr) of rupture of the GUVs induced by the constant tension (σex) due to an external force as a function of σex. We found that larger σex values were required to induce the rupture of GUVs under negative Π compared with GUVs in isotonic solution and quantitatively determined the negative membrane tension induced by Π (σosm) by the difference between these σex values. At small negative Π, the experimental values of negative σosm agree with their theoretical values within experimental error, but as negative Π increases, the deviation increases. Negative tension increased the stability of GUVs because higher tensions were required for GUV rupture, and the rate constant of antimicrobial peptide magainin 2-induced pore formation decreased.

Viscosity-aided electromechanical poration of cells for transfecting molecules

Cell poration technologies offer opportunities not only to understand the activities of biological molecules but also to investigate genetic manipulation possibilities. Unfortunately, transferring large molecules that can carry huge genomic information is challenging. Here, we demonstrate electromechanical poration using a core-shell-structured microbubble generator, consisting of a fine microelectrode covered with a dielectric material. By introducing a microcavity at its tip, we could concentrate the electrical field with the application of electric pulses and generate microbubbles for electromechanical stimulation of cells. Specifically, the technology enables transfection with molecules that are thousands of kDa even into osteoblasts and Chlamydomonas, which are generally considered to be difficult to inject. Notably, we found that the transfection efficiency can be enhanced by adjusting the viscosity of the cell suspension, which was presumably achieved by remodeling of the membrane cytoskeleton. The applicability of the approach to a variety of cell types opens up numerous emerging gene engineering applications.

The Cell Physiome: What Do We Need in a Computational Physiology Framework for Predicting Single-Cell Biology?

Modern biology and biomedicine are undergoing a big data explosion, needing advanced computational algorithms to extract mechanistic insights on the physiological state of living cells. We present the motivation for the Cell Physiome project: a framework and approach for creating, sharing, and using biophysics-based computational models of single-cell physiology. Using examples in calcium signaling, bioenergetics, and endosomal trafficking, we highlight the need for spatially detailed, biophysics-based computational models to uncover new mechanisms underlying cell biology. We review progress and challenges to date toward creating cell physiome models. We then introduce bond graphs as an efficient way to create cell physiome models that integrate chemical, mechanical, electromagnetic, and thermal processes while maintaining mass and energy balance. Bond graphs enhance modularization and reusability of computational models of cells at scale. We conclude with a look forward at steps that will help fully realize this exciting new field of mechanistic biomedical data science. Expected final online publication date for the Annual Review of Biomedical Data Science, Volume 5 is August 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Assessing hypoxic damage to placental trophoblasts by measuring membrane viscosity of extracellular vesicles

INTRODUCTION

As highly sophisticated intercellular communication vehicles in biological systems, extracellular vesicles (EVs) have been investigated as both promising liquid biopsy-based disease biomarkers and drug delivery carriers. Despite tremendous progress in understanding their biological and physiological functions, mechanical characterization of these nanoscale entities remains challenging due to the limited availability of proper techniques. Especially, whether damage to parental cells can be reflected by the mechanical properties of their EVs remains unknown.

METHODS

In this study, we characterized membrane viscosities of different types of EVs collected from primary human trophoblasts (PHTs), including apoptotic bodies, microvesicles and small extracellular vesicles, using fluorescence lifetime imaging microscopy (FLIM). The biochemical origin of EV membrane viscosity was examined by analyzing their phospholipid composition, using mass spectrometry.

RESULTS

We found that different EV types derived from the same cell type exhibit different membrane viscosities. The measured membrane viscosity values are well supported by the lipidomic analysis of the phospholipid compositions. We further demonstrate that the membrane viscosity of microvesicles can faithfully reveal hypoxic injury of the human trophoblasts. More specifically, the membrane of PHT microvesicles released under hypoxic condition is less viscous than its counterpart under standard culture condition, which is supported by the reduction in the phosphatidylethanolamine-to-phosphatidylcholine ratio in PHT microvesicles.

DISCUSSION

Our study suggests that biophysical properties of released trophoblastic microvesicles can reflect cell health. Characterizing EV's membrane viscosity may pave the way for the development of new EV-based clinical applications.

Membrane tension

The cell membrane serves as a barrier that restricts the rate of exchange of diffusible molecules. Tension in the membrane regulates many crucial cell functions involving shape changes and motility, cell signaling, endocytosis, and mechanosensation. Tension reflects the forces contributed by the lipid bilayer, the cytoskeleton, and the extracellular matrix. With a fluid-like bilayer model, membrane tension is presumed uniform and hence propagated instantaneously. In this review, we discuss techniques to measure the mean membrane tension and how to resolve the stresses in different components and consider the role of bilayer heterogeneity.

Effects of lipid heterogeneity on model human brain lipid membranes

Cell membranes naturally contain a heterogeneous lipid distribution. However, homogeneous bilayers are commonly preferred and utilised in computer simulations due to their relative simplicity, and the availability of lipid force field parameters. Recently, experimental lipidomics data for the human brain cell membranes under healthy and Alzheimer's disease (AD) conditions were investigated, since disruption to the lipid composition has been implicated in neurodegenerative disorders, including AD [R. B. Chan et al., J. Biol. Chem., 2012, 287, 2678-2688]. In order to observe the effects of lipid complexity on the various bilayer properties, molecular dynamics simulations were used to study four membranes with increasing heterogeneity: a pure POPC membrane, a POPC and cholesterol membrane in a 1 : 1 ratio (POPC-CHOL), and to our knowledge, the first realistic models of a healthy brain membrane and an Alzheimer's diseased brain membrane. Numerous structural, interfacial, and dynamical properties, including the area per lipid, interdigitation, dipole potential, and lateral diffusion of the two simple models, POPC and POPC-CHOL, were analysed and compared to those of the complex brain models consisting of 27 lipid components. As the membranes gain heterogeneity, a number of alterations were found in the structural and dynamical properties, and more significant differences were observed in the lateral diffusion. Additionally, we observed snorkeling behaviour of the lipid tails that may play a role in the permeation of small molecules across biological membranes. In this work, atomistic description of realistic brain membrane models is provided, which can add insight towards the permeability and transport pathways of small molecules across these membrane barriers.

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.

Amphiphilic copolymers change the nature of the ordered-to-disordered phase transition of lipid membranes from discontinuous to continuous

The phase behaviour as a function of temperature is explored for pure phospholipid (DPPC) and hybrid lipid-polymer (DPPC/Pluronic L64) bilayers with the aid of atomistic MD simulations. The range of the fixed-temperature simulations includes temperatures below and above the known melting temperature (Tm) for DPPC membranes. For the pure lipid bilayer, the main phase transition is discontinuous, as verified by the abrupt changes observed in the membrane structure, elasticity and the lipid diffusivity near the critical temperature Tm, which lies in the region 298.15-303.15 K. A pre-transition step is detected at 298.15 K which has been identified as the ripple phase (Pβ'), where ordered and disordered lipids coexist, causing thickness fluctuations. In the ordered gel phase, the positional ordering as assessed by the lipid radial distribution functions is long-range and some degree of hexagonal packing is measured. The hybrid bilayers on the other hand, transform from a more ordered to a disordered phase in a continuous manner, without finite jumps in their properties. No signs of the ripple phase are identified and the ordered phase exhibits very limited hexagonal packing and some positional ordering that decays fast. The effect of the inserted polymers in the two phases is reversed; at low temperatures, they render the membrane thinner, less cohesive and less ordered compared to the pure one, with the lipids assuming faster diffusion rates, whereas at high temperatures, the polymer interaction with the lipids acts reducing their diffusivity, but also increasing the lipid tail ordering and the membrane stiffness. The ability of the amphiphilic L64 copolymers to change the nature of the main phase transition of lipid membranes and their properties both in the ordered and the disordered phase is of vital importance for the prediction of both the efficacy of hybrid lipid/polymer nanoparticles as drug delivery vehicles as well as their potential adverse implications during interactions with healthy cell membranes.

Tip-Enhanced Infrared Difference-Nanospectroscopy of the Proton Pump Activity of Bacteriorhodopsin in Single Purple Membrane Patches

Photosensitive proteins embedded in the cell membrane (about 5 nm thickness) act as photoactivated proton pumps, ion gates, enzymes, or more generally, as initiators of stimuli for the cell activity. They are composed of a protein backbone and a covalently bound cofactor (e.g. the retinal chromophore in bacteriorhodopsin (BR), channelrhodopsin, and other opsins). The light-induced conformational changes of both the cofactor and the protein are at the basis of the physiological functions of photosensitive proteins. Despite the dramatic development of microscopy techniques, investigating conformational changes of proteins at the membrane monolayer level is still a big challenge. Techniques based on atomic force microscopy (AFM) can detect electric currents through protein monolayers and even molecular binding forces in single-protein molecules but not the conformational changes. For the latter, Fourier-transform infrared spectroscopy (FTIR) using difference-spectroscopy mode is typically employed, but it is performed on macroscopic liquid suspensions or thick films containing large amounts of purified photosensitive proteins. In this work, we develop AFM-assisted, tip-enhanced infrared difference-nanospectroscopy to investigate light-induced conformational changes of the bacteriorhodopsin mutant D96N in single submicrometric native purple membrane patches. We obtain a significant improvement compared with the signal-to-noise ratio of standard IR nanospectroscopy techniques by exploiting the field enhancement in the plasmonic nanogap that forms between a gold-coated AFM probe tip and an ultraflat gold surface, as further supported by electromagnetic and thermal simulations. IR difference-spectra in the 1450-1800 cm^-1 range are recorded from individual patches as thin as 10 nm, with a diameter of less than 500 nm, well beyond the diffraction limit for FTIR microspectroscopy. We find clear spectroscopic evidence of a branching of the photocycle for BR molecules in direct contact with the gold surfaces, with equal amounts of proteins either following the standard proton-pump photocycle or being trapped in an intermediate state not directly contributing to light-induced proton transport. Our results are particularly relevant for BR-based optoelectronic and energy-harvesting devices, where BR molecular monolayers are put in contact with metal surfaces, and, more generally, for AFM-based IR spectroscopy studies of conformational changes of proteins embedded in intrinsically heterogeneous native cell membranes.

Intracellular ion and protein nanoparticle-induced osmotic pressure modify astrocyte swelling and brain edema in response to glutamate stimuli

Intracellular tension activity plays a crucial role in cytotoxic brain edema and astrocyte swelling. Here, a few genetically encoded FRET-based tension probes were designed to detect cytoskeletal structural tension optically, including their magnitude and vectors. The astrocyte swelling resulted in GFAP tension increment, which is associated with the antagonistic effect of inward microfilaments (MFs) and microtubules (MTs) forces. In glutamate-induced astrocyte swelling, GFAP tension rise resulted from outward ion and protein nanoparticle-induced osmotic pressure (PN-OP) increases, where PN-OP could be elicited by MF and MT depolymerization, protein nanoparticle production, and activation of cofilin and stathmin-1. Attenuation of both ion osmotic pressure and PN-OP by drug combinations, together with free-radical scavenger, relieved cerebral edema in vivo. The study suggests that intracellular osmotic pressure (especially PN-OP) has a pivotal role in glutamate-induced astrocyte swelling and brain edema. Recovery of cytoplasmic potential is a promising target to develop new drugs and cure brain edema.

Molecular-Level "Observations" of the Behavior of Gold Nanoparticles in Aqueous Solution and Interacting with a Lipid Bilayer Membrane

We use coarse-grained molecular dynamics simulations to "observe" details of interactions between ligand-covered gold nanoparticles and a lipid bilayer model membrane. In molecular dynamics simulations, one puts the individual atoms and groups of atoms of the physical system to be "observed" into a simulation box, specifies the forms of the potential energies of interactions between them (ultimately quantum based), and lets them individually move classically according to Newton's equations of motion, based on the forces arising from the assumed potential energy forms. The atoms that are chemically bonded to each other stay chemically bonded, following known potentials (force fields) that permit internal degrees of freedom (internal rotation, torsion, vibrations), and the interactions between nonbonded atoms are simplified to Lennard-Jones forms (in our case) and coulombic (where electrical charges are present) in which the parameters are previously optimized to reproduce thermodynamic properties or are based on quantum electronic calculations. The system is started out at a reasonable set of coordinates for all atoms or groups of atoms, and then permitted to develop according to the equations of motion, one small step (usually 10 fs time step) at a time, for millions of steps until the system is at a quasi-equilibrium (usually reached after hundreds of nanoseconds). We then let the system play out its motions further for many nanoseconds to observe the behavior, periodically taking snapshots (saving all positions and energies), and post-processing the snapshots to obtain various average descriptions of the system. Alkanethiols of various lengths serve as examples of hydrophobic ligands and methyl-terminated PEG with various numbers of monomer units serve as examples of hydrophilic ligands. Spherical gold particles of various diameters as well as gold nanorods form the core to which ligands are attached. The nanoparticles are characterized at the molecular level, especially the distributions of ligand configurations and their dependence on ligand length, and surface coverage. Self-assembly of the bilayer from an isotropic solution and observation of membrane properties that correspond well to experimental values validate the simulations. The mechanism of permeation of a gold NP coated with either a hydrophobic or a hydrophilic ligand, and its dependence on surface coverage, ligand length, core diameter, and core shape, is investigated. Lipid response such as lipid flip-flops, lipid extraction, and changes in order parameter of the lipid tails are examined in detail. The mechanism of permeation of a PEGylated nanorod is shown to occur by tilting, lying down, rotating, and straightening up. The nature of the information provided by molecular dynamics simulations permits understanding of the detailed behavior of gold nanoparticles interacting with lipid membranes which in turn helps to understand why some known systems work better than others and aids the design of new particles and improvement of methods for preparing existing ones.

Emerging Role of Plasma Membranes in Vascular Endothelial Mechanosensing

Vascular endothelial cells (ECs) maintain circulatory system homeostasis by changing their functions in response to changes in hemodynamic forces, including shear stress and stretching. However, it is unclear how ECs sense changes in shear stress and stretching and transduce these changes into intracellular biochemical signals. The plasma membranes of ECs have recently been shown to respond to shear stress and stretching differently by rapidly changing their lipid order, fluidity, and cholesterol content. Such changes in the membranes' physical properties trigger the activation of membrane receptors and cell responses specific to each type of force. Artificial lipid-bilayer membranes show similar changes in lipid order in response to shear stress and stretching, indicating that they are physical phenomena rather than biological reactions. These findings suggest that the plasma membranes of ECs act as mechanosensors; in response to mechanical forces, they first alter their physical properties, modifying the conformation and function of membrane proteins, which then activates downstream signaling pathways. This new appreciation of plasma membranes as mechanosensors could help to explain the distinctive features of mechanotransduction in ECs involving shear stress and stretching, which activate a variety of membrane proteins and multiple signal transduction pathways almost simultaneously.

Distribution of mechanical stress in the Escherichia coli cell envelope

The cell envelope in Gram-negative bacteria comprises two distinct membranes with a cell wall between them. There has been a growing interest in understanding the mechanical adaptation of this cell envelope to the osmotic pressure (or turgor pressure), which is generated by the difference in the concentration of solutes between the cytoplasm and the external environment. However, it remains unexplored how the cell wall, the inner membrane (IM), and the outer membrane (OM) effectively protect the cell from this pressure by bearing the resulting surface tension, thus preventing the formation of inner membrane bulges, abnormal cell morphology, spheroplasts and cell lysis. In this study, we have used molecular dynamics (MD) simulations combined with experiments to resolve how and to what extent models of the IM, OM, and cell wall respond to changes in surface tension. We calculated the area compressibility modulus of all three components in simulations from tension-area isotherms. Experiments on monolayers mimicking individual leaflets of the IM and OM were also used to characterize their compressibility. While the membranes become softer as they expand, the cell wall exhibits significant strain stiffening at moderate to high tensions. We integrate these results into a model of the cell envelope in which the OM and cell wall share the tension at low turgor pressure (0.3 atm) but the tension in the cell wall dominates at high values (>1 atm).

Accurate In Silico Modeling of Asymmetric Bilayers Based on Biophysical Principles

Technological advances in the last decade have enabled the study of ever more complex and physiologically relevant model membranes to help dispel the mystery surrounding the role of plasma membrane asymmetry in various cellular processes. The slowly accumulating body of experimental data is fueling renewed interest in and the need for computational methods to support interpretations and address a wide range of problems that are still not amenable to direct experimental study. The specific appeal of molecular dynamics simulations for this purpose lies in their ability to access information at atomic resolution, which is useful for the formulation of testable mechanistic hypotheses. But, the range of questions that can be addressed reliably with such simulations is determined by the appropriate construction and simulation of asymmetric bilayer models. One essential way to achieve this goal is to follow rigorous biophysical criteria and principles. In this context, we show that the requirement for a robust comparison between the properties of simulated asymmetric and symmetric model membranes is for the tension in each bilayer leaflet to be zero. Commonly used methods for constructing asymmetric bilayers, including matching the average areas of the leaflets from the corresponding symmetric systems, do not ensure zero leaflet tension, thus precluding physically realistic changes in the areas of the two leaflets. We present, to our knowledge, a new method for identifying the ideal lipid packing in bilayers with different leaflet compositions that achieves the zero-tension goal, and discuss the basic principles underlying the biophysically correct computational study of asymmetric membranes.

Atomistic-level study of the interactions between hIAPP protofibrils and membranes: Influence of pH and lipid composition

The pathology of type 2 diabetes mellitus is associated with the aggregation of human islet amyloid polypeptide (hIAPP) and aggregation-mediated membrane disruption. The interactions of hIAPP aggregates with lipid membrane, as well as the effects of pH and lipid composition at the atomic level, remain elusive. Herein, using molecular dynamics simulations, we investigate the interactions of hIAPP protofibrillar oligomers with lipids, and the membrane perturbation that they induce, when they are partially inserted in an anionic dipalmitoyl-phosphatidylglycerol (DPPG) membrane or a mixed dipalmitoyl-phosphatidylcholine (DPPC)/DPPG (7:3) lipid bilayer under acidic/neutral pH conditions. We observed that the tilt angles and insertion depths of the hIAPP protofibril are strongly correlated with the pH and lipid composition. At neutral pH, the tilt angle and insertion depth of hIAPP protofibrils at a DPPG bilayer reach ~52° and ~1.62 nm with respect to the membrane surface, while they become ~77° and ~1.75 nm at a mixed DPPC/DPPG membrane. The calculated tilt angle of hIAPP at DPPG membrane is consistent with a recent chiral sum frequency generation spectroscopic study. The acidic pH induces a smaller tilt angle of ~40° and a shallower insertion depth (~1.24 nm) of hIAPP at the DPPG membrane surface, mainly due to protonation of His18 near the turn region. These differences mainly result from a combination of distinct electrostatic, van der Waals, hydrogen bonding and salt-bridge interactions between hIAPP and lipid bilayers. The hIAPP-membrane interaction energy analysis reveals that besides charged residues K1, R11 and H18, aromatic residues Phe15 and Phe23 also exhibit strong interactions with lipid bilayers, revealing the crucial role of aromatic residues in stabilizing the membrane-bound hIAPP protofibrils. hIAPP-membrane interactions disturb the lipid ordering and the local bilayer thickness around the peptides. Our results provide atomic-level information of membrane interaction of hIAPP protofibrils, revealing pH-dependent and membrane-modulated hIAPP aggregation at the early stage.

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.

Osmosensing by the bacterial PhoQ/PhoP two-component system

The PhoQ/PhoP two-component system plays an essential role in the response of enterobacteria to the environment of their mammalian hosts. It is known to sense several stimuli that are potentially associated with the host, including extracellular magnesium limitation, low pH, and the presence of cationic antimicrobial peptides. Here, we show that the PhoQ/PhoP two-component systems of Escherichia coli and Salmonella can also perceive an osmotic upshift, another key stimulus to which bacteria become exposed within the host. In contrast to most previously established stimuli of PhoQ, the detection of osmotic upshift does not require its periplasmic sensor domain. Instead, we show that the activity of PhoQ is affected by the length of the transmembrane (TM) helix as well as by membrane lateral pressure. We therefore propose that osmosensing relies on a conformational change within the TM domain of PhoQ induced by a perturbation in cell membrane thickness and lateral pressure under hyperosmotic conditions. Furthermore, the response mediated by the PhoQ/PhoP two-component system was found to improve bacterial growth recovery under hyperosmotic stress, partly through stabilization of the sigma factor RpoS. Our findings directly link the PhoQ/PhoP two-component system to bacterial osmosensing, suggesting that this system can mediate a concerted response to most of the established host-related cues.

How the Incorporation of Pluronic Block Copolymers Modulates the Response of Lipid Membranes to Mechanical Stress

We employ atomistic molecular dynamics simulations to investigate the effect that the incorporation of the nonionic amphiphilic copolymer known as Pluronic L64 has on the mechanical stability of a DPPC membrane. The simulations reveal that the incorporation of the polymer chains leads to membranes that can sustain increasing mechanical stresses. Analysis of mechanical, structural, and dynamic properties of the membrane shows that the polymer chains interact strongly with the lipids in the vicinity, restraining their mobility and imparting better mechanical stability to the membrane. The hybrid membranes under tension remain thicker, more ordered, and stiffer in comparison to their lipid analogues. Trans-bilayer lipid movements (flip-flop) are observed and appear to be triggered by the presence of the polymer chains. A careful analysis of the pore formation under high tensions reveals two distinctive mechanisms that depend on the distribution of the hydrophilic polymer blocks in the bilayer. Finally, the rate of growth of the formed membrane defects is slowed down in the presence of polymers. These findings show that Pluronic block copolymers could be exploited for the formation of optimized hybrid nanodevices with controlled elastic and dynamic properties.

Tension Independence of Lipid Diffusion and Membrane Viscosity

The diffusion of biomolecules at lipid membranes is governed by the viscosity of the underlying two-dimensionally fluid lipid bilayer. For common three-dimensional fluids, viscosity can be modulated by hydrostatic pressure, and pressure-viscosity data have been measured for decades. Remarkably, the two-dimensional analogue of this relationship, the dependence of molecular mobility on tension, has to the best of our knowledge never been measured for lipid bilayers, limiting our understanding of cellular mechanotransduction as well as the fundamental fluid mechanics of membranes. Here we report both molecular-scale and mesoscopic measures of fluidity in giant lipid vesicles as a function of mechanical tension applied using micropipette aspiration. Both molecular-scale data, from fluorescence recovery after photobleaching, and micron-scale data, from tracking the diffusion of phase-separated domains, show a surprisingly weak dependence of viscosity on tension, in contrast to predictions of recent molecular dynamics simulations, highlighting fundamental gaps in our understanding of membrane fluidity.

Experimental Estimation of Membrane Tension Induced by Osmotic Pressure

Osmotic pressure (Π) induces the stretching of plasma membranes of cells or lipid membranes of vesicles, which plays various roles in physiological functions. However, there have been no experimental estimations of the membrane tension of vesicles upon exposure to Π. In this report, we estimated experimentally the lateral tension of the membranes of giant unilamellar vesicles (GUVs) when they were transferred into a hypotonic solution. First, we investigated the effect of Π on the rate constant, k_p, of constant-tension (σ_ex)-induced rupture of dioleoylphosphatidylcholine (DOPC)-GUVs using the method developed by us recently. We obtained the σ_ex dependence of k_p in GUVs under Π and by comparing this result with that in the absence of Π, we estimated the tension of the membrane due to Π at the swelling equilibrium, σ_osm^eq. Next, we measured the volume change of DOPC-GUVs under small Π. The experimentally obtained values of σ_osm^eq and the volume change agreed with their theoretical values within the limits of the experimental errors. Finally, we investigated the characteristics of the Π-induced pore formation in GUVs. The σ_osm^eq corresponding to the threshold Π at which pore formation is induced is similar to the threshold tension of the σ_ex-induced rupture. The time course of the radius change of GUVs in the Π-induced pore formation depends on the total membrane tension, σ_t; for small σ_t, the radius increased with time to an equilibrium one, which remained constant for a long time until pore formation, but for large σ_t, the radius increased with time and pore formation occurred before the swelling equilibrium was reached. Based on these results, we discussed the σ_osm^eq and the Π-induced pore formation in lipid membranes.

GsMTx4: Mechanism of Inhibiting Mechanosensitive Ion Channels

GsMTx4 is a spider venom peptide that inhibits cationic mechanosensitive channels (MSCs). It has six lysine residues that have been proposed to affect membrane binding. We synthesized six analogs with single lysine-to-glutamate substitutions and tested them against Piezo1 channels in outside-out patches and independently measured lipid binding. Four analogs had ∼20% lower efficacy than the wild-type (WT) peptide. The equilibrium constants calculated from the rates of inhibition and washout did not correlate with the changes in inhibition. The lipid association strength of the WT GsMTx4 and the analogs was determined by tryptophan autofluorescence quenching and isothermal calorimetry with membrane vesicles and showed no significant differences in binding energy. Tryptophan fluorescence-quenching assays showed that both WT and analog peptides bound superficially near the lipid-water interface, although analogs penetrated deeper. Peptide-lipid association, as a function of lipid surface pressure, was investigated in Langmuir monolayers. The peptides occupied a large fraction of the expanded monolayer area, but that fraction was reduced by peptide expulsion as the pressure approached the monolayer-bilayer equivalence pressure. Analogs with compromised efficacy had pressure-area isotherms with steeper slopes in this region, suggesting tighter peptide association. The pressure-dependent redistribution of peptide between "deep" and "shallow" binding modes was supported by molecular dynamics (MD) simulations of the peptide-monolayer system under different area constraints. These data suggest a model placing GsMTx4 at the membrane surface, where it is stabilized by the lysines, and occupying a small fraction of the surface area in unstressed membranes. When applied tension reduces lateral pressure in the lipids, the peptides penetrate deeper acting as "area reservoirs" leading to partial relaxation of the outer monolayer, thereby reducing the effective magnitude of stimulus acting on the MSC gate.

Functionalized lipids and surfactants for specific applications

Synthetic lipids and surfactants that do not exist in biological systems have been used for the last few decades in both basic and applied science. The most notable applications for synthetic lipids and surfactants are drug delivery, gene transfection, as reporting molecules, and as support for structural lipid biology. In this review, we describe the potential of the synergistic combination of computational and experimental methodologies to study the behavior of synthetic lipids and surfactants embedded in lipid membranes and liposomes. We focused on select cases in which molecular dynamics simulations were used to complement experimental studies aiming to understand the structure and properties of new compounds at the atomistic level. We also describe cases in which molecular dynamics simulations were used to design new synthetic lipids and surfactants, as well as emerging fields for the application of these compounds. This article is part of a Special Issue entitled: Biosimulations edited by Ilpo Vattulainen and Tomasz Róg.

Physical Principles of Nanoparticle Cellular Endocytosis

This review article focuses on the physiochemical mechanisms underlying nanoparticle uptake into cells. When nanoparticles are in close vicinity to a cell, the interactions between the nanoparticles and the cell membrane generate forces from different origins. This leads to the membrane wrapping of the nanoparticles followed by cellular uptake. This article discusses how the kinetics, energetics, and forces are related to these interactions and dependent on the size, shape, and stiffness of nanoparticles, the biomechanical properties of the cell membrane, as well as the local environment of the cells. The discussed fundamental principles of the physiochemical causes for nanoparticle-cell interaction may guide new studies of nanoparticle endocytosis and lead to better strategies to design nanoparticle-based approaches for biomedical applications.

The role of substrate topography on the cellular uptake of nanoparticles

Improving targeting efficacy has been a central focus of the studies on nanoparticle (NP)-based drug delivery nanocarriers over the past decades. As cells actively sense and respond to the local physical environments, not only the NP design (e.g., size, shape, ligand density, etc.) but also the cell mechanics (e.g., stiffness, spreading, expressed receptors, etc.) affect the cellular uptake efficiency. While much work has been done to elucidate the roles of NP design for cells seeded on a flat tissue culture surface, how the local physical environments of cells mediate uptake of NPs remains unexplored, despite the widely known effect of local physical environments on cellular responses in vitro and disease states in vivo. Here, we report the active responses of human osteosarcoma cells to fibrous substrate topographies and the subsequent changes in the cellular uptake of NPs. Our experiments demonstrate that surface topography modulates cellular uptake efficacy by mediating cell spreading and membrane mechanics. The findings provide a concrete example of the regulative role of the physical environments of cells on cellular uptake of NPs, therefore advancing the rational design of NPs for enhanced drug delivery in targeted cancer therapy.

Live cell plasma membranes do not exhibit a miscibility phase transition over a wide range of temperatures

Lipid/cholesterol mixtures derived from cell membranes as well as their synthetic reconstitutions exhibit well-defined miscibility phase transitions and critical phenomena near physiological temperatures. This suggests that lipid/cholesterol-mediated phase separation plays a role in the organization of live cell membranes. However, macroscopic lipid-phase separation is not generally observed in cell membranes, and the degree to which properties of isolated lipid mixtures are preserved in the cell membrane remain unknown. A fundamental property of phase transitions is that the variation of tagged particle diffusion with temperature exhibits an abrupt change as the system passes through the transition, even when the two phases are distributed in a nanometer-scale emulsion. We support this using a variety of Monte Carlo and atomistic simulations on model lipid membrane systems. However, temperature-dependent fluorescence correlation spectroscopy of labeled lipids and membrane-anchored proteins in live cell membranes shows a consistently smooth increase in the diffusion coefficient as a function of temperature. We find no evidence of a discrete miscibility phase transition throughout a wide range of temperatures: 14-37 °C. This contrasts the behavior of giant plasma membrane vesicles (GPMVs) blebbed from the same cells, which do exhibit phase transitions and macroscopic phase separation. Fluorescence lifetime analysis of a DiI probe in both cases reveals a significant environmental difference between the live cell and the GPMV. Taken together, these data suggest the live cell membrane may avoid the miscibility phase transition inherent to its lipid constituents by actively regulating physical parameters, such as tension, in the membrane.

Surface-Tension Replica-Exchange Molecular Dynamics Method for Enhanced Sampling of Biological Membrane Systems

Conformational sampling is fundamentally important for simulating complex biomolecular systems. The generalized-ensemble algorithm, especially the temperature replica-exchange molecular dynamics method (T-REMD), is one of the most powerful methods to explore structures of biomolecules such as proteins, nucleic acids, carbohydrates, and also of lipid membranes. T-REMD simulations have focused on soluble proteins rather than membrane proteins or lipid bilayers, because explicit membranes do not keep their structural integrity at high temperature. Here, we propose a new generalized-ensemble algorithm for membrane systems, which we call the surface-tension REMD method. Each replica is simulated in the NPγT ensemble, and surface tensions in a pair of replicas are exchanged at certain intervals to enhance conformational sampling of the target membrane system. We test the method on two biological membrane systems: a fully hydrated DPPC (1,2-dipalmitoyl-sn-glycero-3-phosphatidylcholine) lipid bilayer and a WALP23-POPC (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine) membrane system. During these simulations, a random walk in surface tension space is realized. Large-scale lateral deformation (shrinking and stretching) of the membranes takes place in all of the replicas without collapse of the lipid bilayer structure. There is accelerated lateral diffusion of DPPC lipid molecules compared with conventional MD simulation, and a much wider range of tilt angle of the WALP23 peptide is sampled due to large deformation of the POPC lipid bilayer and through peptide-lipid interactions. Our method could be applicable to a wide variety of biological membrane systems.

Anomalous diffusion of proteins in sheared lipid membranes

We use coarse grained molecular dynamics simulations to investigate diffusion properties of sheared lipid membranes with embedded transmembrane proteins. In membranes without proteins, we find normal in-plane diffusion of lipids in all flow conditions. Protein embedded membranes behave quite differently: by imposing a simple shear flow and sliding the monolayers of the membrane over each other, the motion of protein clusters becomes strongly superdiffusive in the shear direction. In such a circumstance, the subdiffusion regime is predominant perpendicular to the flow. We show that superdiffusion is a result of accelerated chaotic motions of protein-lipid complexes within the membrane voids, which are generated by hydrophobic mismatch or the transport of lipids by proteins.

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.

Substrate stiffness regulates cellular uptake of nanoparticles

Nanoparticle (NP)-bioconjugates hold great promise for more sensitive disease diagnosis and more effective anticancer drug delivery compared with existing approaches. A critical aspect in both applications is cellular internalization of NPs, which is influenced by NP properties and cell surface mechanics. Despite considerable progress in optimization of the NP-bioconjugates for improved targeting, the role of substrate stiffness on cellular uptake has not been investigated. Using polyacrylamide (PA) hydrogels as model substrates with tunable stiffness, we quantified the relationship between substrate stiffness and cellular uptake of fluorescent NPs by bovine aortic endothelial cells (BAECs). We found that a stiffer substrate results in a higher total cellular uptake on a per cell basis, but a lower uptake per unit membrane area. To obtain a mechanistic understanding of the cellular uptake behavior, we developed a thermodynamic model that predicts that membrane spreading area and cell membrane tension are two key factors controlling cellular uptake of NPs, both of which are modulated by substrate stiffness. Our experimental and modeling results not only open up new avenues for engineering NP-based cancer cell targets for more effective in vivo delivery but also contribute an example of how the physical environment dictates cellular behavior and function.

Computer simulation study of nanoparticle interaction with a lipid membrane under mechanical stress

Pore formation of lipid bilayers under mechanical stress is critical to biological processes. A series of coarse grained molecular dynamics simulations of lipid bilayers with carbon nanoparticles different in size have been performed. Surface tension was applied to study the disruption of lipid bilayers by nanoparticles and the formation of pores inside the bilayers. The presence of small nanoparticles enhances the probability of water penetration thus decreasing the membrane rupture tension, while big nanoparticles have the opposite effect. Nanoparticle volume affects bilayer strength indirectly, and particle surface density can complicate the interaction. The structural, dynamic, elastic properties and lateral densities of lipid bilayers with nanoparticles under mechanical stress were analyzed. The results demonstrate the ability of nanoparticles to adjust the structural and dynamic properties of a lipid membrane, and to efficiently regulate the pore formation behavior and hydrophobicity of membranes.

Lipid Bilayers: The Effect of Force Field on Ordering and Dynamics

The sensitivity of the structure and dynamics of a fully hydrated pure bilayer of 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) in molecular dynamics simulations to changes in force-field and simulation parameters has been assessed. Three related force fields (the Gromos 54A7 force field, a Gromos 53A6-derived parameter set and a variant of the Berger parameters) in combination with either particle-mesh Ewald (PME) or a reaction field (RF) were compared. Structural properties such as the area per lipid, carbon-deuterium order parameters, electron density profile and bilayer thicknesses, are reproduced by all the parameter sets within the uncertainty of the available experimental data. However, there are clear differences in the ordering of the glycerol backbone and choline headgroup, and the orientation of the headgroup dipole. In some cases, the degree of ordering was reminiscent of a liquid-ordered phase. It is also shown that, although the lateral diffusion of the lipids in the plane of the bilayer is often used to validate lipid force fields, because of the uncertainty in the experimental measurements and the fact that the lateral diffusion is dependent on the choice of the simulation conditions, it should not be employed as a measure of quality. Finally, the simulations show that the effect of small changes in force-field parameters on the structure and dynamics of a bilayer is more significant than the treatment of the long-range electrostatic interactions using RF or PME. Overall, the Gromos 54A7 best reproduced the range of experimental data examined.

Mechanical properties of lipid bilayers and regulation of mechanosensitive function: from biological to biomimetic channels

Material properties of lipid bilayers, including thickness, intrinsic curvature and compressibility regulate the function of mechanosensitive (MS) channels. This regulation is dependent on phospholipid composition, lateral packing and organization within the membrane. Therefore, a more complete framework to understand the functioning of MS channels requires insights into bilayer structure, thermodynamics and phospholipid structure, as well as lipid-protein interactions. Phospholipids and MS channels interact with each other mainly through electrostatic forces and hydrophobic matching, which are also crucial for antimicrobial peptides. They are excellent models for studying the formation and stabilization of membrane pores. Importantly, they perform equivalent responses as MS channels: (1) tilting in response to tension and (2) dissipation of osmotic gradients. Lessons learned from pore forming peptides could enrich our knowledge of mechanisms of action and evolution of these channels. Here, the current state of the art is presented and general principles of membrane regulation of mechanosensitive function are discussed.

Effect of membrane tension on the physical properties of DOPC lipid bilayer membrane

Molecular dynamics simulations of a dioleoylphosphocholine (DOPC) lipid bilayer were performed to explore its mechanosensitivity. Variations in the bilayer properties, such as area per lipid, volume, thickness, hydration depth (HD), hydration thickness (HT), lateral diffusion coefficient, and changes in lipid structural order were computed in the membrane tension range 0 to 15dyn/cm. We determined that an increase in membrane tension results in a decrease in the bilayer thickness and HD of ~5% and ~5.7% respectively, whereas area per lipid, volume, and HT/HD increased by 6.8%, 2.4%, and 5% respectively. The changes in lipid conformation and orientation were characterized using orientational (S(2)) and deuterium (S(CD)) order parameters. Upon increase of membrane tension both order parameters indicated an increase in lipid disorder by 10-20%, mostly in the tail end region of the hydrophobic chains. The effect of membrane tension on lipid lateral diffusion in the DOPC bilayer was analyzed on three different time scales corresponding to inertial motion, anomalous diffusion and normal diffusion. The results showed that lateral diffusion of lipid molecules is anomalous in nature due to the non-exponential distribution of waiting times. The anomalous and normal diffusion coefficients increased by 20% and 52% when the membrane tension changed from 0 to 15dyn/cm, respectively. In conclusion, our studies showed that membrane tension causes relatively significant changes in the area per lipid, volume, polarity, membrane thickness, and fluidity of the membrane suggesting multiple mechanisms by which mechanical perturbation of the membrane could trigger mechanosensitive response in cells.

Quantitative scheme for full-field polarization rotating fluorescence microscopy using a liquid crystal variable retarder

We present a quantitative scheme for full-field polarization rotating fluorescence microscopy. A quarter-wave plate, in combination with a liquid crystal variable retarder, provides a tunable method to rotate polarization states of light prior to its being coupled into a fluorescence microscope. A calibration of the polarization properties of the incident light is performed in order to correct for elliptical polarization states. This calibration allows the response of the sample to linear polarization states of light to be recovered. Three known polarization states of light can be used to determine the average fluorescent dipole orientations in the presence of a spatially varying dc offset or background polarization-invariant fluorescence signal. To demonstrate the capabilities of this device, we measured a series of full-field fluorescence polarization images from fluorescent analogs incorporated in the lipid membrane of Burkitts lymphoma CA46 cells. The fluorescent lipid-like analogs used in this study are molecules that are labeled by either a DiI (1,1(')-Dioctadecyl 3,3,3('),3(')-Tetramethylindocarbocyanine) fluorophore in its head group or a Bodipy (4,4-difluoro-4-bora-3a,4a-diaza-s-indacene) molecule in its acyl chain. A spatially varying contrast in the normalized amplitude was observed on the cell surface, where the orientation of the DiI molecules is tangential to the cell membrane. The internally labeled cellular structures showed zero response to changes in linear polarization, and the net linear polarization amplitude for these regions was zero. This instrument provides a low cost calibrated method that may be coupled to existing fluorescence microscopes to perform investigations of cellular processes that involve a change in molecular orientations.

Tuning membrane phase separation using nonlipid amphiphiles

Lipid phase separation may be a mechanism by which lipids participate in sorting membrane proteins and facilitate membrane-mediated biochemical signaling in cells. To provide new tools for membrane lipid phase manipulation that avoid direct effects on protein activity and lipid composition, we studied phase separation in binary and ternary lipid mixtures under the influence of three nonlipid amphiphiles, vitamin E (VE), Triton-X (TX)-100, and benzyl alcohol (BA). Mechanisms of additive-induced phase separation were elucidated using coarse-grained molecular dynamics simulations of these additives in a liquid bilayer made from 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) and 1,2-dilinoleoyl-sn-glycero-3-phosphocholine [corrected]. From simulations, the additive's partitioning preference, changes in membrane thickness, and alterations in lipid order were quantified. Simulations showed that VE favored the DPPC phase but partitioned predominantly to the domain boundaries and lowered the tendency for domain formation, and therefore acted as a linactant. This simulated behavior was consistent with experimental observations in which VE promoted lipid mixing and dispersed domains in both gel/liquid and liquid-ordered/liquid-disordered systems. From simulation, BA partitioned predominantly to the DUPC phase, decreased lipid order there, and thinned the membrane. These actions explain why, experimentally, BA promoted phase separation in both binary and ternary lipid mixtures. In contrast, TX, a popular detergent used to isolate raft membranes in cells, exhibited equal preference for both phases, as demonstrated by simulations, but nonetheless, was a strong domain promoter in all lipid mixtures. Further analysis showed that TX increased membrane thickness of the DPPC phase to a greater extent than the DUPC phase and thus increased hydrophobic mismatch, which may explain experimental observation of phase separation in the presence of TX. In summary, these nonlipid amphiphiles provide new tools to tune domain formation in model vesicle systems and could provide the means to form or disperse membrane lipid domains in cells, in addition to the well-known methods involving cholesterol enrichment and sequestration.

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.