Molecular order in cis and trans unsaturated phospholipid bilayers

The Chemical Reactivity of Membrane Lipids

It is well-known that aqueous dispersions of phospholipids spontaneously assemble into bilayer structures. These structures have numerous applications across chemistry and materials science and form the fundamental structural unit of the biological membrane. The particular environment of the lipid bilayer, with a water-poor low dielectric core surrounded by a more polar and better hydrated interfacial region, gives the membrane particular biophysical and physicochemical properties and presents a unique environment for chemical reactions to occur. Many different types of molecule spanning a range of sizes, from dissolved gases through small organics to proteins, are able to interact with membranes and promote chemical changes to lipids that subsequently affect the physicochemical properties of the bilayer. This Review describes the chemical reactivity exhibited by lipids in their membrane form, with an emphasis on conditions where the lipids are well hydrated in the form of bilayers. Key topics include the following: lytic reactions of glyceryl esters, including hydrolysis, aminolysis, and transesterification; oxidation reactions of alkenes in unsaturated fatty acids and sterols, including autoxidation and oxidation by singlet oxygen; reactivity of headgroups, particularly with reactive carbonyl species; and E/Z isomerization of alkenes. The consequences of reactivity for biological activity and biophysical properties are also discussed.

Unsaturated fatty acid synthesis in bacteria: Mechanisms and regulation of canonical and remarkably noncanonical pathways

Unsaturated phospholipid acyl chains are required for membrane function in most bacteria. The double bonds of the cis monoenoic chains arise by two distinct pathways depending on whether oxygen is required. The oxygen-independent pathway (traditionally called the anaerobic pathway) introduces the cis double bond by isomerization of the trans double bond intermediate of the fatty acid elongation cycle. Double bond isomerization occurs at an intermediate chain length (e.g., C10) and the isomerization product is elongated to the C16-C18 chains that become phospholipid monoenoic acyl chains. This pathway was first delineated in Escherichia coli and became the paradigm pathway. However, studies of other bacteria show deviations from this paradigm, the most exceptional being reversal of the fatty acid elongation cycle by a reaction paralleling the initial step in the β-oxidative degradation of fatty acids. In the oxygen-dependent pathway diiron enzymes called desaturases introduce a double bond into a saturated acyl chain by regioselective cis dehydrogenation through activation of molecular oxygen with an active-site diiron cluster. This difficult hydrogen abstraction from a methylene group often occurs at the midpoint of a saturated fatty acyl chain. In bacteria the acyl chain is a phospholipid acyl chain, and the desaturase is membrane bound. Both the oxygen-independent oxygen-dependent pathways are transcriptionally regulated by repressor and activator proteins that respond to small molecule ligands such as acyl-CoAs. However, in Bacillus subtilis the desaturase is synthesized only at low growth temperatures, a process controlled by a signal transduction regulatory pathway dependent on membrane lipid properties.

Functional regulation of aquaporin dynamics by lipid bilayer composition

With the diversity of lipid-protein interactions, any observed membrane protein dynamics or functions directly depend on the lipid bilayer selection. However, the implications of lipid bilayer choice are seldom considered unless characteristic lipid-protein interactions have been previously reported. Using molecular dynamics simulation, we characterize the effects of membrane embedding on plant aquaporin SoPIP2;1, which has no reported high-affinity lipid interactions. The regulatory impacts of a realistic lipid bilayer, and nine different homogeneous bilayers, on varying SoPIP2;1 dynamics are examined. We demonstrate that SoPIP2;1's structure, thermodynamics, kinetics, and water transport are altered as a function of each membrane construct's ensemble properties. Notably, the realistic bilayer provides stabilization of non-functional SoPIP2;1 metastable states. Hydrophobic mismatch and lipid order parameter calculations further explain how lipid ensemble properties manipulate SoPIP2;1 behavior. Our results illustrate the importance of careful bilayer selection when studying membrane proteins. To this end, we advise cautionary measures when performing membrane protein molecular dynamics simulations.

Dynamical Crossover of Interfacial Water upon Melting of a Lipid Bilayer: Influence of Different Parts of the Headgroups

All-atom molecular dynamics simulations of a 1,2-dimyristoyl-sn-glycero-3-phosphocholine bilayer in contact with liquid water were performed at different temperatures ranging from 285 to 320 K. We have investigated the heterogeneity and dynamical transitions in interfacial water as the lipid bilayer undergoes a melting transition. Results are obtained for water at the outer surface of the bilayer and for those buried more deeply in the lipid chains of the bilayer. It is found that lipid bilayer melting influences both the structure and dynamics of interfacial water. The number of interfacial water molecules shows a jump in the melting of the bilayer. The temperature dependence of the diffusivity and orientational relaxation of interfacial water molecules exhibits a dynamical crossover upon melting of the bilayer. The extent of dynamical crossover is found to be rather strong with significant changes in activation barriers for interfacial water around the carbonyl groups, which are deeply buried toward the lipid chains of the bilayer. The dynamical crossover gradually decreases as one moves further away from the outer surface, and it essentially vanishes for water in the region of 5-10 Å from the outer surface. It is found that the lipid melting-induced dynamical crossover of interfacial water is significant only for water that is in close proximity to the bilayer surface or deeply buried into it. The current results reveal that water molecules in different parts of the interface respond differently on melting of the bilayer. The current study also shows that the carbonyl-bound water molecules can play an important role in the phase transition of the bilayer as the temperature is raised through its melting point.

Functional Regulation of Aquaporin Dynamics by Lipid Bilayer Composition

With the diversity of lipid-protein interactions, any observed membrane protein dynamics or functions directly depend on the lipid bilayer selection. However, the implications of lipid bilayer choice are seldom considered unless characteristic lipid-protein interactions have been previously reported. Using molecular dynamics simulation, we characterize the effects of membrane embedding on plant aquaporin SoPIP2;1, which has no reported high-affinity lipid interactions. The regulatory impacts of a realistic lipid bilayer, and nine different homogeneous bilayers, on varying SoPIP2;1 dynamics were examined. We demonstrate that SoPIP2;1s structure, thermodynamics, kinetics, and water transport are altered as a function of each membrane construct's ensemble properties. Notably, the realistic bilayer provides stabilization of non-functional SoPIP2;1 metastable states. Hydrophobic mismatch and lipid order parameter calculations further explain how lipid ensemble properties manipulate SoPIP2;1 behavior. Our results illustrate the importance of careful bilayer selection when studying membrane proteins. To this end, we advise cautionary measures when performing membrane protein molecular dynamics simulations.

Structural Impact of Selected Retinoids on Model Photoreceptor Membranes

Photoreceptor membranes have a unique lipid composition. They contain a high level of polyunsaturated fatty acids including the most unsaturated fatty acid in nature, docosahexaenoic acid (22:6), and are enriched in phosphatidylethanolamines. The phospholipid composition and cholesterol content of the subcellular components of photoreceptor outer segments enables to divide photoreceptor membranes into three types: plasma membranes, young disc membranes, and old disc membranes. A high degree of lipid unsaturation, extended exposure to intensive irradiation, and high respiratory demands make these membranes sensitive to oxidative stress and lipid peroxidation. Moreover, all-trans retinal (AtRAL), which is a photoreactive product of visual pigment bleaching, accumulates transiently inside these membranes, where its concentration may reach a phototoxic level. An elevated concentration of AtRAL leads to accelerated formation and accumulation of bisretinoid condensation products such as A2E or AtRAL dimers. However, a possible structural impact of these retinoids on the photoreceptor-membrane properties has not yet been studied. In this work we focused just on this aspect. The changes induced by retinoids, although noticeable, seem not to be significant enough to be physiologically relevant. This is, however, an positive conclusion because it can be assumed that accumulation of AtRAL in photoreceptor membranes will not affect the transduction of visual signals and will not disturb the interaction of proteins engaged in this process.

Characterizing the Interactions of Cell-Membrane-Disrupting Peptides with Lipid-Functionalized Single-Walled Carbon Nanotubes

Lipid-functionalized single-walled carbon nanotubes (SWNTs) have garnered significant interest for their potential use in a wide range of biomedical applications. In this work, we used molecular dynamics simulations to study the equilibrium properties of SWNTs surrounded by the phosphatidylcholine (POPC) corona phase and their interactions with three cell membrane disruptor peptides: colistin, TAT peptide, and crotamine-derived peptide. Our results show that SWNTs favor asymmetrical positioning within the POPC corona, so that one side of the SWNT, covered by the thinnest part of the corona, comes in contact with charged and polar functional groups of POPC and water. We also observed that colistin and TAT insert deeply into the POPC corona, while crotamine-derived peptide only adsorbs to the corona surface. In separate simulations, we show that three examined peptides exhibit similar insertion and adsorption behaviors when interacting with POPC bilayers, confirming that peptide-induced perturbations to POPC in conjugates and bilayers are similar in nature and magnitude. Furthermore, we observed correlations between the peptide-induced structural perturbations and the near-infrared emission of the lipid-functionalized SWNTs, which suggest that the optical signal of the conjugates transduces the morphological changes in the lipid corona. Overall, our findings indicate that lipid-functionalized SWNTs could serve as simplified cell membrane model systems for prescreening of new antimicrobial compounds that disrupt cell membranes.

Characterizing the Interactions of Cell Membrane-Disrupting Peptides with Lipid-Functionalized Single-Walled Carbon Nanotube Systems for Antimicrobial Screening

Lipid-functionalized single-walled carbon nanotubes (SWNTs) have garnered significant interest for their potential use in a wide range of biomedical applications. In this work, we used molecular dynamics simulations to study the equilibrium properties of SWNTs surrounded by the phosphatidylcholine (POPC) corona phase, and their interactions with three cell membrane disruptor peptides: colistin, TAT peptide, and crotamine-derived peptide. Our results show that SWNTs favor asymmetrical positioning within the POPC corona, so that one side of the SWNT, covered by the thinnest part of the corona, comes in contact with charged and polar functional groups of POPC and water. We also observed that colistin and TAT insert deeply into POPC corona, while crotamine-derived peptide only adsorbs to the corona surface. Compared to crotamine-derived peptide, colistin and TAT also induce larger perturbations in the thinnest region of the corona, by allowing more water molecules to directly contact the SWNT surface. In separate simulations, we show that three examined peptides exhibit similar insertion and adsorption behaviors when interacting with POPC bilayers, confirming that peptide-induced perturbations to POPC in conjugates and bilayers are similar in nature and magnitude. Furthermore, we observed correlations between the peptide-induced structural perturbations and the near-infrared emission of the lipid-functionalized SWNTs, which suggest that the optical signal of the conjugates transduces the morphological changes in the lipid corona. Overall, our findings indicate that lipid-functionalized SWNTs could serve as simplified cell membrane model systems for pre-screening of new antimicrobial compounds that disrupt cell membranes.

Structure of POPC Lipid Bilayers in OPLS3e Force Field

It is crucial for molecular dynamics simulations of biomembranes that the force field parameters give a realistic model of the membrane behavior. In this study, we examined the OPLS3e force field for the carbon-hydrogen order parameters S_CH of POPC (1-palmitoyl-2-oleoylphosphatidylcholine) lipid bilayers at varying hydration conditions and ion concentrations. The results show that OPLS3e behaves similarly to the CHARMM36 force field and relatively accurately follows the experimentally measured S_CH for the lipid headgroup, the glycerol backbone, and the acyl tails. Thus, OPLS3e is a good choice for POPC bilayer simulations under many biologically relevant conditions. The exception are systems with an abundancy of ions, as similarly to most other force fields OPLS3e strongly overestimates the membrane-binding of cations, especially Ca²⁺. This leads to undesirable positive charge of the membrane surface and drastically lowers the concentration of Ca²⁺ in the surrounding solvent, which might cause issues in systems sensitive to correct charge distribution profiles across the membrane.

Effect of encapsulated protein on the dynamics of lipid sponge phase: a neutron spin echo and molecular dynamics simulation study

Lipid membranes are highly mobile systems with hierarchical, time and length scale dependent, collective motions including thickness fluctuations, undulations, and topological membrane changes, which play an important role in membrane interactions. In this work we have characterised the effect of encapsulating two industrially important enzymes, β-galactosidase and aspartic protease, in lipid sponge phase nanoparticles on the dynamics of the lipid membrane using neutron spin echo (NSE) spectroscopy and molecular dynamics (MD) simulations. From NSE, reduced membrane dynamics were observed upon enzyme encapsulation, which were dependent on the enzyme concentration and type. By fitting the intermediate scattering functions (ISFs) with a modified Zilman and Granek model including nanoparticle diffusion, an increase in membrane bending rigidity was observed, with a larger effect for β-galactosidase than aspartic protease at the same concentration. MD simulations for the system with and without aspartic protease showed that the lipids relax more slowly in the system with protein due to the replacement of the lipid carbonyl-water hydrogen bonds with lipid-protein hydrogen bonds. This indicates that the most likely cause of the increase in membrane rigidity observed in the NSE measurements was dehydration of the lipid head groups. The dynamics of the protein itself were also studied, which showed a stable secondary structure of protein over the simulation, indicating no unfolding events occurred.

The interaction of steroids with phospholipid bilayers and membranes

Steroids are critical for various physiological processes and used to treat inflammatory conditions. Steroids act by two distinct pathways. The genomic pathway is initiated by the steroid binding to nuclear receptors while the non-genomic pathway involves plasma membrane receptors. It has been proposed that steroids might also act in a more indirect mechanism by altering biophysical properties of membranes. Yet, little is known about the effect of steroids on membranes, and steroid-membrane interactions are complex and challenging to characterise. The focus of this review is to outline what is currently known about the interactions of steroids with phospholipid bilayers and illustrate the complexity of these systems using cortisone and progesterone as the main examples. The combined findings from current work demonstrate that the hydrophobicity and planarity of the steroid core does not provide a consensus for steroid-membrane interactions. Even small differences in the substituents on the steroid core can result in significant changes in steroid-membrane interactions. Furthermore, steroid-induced changes in phospholipid bilayer properties are often dependent on steroid concentration and lipid composition. This complexity means that currently there is insufficient information to establish a reliable structure-activity relationship to describe the effect of steroids on membrane properties. Future work should address the challenge of connecting the findings from studying the effect of steroids on phospholipid bilayers to cell membranes. Insights from steroid-membrane interactions will benefit our understanding of normal physiology and assist drug development.

Microsecond molecular dynamics studies of cholesterol-mediated myelin sheath degeneration in early Alzheimer's disease

Cholesterol-mediated perturbations of membrane structural integrity are key early events in the molecular pathogenesis of Alzheimer's disease (AD). In AD, protein misfolding (proteopathy) and pro-inflammatory conditions (immunopathy) culminate in neuronal death, a process enabled by altered membrane biophysical properties which render neurons more susceptible to proteopathic and immunopathic cytotoxicities. Since cholesterol is a principal neuronal membrane lipid, normal cholesterol homeostasis is central to membrane health; also, since increased cholesterol composition is especially present in neuronal myelin sheath (i.e. brain "white matter"), recent studies have not surprisingly revealed that white matter atrophy precedes the conventional biomarkers of AD (amyloid plaques, tau tangles). Employing extensive microsecond all-atom molecular dynamics simulations, we investigated biophysical and mechanical properties of myelin sheath membrane as a function of cholesterol mole fraction (χ_CHL). Impaired χ_CHL modulates multiple bilayer properties, including surface area per lipid (APL), chain order, number and mass density profiles, area compressibility and bending moduli, bilayer thickness, lipid tilt angles, H-bonding interactions and tail interdigitation. The increased orientational ordering of both palmitoyl and oleoyl chains in model healthy myelin sheath (HMS) membranes illustrates the condensing effect of cholesterol. With an increase in χ_CHL, number density profiles of water tend to attain bulk water number density more quickly, indicating shrinkage in the interfacial region with increasing χ_CHL. The average tilt value is 11.5° for the C10-C13 angle in cholesterol and 64.2° for the P-N angle in POPC lipids in HMS. These calculations provide a molecular-level understanding of myelin sheath susceptibility to pathology as an early event in the pathogenesis of AD.

Water/oil interfacial tension reduction - an interfacial entropy driven process

The interfacial tension (IFT) of a fluid-fluid interface plays an important role in a wide range of applications and processes. When low IFT is desired, surface active compounds (e.g. surfactants) can be added to the system. Numerous attempts have been made to relate changes in IFT arising from such compounds to the specific nature of the interface. However, the IFT is controlled by an interplay of factors such as temperature and molecular structure of surface-active compounds, which make it difficult to predict IFT as those conditions change. In this study, we present the results from molecular dynamics simulations revealing the specific role surfactants play in IFT. We find that, in addition to reducing direct contact between the two fluids, surfactants serve to increase the disorder at the interface (related to interfacial entropy) and consequently reduce the water/oil IFT, especially when surfactants are present at high surface density. Our results suggest that surfactants that yield more disordered interfacial films (e.g. with flexible and/or unsaturated tails) reduce the water/oil IFT more effectively than surfactants which yield highly ordered interfacial films. Our results shed light on some of the factors that control IFT and could have important practical implications in industrial applications such as the design of cosmetics, food products, and detergents.

Axonal plasma membrane-mediated toxicity of cholesterol in Alzheimer's disease: A microsecond molecular dynamics study

Alzheimer's disease is increasingly being recognized as an immune-mediated disease of brain. Since physiological brain health and brain immune function is dependent upon homeostatic neuronal membrane structure and function, alterations in membrane lipid biochemistry may predispose to disease. Brain is rich in cholesterol, and cholesterol metabolism dysfunction is a known risk factor for AD. Employing extensive microsecond all-atom molecular dynamics simulations, we investigated the properties of model neuronal membranes as a function of cholesterol concentration; phospholipid and phospholipid/cholesterol bilayers were also simulated to compare against available experimental data. Increased cholesterol concentrations compact and stiffen the lipid membrane, reducing permeability while modulating local water densities in the peri-membranous environment. Conversely, lower cholesterol mole fraction yields membranes with increased molecular disorder, enhanced fluidity, higher molecular tilting, and augmented interdigitation between bilayer leaflet lipids. Our findings provide a molecular insight on effect of cholesterol composition on various biochemical processes occurring at neuronal axon plasma membrane. These calculations also endeavor to establish a membrane-based link between cholesterol as an AD risk factor and possible AD pathology.

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-dynamics-simulation-guided membrane engineering allows the increase of membrane fatty acid chain length in Saccharomyces cerevisiae

The use of lignocellulosic-based fermentation media will be a necessary part of the transition to a circular bio-economy. These media contain many inhibitors to microbial growth, including acetic acid. Under industrially relevant conditions, acetic acid enters the cell predominantly through passive diffusion across the plasma membrane. The lipid composition of the membrane determines the rate of uptake of acetic acid, and thicker, more rigid membranes impede passive diffusion. We hypothesized that the elongation of glycerophospholipid fatty acids would lead to thicker and more rigid membranes, reducing the influx of acetic acid. Molecular dynamics simulations were used to predict the changes in membrane properties. Heterologous expression of Arabidopsis thaliana genes fatty acid elongase 1 (FAE1) and glycerol-3-phosphate acyltransferase 5 (GPAT5) increased the average fatty acid chain length. However, this did not lead to a reduction in the net uptake rate of acetic acid. Despite successful strain engineering, the net uptake rate of acetic acid did not decrease. We suggest that changes in the relative abundance of certain membrane lipid headgroups could mitigate the effect of longer fatty acid chains, resulting in a higher net uptake rate of acetic acid.

Long-Chain Fatty Acids in Bones and Their Link to Submicroscopic Vascularization Network: NMR Assignment and Relaxation Studies under Magic Angle Spinning Conditions in Intramuscular Bones of Atlantic Herring Fish

The long-lasting proton signals in bones are identified as long-chain fatty acids, including saturated, mono-, and di-unsaturated fatty acids, with direct nuclear magnetic resonance evidence. We used intramuscular bones from Atlantic Herring fish to avoid interference from lipid-rich marrows. The key is to recognize that these signals are from mobile phase materials and study them with J-coupled correlation spectroscopies under magic angle spinning conditions. We kept extensive ¹H-spin-echo records that allowed us to examine the effect of magic angle spinning on the transverse relaxation time of water and lipids over time. While it is impossible to distinguish based on chemical shifts, the relaxation data suggest that the signals are more consistent with the interpretation of phospholipid membranes than triglycerides in lipid droplets. In particular, the simultaneous T₂ changes in water and lipids suggest that the centrifugal impact of magic angle spinning alters the lipid's structure in very tight spaces. Additional evidence of phospholipid membranes came from the choline-γ resonance at 3.2 ppm in fresh samples, which disappears with magic angle spinning. Thus, the fatty acid signals are at least partially from membrane bilayer structures, and we propose that they are linked to the submicroscopic vascularization channels similar to the dense canaliculi network in mammalian bones. Our detection of phospholipids from bones depended critically on two factors: (1) the elimination of the overwhelming triglyceride signals from marrows and (2) the preservation of water that biomembranes require. The relaxation data reveal aspects of lipid fluidity that have not been elucidated by previous order parameter studies on model membranes. Relaxation times have long been considered difficult to interpret. A robust and renewed understanding may be beneficial.

Tauroursodeoxycholic Acid (TUDCA)-Lipid Interactions and Antioxidant Properties of TUDCA Studied in Model of Photoreceptor Membranes

Tauroursodeoxycholic acid (TUDCA), a hydrophilic bile acid containing taurine conjugated with the ursodeoxycholic acid (UDCA), has been known and used from ancient times as a therapeutic compound in traditional Chinese medicine. TUDCA has recently been gaining significant interest as a neuroprotective agent, also exploited in the visual disorders. Among several mechanisms of TUDCA's protective action, its antioxidant activity and stabilizing effect on mitochondrial and plasma membranes are considered. In this work we investigated antioxidant activity of TUDCA and its impact on structural properties of model membranes of different composition using electron paramagnetic resonance spectroscopy and the spin labeling technique. Localization of TUDCA molecules in a pure POPC bilayer has been studied using a molecular dynamics simulation (MD). The obtained results indicate that TUDCA is not an efficient singlet oxygen (1O2 (1Δg)) quencher, and the determined rate constant of its interaction with 1O2 (1Δg) is only 1.9 × 105 M-1s-1. However, in lipid oxidation process induced by a Fenton reaction, TUDCA reveals substantial antioxidant activity significantly decreasing the rate of oxygen consumption in the system studied. In addition, TUDCA induces slight, but noticeable changes in the polarity and fluidity of the investigated model membranes. The results of performed MD simulation correspond very well with the experimental results.

Effects of Conformational Variation on Structural Insights from Solution-Phase Surface-Enhanced Raman Spectroscopy

Surface-enhanced Raman scattering (SERS) spectra contain information on the chemical structure on nanoparticle surfaces through the position and alignment of molecules with the electromagnetic near field. Time-dependent density functional theory (TDDFT) can provide the Raman tensors needed for a detailed interpretation of SERS spectra. Here, the impact of molecular conformations on SERS spectra is considered. TDDFT calculations of the surfactant cetyltrimethylammonium bromide with five conformers produced more accurate unenhanced Raman spectra than a simple all-trans structure. The calculations and measurements also demonstrated a loss of structural information in the CH2/CH3 scissor vibration band at 1450 cm-1 in the SERS spectra. To study lipid bilayers, TDDFT calculations on conformers of methyl phosphorylcholine and cis-5-decene served as models for the symmetric choline stretch in the lipid headgroup and the C═C stretch in the acyl chains of 1,2-oleoyl-glycero-3-phosphocholine. Conformer considerations enabled a measurement of the distribution of double-bond orientations with an order parameter of SC═C = 0.53.

Effect of TMAO on the Structure and Phase Transition of Lipid Membranes: Potential Role of TMAO in Stabilizing Cell Membranes under Osmotic Stress

Extremophiles adopt strategies to deal with different environmental stresses, some of which are severely damaging to their cell membrane. To combat high osmotic stress, deep-sea organisms synthesize osmolytes, small polar organic molecules, like trimethylamine-N-oxide (TMAO), and incorporate them in the cell. TMAO is known to protect cells from high osmotic or hydrostatic pressure. Several experimental and simulation studies have revealed the roles of such osmolytes on stabilizing proteins. In contrast, the effect of osmolytes on the lipid membrane is poorly understood and broadly debated. A recent experiment has found strong evidence of the possible role of TMAO in stabilizing lipid membranes. Using the molecular dynamics (MD) simulation technique, we have demonstrated the effect of TMAO on two saturated fully hydrated lipid membranes in their fluid and gel phases. We have captured the impact of TMAO's concentration on the membrane's structural properties along with the fluid/gel phase transition temperatures. On increasing the concentration of TMAO, we see a substantial increase in the packing density of the membrane (estimated by area, thickness, and volume) and enhancement in the orientational order of lipid molecules. Having repulsive interaction with the lipid head group, the TMAO molecules are expelled away from the membrane surface, which induces dehydration of the lipid head groups, increasing the packing density. The addition of TMAO also increases the fluid/gel phase transition temperature of the membrane. All of these results are in close agreement with the experimental observations. This study, therefore, provides a molecular-level understanding of how TMAO can influence the cell membrane of deep-sea organisms and help in combating the osmotic stress condition.

Multivalent lipid targeting by the calcium-independent C2A domain of synaptotagmin-like protein 4/granuphilin

Synaptotagmin-like protein 4 (Slp-4), also known as granuphilin, is a Rab effector responsible for docking secretory vesicles to the plasma membrane before exocytosis. Slp-4 binds vesicular Rab proteins via an N-terminal Slp homology domain, interacts with plasma membrane SNARE complex proteins via a central linker region, and contains tandem C-terminal C2 domains (C2A and C2B) with affinity for phosphatidylinositol-(4,5)-bisphosphate (PIP₂). The Slp-4 C2A domain binds with low nanomolar apparent affinity to PIP₂ in lipid vesicles that also contain background anionic lipids such as phosphatidylserine but much weaker when either the background anionic lipids or PIP₂ is removed. Through computational and experimental approaches, we show that this high-affinity membrane binding arises from concerted interaction at multiple sites on the C2A domain. In addition to a conserved PIP₂-selective lysine cluster, a larger cationic surface surrounding the cluster contributes substantially to the affinity for physiologically relevant lipid compositions. Although the K398A mutation in the lysine cluster blocks PIP₂ binding, this mutated protein domain retains the ability to bind physiological membranes in both a liposome-binding assay and MIN6 cells. Molecular dynamics simulations indicate several conformationally flexible loops that contribute to the nonspecific cationic surface. We also identify and characterize a covalently modified variant that arises through reactivity of the PIP₂-binding lysine cluster with endogenous bacterial compounds and binds weakly to membranes. Overall, multivalent lipid binding by the Slp-4 C2A domain provides selective recognition and high-affinity docking of large dense core secretory vesicles to the plasma membrane.

Optimization of Slipids Force Field Parameters Describing Headgroups of Phospholipids

The molecular mechanics force field Slipids developed in a series of works by Jämbeck and Lyubartsev (J. Phys. Chem. B 2012, 116, 3164-3179; J. Chem. Theory Comput. 2012, 8, 2938-2948) generally provides a good description of various lipid bilayer systems. However, it was also found that order parameters of C-H bonds in the glycerol moiety of the phosphatidylcholine headgroup deviate significantly from NMR results. In this work, the dihedral force field parameters have been reparameterized in order to improve the agreement with experiment. For this purpose, we have computed energies for a large amount of lipid headgroup conformations using density functional theory on the B3P86/cc-pvqz level and optimized dihedral angle parameters simultaneously to provide the best fit to the quantum chemical energies. The new parameter set was validated for three lipid bilayer systems against a number of experimental properties including order parameters, area per lipid, scattering form factors, bilayer thickness, area compressibility and lateral diffusion coefficients. In addition, the order parameter dependence on cholesterol content in the POPC bilayer was investigated. It is shown that the new force field significantly improves agreement with the experimental order parameters for the lipid headgroup while keeping good agreement with other experimentally measured properties.

Side-by-side comparison of Notch- and C83 binding to γ-secretase in a complete membrane model at physiological temperature

γ-Secretase cleaves the C99 fragment of the amyloid precursor protein, leading to formation of aggregated β-amyloid peptide central to Alzheimer's disease, and Notch, essential for cell regulation. Recent cryogenic electron microscopy (cryo-EM) structures indicate major changes upon substrate binding, a β-sheet recognition motif, and a possible helix unwinding to expose peptide bonds towards nucleophilic attack. Here we report side-by-side comparison of the 303 K dynamics of the two proteins in realistic membranes using molecular dynamics simulations. Our ensembles agree with the cryo-EM data (full-protein Cα-RMSD = 1.62–2.19 Å) but reveal distinct presenilin helix conformation states and thermal β-strand to coil transitions of C83 and Notch100. We identify distinct 303 K hydrogen bond dynamics and water accessibility of the catalytic sites. The RKRR motif (1758–1761) contributes significantly to Notch binding and serves as a “membrane anchor” that prevents Notch displacement. Water that transiently hydrogen bonds to G1753 and V1754 probably represents the catalytic nucleophile. At 303 K, Notch and C83 binding induce different conformation states, with Notch mostly present in a closed state with shorter Asp–Asp distance. This may explain the different outcome of Notch and C99 cleavage, as the latter is more imprecise with many products. Our identified conformation states may aid efforts to develop conformation-selective drugs that target C99 and Notch cleavage differently, e.g. Notch-sparing γ-secretase modulators.

Distinct membrane dynamics and conformations of C83- and Notch-bound γ-secretase may aid the development of Notch-sparing treatments of Alzheimer's disease.

Update of the CHARMM36 United Atom Chain Model for Hydrocarbons and Phospholipids

Accurate lipid force field (FF) parameters used in molecular dynamics (MD) simulations are crucial for understanding the properties of lipid-containing systems and biological processes related to lipids. The last update of the CHARMM36 united atom chain model (C36UA) was in 2013 [Lee, S. J. Phys. Chem. B 2014, 118, 547 556]; it utilized CHARMM36 (C36) lipid FF parameters for headgroups and OPLS-UA Lennard-Jones (LJ) parameters for tails. Simulations with the FF were able to reproduce many experimental observables of lipid bilayers accurately, but to be more applicable for a wide range of lipids, additional FF parameter optimization was needed. In this work, we present an update of the model, named C36UAr. The parameterization included the LJ parameters for hydrocarbons and related dihedrals. Bulk liquid properties (density, heat of vaporization, isothermal compressibility, and diffusion constant) of model compounds were used as targets for the LJ parameter fitting, and dihedrals were fit to either quantum mechanical (QM) or potential of mean force (PMF) calculations using C36. Thermodynamic reweighting was used to further improve the parameters. Bilayer simulations of various lipid headgroups (phosphatidylcholine, phosphatidylethanolamine, and phosphatidylglycerol) and tails (saturated, monounsaturated, and polyunsaturated) were performed to validate the model, and significant improvements were seen in bilayer properties, including surface area, membrane thicknesses, NMR deuterium order parameters, and density profiles. C36UAr was also compared to the hydrogen mass repartitioning (HMR) method. The high accuracy and competitive efficiency shown in this study make C36UAr one of the best choices for studies of membrane structure and membrane-associated proteins.

The Effect of Force-Field Parameters on Cytochrome P450-Membrane Interactions: Structure and Dynamics

The simulation of membrane proteins requires compatible protein and lipid force fields that reproduce the properties of both the protein and the lipid bilayer. Cytochrome P450 enzymes are bitopic membrane proteins with a transmembrane helical anchor and a large cytosolic globular domain that dips into the membrane. As such, they are representative and challenging examples of membrane proteins for simulations, displaying features of both peripheral and integral membrane proteins. We performed molecular dynamics simulations of three cytochrome P450 isoforms (2C9, 2C19 and 1A1) in a 2-oleoyl-1-palmitoyl-sn-glycerol-3-phosphocholine bilayer using two AMBER force field combinations: GAFF-LIPID with ff99SB for the protein, and LIPID14 with ff14SB for the protein. Comparison of the structural and dynamic properties of the proteins, the lipids and the protein-membrane interactions shows differing sensitivity of the cytochrome P450 isoforms to the choice of force field, with generally better agreement with experiment for the LIPID14 + ff14SB combination.

A novel synthesis of trans-unsaturated fatty acids by the Gram-positive commensal bacterium Enterococcus faecalis FA2-2

A key mechanism of Pseudomonas spp. adaptation to environmental stressors is their ability to convert the cis-unsaturated fatty acids of the membrane lipids to their trans-isomers to rigidify the membrane and thereby resist stresses. Although this Cti-catalyzed enzymatic isomerization has been well investigated in the P. putida paradigm, several bacterial species have been found to produce trans-unsaturated fatty acids. Although cti orthologs have only been reported in Gram-negative bacteria, we report that E. faecalis FA2-2 cultures synthesize trans-unsaturated fatty acids during growth by a mechanism similar of P. putida. Although the role of trans-unsaturated fatty acids (trans-UFAs) in E. faecalis remains obscure, our results indicate that organic solvents, as well as the membrane altering antibiotic, daptomycin, had no effect on trans-UFA formation in E. faecalis FA2-2. Moreover trans-UFA production in E. faecalis FA2-2 membranes was constant in oxidative stress conditions or when metal chelator EDTA was added, raising the question about the role of heme domain in cis-trans isomerization in E. faecalis FA2-2. Although growth temperature and growth phase had significant effects on cis-trans isomerization, the bulk physical properties of the membranes seems unlikely to be altered by the low levels of trans-UFA. Hence, any effects seems likely to be on membrane proteins and membrane enzyme activities. We also report investigations of cti gene distribution in bacteria was and suggest the distribution to be triggered by habitat population associations. Three major Cti clusters were defined, corresponding to Pseudomonas, Pseudoalteromonas and Vibrio Cti proteins.

Coping with arsenic stress: Adaptations of arsenite-oxidizing bacterial membrane lipids to increasing arsenic levels

Elevated levels of arsenic (As) in aquifers of South East Asia have caused diverse health problems affecting millions of people who drink As-rich groundwater and consume various contaminated agriculture products. The biogeochemical cycling and mobilization/immobilization of As from its mineral-bound phase is controlled by pH, oxic/anoxic conditions, and different microbial processes. The increased As flux generated from ongoing biogeochemical processes in the subsurface in turn affects the in situ microbial communities. This study analyzes how the indigenous arsenite-oxidizing bacteria combat As stress by various biophysical alterations and self-adaptation mechanisms. Fifteen arsenite-oxidizing bacterial strains were isolated and identified using a polyphasic approach. The bacterial strains isolated from these aquifers belong predominantly to arsenite-oxidizing bacterial groups. Of these, the membrane-bound phospholipid fatty acids (PLFA) were characterized in seven selected bacterial isolates grown at different concentrations of As(III) in the medium. One of the significant findings of this study is how the increase in external stress can induce alteration of membrane PLFAs. The change in fatty acid saturation and alteration of their steric conformation suggests alteration of membrane fluidity due to change in As-related stress. However, different bacterial groups can have different degrees of alteration that can affect sustainability in As-rich aquifers of the Bengal Delta Plain.

Simulations of lipid bilayers using the CHARMM36 force field with the TIP3P-FB and TIP4P-FB water models

The CHARMM36 force field for lipids is widely used in simulations of lipid bilayers. The CHARMM family of force fields were developed for use with the mTIP3P water model. This water model has an anomalously high dielectric constant and low viscosity, which limits its accuracy in the calculation of quantities like permeability coefficients. The TIP3P-FB and TIP4P-FB water models are more accurate in terms of the dielectric constant and transport properties, which could allow more accurate simulations of systems containing water and lipids. To test whether the CHARMM36 lipid force field is compatible with the TIP3P-FB and TIP4P-FB water models, we have performed simulations of 1,2-dipalmitoyl-sn-glycero-3-phosphocholine and 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine bilayers. The calculated headgroup area, compressibility, order parameters, and X-ray form factors are in good agreement with the experimental values, indicating that these improved water models can be used with the CHARMM36 lipid force field without modification when calculating membrane physical properties. The water permeability predicted by these models is significantly different; the mTIP3P-model diffusion in solution and at the lipid-water interface is anomalously fast due to the spuriously low viscosity of mTIP3P-model water, but the potential of mean force of permeation is higher for the TIP3P-FB and TIP4P-FB models due to their high excess chemical potentials. As a result, the rates of water permeation calculated the FB water models are slower than the experimental value by a factor of 15-17, while simulations with the mTIP3P model only underestimate the water permeability by a factor of 3.

All n-3 PUFA are not the same: MD simulations reveal differences in membrane organization for EPA, DHA and DPA

Eicosapentaenoic (EPA, 20:5), docosahexaenoic (DHA, 22:6) and docosapentaenoic (DPA, 22:5) acids are omega-3 polyunsaturated fatty acids (n-3 PUFA) obtained from dietary consumption of fish oils that potentially alleviate the symptoms of a range of chronic diseases. We focus here on the plasma membrane as a site of action and investigate how they affect molecular organization when taken up into a phospholipid. All atom MD simulations were performed to compare 1-stearoyl-2-eicosapentaenoylphosphatylcholine (EPA-PC, 18:0-20:5PC), 1-stearoyl-2-docosahexaenoylphosphatylcholine (DHA-PC, 18:0-22:6PC), 1-stearoyl-2-docosapentaenoylphosphatylcholine (DPA-PC, 18:0-22:5PC) and, as a monounsaturated control, 1-stearoyl-2-oleoylphosphatidylcholine (OA-PC, 18:0-18:1PC) bilayers. They were run in the absence and presence of 20mol% cholesterol. Multiple double bonds confer high disorder on all three n-3 PUFA. The different number of double bonds and chain length for each n-3 PUFA moderates the reduction in membrane order exerted (compared to OA-PC, S¯CD=0.152). EPA-PC (S¯CD=0.131) is most disordered, while DPA-PC (S¯CD=0.140) is least disordered. DHA-PC (S¯CD=0.139) is, within uncertainty, the same as DPA-PC. Following the addition of cholesterol, order in EPA-PC (S¯CD=0.169), DHA-PC (S¯CD=0.178) and DPA-PC (S¯CD=0.182) is increased less than in OA-PC (S¯CD=0.214). The high disorder of n-3 PUFA is responsible, preventing the n-3 PUFA-containing phospholipids from packing as close to the rigid sterol as the monounsaturated control. Our findings establish that EPA, DHA and DPA are not equivalent in their interactions within membranes, which possibly contributes to differences in clinical efficacy.

Immediate response mechanisms of Gram-negative solvent-tolerant bacteria to cope with environmental stress: cis-trans isomerization of unsaturated fatty acids and outer membrane vesicle secretion

Bacteria have evolved an array of adaptive mechanisms enabling them to survive and grow in the presence of different environmental stresses. These mechanisms include either modifications of the membrane or changes in the overall energy status, cell morphology, and cell surface properties. Long-term adaptations are dependent on transcriptional regulation, the induction of anabolic pathways, and cell growth. However, to survive sudden environmental changes, bacterial short-term responses are essential to keep the cells alive after the occurrence of an environmental stress factor such as heat shock or the presence of toxic organic solvents. Thus far, two main short-term responses are known. On the one hand, a fast isomerization of cis into trans unsaturated fatty leads to a quick rigidification of the cell membrane, a mechanism known in some genera of Gram-negative bacteria. On the other hand, a fast, effective, and ubiquitously present countermeasure is the release of outer membrane vesicles (OMVs) from the cell surface leading to a rapid increase in cell surface hydrophobicity and finally to the formation of cell aggregates and biofilms. These immediate response mechanisms just allow the bacteria to stay physiologically active and to employ long-term responses to assure viability upon changing environmental conditions. Here, we provide insight into the two aforementioned rapid adaptive mechanisms affecting ultimately the cell envelope of Gram-negative bacteria.

Acemetacin–phosphatidylcholine interactions are determined by the drug ionization state

Complementary biophysical techniques depicted the differential effects of acemetacin ionic forms on phosphatidylcholine bilayers.

Effect of ergosterol on the fungal membrane properties. All-atom and coarse-grained molecular dynamics study

Cell membranes are complex multicomponent systems consisting of thousands of different lipids with numerous embedded membrane proteins and many types of sterols. We used all-atom and coarse-grained molecular dynamics simulations to study the structural and dynamical properties of phospholipid bilayers containing four types of phospholipids and different amount of ergosterol, main sterol component in the fungal membranes. To characterize the influence of ergosterol on the membrane properties we analyzed the surface area per lipid, bilayer thickness, area compressibility modulus, mass density profiles, deuterium order parameters, and lateral diffusion coefficients. The presence of ergosterol induces the ordering of lipids leading to their denser packing, to reducing the lateral diffusion of lipids and lipid surface area, to increasing the thickness of bilayer and compressibility modulus. In addition, we evaluated each calculated property between the two simulation methods.

Molecular Dynamics Simulations of Ceramide and Ceramide-Phosphatidylcholine Bilayers

Recent studies in lipid raft formation and stratum corneum permeability have focused on the role of ceramides (CER). In this study, we use the all-atom CHARMM36 (C36) force field to simulate bilayers using N-palmitoylsphingosine (CER16) or α-hydroxy-N-stearoyl phytosphingosine (CER[AP]) in 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) or 1-palmitoyl-2-oleoylphosphatidylcholine (POPC), which serve as general membrane models. Conditions are replicated from experimental studies for comparison purposes, and concentration (X_CER) is varied to probe the effect of CER on these systems. Comparisons with experiment based on deuterium order parameters and bilayer thickness demonstrate good agreement, thus supporting further use of the C36 force field. CER concentration is shown to have a profound effect on nearly all membrane properties including surface area per lipid, chain order and tilt, area compressibility moduli, bilayer thickness, hydrogen bonding, and lipid clustering. Hydrogen bonding in particular can significantly affect other membrane properties and can even encourage transition to a gel phase. Despite CER's tendency to condense the membrane, an expansion of CER lipids with increasing X_CER is possible depending on how the balance between various hydrogen-bond pairs and lipid clustering is perturbed. Based on gel phase transitions, support is given for phytosphingosine's role as a hydrogen-bond bridge between sphingosine ordered domains in the stratum corneum.

On the Calculation of Acyl Chain Order Parameters from Lipid Simulations

For molecular dynamics simulations of biological membrane systems to live up to the potential of providing accurate atomic level detail into membrane properties and functions, it is essential that the force fields used to model such systems are as accurate as possible. One membrane property that is often used to assess force field accuracy is the carbon-hydrogen (or carbon-deuterium) order parameters of the lipid tails, which can be accurately measured using experimental NMR techniques. There are a variety of analysis tools available to calculate these order parameters from simulations and it is essential that these computational tools work correctly to ensure the accurate assessment of the simulation force fields. In this work we compare many of these computational tools for calculating the order parameters of POPC membranes. While tools that work on all-atom systems and tools that work on saturated lipid tails in general work extremely well, we demonstrate that the majority of the tested tools that calculate the order parameters for unsaturated united-atom lipid tails do so incorrectly. We identify tools that do perform accurate calculations and include one such program with this work, enabling rapid and accurate calculation of united-atom lipid order parameters. Furthermore, we discuss cases in which it is nontrivial to appropriately predict the unsaturated carbon order parameters in united-atom systems. Finally, we examine order parameter splitting for carbon 2 in sn-2 lipid chains, demonstrating substantial deviations from experimental values in several all-atom and united-atom lipid force fields.

Concepts and Methods of Solid-State NMR Spectroscopy Applied to Biomembranes

Concepts of solid-state NMR spectroscopy and applications to fluid membranes are reviewed in this paper. Membrane lipids with ²H-labeled acyl chains or polar head groups are studied using ²H NMR to yield knowledge of their atomistic structures in relation to equilibrium properties. This review demonstrates the principles and applications of solid-state NMR by unifying dipolar and quadrupolar interactions and highlights the unique features offered by solid-state ²H NMR with experimental illustrations. For randomly oriented multilamellar lipids or aligned membranes, solid-state ²H NMR enables direct measurement of residual quadrupolar couplings (RQCs) due to individual C-²H-labeled segments. The distribution of RQC values gives nearly complete profiles of the segmental order parameters S_CD⁽ⁱ⁾ as a function of acyl segment position (i). Alternatively, one can measure residual dipolar couplings (RDCs) for natural abundance lipid samples to obtain segmental S_CH order parameters. A theoretical mean-torque model provides acyl-packing profiles representing the cumulative chain extension along the normal to the aqueous interface. Equilibrium structural properties of fluid bilayers and various thermodynamic quantities can then be calculated, which describe the interactions with cholesterol, detergents, peptides, and integral membrane proteins and formation of lipid rafts. One can also obtain direct information for membrane-bound peptides or proteins by measuring RDCs using magic-angle spinning (MAS) in combination with dipolar recoupling methods. Solid-state NMR methods have been extensively applied to characterize model membranes and membrane-bound peptides and proteins, giving unique information on their conformations, orientations, and interactions in the natural liquid-crystalline state.

Drude Polarizable Force Field for Molecular Dynamics Simulations of Saturated and Unsaturated Zwitterionic Lipids

Additive force fields are designed to account for induced electronic polarization in a mean-field average way, using effective empirical fixed charges. The limitation of this approximation is cause for serious concerns, particularly in the case of lipid membranes, where the molecular environment undergoes dramatic variations over microscopic length scales. A polarizable force field based on the classical Drude oscillator offers a practical and computationally efficient framework for an improved representation of electrostatic interactions in molecular simulations. Building on the first-generation Drude polarizable force field for the dipalmitoylphosphatidylcholine 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) molecule, the present effort was undertaken to improve this initial model and expand the force field to a wider range of phospholipid molecules. New lipids parametrized include dimyristoylphosphatidylcholine (DMPC), dilauroylphosphatidylcholine (DLPC), 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC), 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC), dipalmitoylphosphatidylethanolamine (DPPE), 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoethanolamine (POPE), and 1,2-dioleoyl-sn-glycero-3-phosphoethanolamine (DOPE). The iterative optimization protocol employed in this effort led to lipid models that achieve a good balance between reproducing quantum mechanical data on model compound representative of phospholipids and reproducing a range of experimental condensed phase properties of bilayers. A parametrization strategy based on a restrained ensemble-maximum entropy methodology was used to help accurately match the experimental NMR order parameters in the polar headgroup region. All the parameters were developed to be compatible with the remainder of the Drude polarizable force field, which includes water, ions, proteins, DNA, and selected carbohydrates.

Unperturbed hydrocarbon chains and liquid phase bilayer lipid chains: a computer simulation study

In this work, the properties of saturated and unsaturated fatty acid acyl chains 16:0, 18:0, 18:1(n-9)cis, 18:2(n-6)cis, 18:3(n-3)cis, 18:4(n-3)cis, 18:5(n-3)cis, 20:4(n-6)cis, 20:5(n-3)cis and 22:6(n-3)cis in a bilayer liquid crystalline state and similar hydrocarbon chains (with CH[Formula: see text] terminal groups instead of C=O groups) in the unperturbed state characterised by a lack of long-range interaction were investigated. The unperturbed hydrocarbon chains were modelled by Monte Carlo simulations at temperature [Formula: see text] K; sixteen fully hydrated homogeneous liquid crystalline phosphatidylcholine bilayers containing these chains were studied by molecular dynamics simulations at the same temperature. To eliminate effects of the simulation parameters, the molecular dynamics and Monte Carlo simulations were carried out using the same structural data and force field coefficients. From these computer simulations, the average distances between terminal carbon atoms of the chains (end-to-end distances) were calculated and compared. The trends in the end-to-end distances obtained for the unperturbed chains were found to be qualitatively similar to those obtained for the same lipid chains in the bilayers. So, for understanding of a number of processes in biological membranes (e.g., changes in fatty acid composition caused by environmental changes such as temperature and pressure), it is possible to use, at least as a first approximation, the relationships between the structure and properties for unperturbed or isolated hydrocarbon chains.

Molecular dynamics simulations of T-2410 and T-2429 HIV fusion inhibitors interacting with model membranes: Insight into peptide behavior, structure and dynamics

T-2410 and T-2429 are HIV fusion inhibitor peptides (FI) designed to present a higher efficiency even against HIV strains that developed resistance against other FIs. Similar peptides were shown to interact with model membranes both in the liquid disordered phase and in the liquid ordered state. Those results indicated that such interaction is important to function and could be correlated with their effectiveness. Extensive molecular dynamics simulations were carried out to investigate the interactions between both T-2410 and T-2429 with bilayers of pure 1-palmitoyl-2-oleoyl-phosphatidylcholine (POPC) and a mixture of POPC/cholesterol (Chol) (1:1). It was observed that both peptides interact strongly with both membrane systems, especially with the POPC/Chol systems, where these peptides show the highest number of H-bonds observed so far. T-2410 and T-2429 showed higher extent of interaction with bilayers when compared to T-20 or T-1249 in previous studies. This is most notable in POPC/Chol membranes where, although able to form H-bonds with Chol, they do so to a lesser extent than T-1249 does, the latter being the only FI peptide so far that was observed to form H-bonds with Chol. This behavior suggests that interaction of FI peptides with rigid Chol rich membranes may not be as dependent from peptide/Chol H-bond formation as previous results of T-1249 behavior led to believe. As in other similar peptides, the higher ability to interact with membranes shown by T-2410 and T2429 is probably correlated with its higher inhibitory efficiency.

Effect of cholesterol on the molecular structure and transitions in a clinical-grade lung surfactant extract

The lipid-protein film covering the interface of the lung alveolar in mammals is vital for proper lung function and its deficiency is related to a range of diseases. Here we present a molecular-level characterization of a clinical-grade porcine lung surfactant extract using a multitechnique approach consisting of [Formula: see text]-[Formula: see text] solid-state nuclear magnetic spectroscopy, small- and wide-angle X-ray scattering, and mass spectrometry. The detailed characterization presented for reconstituted membranes of a lung extract demonstrates that the molecular structure of lung surfactant strongly depends on the concentration of cholesterol. If cholesterol makes up about 11% of the total dry weight of lung surfactant, the surfactant extract adopts a single liquid-ordered lamellar phase, [Formula: see text], at physiological temperatures. This [Formula: see text] phase gradually changes into a liquid-disordered lamellar phase, [Formula: see text], when the temperature is increased by a few degrees. In the absence of cholesterol the system segregates into one lamellar gel phase and one [Formula: see text] phase. Remarkably, it was possible to measure a large set of order parameter magnitudes [Formula: see text] from the liquid-disordered and -ordered lamellar phases and assign them to specific C-H bonds of the phospholipids in the biological extract with no use of isotopic labeling. These findings with molecular details on lung surfactant mixtures together with the presented NMR methodology may guide further development of pulmonary surfactant pharmaceuticals that better mimic the physiological self-assembly compositions for treatment of pathological states such as respiratory distress syndrome.

Assessing gastric toxicity of xanthone derivatives of anti-inflammatory activity using simulation and experimental approaches

Xanthones are tricyclic compounds of natural or synthetic origin exhibiting a broad spectrum of therapeutic activities. Three synthetic xanthone derivatives (KS1, KS2, and KS3) with properties typical for nonsteroidal anti-inflammatory drugs (NSAID) were objects of the presented model study. NSAIDs are in common use however; several of them exhibit gastric toxicity predominantly resulting from their direct interactions with the outermost lipid layer of the gastric mucosa that impair its hydrophobic barrier property. Among the studied xanthones, gastric toxicity of only KS2 has been determined in previous pharmacological studies, and it is low. In this study, carried out using X-ray diffraction and computer simulation, a palmitoyloleoylphosphatidylcholine-cholesterol bilayer (POPC-Chol) was used as a model of a hydrophobic layer of lipids protecting gastric mucosa as POPC and Chol are the main lipids in human mucus. X-ray diffraction data were used to validate the computer model. The aim of the study was to assess potential gastric toxicity of the xanthones by analysing their atomic level interactions with lipids, ions, and water in the lipid bilayer and their effect on the bilayer physicochemical properties. The results show that xanthones have small effect on the bilayer properties except for its rigidity whereas their interactions with water, ions, and lipids depend on their protonation state and for a given state, are similar for all the xanthones. As gastric toxicity of KS2 is low, based on MD simulations one can predict that toxicity of KS1 and KS3 is also low.