Sensing hydration and behavior of pyrene in POPC and POPC/cholesterol bilayers: A molecular dynamics study

Effects of Nitro-Oxidative Stress on Biomolecules: Part 1-Non-Reactive Molecular Dynamics Simulations

Plasma medicine, or the biomedical application of cold atmospheric plasma (CAP), is an expanding field within plasma research. CAP has demonstrated remarkable versatility in diverse biological applications, including cancer treatment, wound healing, microorganism inactivation, and skin disease therapy. However, the precise mechanisms underlying the effects of CAP remain incompletely understood. The therapeutic effects of CAP are largely attributed to the generation of reactive oxygen and nitrogen species (RONS), which play a crucial role in the biological responses induced by CAP. Specifically, RONS produced during CAP treatment have the ability to chemically modify cell membranes and membrane proteins, causing nitro-oxidative stress, thereby leading to changes in membrane permeability and disruption of cellular processes. To gain atomic-level insights into these interactions, non-reactive molecular dynamics (MD) simulations have emerged as a valuable tool. These simulations facilitate the examination of larger-scale system dynamics, including protein-protein and protein-membrane interactions. In this comprehensive review, we focus on the applications of non-reactive MD simulations in studying the effects of CAP on cellular components and interactions at the atomic level, providing a detailed overview of the potential of CAP in medicine. We also review the results of other MD studies that are not related to plasma medicine but explore the effects of nitro-oxidative stress on cellular components and are therefore important for a broader understanding of the underlying processes.

Fluorescent Probes cis- and trans-Parinaric Acids in Fluid and Gel Lipid Bilayers: A Molecular Dynamics Study

Fluorescence probes are indispensable tools in biochemical and biophysical membrane studies. Most of them possess extrinsic fluorophores, which often constitute a source of uncertainty and potential perturbation to the host system. In this regard, the few available intrinsically fluorescent membrane probes acquire increased importance. Among them, cis- and trans-parinaric acids (c-PnA and t-PnA, respectively) stand out as probes of membrane order and dynamics. These two compounds are long-chained fatty acids, differing solely in the configurations of two double bonds of their conjugated tetraene fluorophore. In this work, we employed all-atom and coarse-grained molecular dynamics simulations to study the behavior of c-PnA and t-PnA in lipid bilayers of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) and 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC), representative of the liquid disordered and solid ordered lipid phases, respectively. All-atom simulations indicate that the two probes show similar location and orientation in the simulated systems, with the carboxylate facing the water/lipid interface and the tail spanning the membrane leaflet. The two probes establish interactions with the solvent and lipids to a similar degree in POPC. However, the almost linear t-PnA molecules have tighter lipid packing around them, especially in DPPC, where they also interact more with positively charged lipid choline groups. Probably for these reasons, while both probes show similar partition (assessed from computed free energy profiles across bilayers) to POPC, t-PnA clearly partitions more extensively than c-PnA to the gel phase. t-PnA also displays more hindered fluorophore rotation, especially in DPPC. Our results agree very well with experimental fluorescence data from the literature and allow deeper understanding of the behavior of these two reporters of membrane organization.

Effect of cholesterol on the membrane partitioning dynamics of hepatitis A virus-2B peptide

Understanding viral peptide detection and partitioning and the subsequent host membrane composition-based response is essential for gaining insights into the viral mechanism. Here, we probe the crucial role of the presence of membrane lipid packing defects, depending on the membrane composition, in allowing the viral peptide belonging to C-terminal Hepatitis A Virus-2B (HAV-2B) to detect, attach and subsequently partition into host cell membrane mimics. Using molecular dynamics simulations, we conclusively show that the hydrophobic residues in the viral peptide detect transiently present lipid packing defects, insert themselves into such defects, form anchor points and facilitate the partitioning of the peptide, thereby inducing membrane disruption. We also show that the presence of cholesterol significantly alters such lipid packing defects, both in size and in number, thus mitigating the partitioning of the membrane active viral peptide into cholesterol-rich membranes. Our results are in excellent agreement with previously published experimental data and further explain the role of lipid defects in understanding such data. These results show differential ways in which the presence and absence of cholesterol can alter the permeability of the host membranes to the membrane active peptide component of HAV-2B virus, via lipid packing defects, and can possibly be a part of the general membrane detection mechanism for viroporins.

Modulation of Phospholipid Bilayer Properties by Simvastatin

Simvastatin (Zocor) is one of the most prescribed drugs for reducing high cholesterol. Although simvastatin is ingested in its inactive lactone form, it is converted to its active dihydroxyheptanoate form by carboxylesterases in the liver. The dihydroxyheptanoate form can also be converted back to its original lactone form. Unfortunately, some of the side effects associated with the intake of simvastatin and other lipophilic statins at higher doses include statin-associated myopathy (SAM) and, in more severe cases, kidney failure. While the cause of SAM is unknown, it is hypothesized that these side effects are dependent on the localization of statins in lipid bilayers and their impact on bilayer properties. In this work, we carry out all-atom molecular dynamics simulations on both the lactone and dihydroxyheptanoate forms of simvastatin (termed "SN" and "SA", respectively) with a pure 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) lipid bilayer and a POPC/cholesterol (30 mol %) binary mixture as membrane models. Additional simulations were carried out with multiple simvastatin molecules to mimic in vitro conditions that produced pleiotropic effects. Both SN and SA spontaneously diffused into the lipid bilayer, and a longer simulation time of 4 μs was needed for the complete incorporation of multiple SAs into the bilayer. By constructing potential mean force and electron density profiles, we find that SN localizes deeper within the hydrophobic interior of the bilayer and that SA has a greater tendency to form hydrogen-bonding interactions with neighboring water molecules and lipid headgroups. For the pure POPC bilayer, both SN and SA increase membrane order, while membrane fluidity increases for the POPC/cholesterol bilayer.

Effects of cholesterol on chlorzoxazone translocation across POPC bilayer

Cholesterol plays a crucial role in modulating the physicochemical properties of membranes, thus influencing the membrane transport of drugs. In this paper, the effects caused by cholesterol on the membrane transport of chlorzoxazone (CZX), a centrally acting muscle relaxant drug, were probed through molecular dynamics simulations. POPC was selected as the model lipid, and three different cholesterol concentrations (0%, 20%, and 50% CHOL) were considered. The outcomes reveal that the area per lipid of POPC decreases and the order parameter increases with enhanced concentration of CHOL. CZX prefers to localize at the interface between the headgroup region and the hydrophobic tail region of POPC, and the main energy barrier occurs in the hydrophobic region. The impact of CHOL on the free energy profile is correlated with concentration: low concentration facilitates CZX permeation, while high concentration hinders CZX permeation. Our findings coincide with experimental results, enhancing the mechanism understanding of how drug molecules are transported through membranes in the presence of CHOL. • The effects caused by cholesterol (CHOL) on the membrane transport of chlorzoxazone (CZX) were studied. • Low CHOL concentration facilitates CZX permeation, while high concentration hinders CZX permeation. • Our findings improve the mechanism understanding of CHOL effects on CZX translocation across membrane.

Oxidation degree of a cell membrane model and its response to structural changes, a coarse-grained molecular dynamics approach

Oxidative stress plays an essential role in the regulation of vital processes in living organisms. Reactive oxygen species can react chemically with the constituents of the cells leading to irreversible damage. The first structure of the cell in contact with the environment that surrounds it is the membrane, which protects it and allows the exchange of substances. Some signals manifest when the components of a bilayer are undergoing oxidation, like an increase in the lipid area, decrease in the thickness of the bilayer, and exchange of the oxidized groups toward the bilayer surface. In this investigation, a molecular dynamics simulation was done on a set of Dioleoylphosphatidylcholine membranes with different percentage of oxidized lipids, in order to observe the effect of the oxidation degree on the membrane structure. It was found that, as higher the concentration of oxidized lipids is, the larger the damage of the membrane. This is reflected in the increase in the lipid area and the decrease in the thickness and membrane packing. Also, it was observed that hydrophobicity inside the membrane decreases as the oxidation percentage increases.Communicated by Ramaswamy H. Sarma.

Time-Dependent Lipid Dynamics, Organization and Peptide-Lipid Interaction in Phospholipid Bilayers with Incorporated β-Amyloid Oligomers

Nonfibrillar β-amyloid (Aβ) oligomers are considered as major neurotoxic species in the pathology of Alzheimer's disease. The presence of Aβ oligomers was shown to cause membrane disruptions in a broad range of model systems. However, the molecular basis of such a disruption process remains unknown. We previously demonstrated that membrane-incorporated 40-residue Aβ (Aβ₄₀) oligomers could form coaggregates with phospholipids. This process occurred more rapidly than the fibrillization of Aβ₄₀ and led to more severe membrane disruption. The present study probes the time-dependent changes in lipid dynamics, bilayer structures, and peptide-lipid interactions along the time course of the oligomer-induced membrane disruption, using solid-state NMR spectroscopy. Our results suggest the presence of certain intermediate states with phospholipid molecules entering the C-terminal hydrogen-bonding networks of the Aβ₄₀ oligomeric cores. This work provides insights on the molecular mechanisms of Aβ₄₀-oligomer-induced membrane disruption.

Affinity for the Interface Underpins Potency of Antibodies Operating In Membrane Environments

The contribution of membrane interfacial interactions to recognition of membrane-embedded antigens by antibodies is currently unclear. This report demonstrates the optimization of this type of antibodies via chemical modification of regions near the membrane but not directly involved in the recognition of the epitope. Using the HIV-1 antibody 10E8 as a model, linear and polycyclic synthetic aromatic compounds are introduced at selected sites. Molecular dynamics simulations predict the favorable interactions of these synthetic compounds with the viral lipid membrane, where the epitope of the HIV-1 glycoprotein Env is located. Chemical modification of 10E8 with aromatic acetamides facilitates the productive and specific recognition of the native antigen, partially buried in the crowded environment of the viral membrane, resulting in a dramatic increase of its capacity to block viral infection. These observations support the harnessing of interfacial affinity through site-selective chemical modification to optimize the function of antibodies that target membrane-proximal epitopes.

The Secret Lives of Fluorescent Membrane Probes as Revealed by Molecular Dynamics Simulations

Fluorescent probes have been employed for more than half a century to study the structure and dynamics of model and biological membranes, using spectroscopic and/or microscopic experimental approaches. While their utilization has led to tremendous progress in our knowledge of membrane biophysics and physiology, in some respects the behavior of bilayer-inserted membrane probes has long remained inscrutable. The location, orientation and interaction of fluorophores with lipid and/or water molecules are often not well known, and they are crucial for understanding what the probe is actually reporting. Moreover, because the probe is an extraneous inclusion, it may perturb the properties of the host membrane system, altering the very properties it is supposed to measure. For these reasons, the need for independent methodologies to assess the behavior of bilayer-inserted fluorescence probes has been recognized for a long time. Because of recent improvements in computational tools, molecular dynamics (MD) simulations have become a popular means of obtaining this important information. The present review addresses MD studies of all major classes of fluorescent membrane probes, focusing in the period between 2011 and 2020, during which such work has undergone a dramatic surge in both the number of studies and the variety of probes and properties accessed.

Raman Spectroscopy Study of Curvature-Mediated Lipid Packing and Sorting in Single Lipid Vesicles

Cellular plasma membrane deformability and stability is important in a range of biological processes. Changes in local curvature of the membrane affect the lateral movement of lipids, affecting the biophysical properties of the membrane. An integrated holographic optical tweezers and Raman microscope was used to investigate the effect of curvature gradients induced by optically stretching individual giant unilamellar vesicles (GUVs) on lipid packing and lateral segregation of cholesterol in the bilayer. The spatially resolved Raman analysis enabled detection of induced phase separation and changes in lipid ordering in individual GUVs. Using deuterated cholesterol, the changes in lipid ordering and phase separation were linked to lateral sorting of cholesterol in the stretched GUVs. Stretching the GUVs in the range of elongation factors 1-1.3 led to an overall decrease in cholesterol concentration at the edges compared to the center of stretched GUVs. The Raman spectroscopy results were consistent with a model of the bilayer accounting for cholesterol sorting in both bilayer leaflets, with a compositional asymmetry of 0.63 ± 0.04 in favor of the outer leaflet. The results demonstrate the potential of the integrated holographic optical tweezers-Raman technique to induce deformations to individual lipid vesicles and to simultaneously provide quantitative and spatially resolved molecular information. Future studies can extend to include more realistic models of cell membranes and potentially live cells.

Physical properties of model biological lipid bilayers: insights from all-atom molecular dynamics simulations

The physical properties of lipid bilayers are sensitive to the specific type and composition of the lipids that make up the many different types of cell membranes. Studying model bilayers of representative heterogeneous compositions can provide key insights into membrane functionality. In this work, we use atomistic molecular dynamics simulations to characterize key properties in a number of bilayer membranes of varying composition. We first examine basic properties, such as lipid area, volume, and bilayer thickness, of simple, homogeneous bilayers comprising several lipid types, which are prevalent in biological membranes. Such lipids are then used in simulations of heterogeneous systems representative of bacterial, mammalian, and cancer membranes. Our analysis is especially focused on depth-dependent, transmembrane profiles; in particular, we calculate lateral pressure and dipole potential profiles, two fundamental properties which play key roles in a large number of biological functions.

Interaction of NBD-labelled fatty amines with liquid-ordered membranes: a combined molecular dynamics simulation and fluorescence spectroscopy study

A complete homologous series of fluorescent 7-nitrobenz-2-oxa-1,3-diazol-4-yl-(NBD) labelled fatty amines of varying alkyl chain lengths, NBD-Cn, inserted in 1-palmitoyl, 2-oleoyl-sn-glycero-3-phosphocholine (POPC) or N-palmitoyl sphingomyelin (SpM) bilayers, with 50 mol% and 40 mol% cholesterol (Chol), respectively, was studied using atomistic molecular dynamics simulations. For all amphiphiles in both bilayers, the NBD fluorophore locates at the interface, in a more external position than that previously observed for pure POPC bilayers. This shallower location of the NBD group agrees with the lower fluorescent quantum yield, shorter fluorescence lifetime, and higher ionisation constants (smaller pKa) determined experimentally. The more external location is also consistent with the changes measured in steady-state fluorescence anisotropy from POPC to POPC/Chol (1 : 1) vesicles. Accordingly, the equilibrium location of the NBD group within the various bilayers is mainly dictated by bilayer compositions, and is mostly unaffected by the length of the attached alkyl chain. Similarly to the behaviour observed in POPC bilayers, the longer-chained NBD-Cn amphiphiles show significant mass density near the mixed bilayers' midplanes, and the alkyl chains of the longer derivatives, mainly NBD-C16, penetrate the opposite bilayer leaflet to some extent. However, this effect is quantitatively less pronounced in these ordered bilayers than in POPC. Similarly to POPC bilayers, the effects of these amphiphiles on the structure and dynamics of the host lipid were found to be relatively mild, in comparison with acyl-chain phospholipid analogues.

Probing Microscopic Orientation in Membranes by Linear Dichroism

The cell membrane is an ordered environment, which anisotropically affects the structure and interactions of all of its molecules. Monitoring membrane orientation at a local level is rather challenging but could reward crucial information on protein conformation and interactions in the lipid bilayer. We monitored local lipid ordering changes upon varying the cholesterol concentration using polarized light spectroscopy and pyrene as a membrane probe. Pyrene, with a shape intermediate between a disc and a rod, can detect microscopic orientation variations at the level of its size. The global membrane orientation was determined using curcumin, a probe with nonoverlapping absorption relative to that of pyrene. While the macroscopic orientation of a liquid-phase bilayer decreases with increasing cholesterol concentration, the local orientation is improved. Pyrene is found to be sensitive to the local effects induced by cholesterol and temperature on the bilayer. Disentangling local and global orientation effects in membranes could provide new insights into functionally significant interactions of membrane proteins.

Two-Dimensional Potentials of Mean Force of Nile Red in Intact and Damaged Model Bilayers. Application to Calculations of Fluorescence Spectra

Fluorescent dyes revolutionized and expanded our understanding of biological membranes. The interpretation of experimental fluorescence data in terms of membrane structure, however, requires detailed information about the molecular environment of the dyes. Nile red is a fluorescent molecule whose excitation and emission maxima depend on the polarity of the solvent. It is mainly used as a probe to study lipid microenvironments, for example in imaging the progression of damage to the myelin sheath in multiple sclerosis. In this study, we determine the position and orientation of Nile red in lipid bilayers by calculating two-dimensional Potential of Mean Force (2D-PMF) profiles in a defect-free 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) bilayer and in damaged bilayers containing two mixtures of the oxidized lipid 1-palmitoyl-2-(9'-oxo-nonanoyl)-sn-glycero-3-phosphocholine and POPC. From 2D-PMF simulations we obtain positions and orientations of Nile Red corresponding to the minimum on the binding free energy surface in three different membrane environments with increasing amounts of water, mimicking damage in biological tissue. Using representative snapshots from the simulations, we use combined quantum mechanical/molecular mechanical (QM/MM) models to calculate the emission spectrum of Nile red as a function of its local solvation environment. The results of QM and QM/MM computations are in qualitative agreement with the experimentally observed shift in fluorescence for the dye moving from aqueous solution to the more hydrophobic environment of the lipid interiors. The range of the conformation dependent values of the computed absorption-emission spectra and the lack of solvent relaxation effects in the QM/MM calculations made it challenging to delineate specific differences between the intact and damaged bilayers.

From molecular modelling to photophysics of neutral oligo- and polyfluorenes incorporated into phospholipid bilayers

The combination of various experimental techniques with theoretical simulations has allowed elucidation of the mode of incorporation of fluorene based derivatives into phospholipid bilayers. Molecular dynamics (MD) simulations on a fully hydrated 1,2-dimyristoyl-sn-glycero-3-phosphatidylcholine (DMPC) bilayer, with benzene (B), biphenyl (BP), fluorene (F) and tri-(9,9-di-n-octylfluorenyl-2,7-diyl), TF, have provided insights into the topography of these molecules when they are present in the phospholipid bilayer, and suggest marked differences between the behavior of the small molecules and the oligomer. Further information on the interaction of neutral fluorenes within the phospholipid bilayer was obtained by an infrared (IR) spectroscopic study of films of DMPC and of the phospholipid with PFO deuterated specifically on its alkyl chains (DMPC-PFO-d34). This was complemented by measurements of the effect of F, TF and two neutral polymers: polyfluorene poly(9,9-di-n-octylfluorenyl-2,7-diyl), PFO, and poly(9,9-di-n-dodecylfluorenyl-2,7-diyl), PFD, on the phospholipid phase transition temperature using differential scanning calorimetry (DSC). Changes in liposome size upon addition of F and PFO were followed by dynamic light scattering. In addition, the spectroscopic properties of F, TF, PFO and PFD solubilised in DMPC liposomes (absorption, steady-state and time-resolved fluorescence) were compared with those of the same probes in typical organic solvents (chloroform, cyclohexane and ethanol). Combining the insight from MD simulations with the results at the molecular level from the various experimental techniques suggests that while the small molecules have a tendency to be located in the phospholipid head group region, the polymers are incorporated within the lipid bilayers, with the backbone predominantly orthogonal to the phospholipid alkyl chains and with interdigitation of them and the PFO alkyl chains.

Interactions of Lipid Membranes with Fibrillar Protein Aggregates

Amyloid fibrils are an intriguing class of protein aggregates with distinct physicochemical, structural and morphological properties. They display peculiar membrane-binding behavior, thus adding complexity to the problem of protein-lipid interactions. The consensus that emerged during the past decade is that amyloid cytotoxicity arises from a continuum of cross-β-sheet assemblies including mature fibrils. Based on literature survey and our own data, in this chapter we address several aspects of fibril-lipid interactions, including (i) the effects of amyloid assemblies on molecular organization of lipid bilayer; (ii) competition between fibrillar and monomeric membrane-associating proteins for binding to the lipid surface; and (iii) the effects of lipids on the structural morphology of fibrillar aggregates. To illustrate some of the processes occurring in fibril-lipid systems, we present and analyze fluorescence data reporting on lipid bilayer interactions with fibrillar lysozyme and with the N-terminal 83-residue fragment of amyloidogenic mutant apolipoprotein A-I, 1-83/G26R/W@8. The results help understand possible mechanisms of interaction and mutual remodeling of amyloid fibers and lipid membranes, which may contribute to amyloid cytotoxicity.

Cholesterol, sphingolipids, and glycolipids: what do we know about their role in raft-like membranes?

Lipids rafts are considered to be functional nanoscale membrane domains enriched in cholesterol and sphingolipids, characteristic in particular of the external leaflet of cell membranes. Lipids, together with membrane-associated proteins, are therefore considered to form nanoscale units with potential specific functions. Although the understanding of the structure of rafts in living cells is quite limited, the possible functions of rafts are widely discussed in the literature, highlighting their importance in cellular functions. In this review, we discuss the understanding of rafts that has emerged based on recent atomistic and coarse-grained molecular dynamics simulation studies on the key lipid raft components, which include cholesterol, sphingolipids, glycolipids, and the proteins interacting with these classes of lipids. The simulation results are compared to experiments when possible.

Molecular dynamics simulations of depth distribution of spin-labeled phospholipids within lipid bilayer

Spin-labeled lipids are commonly used as fluorescence quenchers in studies of membrane penetration of dye-labeled proteins and peptides using depth-dependent quenching. Accurate calculations of depth of the fluorophore rely on the use of several spin labels placed in the membrane at various positions. The depth of the quenchers (spin probes) has to be determined independently; however, experimental determination of transverse distributions of spin probe depths is difficult. In this Article, we use molecular dynamics (MD) simulations to study the membrane behavior and depth distributions of spin-labeled phospholipids in a 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) bilayer. To probe different depths within the bilayer, a series containing five Doxyl-labeled lipids (n-Doxyl PC) has been studied, in which a spin moiety was covalently attached to nth carbon atoms (where n = 5, 7, 10, 12, and 14) of the sn-2 stearoyl chain of the host phospholipid. Our results demonstrate that the chain-attached spin labels are broadly distributed across the model membrane and their environment is characterized by a high degree of mobility and structural heterogeneity. Despite the high thermal disorder, the depth distributions of the Doxyl labels were found to correlate well with their attachment positions, indicating that the distribution of the spin label within the model membrane is dictated by the depth of the nth lipid carbon atom and not by intrinsic properties of the label. In contrast, a much broader and heterogeneous distribution was observed for a headgroup-attached Tempo spin label of Tempo-PC lipids. MD simulations reveal that, due to the hydrophobic nature, a Tempo moiety favors partitioning from the headgroup region deeper into the membrane. Depending on the concentration of Tempo-PC lipids, the probable depth of the Tempo moiety could span a range from 14.4 to 18.2 Å from the membrane center. Comparison of the MD-estimated immersion depths of Tempo and n-Doxyl labels with their suggested experimental depth positions allows us to review critically the possible sources of error in depth-dependent fluorescence quenching studies.