DPPC-cholesterol phase diagram using coarse-grained Molecular Dynamics simulations

Effect of cholesterol on nanoparticle translocation across a lipid bilayer

Nanoparticles (NPs) have attracted significant attention as carriers for the delivery of drugs, genes, and macromolecules for biomedical and therapeutic applications. These technologies require NPs to be delivered to the interior of the cell. However, this translocation is unlikely because of the presence of a cell membrane composed of phospholipids, cholesterol, proteins, and glycans. The cell membrane composition can influence its rigidity; thus, membrane composition is a crucial factor in determining the translocation of NPs across the cell membrane. Here, we focus on cholesterol, which is an essential component of biological cell membranes, and investigate NP translocation across membranes containing cholesterol under an applied electric field using a coarse-grained molecular dynamics simulation. We found that NPs could translocate directly across cholesterol-containing membranes without irreversible membrane disruption. This unique translocation was induced by two key phenomena. Before NP translocation, a phospholipid-rich/cholesterol-poor domain was formed at the NP-membrane contact interface. Second, a smaller transmembrane pore was formed in the cholesterol-containing membrane during membrane crossing of the NP. Our findings imply that the delivery of NPs to the cell interior across the cholesterol-containing membrane can be achieved by appropriately controlling the strength of the applied electric field, depending on the cholesterol content in the membrane.

Influence of the Water Model on the Structure and Interactions of the GPR40 Protein with the Lipid Membrane and the Solvent: Rigid versus Flexible Water Models

G protein-coupled receptors (GPCR) are responsible for modulating various physiological functions and are thus related to the pathophysiology of different diseases. Being potential therapeutic targets, multiple computational methodologies have been developed to analyze their behavior and interactions with other species. The solvent, on the other hand, has received much less attention. In this work, we analyzed the effect of four explicit water models on the structure and interactions of the GPR40 receptor in its apo form. We employed the rigid SPC/E and TIP4P models, and their flexible versions, the FBA/ϵ and TIP4P/ϵflex. We explored the structural changes and their correlation with some bulk dynamic properties of water. Our results showed an adverse effect on the conservation of the secondary structure of the receptor with all the models due to the breaking of the intramolecular hydrogen bond network, being more evident for the TIP4P models. Notably, all four models brought the receptor to states similar to the active one, modifying the intracellular part of the TM5 and TM6 domains in a "hinge" type movement, allowing the opening of the structure. Regarding the dynamic properties, the rigid models showed results comparable to those obtained in other studies on membrane systems. However, flexible models exhibit disparities in the molecular representation of systems. Surprisingly, the FBA/ϵ model improves the molecular picture of several properties, even though their agreement with bulk diffusion is poorer. These findings reinforce our idea that exploring other water models or improving the current ones, to better represent the membrane interface, can lead to a positive impact on the description of the signal transduction mechanisms and the search of new drugs by targeting these receptors.

Manufacturing process of liposomal Formation: A coarse-grained molecular dynamics simulation

A method of producing liposomes has been previously developed using a continuous manufacturing technology that involves a co-axial turbulent jet in co-flow. In this study, coarse-grained molecular dynamics (CG-MD) simulations were used to gain a deeper understanding of how the self-assembly process of liposomes is affected by the material attributes (such as the concentration of ethanol) and the process parameters (such as temperature), while also providing detailed information on a nano-scale molecular level. Specifically, the CG-MD simulations yield a comprehensive internal view of the structure and formation mechanisms of liposomes containing DPPC, DPPG, and cholesterol molecules. The importance of this work is that structural details on the molecular level are proposed, and such detail is not possible to obtain through experimental studies alone. The assessment of structural properties, including the area per lipid, diffusion coefficient, and order parameters, indicated that a thicker bilayer was observed at higher ethanol concentrations, while a thinner bilayer was present at higher temperatures. These conditions led to more water penetrating the interior of the bilayer and an unstable structure, as indicated by a larger contact area between lipids and water, and a higher coefficient of lipid lateral diffusion. However, stable liposomes were found through these evaluations at lower ethanol concentrations and/or lower process temperatures. Furthermore, the CG-MD model was further compared and validated with experimental and computational data including liposomal bilayer thickness and area per lipid measurements.

Martini 3 Force Field Parameters for Protein Lipidation Post-Translational Modifications

Protein lipidations are vital co/post-translational modifications that tether lipid tails to specific protein amino acids, allowing them to anchor to biological membranes, switch their subcellular localization, and modulate association with other proteins. Such lipidations are thus crucial for multiple biological processes including signal transduction, protein trafficking, and membrane localization and are implicated in various diseases as well. Examples of lipid-anchored proteins include the Ras family of proteins that undergo farnesylation; actin and gelsolin that are myristoylated; phospholipase D that is palmitoylated; glycosylphosphatidylinositol-anchored proteins; and others. Here, we develop parameters for cysteine-targeting farnesylation, geranylgeranylation, and palmitoylation, as well as glycine-targeting myristoylation for the latest version of the Martini 3 coarse-grained force field. The parameters are developed using the CHARMM36m all-atom force field parameters as reference. The behavior of the coarse-grained models is consistent with that of the all-atom force field for all lipidations and reproduces key dynamical and structural features of lipid-anchored peptides, such as the solvent-accessible surface area, bilayer penetration depth, and representative conformations of the anchors. The parameters are also validated in simulations of the lipid-anchored peripheral membrane proteins Rheb and Arf1, after comparison with independent all-atom simulations. The parameters, along with mapping schemes for the popular martinize2 tool, are available for download at 10.5281/zenodo.7849262 and also as supporting information.

Anomalous lateral diffusion of lipids during the fluid/gel phase transition of a lipid membrane

A lipid membrane undergoes a phase transition from fluid to gel phase upon changing external thermodynamic conditions, such as decreasing temperature and increasing pressure. Extremophilic organisms face the challenge of preventing this deleterious phase transition. The main focus of their adaptive strategy is to facilitate effective temperature sensing through sensor proteins, relying on the drastic changes in packing density and membrane fluidity during the phase transition. Although the changes in packing density parameters due to the fluid/gel phase transition are studied in detail, the impact on membrane fluidity is less explored in the literature. Understanding the lateral diffusive dynamics of lipids in response to temperature, particularly during the fluid/gel phase transition, is albeit crucial. Here we have simulated the phase transition of a single component lipid membrane composed of dipalmitoylphosphatidylcholine (DPPC) lipids using a coarse-grained (CG) model and studied the changes of the structural and dynamical properties. It is observed that near the phase transition point, both fluid and gel phase domains coexist together. The dynamics remains highly non-Gaussian for a long time even when the mean square displacement reaches the Fickian regime at a much earlier time. This Fickian yet non-Gaussian diffusion (FnGD) is a characteristic of a highly heterogeneous system, previously observed for the lateral diffusion of lipids in raft mimetic membranes having liquid-ordered and liquid-disordered phases co-existing together. We have analyzed the molecular trajectories and calculated the jump-diffusion of the lipids, stemming from sudden jump translations, using a translational jump-diffusion (TJD) approach. An overwhelming contribution of the jump-diffusion of the lipids is observed suggesting anomalous diffusion of lipids during fluid/gel phase transition of the membrane. These results are important in unravelling the intricate nature of lipid diffusion during the phase transition of the membrane and open up a new possibility of investigating the most significant change of membrane properties during phase transition, which can be effectively sensed by proteins.

Capric Acid and Myristic Acid Permeability Enhancers in Curved Liposome Membranes

Liposomes are considered as advanced drug delivery systems for cancer treatment. A generation of pH-sensitive liposomes is being developed that use fatty acids (FAs) as a trigger for drug release in tumor tissues. However, FAs are also known to enhance permeability, and it is unclear whether FAs in liposomes may cause drug leakage or premature drug release. The passive permeability of the drug through the membrane of the liposome is thus a crucial factor for timely drug delivery. To investigate how the curvature and lipid composition of liposomes affect their passive permeability, coarse-grained molecular dynamics were performed. The permeability was determined with a counting method. Flat bilayers and three liposomes with varying diameters were studied, which had varying lipid compositions of dipalmitoylphosphatidylcholine, cholesterol, and deprotonated or neutral saturated FAs. The investigated permeants were water and two other small permeants, which have different free energy profiles (solubility) across the membrane. First, for the curvature effect, our results showed that curvature increases the water permeability by reducing the membrane thickness. The permeability increase for water is about a factor of 1.7 for the most curved membranes. However, a high curvature decreases permeability for permeants with free energy profiles that are a mix of wells and barriers in the headgroup region of the membrane. Importantly, the type of experimental setup is expected to play a dominant role in the permeability value, i.e., whether permeants are escaping or entering the liposomes. Second, for the composition effect, FAs decrease both the area per lipid (APL) and the membrane thickness, resulting in permeability increases of up to 55%. Cholesterol has a similar effect on the APL but has the opposite impact on membrane thickness and permeability. Therefore, FAs and cholesterol have opposing effects on permeability, with cholesterol's effect being slightly stronger in our simulated bilayers. As all permeability values were well within a factor of 2, and with liposomes usually being larger and less curved in experimental applications, it can be concluded that the passive drug release from a pH-sensitive liposome does not seem to be significantly affected by the presence of FAs.

Understanding the Mechanical Properties of Ultradeformable Liposomes Using Molecular Dynamics Simulations

Improving drug delivery efficiency to solid tumor sites is a central challenge in anticancer therapeutic research. Our previous experimental study (Guo et al., Nat. Commun. 2018, 9, 130) showed that soft, elastic liposomes had increased uptake and accumulation in cancer cells and tumors in vitro and in vivo respectively, relative to rigid particles. As a first step toward understanding how liposomes' molecular structure and composition modulates their elasticity, we performed all-atom and coarse-grained classical molecular dynamics (MD) simulations of lipid bilayers formed by mixing a long-tailed unsaturated phospholipid with a short-tailed saturated lipid with the same headgroup. The former types of phospholipids considered were 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC) and 1,2-dipalmitoleoyl-sn-glycero-3-phosphocholine (termed here DPMPC). The shorter saturated lipids examined were 1,2-diheptanoyl-sn-glycero-3-phosphocholine (DHPC), 1,2-didecanoyl-sn-glycero-3-phosphocholine (DDPC), 1,2-dilauroyl-sn-glycero-3-phosphocholine (DLPC), and 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC). Several lipid concentrations and surface tensions were considered. Our results show that DOPC or DPMPC systems having 25-35 mol % of the shortest lipids DHPC or DDPC are the least rigid, having area compressibility moduli KA that are ∼10% smaller than the values observed in pure DOPC or DPMPC bilayers. These results agree with experimental measurements of the stretching modulus and lysis tension in liposomes with the same compositions. These mixed systems also have lower areas per lipid and form more uneven x-y interfaces with water, the tails of both primary and secondary lipids are more disordered, and the terminal methyl groups in the tails of the long lipid DOPC or DPMPC wriggle more in the vertical direction, compared to pure DOPC or DPMPC bilayers or their mixtures with the longer saturated lipid DLPC or DMPC. These observations confirm our hypothesis that adding increasing concentrations of the short unsaturated lipid DHPC or DDPC to DOPC or DPMPC bilayers alters lipid packing and thus makes the resulting liposomes more elastic and less rigid. No formation of lipid nanodomains was noted in our simulations, and no clear trends were observed in the lateral diffusivities of the lipids as the concentration, type of secondary lipid, and surface tension were varied.

Membrane homeostasis beyond fluidity: control of membrane compressibility

Biomembranes are complex materials composed of lipids and proteins that compartmentalize biochemistry. They are actively remodeled in response to physical and metabolic cues, as well as during cell differentiation and stress. The concept of homeoviscous adaptation has become a textbook example of membrane responsiveness. Here, we discuss limitations and common misconceptions revolving around it. By highlighting key moments in the life cycle of a transmembrane protein, we illustrate that membrane thickness and a finely regulated membrane compressibility are crucial to facilitate proper membrane protein insertion, function, sorting, and inheritance. We propose that the unfolded protein response (UPR) provides a mechanism for endoplasmic reticulum (ER) membrane homeostasis by sensing aberrant transverse membrane stiffening and triggering adaptive responses that re-establish membrane compressibility.

Thylakoid Composition Facilitates Chlorophyll a Dimerization through Stronger Interlipid Interactions

Plant thylakoid membrane serves as a crucial matrix for the aggregation of chlororophyll a (CLA) pigments, essential for light harvesting. To understand the role of lipid compositions in the stability of CLA aggregates, dimerization of chlorophyll a molecules (CLA) is studied in the presence of the thylakoid and the bilayers comprising either the least or the highest unsaturated lipids by using coarse-grained molecular dynamics simulations. The thylakoid membrane enhances the stability of the CLA dimer compared with other membranes due to very strong lipid-lipid interactions. The thylakoid exhibits a distinct distribution of lipids around the CLA dimer. Less unsaturated lipids reside in close proximity to the dimer, promoting increased order and efficient packing. Conversely, higher unsaturated lipids are depleted from the dimer, imparting flexibility to the membrane. The combination of tight packing near the dimer and membrane flexibility away from the dimer enhances the stability of the dimer in the thylakoid membrane. Our results suggest that lipid mixing, rather than lipid unsaturation, plays a critical role in facilitating CLA dimerization by modulating the membrane microenvironment through stronger lipid-lipid interactions. These insights will be useful in understanding how lipid compositions affect efficient light absorption and energy transfer during photosynthesis in the future.

Martini 3 Coarse-Grained Force Field for Cholesterol

Cholesterol plays a crucial role in biomembranes by regulating various properties, such as fluidity, rigidity, permeability, and organization of lipid bilayers. The latest version of the Martini model, Martini 3, offers significant improvements in interaction balance, molecular packing, and inclusion of new bead types and sizes. However, the release of the new model resulted in the need to reparameterize many core molecules, including cholesterol. Here, we describe the development and validation of a Martini 3 cholesterol model, addressing issues related to its bonded setup, shape, volume, and hydrophobicity. The proposed model mitigates some limitations of its Martini 2 predecessor while maintaining or improving the overall behavior.

Box-Behnken design of thermo-responsive nano-liposomes loaded with a platinum(iv) anticancer complex: evaluation of cytotoxicity and apoptotic pathways in triple negative breast cancer cells

Herein, thermo-responsive liposomes (TLs) loaded with Asp (Asp/TLs) were produced by self-assembling DPPC, DSPE-PEG2000, and cholesterol. The preparation variables were optimized using the Box-Behnken design (BBD). The optimized Asp/TLs exhibited an average particle size of 114.05 ± 1.56 nm, PDI of 0.15 ± 0.015, zeta potential of -15.24 ± 0.65 mV, and entrapment efficiency (EE%) of 84.08 ± 2.75%. In addition, under physiological conditions, Asp/TLs showed spherical shape, outstanding stability and thermo-triggered the release of Asp at 38 °C, reaching the maximum Asp release at 40 °C. The MTT assay showed that the optimal Asp/TLs exhibited the highest cytotoxic activity upon exposure to mild hyperthermia (40 °C) against the invasive triple-negative breast cancer cell line (MDA-MB-231) when compared to other preparations. The IC50 of Asp/TLs (40 °C) was estimated at 0.9 μg mL-1, while that of free Asp (40 °C) was 3.83 μg mL-1. As such, the optimal Asp/TLs were shown to increase the cytotoxic activity of Asp by 4-fold upon exposure to mild hyperthermia. The IC50 values of Asp and Asp/TLs without exposure to 40 °C were 6.6 μg mL-1 and 186 μg mL-1, respectively. This indicated that Asp was released only when placed at 40 °C. The apoptosis assay revealed that Asp/TLs (40 °C) caused a remarkable increase in the percentage of cell population among both the late apoptosis and necrosis quartiles, as well as a significant decline in the viable cell quartile (P ≤ 0.001) when compared to Asp (40 °C). Asp/TLs (40 °C) and Asp (40 °C) could stimulate the intrinsic apoptosis pathway by upregulating the apoptotic genes Bak and Bax, while downregulating the anti-apoptotic genes, BCL-xL and BCL-2. The free Asp (40 °C) increased the gene expression of Bak and Bax by 4.4- and 5.2-folds, while reducing the expression of BCL-xL and BCL-2 by 50% and 73%, respectively. The optimal Asp TLs (40 °C) manifested more potent effects as demonstrated by the upregulation of Bak, Bax, and P53 by 5.6-, 7.2-, and 1.3-folds, as well as the downregulation of BCL-xL and BCL-2 by 70% and 85%, respectively. As such, the optimal Asp TLs (40 °C) treatment displayed the most potent cytotoxic profile and induced both apoptosis and necrosis in MDA-MB-231.

Unraveling the Phase Behavior, Mechanical Stability, and Protein Reconstitution Properties of Polymer-Lipid Hybrid Vesicles

Hybrid vesicles consisting of natural phospholipids and synthetic amphiphilic copolymers have shown remarkable material properties and potential for biotechnology, combining the robustness of polymers with the biocompatibility of phospholipid membranes. To predict and optimize the mixing behavior of lipids and copolymers, as well as understand the interaction between the hybrid membrane and macromolecules like membrane proteins, a comprehensive understanding at the molecular level is essential. This can be achieved by a combination of molecular dynamics simulations and experiments. Here, simulations of POPC and PBD22-b-PEO14 hybrid membranes are shown, uncovering different copolymer configurations depending on the polymer-to-lipid ratio. High polymer concentrations created thicker membranes with an extended polymer conformation, while high lipid content led to the collapse of the polymer chain. High concentrations of polymer were further correlated with a decreased area compression modulus and altered lateral pressure profiles, hypothesized to result in the experimentally observed improvement in membrane protein reconstitution and resistance toward destabilization by detergents. Finally, simulations of a WALP peptide embedded in the bilayer showed that only membranes with up to 50% polymer content favored a transmembrane configuration. These simulations correlate with previous and new experimental results and provide a deeper understanding of the properties of lipid-copolymer hybrid membranes.

Viscosity-Sensitive Membrane Dyes as Tools To Estimate the Crystalline Structure of Lipid Bilayers

Lipid membranes are crucial for cellular integrity and regulation, and tight control of their structural and mechanical properties is vital to ensure that they function properly. Fluorescent probes sensitive to the membrane's microenvironment are useful for investigating lipid membrane properties; however, there is currently a lack of quantitative correlation between the exact parameters of lipid organization and a readout from these dyes. Here, we investigate this relationship for "molecular rotors", or microviscosity sensors, by simultaneously measuring their fluorescence lifetime to determine the membrane viscosity, while using X-ray diffraction to determine the membrane's structural properties. Our results reveal a phase-dependent correlation between the membrane's structural parameters and mechanical properties measured by a BODIPY-based molecular rotor, giving excellent predictive power for the structural descriptors of the lipid bilayer. We also demonstrate that differences in membrane thickness between different lipid phases are not a prerequisite for the formation of lipid microdomains and that this requirement can be disrupted by the presence of line-active molecules. Our results underpin the use of membrane-sensitive dyes as reporters of the structure of lipid membranes.

Interdependence of cholesterol distribution and conformational order in lipid bilayers

We show, via molecular simulations, that not only does cholesterol induce a lipid order, but the lipid order also enhances cholesterol localization within the lipid leaflets. Therefore, there is a strong interdependence between these two phenomena. In the ordered phase, cholesterol molecules are predominantly present in the bilayer leaflets and orient themselves parallel to the bilayer normal. In the disordered phase, cholesterol molecules are mainly present near the center of the bilayer at the midplane region and are oriented orthogonal to the bilayer normal. At the melting temperature of the lipid bilayers, cholesterol concentration in the leaflets and the bilayer midplane is equal. This result suggests that the localization of cholesterol in the lipid bilayers is mainly dictated by the degree of ordering of the lipid bilayer. We validate our findings on 18 different lipid bilayer systems, obtained from three different phospholipid bilayers with varying concentrations of cholesterol. To cover a large temperature range in simulations, we employ the Dry Martini force field. We demonstrate that the Dry and the Wet Martini (with polarizable water) force fields produce comparable results.

CHARMM-GUI Membrane Builder: Past, Current, and Future Developments and Applications

Molecular dynamics simulations of membranes and membrane proteins serve as computational microscopes, revealing coordinated events at the membrane interface. As G protein-coupled receptors, ion channels, transporters, and membrane-bound enzymes are important drug targets, understanding their drug binding and action mechanisms in a realistic membrane becomes critical. Advances in materials science and physical chemistry further demand an atomistic understanding of lipid domains and interactions between materials and membranes. Despite a wide range of membrane simulation studies, generating a complex membrane assembly remains challenging. Here, we review the capability of CHARMM-GUI Membrane Builder in the context of emerging research demands, as well as the application examples from the CHARMM-GUI user community, including membrane biophysics, membrane protein drug-binding and dynamics, protein-lipid interactions, and nano-bio interface. We also provide our perspective on future Membrane Builder development.

The impact of lipid polyunsaturation on the physical and mechanical properties of lipid membranes

The lipid composition of cellular membranes and the balance between the different lipid components can be impacted by aging, certain pathologies, specific diets and other factors. This is the case in a subgroup of individuals with psychiatric disorders, such as schizophrenia, where cell membranes of patients have been shown to be deprived in polyunsaturated fatty acids (PUFAs), not only in brain areas where the target receptors are expressed but also in peripheral tissues. This PUFA deprivation thus represents a biomarker of such disorders that might impact not only the interaction of antipsychotic medications with these membranes but also the activation and signaling of the targeted receptors embedded in the lipid membrane. Therefore, it is crucial to understand how PUFAs levels alterations modulate the different physical properties of membranes. In this paper, several biophysical approaches were combined (Laurdan fluorescence spectroscopy, atomic force microscopy, differential scanning calorimetry, molecular modeling) to characterize membrane properties such as fluidity, elasticity and thickness in PUFA-enriched cell membranes and lipid model systems reflecting the PUFA imbalance observed in some diseases. The impact of both the number of unsaturations and their position along the chain on the above properties was investigated. Briefly, data revealed that PUFA presence in membranes increases membrane fluidity, elasticity and flexibility and decreases its thickness and order parameter. Both the level of unsaturation and their position affect these membrane properties.

The effects of serum albumin pre-adsorption of nanoparticles on protein corona and membrane interaction: A molecular simulation study

As a platform to deliver imaging and therapeutic agents to targeted sites in vivo, nanoparticles (NPs) have widespread applications in diagnosis and treatment of cancer. However, the poor in vivo delivery efficiency of nanoparticles limits its potential for further application. Once enter the physiological environment, nanoparticles immediately interact with proteins and form protein corona, which changes the physicochemical properties of nanoparticle surface and further affects their transport. In this study, we performed molecular dynamics simulations to study the adsorption mechanism of nanoparticles with various surface modifications and different proteins (e.g., human serum albumin, complement protein C3b), and their interactions with cell membrane. The results show that protein human serum albumin prefers to interact with hydrophobic and positively charged nanoparticles, while the protein C3b prefers the hydrophobic and charged nanoparticles. The pre-adsorption of human serum albumin on the nanoparticle surface obviously decreases the interaction of nanoparticle with C3b. Furthermore, the high amount of protein pre-adsorption could decrease the probability of nanoparticle-membrane interaction. These results indicate that appropriate modification of nanoparticles with protein provides nanoparticles with better capability of targeting, which could be used to guide nanoparticle design and improve transport efficiency.

Multiscale compression-induced restructuring of stacked lipid bilayers: From buckling delamination to molecular packing

Lipid membranes in nature adapt and reconfigure to changes in composition, temperature, humidity, and mechanics. For instance, the oscillating mechanical forces on lung cells and alveoli influence membrane synthesis and structure during breathing. However, despite advances in the understanding of lipid membrane phase behavior and mechanics of tissue, there is a critical knowledge gap regarding the response of lipid membranes to micromechanical forces. Most studies of lipid membrane mechanics use supported lipid bilayer systems missing the structural complexity of pulmonary lipids in alveolar membranes comprising multi-bilayer interconnected stacks. Here, we elucidate the collective response of the major component of pulmonary lipids to strain in the form of multi-bilayer stacks supported on flexible elastomer substrates. We utilize X-ray diffraction, scanning probe microscopy, confocal microscopy, and molecular dynamics simulation to show that lipid multilayered films both in gel and fluid states evolve structurally and mechanically in response to compression at multiple length scales. Specifically, compression leads to increased disorder of lipid alkyl chains comparable to the effect of cholesterol on gel phases as a direct result of the formation of nanoscale undulations in the lipid multilayers, also inducing buckling delamination and enhancing multi-bilayer alignment. We propose this cooperative short- and long-range reconfiguration of lipid multilayered films under compression constitutes a mechanism to accommodate stress and substrate topography. Our work raises fundamental insights regarding the adaptability of complex lipid membranes to mechanical stimuli. This is critical to several technologies requiring mechanically reconfigurable surfaces such as the development of electronic devices interfacing biological materials.

Molecular dynamics insights on temperature and pressure effects on electroporation

Electroporation is a cell-level phenomenon caused by an ionic imbalance in the membrane, being of great relevance in various fields of knowledge. A dependence of the pore formation kinetics on the environmental conditions (temperature and pressure) of the cell membrane has already been reported, but further clarification regarding how these variables affect the pore formation/resealing dynamics and the transport of molecules through the membrane is still lacking. The objective of the present study was to investigate the temperature (288-348 K) and pressure (1-5000 atm) effects on the electroporation kinetics using coarse-grained molecular dynamics simulations. Results shown that the time for pore formation and resealing increased with pressure and decreased with temperature, whereas the maximum pore radius increased with temperature and decreased with pressure. This behavior influenced the ion migration through the bilayer, and the higher ionic mobility was obtained in the 288 K/1000 atm simulations, i.e., a combination of low temperature and (not excessively) high pressure. These results were used to discuss some experimental observations regarding the extraction of intracellular compounds applying this technique. This study contributes to a better understanding of electroporation under different thermodynamic conditions and to an optimal selection of processing parameters in practical applications which exploit this phenomenon.

Formation of lipid raft nanodomains in homogeneous ternary lipid mixture of POPC/DPSM/cholesterol: Theoretical insights

Lipid rafts, in biological membranes, are cholesterol-rich nanodomains that regulate many protein activities and cellular processes. Understanding the formation of the lipid-raft nanodomains helps us elucidate many complex interactions in the cell. In this study, the formation of lipid-raft nanodomains in a ternary palmitoyl-oleoyl-phosphatidylcholine/stearoyl-sphingomyelin/cholesterol (POPC/DPSM/Chol) lipid mixture, the most realistic surrogate model for biological membranes, has been successfully observed for the first time in-silico using microsecond timescale molecular dynamics simulations. The model reveals the formation of cholesterol-induced nanodomains with raft-like characteristics and their underlying mechanism: the cholesterol molecules segregate themselves into cholesterol nanodomains and then enrich the cholesterol-rich domain with sphingomyelin molecules to form a lipid raft thanks to the weak bonding of cholesterol with sphingomyelin. Besides, it is found that the increase in cholesterol concentration enhances the biophysical properties (e.g., bilayer thickness, area per lipid headgroup, and order parameter) of the lipid raft nanodomains. Such findings suggest that the POPC/DPSM/Chol bilayer is a suitable model to fundamentally extend the nanodomain evolution to investigate their lifetime and kinetics as well as to study protein-lipid interaction, protein-protein interaction, and selection of therapeutic molecules in the presence of lipid rafts.

Homeoviscous Adaptation of the Lipid Membrane of a Soil Bacterium Surviving under Diurnal Temperature Variation: A Molecular Simulation Perspective

A recent experiment has reported the lipidome remodeling of a soil-based plant-associated bacterium Methylobacterium extorquens due to diurnal temperature variations. The key adaptation strategy is the headgroup-specific remodeling of the acyl chain. To understand the idiosyncratic adaptation at the molecular level, we simulate the model membrane of the same bacterium using the reported lipidome compositions at four different experimental temperatures. We investigate the temperature-dependent packing density and fluidity of the membrane, the constancy of which is key to the homeoviscous adaptation. The results show that complex lipidome remodeling approximately preserves membrane properties under heat and cold stress. The headgroup-specific remodeling of the acyl chain serves to fine-tune the packing density and fluidity of the membrane at different temperatures. While lipids with strongly interacting headgroups are more abundant at higher temperatures, the lipidome is more dominated by lipids with weaker interacting headgroups at lower temperatures. This adaptation alleviates lipid membrane disruption caused by heat and cold stress. This study provides a molecular picture of the homeoviscous adaptation of the realistic lipid membrane of a soil-based bacterium.

How Does Temperature Affect the Dynamics of SARS-CoV-2 M Proteins? Insights from Molecular Dynamics Simulations

Enveloped viruses, in general, have several transmembrane proteins and glycoproteins, which assist the virus in entry and attachment onto the host cells. These proteins also play a significant role in determining the shape and size of the newly formed virus particles. The lipid membrane and the embedded proteins affect each other in non-trivial ways during the course of the viral life cycle. Unraveling the nature of the protein-protein and protein-lipid interactions, under various environmental and physiological conditions, could therefore prove to be crucial in development of therapeutics. Here, we study the M protein of SARS-CoV-2 to understand the effect of temperature on the properties of the protein-membrane system. The membrane-embedded dimeric M proteins were studied using atomistic and coarse-grained molecular dynamics simulations at temperatures ranging between 10 and 50 °C. While temperature-induced fluctuations are expected to be monotonic, we observe a steady rise in the protein dynamics up to 40 °C, beyond which it surprisingly reverts back to the low-temperature behavior. Detailed investigation reveals disordering of the membrane lipids in the presence of the protein, which induces additional curvature around the transmembrane region. Coarse-grained simulations indicate temperature-dependent aggregation of M protein dimers. Our study clearly indicates that the dynamics of membrane lipids and integral M protein of SARS-CoV-2 enables it to better associate and aggregate only at a certain temperature range (i.e., ~ 30-40 °C). This can have important implications in the protein aggregation and subsequent viral budding/fission processes.

Curvature-based sorting of eight lipid types in asymmetric buckled plasma membrane models

Curvature is a fundamental property of biological membranes and has essential roles in cellular function. Bending of membranes can be induced by their lipid and protein compositions, as well as peripheral proteins, such as those that make up the cytoskeleton. An important aspect of membrane function is the grouping of lipid species into microdomains, or rafts, which serve as platforms for specific biochemical processes. The fluid mosaic model of membranes has evolved to recognize the importance of curvature and leaflet asymmetry, and there are efforts towards evaluating their functional roles. This work investigates the effect of curvature on the sorting of lipids in buckled asymmetric bilayers containing eight lipid types, approximating an average mammalian plasma membrane, through coarse-grained (CG) molecular dynamics (MD) simulations with the Martini force field. The simulations reveal that i) leaflet compositional asymmetry can induce curvature asymmetry, ii) lipids are sorted by curvature to different extents, and iii) curvature-based partitioning trends show moderate to strong correlations with lipid molecular volumes and head to tail bead ratios, respectively. The findings provide unique insights into the role of curvature in membrane organization, and the curvature-based sorting trends should be useful references for later investigations, and potentially interpreting the functional roles of specific lipids.

Effect of the cholesterol on electroporation of planar lipid bilayer

Electroporation threshold depends on the membrane composition, with cholesterol being one of its key components already studied in the past, but the results were inconclusive. The aim of our study was to determine behaviour of planar lipid bilayers with varying cholesterol concentrations under electric field. This would give us a better insight into cholesterol effect on membrane properties during electroporation process, since cholesterol is one of the major components of biological membranes and plays a crucial role in membrane organisation, dynamics, and function. Planar lipid bilayers were prepared from phosphatidylcholine lipids with 0, 20, 30, 50 and 80 mol% cholesterol. Capacitance was measured using the discharge method. Results show no statistical difference of cBLM between the cholesterol concentrations. Breakdown voltage Ubr of planar lipid bilayers was measured by means of linear rising voltage with seven different slopes. Obtained results were fitted to a strength-duration curve, where parameter Ubrmin represents minimal breakdown voltage, and parameter τRC represents the inclination of the strength-duration curve. Adding cholesterol to planar lipid bilayer gradually increased its Ubrmin until 50 mol% cholesterol concentration. Afterwards at 80 mol% Ubrmin does not further increase, in fact it reduces by 20% of the Ubrmin at 50 mol% cholesterol concentration.

Effect of Ionic Strength on Ibuprofenate Adsorption on a Lipid Bilayer of Dipalmitoylphosphatidylcholine from Molecular Dynamics Simulations

In this work, the free energy change in the process of transferring ibuprofenate from the bulk solution to the center of a model of the dipalmitoylphosphatidylcholine bilayer at different NaCl concentrations was calculated. Two minima were found in the free energy profile: a local minimum, located in the vicinity of the membrane, and the global free energy minimum, found near the headgroup region. The downward shift of free energy minima with increasing NaCl concentration is consistent with the results of previous works. Conversely, the upward shift of the free energy maximum with increasing ionic strength is due to the competition of sodium ions and lipids molecules to coordinate with ibuprofenate and neutralize its charge. In addition, normal molecular dynamics simulations were performed to study the effects of the ibuprofenate on the lipid bilayer and in the presence of a high ibuprofenate concentration. The effect of ionic strength on the properties of the lipid bilayer and on lipid-drug interactions was analyzed. The area per lipid shrinking with increasing ionic strength, volume of lipids, and thickness of the bilayer is consistent with the experimental results. At a very high ibuprofenate concentration, the lipid bilayer dehydrates, and it consequently transforms into the gel phase, thus blocking the permeation.

The Two Faces of the Liquid Ordered Phase

Coexisting liquid ordered (L_o) and liquid disordered (L_d) lipid phases in synthetic and plasma membrane-derived vesicles are commonly used to model the heterogeneity of biological membranes, including their putative ordered rafts. However, raft-associated proteins exclusively partition to the L_d and not the L_o phase in these model systems. We believe that the difference stems from the different microscopic structures of the lipid rafts at physiological temperature and the L_o phase studied at room temperature. To probe this structural diversity across temperatures, we performed atomistic molecular dynamics simulations, differential scanning calorimetry, and fluorescence spectroscopy on L_o phase membranes. Our results suggest that raft-associated proteins are excluded from the L_o phase at room temperature due to the presence of a stiff, hexagonally packed lipid structure. This structure melts upon heating, which could lead to the preferential solvation of proteins by order-preferring lipids. This structural transition is manifested as a subtle crossover in membrane properties; yet, both temperature regimes still fulfill the definition of the L_o phase. We postulate that in the compositionally complex plasma membrane and in vesicles derived therefrom, both molecular structures can be present depending on the local lipid composition. These structural differences must be taken into account when using synthetic or plasma membrane-derived vesicles as a model for cellular membrane heterogeneity below the physiological temperature.

Role of cholesterol flip-flop in oxidized lipid bilayers

We performed a series of molecular dynamics simulations of cholesterol (Chol) in nonoxidized 1-palmitoyl-2-linoleoyl-sn-glycero-3-phosphatidylcholine (PLPC) bilayer and in binary mixtures of PLPC-oxidized-lipid-bilayers with 0-50% Chol concentration and oxidized lipids with hydroperoxide and aldehyde oxidized functional groups. From the 60 unbiased molecular dynamics simulations (total of 161 μs), we found that Chol inhibited pore formation in the aldehyde-containing oxidized lipid bilayers at concentrations greater than 11%. For both pure PLPC bilayer and bilayers with hydroperoxide lipids, no pores were observed at any Chol concentration. Furthermore, increasing cholesterol concentration led to a change of phase state from the liquid-disordered to the liquid-ordered phase. This condensing effect of Chol was observed in all systems. Data analysis shows that the addition of Chol results in an increase in bilayer thickness. Interestingly, we observed Chol flip-flop only in the aldehyde-containing lipid bilayer but neither in the PLPC nor the hydroperoxide bilayers. Umbrella-sampling simulations were performed to calculate the translocation free energies and the Chol flip-flop rates. The results show that Chol's flip-flop rate depends on the lipid bilayer type, and the highest rate are found in aldehyde bilayers. As the main finding, we shown that Chol stabilizes the oxidized lipid bilayer by confining the distribution of the oxidized functional groups.

Insight into the Exosomal Membrane: From Viewpoints of Membrane Fluidity and Polarity

Numerous research studies have been done for exosomes, particularly focusing on membrane proteins and included nucleic acids, and the volume of the knowledge about the lipids in the exosomal membrane has been increasing. However, the dynamic property of the exosomal membrane is hardly studied. By employing milk exosome as an example, herein the exosomal membrane was characterized focusing on the membrane fluidity and polarity. The lipid composition and phase state of milk exosome (exosome from bovine milk) were estimated. The milk exosome contained enriched Chol (43.6 mol % in total lipid extracts), which made the membrane in the liquid-ordered (l_o) phase by interacting with phospholipids. To suggest a model of exosomal vesicle cargo, the liposome compositions that mimic milk exosome were studied: liposomes were made of cholesterol (Chol), milk sphingomyelin (milk SM), and 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC). By using fluorescent probes 1,6-diphenyl-1,3,5-hexatriene and 6-dodecanoyl-2-dimethylaminonaphthalene, the microenvironments of submicron-sized membranes of exosome and model liposomes were investigated. The membrane fluidity of milk exosome was slightly higher than those of Chol/milk SM/POPC liposomes with a similar content of Chol, suggesting the presence of enriched unsaturated lipids. The most purposeful membrane property was obtained by the liposome composition of Chol/milk SM/POPC = 40/15/45. From the above, it is concluded that Chol is a fundamental component of the milk exosomal membrane to construct the enriched l_o phase, which could increase the membrane rigidity and contribute to the function of exosome.

Unravelling the Effects of Cholesterol on the Second-Order Nonlinear Optical Responses of Di-8-ANEPPS Dye Embedded in Phosphatidylcholine Lipid Bilayers

Cholesterol is known for its role in maintaining the correct fluidity and rigidity of the animals cell membranes and thus their functions. Assessing the content and the role of cholesterol in lipid bilayers is therefore of crucial importance for a deeper understanding and control of membrane functioning. In this computational work, we investigate bilayers built from three types of glycerophospholipid phosphatidylcholine (PC) lipids, namely dipalmitoylphosphatidylcholine (DPPC), 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC), and dioleoylphosphatidylcholine (DOPC), and containing different amounts of cholesterol by analyzing the second-harmonic generation (SHG) nonlinear optical (NLO) response of a probe molecule, di-8-ANEPPS, inserted into the membranes. This molecular property presents the advantage to be specific to interfacial regions such as lipid bilayers. To unravel these effects, Molecular Dynamics (MD) simulations have been performed on both DPPC and DOPC lipids by varying the cholesterol mole fraction (from 0 to 0.66), while POPC was only considered as a pure bilayer. In the case of the structural properties of the bilayers, all the analyses converge toward the same conclusion: as the mole fraction of cholesterol increases, the systems become more rigid, confirming the condensing effect of cholesterol. In addition, the chromophore is progressively more aligned with respect to the normal to the bilayer. On the contrary, addition of unsaturation disorders the lipid bilayers, with barely no impact on the alignment of the chromophore. Then, using the frames obtained from the MD simulations, the first hyperpolarizability β of the dye in its environment has been computed at the TDDFT level. On the one hand, the addition of cholesterol induces a progressive increase of the diagonal component the β tensor parallel to the bilayer normal. On the other hand, larger β values have been calculated for the unsaturated than for the saturated lipid systems. In summary, this study illustrates the relationship between the composition and structure of the bilayers and the NLO responses of the embedded dye.

Molecular Mechanism of Ultrasound-Induced Structural Defects in Liposomes: A Nonequilibrium Molecular Dynamics Simulation Study

The use of ultrasound in combination with liposomes is a promising approach to improve drug delivery. To achieve an optimal drug release rate, it is important to understand how ultrasound induces pathways on the liposome surface where drugs can be released from the liposome. To this end, we carry out large-scale ultrasound-induced molecular dynamics simulations for three single lipid component liposomes formed from the commonly used phospholipids: 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC), 1,2-dipalmitoylphosphatidylcholine (DPPC), or phosphatidylcholine (POPC). The results show that ultrasound induces the detachment of two leaflets of the DOPC surface, suggesting that the drug release pathway may be through the low lipid packing areas on the stretched surface. In contrast, ultrasound induces pore formation on the surface of DPPC and DOPC, where drugs could escape from the liposomes. While the leaflet detachment and transient pore formation are the mechanisms of DOPC and DPPC, respectively, in both liquid-ordered and liquid-disordered phases, the leaflet detachment mechanism is switched to the transient pore formation mechanism on going from the liquid-ordered phase to the liquid-disordered phase in the POPC liposome. By adding 30% mol cholesterol, the leaflet detachment mechanism is observed in all liposomes. We found that the molecular origin that determines a mechanism is the competition between the intraleaflet and interleaflet interacting energy of lipids. The connection to experimental and theoretical modeling is discussed in some detail.

Phase Transitions and Structural Changes in DPPC Liposomes Induced by a 1-Carba-Alpha-Tocopherol Analogue

Steady-state emission spectroscopy of 1-anilino-8- naphthalene sulfonate (ANS) and 1,6-diphenyl-1,3,5-hexatriene (DPH), fluorescence anisotropy, and DSC methods were used to characterize the interactions of the newly synthesized 1-carba-alpha-tocopherol (CT) with a 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) membrane. The DSC results showed significant perturbations in the DPPC structure for CT concentrations as low as 2 mol%. The main phase transition peak was broadened and shifted to lower temperatures in a concentration-dependent manner, and pretransition was abolished. Increasing CT concentrations induced the formation of new phases in the DPPC structure, leading to melting at lower temperatures and, finally, disruption of the ordered DPPC structure. Hydration and structural changes of the DPPC liposomes using ANS and DPH fluorescent probes, which are selectively located at different places in the bilayer, were studied. With the increased concentration of CT molecules in the DPPC liposomes, structural changes with the simultaneous formation of different phases of such mixture were observed. Temperature studies of such mixtures revealed a decrease in the temperature of the main phase transition and fluidization at decreasing temperatures related to increasing hydration in the bilayer. Contour plots obtained from concentration-temperature data with fluorescent probes allowed for identification of different phases, such as gel, ordered liquid, disordered liquid, and liquid crystalline phases. The CT molecule with a modified chromanol ring embedded in the bilayer led to H-bonding interactions, expelling water molecules from the interphase, thus introducing disorder and structural changes to the highly ordered gel phase.

Molecular Dynamics of Janus Nanodimers Dispersed in Lamellar Phases of a Block Copolymer

We investigate structural and dynamical properties of Janus nanodimers (NDs) dispersed in lamellar phases of a diblock copolymer. By performing molecular dynamics simulations, we show that an accurate tuning of the interactions between NDs and copolymer blocks can lead to a close control of NDs' space distribution and orientation. In particular, NDs are preferentially found within the lamellae if enthalpy-driven forces offset their entropic counterpart. By contrast, when enthalpy-driven forces are not significant, the distribution of NDs, preferentially observed within the inter-lamellar spacing, is mostly driven by excluded-volume effects. Not only does the degree of affinity between host and guest species drive the NDs' distribution in the polymer matrix, but it also determines their space orientation. In turn, these key structural properties influence the long-time dynamics and the ability of NDs to diffuse through the polymer matrix.

Evaluating the Efficiency of the Martini Force Field to Study Protein Dimerization in Aqueous and Membrane Environments

Protein-protein complex assembly is one of the major drivers of biological response. Understanding the mechanisms of protein oligomerization/dimerization would allow one to elucidate how these complexes participate in biological activities and could ultimately lead to new approaches in designing novel therapeutic agents. However, determining the exact association pathways and structures of such complexes remains a challenge. Here, we use parallel tempering metadynamics simulations in the well-tempered ensemble to evaluate the performance of Martini 2.2P and Martini open-beta 3 (Martini 3) force fields in reproducing the structure and energetics of the dimerization process of membrane proteins and proteins in an aqueous solution in reasonable accuracy and throughput. We find that Martini 2.2P systematically overestimates the free energy of association by estimating large barriers in distinct areas, which likely leads to overaggregation when multiple monomers are present. In comparison, the less viscous Martini 3 results in a systematic underestimation of the free energy of association for proteins in solution, while it performs well in describing the association of membrane proteins. In all cases, the near-native dimer complexes are identified as minima in the free energy surface albeit not always as the lowest minima. In the case of Martini 3, we find that the spurious supramolecular protein aggregation present in Martini 2.2P multimer simulations is alleviated and thus this force field may be more suitable for the study of protein oligomerization. We propose that the use of enhanced sampling simulations with a refined coarse-grained force field and appropriately defined collective variables is a robust approach for studying the protein dimerization process, although one should be cautious of the ranking of energy minima.

The effect of polydispersity, shape fluctuations and curvature on small unilamellar vesicle small-angle X-ray scattering curves

It is a challenge to distinguish the effect of shape fluctuations and size polydispersity on experimental small-angle X-ray scattering curves of small unilamellar vesicles. Here it is shown that both effects have distinguishable spectral patterns, and an efficient simulation tool is presented for simulating and analysing experimental data. The importance of curvature-induced electron-density profile asymmetry for estimating the vesicle size from SAXS scattering curves is also demonstrated.

Small unilamellar vesicles (20–100 nm diameter) are model systems for strongly curved lipid membranes, in particular for cell organelles. Routinely, small-angle X-ray scattering (SAXS) is employed to study their size and electron-density profile (EDP). Current SAXS analysis of small unilamellar vesicles (SUVs) often employs a factorization into the structure factor (vesicle shape) and the form factor (lipid bilayer electron-density profile) and invokes additional idealizations: (i) an effective polydispersity distribution of vesicle radii, (ii) a spherical vesicle shape and (iii) an approximate account of membrane asymmetry, a feature particularly relevant for strongly curved membranes. These idealizations do not account for thermal shape fluctuations and also break down for strong salt- or protein-induced deformations, as well as vesicle adhesion and fusion, which complicate the analysis of the lipid bilayer structure. Presented here are simulations of SAXS curves of SUVs with experimentally relevant size, shape and EDPs of the curved bilayer, inferred from coarse-grained simulations and elasticity considerations, to quantify the effects of size polydispersity, thermal fluctuations of the SUV shape and membrane asymmetry. It is observed that the factorization approximation of the scattering intensity holds even for small vesicle radii (∼30 nm). However, the simulations show that, for very small vesicles, a curvature-induced asymmetry arises in the EDP, with sizeable effects on the SAXS curve. It is also demonstrated that thermal fluctuations in shape and the size polydispersity have distinguishable signatures in the SAXS intensity. Polydispersity gives rise to low-q features, whereas thermal fluctuations predominantly affect the scattering at larger q, related to membrane bending rigidity. Finally, it is shown that simulation of fluctuating vesicle ensembles can be used for analysis of experimental SAXS curves.

Microscopic Interactions of Melatonin, Serotonin and Tryptophan with Zwitterionic Phospholipid Membranes

The interactions at the atomic level between small molecules and the main components of cellular plasma membranes are crucial for elucidating the mechanisms allowing for the entrance of such small species inside the cell. We have performed molecular dynamics and metadynamics simulations of tryptophan, serotonin, and melatonin at the interface of zwitterionic phospholipid bilayers. In this work, we will review recent computer simulation developments and report microscopic properties, such as the area per lipid and thickness of the membranes, atomic radial distribution functions, angular orientations, and free energy landscapes of small molecule binding to the membrane. Cholesterol affects the behaviour of the small molecules, which are mainly buried in the interfacial regions. We have observed a competition between the binding of small molecules to phospholipids and cholesterol through lipidic hydrogen-bonds. Free energy barriers that are associated to translational and orientational changes of melatonin have been found to be between 10-20 kJ/mol for distances of 1 nm between melatonin and the center of the membrane. Corresponding barriers for tryptophan and serotonin that are obtained from reversible work methods are of the order of 10 kJ/mol and reveal strong hydrogen bonding between such species and specific phospholipid sites. The diffusion of tryptophan and melatonin is of the order of 10-7 cm2/s for the cholesterol-free and cholesterol-rich setups.

Tailoring patchy nanoparticle design to modulate serum albumin adsorption and membrane interaction

When nanoparticles (NPs) enter into the biological system, a wide range of proteins will coat on their surfaces forming protein corona, which changes the initial synthetic characteristics of NPs to the biological identity, resulting in the loss of their targets or specially designed properties. Although pre-coating with proteins would reduce the protein corona formation, they may diminish the targeting moieties in the transport process. Patchy NPs can offer unique advantages of asymmetry, heterogeneity, and multi-functions. This has inspired us to use the asymmetry to realize the versatility of NPs, to accommodate stealth and targeting functions. In this study, we performed molecular dynamics simulations to investigate the adsorption mechanism between patchy NPs and human serum albumin, and the interaction mechanism between NP-HSA and the membrane. The results show that there is a high probability for HSA to interact with the hydrophobic, or charged brushes of patchy NPs. The adsorption sites, as calculated through the contact probability between NPs and the residues, depend on the NP surface properties. Furthermore, the HSA adsorption on NPs could improve the NP-membrane interaction. The simulation results provide deep understanding of the NP interaction mechanism, which would help the NP design for their biomedical applications.

Comparison of umbrella sampling and steered molecular dynamics methods for computing free energy profiles of aromatic substrates through phospholipid bilayers

Understanding the permeation of molecules through lipid membranes is fundamental for predicting the cellular uptake of solutes and drug delivery mechanisms. In molecular simulations, the usual approach is to compute the free energy (FE) profile of a molecule across a model lipid bilayer, which can then be used to estimate the permeability of the molecule. Umbrella Sampling (US), which involves carrying out a series of biased simulations along a defined reaction coordinate (usually the bilayer normal direction), is a popular method for the computation of such FE profiles. However, US can be challenging to implement because the results are dependent on the strength of the biasing potential and the spacing of windows along the reaction coordinate, which, in practice, are usually optimized by an inefficient trial and error approach. The Steered Molecular Dynamics implementation of the Jarzynski Equality (JE-SMD) has been identified as an alternative to equilibrium sampling methods for measuring the FE change across a reaction coordinate. In the JE-SMD approach, equilibrium FE values are evaluated from the average of rapid non-equilibrium trajectories, thus avoiding the practical issues that come with US. Here, we use three different corrections of the JE-SMD method to calculate the FE change for the translocation of two aromatic substrates, phenylalanine and toluene, across a lipid bilayer and compare the accuracy and computational efficiency of these approaches to the results obtained using US. We show evidence that when computing the free energy profile, the JE-SMD approach suffers from insufficient sampling convergence of the bilayer environment and is dependent on the characteristic of the aromatic substrate itself. We deduce that, despite its drawbacks, US remains the more viable approach of the two for computing the FE profile.

Native-like membrane models of E. coli polar lipid extract shed light on the importance of lipid composition complexity

BACKGROUND

Lipid-protein interactions stabilize protein oligomers, shape their structure, and modulate their function. Whereas in vitro experiments already account for the functional importance of lipids by using natural lipid extracts, in silico methods lack behind by embedding proteins in single component lipid bilayers. However, to accurately complement in vitro experiments with molecular details at very high spatio-temporal resolution, molecular dynamics simulations have to be performed in natural(-like) lipid environments.

RESULTS

To enable more accurate MD simulations, we have prepared four membrane models of E. coli polar lipid extract, a typical model organism, each at all-atom (CHARMM36) and coarse-grained (Martini3) representations. These models contain all main lipid headgroup types of the E. coli inner membrane, i.e., phosphatidylethanolamines, phosphatidylglycerols, and cardiolipins, symmetrically distributed between the membrane leaflets. The lipid tail (un)saturation and propanylation stereochemistry represent the bacterial lipid tail composition of E. coli grown at 37∘C until 3/4 of the log growth phase. The comparison of the Simple three lipid component models to the complex 14-lipid component model Avanti over a broad range of physiologically relevant temperatures revealed that the balance of lipid tail unsaturation and propanylation in different positions and inclusion of lipid tails of various length maintain realistic values for lipid mobility, membrane area compressibility, lipid ordering, lipid volume and area, and the bilayer thickness. The only Simple model that was able to satisfactory reproduce most of the structural properties of the complex Avanti model showed worse agreement of the activation energy of basal water permeation with the here performed measurements. The Martini3 models reflect extremely well both experimental and atomistic behavior of the E. coli polar lipid extract membranes. Aquaporin-1 embedded in our native(-like) membranes causes partial lipid ordering and membrane thinning in its vicinity. Moreover, aquaporin-1 attracts and temporarily binds negatively charged lipids, mainly cardiolipins, with a distinct cardiolipin binding site in the crevice at the contact site between two monomers, most probably stabilizing the tetrameric protein assembly.

CONCLUSIONS

The here prepared and validated membrane models of E. coli polar lipids extract revealed that lipid tail complexity, in terms of double bond and cyclopropane location and varying lipid tail length, is key to stabilize membrane properties over a broad temperature range. In addition, they build a solid basis for manifold future simulation studies on more realistic lipid membranes bridging the gap between simulations and experiments.

Interplay of cholesterol, membrane bilayers and the AT1R: A cholesterol consensus motif on AT1R is revealed

Hypertension, mediated by the Angiotensin II receptor type 1 (AT1R), is still the major cause of premature death despite the discovery of novel therapeutics, highlighting the importance of an in depth understanding of the drug-AT1R recognition mechanisms coupled with the impact of the membrane environment on the interaction of drugs with AT1R. Herein, we examine the interplay of cholesterol-lipid-candesartan and the AT1R using Molecular Dynamics simulations of a model membrane consisting of 60:40 mol%. DPPC:cholesterol, candesartan and the AT1R, mimicking the physiological cholesterol concentration in sarcolemma membranes. The simulations of the model membrane of 60:40 mol%. DPPC:cholesterol were further validated using DOSY NMR experiments. Interestingly, our results suggest a significant role of cholesterol in the AT1R function imposed through a Cholesterol Consensus Motif (CCM) in the receptor, which could be crucial in the drug binding process. Candesartan diffusion towards AT1R through incorporation into lipid bilayers, appears to be retarded by the presence of cholesterol. However, its direct approach towards AT1R may be facilitated through the mobility induced on the N-terminus by the cholesterol binding on the CCM these novel insights could pave the way towards the development of more potent pharmaceutical agents to combat hypertension more effectively.

Coarse-grained implicit solvent lipid force field with a compatible resolution to the Cα protein representation

Biological membranes have been prominent targets for coarse-grained (CG) molecular dynamics simulations. While minimal CG lipid models with three beads per lipid and quantitative CG lipid models with >10 beads per lipid have been well studied, in between them, CG lipid models with a compatible resolution to residue-level CG protein models are much less developed. Here, we extended a previously developed three-bead lipid model into a five-bead model and parameterized it for two phospholipids, POPC (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine) and DPPC (1,2-dipalmitoyl-sn-glycero-3-phosphatidylcholine). The developed model, iSoLF, reproduced the area per lipid, hydrophobic thickness, and phase behaviors of the target phospholipid bilayer membranes at the physiological temperature. The model POPC and DPPC membranes were in liquid and gel phases, respectively, in accordance with experiments. We further examined the spontaneous formation of a membrane bilayer, the temperature dependence of physical properties, the vesicle dynamics, and the POPC/DPPC two-component membrane dynamics of the CG lipid model, showing some promise. Once combined with standard Cα protein models, the iSoLF model will be a powerful tool to simulate large biological membrane systems made of lipids and proteins.

Interleaflet Decoupling in a Lipid Bilayer at Excess Cholesterol Probed by Spectroscopic Ellipsometry and Simulations

Artificial lipid membranes are often investigated as a replica of the cell membrane in the form of supported lipid bilayers (SLBs). In SLBs, the phase state of a lipid bilayer strongly depends on the presence of molecules such as cholesterol, ceramide, and physical parameters such as temperature. Cholesterol is a key molecule of biological membranes and it exerts condensing effect on lipid bilayers. In this paper, we demonstrate the influence of excess cholesterol content on a supported lipid bilayer of 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC) (fluid-phase) using spectroscopic ellipsometry (SE) and coarse-grained (CG) molecular dynamics (MD) simulations. The results show the condensation effect due to cholesterol addition up to 30% and interleaflet decoupling at excess cholesterol beyond 30%. SE results show the separation of individual leaflets of the bilayer and influence of cholesterol on the biophysical properties such as thickness and optical index. CG simulations were performed at different ratios of DOPC:cholesterol mixtures to explore cholesterol-driven bilayer properties and stability. The simulations displayed the accumulation of cholesterol molecules at the interface of the lower and upper leaflets of the bilayer, thus leading to undulations in the bilayer. This work reports the successful application of SE technique to study lipid-cholesterol interactions for the first time.

Structure of supported DPPC/cholesterol bilayers studied via X-ray reflectivity

The electron density profile of bilayers of DPPC/cholesterol mixtures supported on semiconductor grade silicon substrates were studied with the objective of determining how the proximity of a solid interface modifies the phase diagram of mixed bilayers. The bilayers were studied in situ immersed in water via synchrotron X-ray reflectivity (XRR). Measurements were performed as a function of temperature through the main phase transition and cholesterol mole fractions up to 40%. Analysis of XRR yields the bilayer thickness, roughness and leaflet asymmetry. We find that the structure of the pure DPPC bilayers in the gel phase is in agreement with previous X-ray measurements of supported bilayers deposited via vesicle fusion and multilamellar vesicles but show more clearly defined features than measurements made on films formed using Langmuir-Blodget Langmuir-Shaffer (LB) deposition. Examination of bilayer thickness vs. temperature shows that the melting temperature for supported bilayers is shifted upwards by approximately 4 °C relative to multilamellar vesicles and that the melting temperature decreases with increasing cholesterol content up to 20%. For pure DPPC bilayers the leaflets melt in two stages with the distal leaflet melting first. For cholesterol concentrations of 10% and 20% there is no clear indication of separate melting. For 33% and 40% cholesterol content no clear transition is seen in the bilayer thickness, but an abrupt change in roughness indicates possible microdomain formation in the 40% cholesterol sample.

Insights into the interactions of bisphenol and phthalate compounds with unamended and carnitine-amended montmorillonite clays

Montmorillonite clays could be promising sorbents to mitigate toxic compound exposures. Bisphenols A (BPA) and S (BPS) as well as phthalates, dibutyl phthalate (DBP) and di-2-ethylhexyl phthalate (DEHP), are ubiquitous environmental contaminants linked to adverse health effects. Here, we combined computational and experimental methods to investigate the ability of montmorillonite clays to sorb these compounds. Molecular dynamics simulations predicted that parent, unamended, clay has higher binding propensity for BPA and BPS than for DBP and DEHP; carnitine-amended clay improved BPA and BPS binding, through carnitine simultaneously anchoring to the clay through its quaternary ammonium cation and forming hydrogen bonds with BPA and BPS. Experimental isothermal analysis confirmed that carnitine-amended clay has enhanced BPA binding capacity, affinity and enthalpy. Our studies demonstrate how computational and experimental methods, combined, can characterize clay binding and sorption of toxic compounds, paving the way for future investigation of clays to reduce BPA and BPS exposure.

Regulating Lipid Composition Rationalizes Acyl Tail Saturation Homeostasis in Ectotherms

Cell membranes mainly consist of lipid bilayers with an actively regulated composition. The underlying processes are still poorly understood, in particular, how the hundreds of components are controlled. Cholesterol has been found to correlate with phospholipid saturation for reasons that remain unclear. To better understand the link between cell membrane regulation and chemical composition, we establish a computational framework based on chemical reaction networks, resulting in multiple semigrand canonical ensembles. By running computer simulations, we show that regulating the chemical potential of lipid species is sufficient to reproduce the experimentally observed increase in acyl tail saturation with added cholesterol. Our model proposes a different picture of lipid regulation in which components can be regulated passively instead of actively. In this picture, phospholipid acyl tail composition naturally adapts to added molecules such as cholesterol or proteins. A comparison between regulated membranes with commonly studied ternary model membranes shows stark differences: for instance, correlation lengths and viscosities observed are independent of lipid chemical affinity.

Principal component analysis highlights the influence of temperature, curvature and cholesterol on conformational dynamics of lipids

Membrane lipids are inherently highly dynamic molecules. Currently, it is difficult to probe the structures of individual lipids experimentally at the timescales corresponding to atomic motions, and consequently molecular dynamics simulations are used widely. In our previous work, we have introduced the principal component analysis (PCA) as a convenient framework for comprehensive quantitative description of lipid motions. Here, we present a newly developed open source script, PCAlipids, which automates the analysis and allows us to refine the approach and test its limitations. We use PCAlipids to determine the influence of temperature, cholesterol and curvature on individual lipids, and show that the most prominent lipid tail scissoring motion is strongly affected by these factors and allows tracking of phase transition. Addition of cholesterol affects the conformations and selectively changes the dynamics of lipid molecules, impacting the large-amplitude motions. Introduction of curvature biases the conformational ensembles towards more extended structures. We hope that the developed approach will be useful for understanding the molecular basis of different processes occurring in lipid membrane systems and will stimulate development of complementary experimental techniques probing the conformations of individual lipid molecules.