What Can We Learn about Cholesterol's Transmembrane Distribution Based on Cholesterol-Induced Changes in Membrane Dipole Potential?

Polarizable atomic multipole-based force field for cholesterol

Cholesterol is one of the essential component of lipid in membrane. We present a polarizable atomic multipole force field (FF) for the molecular dynamic simulation of cholesterol. The FF building process follows the computational framework as the atomic multipole optimized energetics for biomolecular applications (AMOEBA) model. In this framework, the electronics parameters, including atomic monopole moments, dipole moments, and quadrupole moments calculated from ab initio calculations in the gas phase, are applied to represent the charge distribution. Furthermore, the many-body polarization is modeled by following the same pattern of distributed atomic polarizabilities. Then, the bilayers composed of two typical phospholipid molecules, 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC), and 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC), in a range of different cholesterol concentrations are built and implemented by molecular dynamics (MD) simulations based on the proposed polarizable FF. The simulation results are statistically analyzed to validate the feasibility of the proposed FF. The structural properties of the bilayers are calculated to compare with the related experimental values. The MD values show the same trend of experimental values changes.Communicated by Ramaswamy H. Sarma.

Peptide Power: Mechanistic Insights into the Effect of Mitochondria-Targeted Tetrapeptides on Membrane Electrostatics from Molecular Simulations

Mitochondrial dysfunction is implicated in nine of the ten leading causes of death in the US, yet there are no FDA-approved therapeutics to treat it. Synthetic mitochondria-targeted peptides (MTPs), including the lead compound SS-31, offer promise, as they have been shown to restore healthy mitochondrial function and treat a variety of common diseases. At the cellular level, research has shown that MTPs accumulate strongly at the inner mitochondrial membrane (IMM), slow energy sinks (e.g., proton leaks), and improve ATP production. Modulation of electrostatic fields around the IMM has been implicated as a key aspect in the mechanism of action (MoA) of these peptides; however, molecular and mechanistic details have remained elusive. In this study, we employed all-atom molecular dynamics simulations (MD) to investigate the interactions of four MTPs with lipid bilayers and calculate their effect on structural and electrostatic properties. In agreement with previous experimental findings, we observed the modulation of the membrane surface and dipole potentials by MTPs. The simulations reveal that the MTPs achieve a reduction in the dipole potential by acting to disorder both lipid head groups and water layers proximal to the bilayer surface. We also find that MTPs decrease the bilayer thickness and increase the membrane's capacitance. These changes suggest that MTPs may enhance how much potential energy can be stored across the IMM at a given transmembrane potential difference. The MTPs also displace cations away from the bilayer surface, modulating the surface potential and offering an alternative mechanism for how these MTPs reduce mitochondrial energy sinks like proton leaks and mitigate Ca2+ accumulation stress. In conclusion, this study highlights the therapeutic potential of MTPs and underlines how interactions of MTPs with lipid bilayers serve as a fundamental component of their MoA.

Membrane Domain Anti-Registration Induces an Intrinsic Transmembrane Potential

Plasma membrane segregation into various nanoscale membrane domains is driven by distinct interactions between diverse lipids and proteins. Among them, liquid-ordered (Lo) membrane domains are defined as "lipid rafts" and liquid-disordered (Ld) ones as "lipid non-rafts". Using model membrane systems, both intra-leaflet and inter-leaflet dynamics of these membrane domains are widely studied. Nevertheless, the biological impact of the latter, which is accompanied by membrane domain registration/anti-registration, is far from clear. Hence, in this work, we studied the biological relevance of the membrane domain anti-registration using both all-atom molecular dynamics (MD) simulations and confocal fluorescence microscopy. All-atom MD simulations suggested an intrinsic transmembrane potential for the case of the membrane anti-registration (Lo/Ld). Meanwhile, confocal fluorescence microscopy experiments of HeLa and 293T cell lines indicated that membrane cholesterol depletion could significantly alter the transmembrane potential of cells. Considering differences in the cholesterol content between Lo and Ld membrane domains, our confocal fluorescence microscopy experiments are consistent with our all-atom MD simulations. In short, membrane domain anti-registration induces local membrane asymmetry and, thus, an intrinsic transmembrane potential.

Membrane Dipole Potential: An Emerging Approach to Explore Membrane Organization and Function

Biological membranes are complex organized molecular assemblies of lipids and proteins that provide cells and membrane-bound intracellular organelles their individual identities by morphological compartmentalization. Membrane dipole potential originates from the electrostatic potential difference within the membrane due to the nonrandom arrangement (orientation) of amphiphile and solvent (water) dipoles at the membrane interface. In this Feature Article, we will focus on the measurement of dipole potential using electrochromic fluorescent probes and highlight interesting applications. In addition, we will focus on ratiometric fluorescence microscopic imaging technique to measure dipole potential in cellular membranes, a technique that can be used to address novel problems in cell biology which are otherwise difficult to address using available approaches. We envision that membrane dipole potential could turn out to be a convenient tool in exploring the complex interplay between membrane lipids and proteins and could provide novel insights in membrane organization and function.

Label-Free Imaging of Membrane Potentials by Intramembrane Field Modulation, Assessed by Second Harmonic Generation Microscopy

Optical interrogation of cellular electrical activity has proven itself essential for understanding cellular function and communication in complex networks. Voltage-sensitive dyes are important tools for assessing excitability but these highly lipophilic sensors may affect cellular function. Label-free techniques offer a major advantage as they eliminate the need for these external probes. In this work, it is shown that endogenous second-harmonic generation (SHG) from live cells is highly sensitive to changes in transmembrane potential (TMP). Simultaneous electrophysiological control of a living human embryonic kidney (HEK293T) cell, through a whole-cell voltage-clamp reveals a linear relation between the SHG intensity and membrane voltage. The results suggest that due to the high ionic strengths and fast optical response of biofluids, membrane hydration is not the main contributor to the observed field sensitivity. A conceptual framework is further provided that indicates that the SHG voltage sensitivity reflects the electric field within the biological asymmetric lipid bilayer owing to a nonzero χeff(2) tensor. Changing the TMP without surface modifications such as electrolyte screening offers high optical sensitivity to membrane voltage (≈40% per 100 mV), indicating the power of SHG for label-free read-out. These results hold promise for the design of a non-invasive label-free read-out tool for electrogenic cells.

Phosphatidylserine Exposed Lipid Bilayer Models for Understanding Cancer Cell Selectivity of Natural Compounds: A Molecular Dynamics Simulation Study

Development of drugs that are selectively toxic to cancer cells and safe to normal cells is crucial in cancer treatment. Evaluation of membrane permeability is a key metric for successful drug development. In this study, we have used in silico molecular models of lipid bilayers to explore the effect of phosphatidylserine (PS) exposure in cancer cells on membrane permeation of natural compounds Withaferin A (Wi-A), Withanone (Wi-N), Caffeic Acid Phenethyl Ester (CAPE) and Artepillin C (ARC). Molecular dynamics simulations were performed to compute permeability coefficients. The results indicated that the exposure of PS in cancer cell membranes facilitated the permeation of Wi-A, Wi-N and CAPE through a cancer cell membrane when compared to a normal cell membrane. In the case of ARC, PS exposure did not have a notable influence on its permeability coefficient. The presented data demonstrated the potential of PS exposure-based models for studying cancer cell selectivity of drugs.

Model Membrane Systems Used to Study Plasma Membrane Lipid Asymmetry

It is well known that the lipid distribution in the bilayer leaflets of mammalian plasma membranes (PMs) is not symmetric. Despite this, model membrane studies have largely relied on chemically symmetric model membranes for the study of lipid–lipid and lipid–protein interactions. This is primarily due to the difficulty in preparing stable, asymmetric model membranes that are amenable to biophysical studies. However, in the last 20 years, efforts have been made in producing more biologically faithful model membranes. Here, we review several recently developed experimental and computational techniques for the robust generation of asymmetric model membranes and highlight a new and particularly promising technique to study membrane asymmetry.

How to control interactions of cellulose-based biomaterials with skin: the role of acidity in the contact area

Being able to control the interactions of biomaterials with living tissues and skin is highly desirable for many biomedical applications. This is particularly the case for cellulose-based materials which provide one of the most versatile platforms for tissue engineering due to their strength, biocompatibility and abundance. Achieving such control, however, requires detailed molecular-level knowledge of the dominant interaction mechanisms. Here, we employed both biased and unbiased atomic-scale molecular dynamics simulations to explore how cellulose crystals interact with model stratum corneum bilayers, ternary mixtures of ceramides, cholesterol, and free fatty acids. Our findings show that acidity in the contact area directly affects binding between cellulose and the stratum corneum bilayer: Protonation of free fatty acids in the bilayer promotes attractive cellulose-bilayer interactions. We identified two major factors that control the cellulose-skin interactions: (i) the electrostatic repulsion between a cellulose crystal and the charged (anionic due to deprotonated fatty acids) surface of a stratum corneum bilayer and (ii) the cellulose-stratum corneum hydrogen bonding. When less than half of the fatty acids in the bilayer are protonated, the first factor dominates and there is no binding to skin. At a larger degree of fatty acid protonation the cellulose-stratum corneum hydrogen bonding prevails yielding a tight binding. Remarkably, we found that ceramide molecules are the key component in hydrogen bonding with cellulose. Overall, our findings highlight the critical role of fatty acid protonation in biomaterial-stratum corneum interactions and can be used for optimizing the surface properties of cellulose-based materials aimed at biomedical applications such as wound dressings.

Transmembrane potential of physiologically relevant model membranes: Effects of membrane asymmetry

Transmembrane potential difference (V_m) plays important roles in regulating various biological processes. At the macro level, V_m can be experimentally measured or calculated using the Nernst or Goldman-Hodgkin-Katz equation. However, the atomic details responsible for its generation and impact on protein and lipid dynamics still need to be further elucidated. In this work, we performed a series of all-atom molecular dynamics (MD) simulations of symmetric model membranes of various lipid compositions and cation contents to evaluate the relationship between membrane asymmetry and V_m. Specifically, we studied the impact of the asymmetric distribution of POPS (1-palmitoyl-2-oleoyl-sn-glycero-3-phospho-l-serine), PIP2 (phosphatidylinositol 4,5-bisphosphate), as well as Na⁺ and K⁺ on V_m using atomically detailed MD simulations of symmetric model membranes. The results suggest that, for an asymmetric POPC-POPC/POPS bilayer in the presence of NaCl, the presence of the monovalent anionic lipid POPS in the inner leaflet polarizes the membrane (ΔV_m < 0). Intriguingly, replacing a third of the POPS lipids by the polyvalent anionic signaling lipid PIP2 counteracts this effect, resulting in a smaller negative membrane potential. We also found that replacing Na⁺ ions in the inner region by K⁺ depolarizes the membrane (ΔV_m > 0). These divergent effects arise from variations in the strength of cation-lipid interactions and are correlated with changes in lipid chain order and head-group orientation.

Influence of Membrane Phase on the Optical Properties of DPH

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

Multiscale Simulations of Biological Membranes: The Challenge To Understand Biological Phenomena in a Living Substance

Biological membranes are tricky to investigate. They are complex in terms of molecular composition and structure, functional over a wide range of time scales, and characterized by nonequilibrium conditions. Because of all of these features, simulations are a great technique to study biomembrane behavior. A significant part of the functional processes in biological membranes takes place at the molecular level; thus computer simulations are the method of choice to explore how their properties emerge from specific molecular features and how the interplay among the numerous molecules gives rise to function over spatial and time scales larger than the molecular ones. In this review, we focus on this broad theme. We discuss the current state-of-the-art of biomembrane simulations that, until now, have largely focused on a rather narrow picture of the complexity of the membranes. Given this, we also discuss the challenges that we should unravel in the foreseeable future. Numerous features such as the actin-cytoskeleton network, the glycocalyx network, and nonequilibrium transport under ATP-driven conditions have so far received very little attention; however, the potential of simulations to solve them would be exceptionally high. A major milestone for this research would be that one day we could say that computer simulations genuinely research biological membranes, not just lipid bilayers.

Rapid diffusion of cholesterol along polyunsaturated membranes via deep dives

Cholesterol regulates the function of membrane proteins either via direct or membrane-mediated effects. Therefore, its ready availability is crucial for many protein-governed cellular processes. Recent studies suggest that cholesterol can partition to the core of polyunsaturated membranes, where cholesterol binding sites of many membrane proteins are also located. This core region is characterized by a lower viscosity. Therefore, we hypothesized that cholesterol partitioning into the membrane interior increases the rate of its diffusion in polyunsaturated membrane environments. We studied the behavior of cholesterol in membranes with increasing level of lipid chain unsaturation using a combination of atomistic and coarse-grained molecular dynamics simulations. Our simulations suggest a strong correlation between entropy-driven enhanced cholesterol partitioning to the membrane core and its faster lateral diffusion, which indicates that the less viscous membrane core indeed provides an efficient means for cholesterol movement in polyunsaturated membrane environments.

Microscopic view of lipids and their diverse biological functions

Biological membranes and their diverse lipid constituents play key roles in a broad spectrum of cellular and physiological processes. Characterization of membrane-associated phenomena at a microscopic level is therefore essential to our fundamental understanding of such processes. Due to the semi-fluid and dynamic nature of lipid bilayers, and their complex compositions, detailed characterization of biological membranes at an atomic scale has been refractory to experimental approaches. Computational modeling and simulation offer a highly complementary toolset with sufficient spatial and temporal resolutions to fill this gap. Here, we review recent molecular dynamics studies focusing on the diversity of lipid composition of biological membranes, or aiming at the characterization of lipid-protein interaction, with the overall goal of dissecting how lipids impact biological roles of the cellular membranes.

Transverse distribution of plasma membrane bilayer cholesterol: Picking sides

The transverse asymmetry (sidedness) of phospholipids in plasma membrane bilayers is well characterized, distinctive, actively maintained and functionally important. In contrast, numerous studies using a variety of techniques have concluded that plasma membrane bilayer cholesterol is either mostly in the outer leaflet or the inner leaflet or is fairly evenly distributed. Sterols might simply partition according to their differing affinities for the asymmetrically disposed phospholipids, but some studies have proposed that it is actively transported to the outer leaflet. Other work suggests that the sterol is enriched in the inner leaflet, driven by either positive interactions with the phosphatidylethanolamine on that side or by its exclusion from the outer leaflet by the long chain sphingomyelin molecules therein. This uncertainty raises three questions: is plasma membrane cholesterol sidedness fixed in a given cell or cell type; is it generally the same among mammalian species; and does it serve specific physiological functions? This review grapples with these issues.

Interleaflet Coupling, Pinning, and Leaflet Asymmetry-Major Players in Plasma Membrane Nanodomain Formation

The plasma membrane has a highly asymmetric distribution of lipids and contains dynamic nanodomains many of which are liquid entities surrounded by a second, slightly different, liquid environment. Contributing to the dynamics is a continuous repartitioning of components between the two types of liquids and transient links between lipids and proteins, both to extracellular matrix and cytoplasmic components, that temporarily pin membrane constituents. This make plasma membrane nanodomains exceptionally challenging to study and much of what is known about membrane domains has been deduced from studies on model membranes at equilibrium. However, living cells are by definition not at equilibrium and lipids are distributed asymmetrically with inositol phospholipids, phosphatidylethanolamines and phosphatidylserines confined mostly to the inner leaflet and glyco- and sphingolipids to the outer leaflet. Moreover, each phospholipid group encompasses a wealth of species with different acyl chain combinations whose lateral distribution is heterogeneous. It is becoming increasingly clear that asymmetry and pinning play important roles in plasma membrane nanodomain formation and coupling between the two lipid monolayers. How asymmetry, pinning, and interdigitation contribute to the plasma membrane organization is only beginning to be unraveled and here we discuss their roles and interdependence.

Interleaflet Coupling, Pinning, and Leaflet Asymmetry-Major Players in Plasma Membrane Nanodomain Formation

The plasma membrane has a highly asymmetric distribution of lipids and contains dynamic nanodomains many of which are liquid entities surrounded by a second, slightly different, liquid environment. Contributing to the dynamics is a continuous repartitioning of components between the two types of liquids and transient links between lipids and proteins, both to extracellular matrix and cytoplasmic components, that temporarily pin membrane constituents. This make plasma membrane nanodomains exceptionally challenging to study and much of what is known about membrane domains has been deduced from studies on model membranes at equilibrium. However, living cells are by definition not at equilibrium and lipids are distributed asymmetrically with inositol phospholipids, phosphatidylethanolamines and phosphatidylserines confined mostly to the inner leaflet and glyco- and sphingolipids to the outer leaflet. Moreover, each phospholipid group encompasses a wealth of species with different acyl chain combinations whose lateral distribution is heterogeneous. It is becoming increasingly clear that asymmetry and pinning play important roles in plasma membrane nanodomain formation and coupling between the two lipid monolayers. How asymmetry, pinning, and interdigitation contribute to the plasma membrane organization is only beginning to be unraveled and here we discuss their roles and interdependence.