Chloroform alters interleaflet coupling in lipid bilayers: an entropic mechanism

Mechanistic Understanding from Molecular Dynamics in Pharmaceutical Research 2: Lipid Membrane in Drug Design

We review the use of molecular dynamics (MD) simulation as a drug design tool in the context of the role that the lipid membrane can play in drug action, i.e., the interaction between candidate drug molecules and lipid membranes. In the standard "lock and key" paradigm, only the interaction between the drug and a specific active site of a specific protein is considered; the environment in which the drug acts is, from a biophysical perspective, far more complex than this. The possible mechanisms though which a drug can be designed to tinker with physiological processes are significantly broader than merely fitting to a single active site of a single protein. In this paper, we focus on the role of the lipid membrane, arguably the most important element outside the proteins themselves, as a case study. We discuss work that has been carried out, using MD simulation, concerning the transfection of drugs through membranes that act as biological barriers in the path of the drugs, the behavior of drug molecules within membranes, how their collective behavior can affect the structure and properties of the membrane and, finally, the role lipid membranes, to which the vast majority of drug target proteins are associated, can play in mediating the interaction between drug and target protein. This review paper is the second in a two-part series covering MD simulation as a tool in pharmaceutical research; both are designed as pedagogical review papers aimed at both pharmaceutical scientists interested in exploring how the tool of MD simulation can be applied to their research and computational scientists interested in exploring the possibility of a pharmaceutical context for their research.

Cholesterol sequestration by xenon nano bubbles leads to lipid raft destabilization

Combined coarse-grained (CG) and atomistic molecular dynamics (MD) simulations were performed to study the interactions of xenon with model lipid rafts consisting of 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC), 1,2-dilauroyl-sn-glycero-3-phosphocholine (DLPC) and cholesterol (Chol). At a concentration of 2 Xe/lipid we observed an unexpected result: spontaneous nucleation of Xe nano bubbles which rapidly plunged into the bilayer. In this process Chol, essential for raft stabilization, was pulled out from the raft into the hydrophobic zone. When concentration was further increased (3 Xe/lipid), the bubbles increase in size and disrupted both the membrane and raft. We computed the radial distribution functions, pair-wise potentials, second virial coefficients and Schlitter entropy to scrutinize the nature of the interactions. Our findings, concurring with a recent report on the origin of general anaesthesia (M. A. Pavel, E. N. Petersen, H. Wang, R. A. Lerner and S. B. Hansen, Proc. Natl. Acad. Sci. U. S. A., 2020, 117(24), 13757-13766), suggest that the well-known anaesthetic effect of Xe could be mediated by sequestration of Chol, which, in turn, compromises the stability of rafts where specialized proteins needed to produce the nervous signal are anchored.

Interaction of Graphene Nanoparticles and Lipid Membranes Displaying Different Liquid Orderings: A Molecular Dynamics Study

Understanding the effects of graphene-based nanomaterials on lipid membranes is fundamental to determine their environmental impact and the efficiency of their biomedical use. By means of molecular dynamics simulations of simple model lipid bilayers, we analyze in detail the different interaction modes. We have studied bilayers consisting of lipid species (including cholesterol) which display different internal liquid orderings. Nanometric graphene layers can be transiently adsorbed onto the lipid membrane and/or inserted in its hydrophobic region. Once inserted, graphene nanometric flakes display a diffusive dynamics in the membrane plane, they adopt diverse orientations depending on their size and oxidation degree, and they show a particular aversion to be placed close to cholesterol molecules in the membrane. Addition of graphene to phase-segregated ternary membranes is also investigated in the context of the lipid raft model for the lipid organization of biological membranes. Our simulation results show that graphene layers can be inserted indistinctly in the ordered and disordered regions. Once inserted, nanometric flakes migrate to disordered and cholesterol-poor lipid phases.

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.

Effects of Passive Phospholipid Flip-Flop and Asymmetric External Fields on Bilayer Phase Equilibria

Compositional asymmetry between the leaflets of bilayer membranes modifies their phase behavior and is thought to influence other important features such as mechanical properties and protein activity. We address here how phase behavior is affected by passive phospholipid flip-flop, such that the compositional asymmetry is not fixed. We predict transitions from "pre-flip-flop" behavior to a restricted set of phase equilibria that can persist in the presence of passive flip-flop. Surprisingly, such states are not necessarily symmetric. We further account for external symmetry breaking, such as a preferential substrate interaction, and show how this can stabilize strongly asymmetric equilibrium states. Our theory explains several experimental observations of flip-flop-mediated changes in phase behavior and shows how domain formation and compositional asymmetry can be controlled in concert, by manipulating passive flip-flop rates and applying external fields.

Effects of fullerene on lipid bilayers displaying different liquid ordering: a coarse-grained molecular dynamics study

BACKGROUND

The toxic effects and environmental impact of nanomaterials, and in particular of Fullerene particles, are matters of serious concern. It has been reported that fullerene molecules enter the cell membrane and occupy its hydrophobic region. Understanding the effects of carbon-based nanoparticles on biological membranes is therefore of critical importance to determine their exposure risks.

METHODS

We report on a systematic coarse-grained molecular dynamics study of the interaction of fullerene molecules with simple model cell membranes. We have analyzed bilayers consisting of lipid species with different degrees of unsaturation and a variety of cholesterol fractions. Addition of fullerene particles to phase-segregated ternary membranes is also investigated in the context of the lipid raft model for the organization of the cell membrane.

RESULTS

Fullerene addition to lipid membranes modifies their structural properties like thickness, area and internal ordering of the lipid species, as well as dynamical aspects such as molecular diffusion and cholesterol flip-flop. Interestingly, we show that phase-segregating ternary lipid membranes accumulate fullerene molecules preferentially in the liquid-disordered domains promoting phase-segregation and domain alignment across the membrane.

CONCLUSIONS

Lipid membrane internal ordering determines the behavior and distribution of fullerene particle, and this, in turn, determines the influence of fullerene on the membrane. Lipid membranes are good solvents of fullerene molecules, and in particular those with low internal ordering.

GENERAL SIGNIFICANCE

Preference of fullerene molecules to be dissolved in the more disordered hydrophobic regions of a lipid bilayer and the consequent alteration of its phase behavior may have important consequences on the activity of biological cell membranes and on the bioconcentration of fullerene in living organisms.

Alteration of interleaflet coupling due to compounds displaying rapid translocation in lipid membranes

The spatial coincidence of lipid domains at both layers of the cell membrane is expected to play an important role in many cellular functions. Competition between the surface interleaflet tension and a line hydrophobic mismatch penalty are conjectured to determine the transversal behavior of laterally heterogeneous lipid membranes. Here, by a combination of molecular dynamics simulations, a continuum field theory and kinetic equations, I demonstrate that the presence of small, rapidly translocating molecules residing in the lipid bilayer may alter its transversal behavior by favoring the spatial coincidence of similar lipid phases.

Roles of Interleaflet Coupling and Hydrophobic Mismatch in Lipid Membrane Phase-Separation Kinetics

Characterizing the nanoscale dynamic organization within lipid bilayer membranes is central to our understanding of cell membranes at a molecular level. We investigate phase separation and communication across leaflets in ternary lipid bilayers, including saturated lipids with between 12 and 20 carbons per tail. Coarse-grained molecular dynamics simulations reveal a novel two-step kinetics due to hydrophobic mismatch, in which the initial response of the apposed leaflets upon quenching is to increase local asymmetry (antiregistration), followed by dominance of symmetry (registration) as the bilayer equilibrates. Antiregistration can become thermodynamically preferred if domain size is restricted below ∼20 nm, with implications for the symmetry of rafts and nanoclusters in cell membranes, which have similar reported sizes. We relate our findings to theory derived from a semimicroscopic model in which the leaflets experience a “direct” area-dependent coupling, and an “indirect” coupling that arises from hydrophobic mismatch and is most important at domain boundaries. Registered phases differ in composition from antiregistered phases, consistent with a direct coupling between the leaflets. Increased hydrophobic mismatch purifies the phases, suggesting that it contributes to the molecule-level lipid immiscibility. Our results demonstrate an interplay of competing interleaflet couplings that affect phase compositions and kinetics, and lead to a length scale that can influence lateral and transverse bilayer organization within cells.

Nucleation of symmetric domains in the coupled leaflets of a bilayer

We study the kinetics governing the attainment of inter-leaflet domain symmetry in a phase-separating amphiphilic bilayer. "Indirect" inter-leaflet coupling via hydrophobic mismatch can induce an instability towards a metastable pattern of locally asymmetric domains upon quenching from high temperature. This necessitates a nucleation step to form the conventional symmetric domains, which are favoured by a "direct" inter-leaflet coupling. We model the energetics for a symmetric domain to nucleate from the metastable state, and find that an interplay between hydrophobic mismatch and thickness stretching/compression causes the effective hydrophobic mismatch, and thus line tension, to depend on domain size. This leads to strong departure from classical nucleation theory. We speculate on implications for cell membrane rafts or clusters, whose size may be of similar magnitude to estimated critical radii for domain symmetry.

Kinetics of symmetry and asymmetry in a phase-separating bilayer membrane

We simulate a phase-separating bilayer in which the leaflets experience a direct coupling favoring local compositional symmetry ("registered" bilayer phases), and an indirect coupling due to hydrophobic mismatch that favors strong local asymmetry ("antiregistered" bilayer phases). For wide ranges of overall leaflet compositions, multiple competing states are possible. For estimated physical parameters, a quenched bilayer may first evolve toward a metastable state more asymmetric than if the leaflets were uncorrelated; subsequently, it must nucleate to reach its equilibrium, more symmetric, state. These phase-transition kinetics exhibit characteristic signatures through which fundamental and opposing interleaflet interactions may be probed. We emphasize how bilayer phase diagrams with a separate axis for each leaflet can account for overall and local symmetry or asymmetry, and capture a range of observations in the experiment and simulation literature.