
==== Front
RSC Adv
RSC Adv
RA
RSCACL
RSC Advances
2046-2069
The Royal Society of Chemistry

d4ra04846f
10.1039/d4ra04846f
Chemistry
Cholesterol inhibits oxygen permeation through biological membranes: mechanism against double-bond peroxidation†‡
† Videos of POPC bilayers with O2 molecules both with and without cholesterol are available at DOI: https://doi.org/10.5281/zenodo.12570598.

‡ Electronic supplementary information (ESI) available: The trajectories of 16 oxygen (O2) in lipid bilayers (Fig. S1–S6) and the density profiles of lipid and cholesterol in bilayers (Fig. S7) are provided. See DOI: https://doi.org/10.1039/d4ra04846f

Boonnoy Phansiri ab
Janlad Minchakarn ab
https://orcid.org/0000-0003-0887-5074
Bagheri Behnaz cd
Dias Cristiano e
https://orcid.org/0000-0002-8626-3033
G-4531-2010
Karttunen Mikko fg
https://orcid.org/0000-0002-3651-9870
Wong-ekkabut Jirasak ab
a Department of Physics, Faculty of Science, Kasetsart University 50 Ngamwongwan Rd, Chatuchak Bangkok 10900 Thailand jirasak.w@ku.th

b Computational Biomodelling Laboratory for Agricultural Science and Technology (CBLAST), Faculty of Science, Kasetsart University 50 Ngamwongwan Rd, Chatuchak Bangkok 10900 Thailand
c Department of Applied Physics and Science Education, Technical University of Eindhoven PO Box 513 5600 MB Eindhoven The Netherlands b.bagheri@tue.nl

d Institute for Complex Molecular Systems PO Box 513 5600 MB Eindhoven The Netherlands
e Department of Physics, New Jersey Institute of Technology Newark New Jersey 07102-1982 USA cld@njit.edu

f Department of Chemistry, Western University 1151 Richmond Street London Ontario N6A 5B7 Canada mkarttu@uwo.ca

g Department of Physics and Astronomy, Western University 1151 Richmond Street London Ontario N6A 3K7 Canada
13 9 2024
12 9 2024
13 9 2024
14 40 2911329121
4 7 2024
23 8 2024
This journal is © The Royal Society of Chemistry
2024
The Royal Society of Chemistry
https://creativecommons.org/licenses/by-nc/3.0/ This article is licensed under a Creative Commons Attribution-Non Commercial 3.0 Unported Licence. You can use material from this article in other publications without requesting further permissions from the RSC, provided that the correct acknowledgement is given and it is not used for commercial purposes.
The presence of oxygen molecules (O2) in biological membranes promotes lipid peroxidation of phospholipids with unsaturated acyl chains. On the other hand, cholesterol is considered to be an antioxidant molecule as it has a significant barrier effect on the permeation of O2 across membranes. However, a comprehensive explanation of how cholesterol affects the distribution and diffusion of O2 within lipid bilayers is yet to be established. In this study, we investigated the interaction of oxygen molecules with polyunsaturated lipid bilayers using molecular dynamics (MD) simulations. The degree of lipid unsaturation and the concentration of cholesterol were varied to study the permeation of O2. The free energy profile of O2 diffusing from the water phase to the lipid bilayer was calculated using biased umbrella MD simulations. The results show that O2 passively translocates into the membrane without changing the physical properties of the bilayer. Interestingly, in the unsaturated lipid bilayers the presence of cholesterol led to a significantly decreased permeation of O2 and an increase in the lipid chain order. Our results indicate that the hydroxyl groups of cholesterol strongly interact with the O2 molecules effectively inhibiting interactions between the oxygens and the double bonds in unsaturated lipid tails. In addition, a linear relationship between permeation and the ratio of membrane thickness and area per lipid was found. These insights can help our understanding of how the degree of unsaturation in a lipid tail and cholesterol affect lipid peroxidation at the molecular level.

Cholesterol decreases oxygen permeability across the unsaturated lipid bilayer.

National Research Council of Thailand 10.13039/501100004704 N41A640080 Kasetsart University Research and Development Institute 10.13039/501100005621 FF(KU)33.67 Natural Sciences and Engineering Research Council of Canada 10.13039/501100000038 Unassigned Canada Research Chairs 10.13039/501100001804 Unassigned Alliance de recherche numérique du Canada 10.13039/501100021202 Unassigned National Institute of General Medical Sciences 10.13039/100000057 1R15GM148982-01 pubstatusPaginated Article
==== Body
pmc1 Introduction

Biological membranes separate internal and external environments of cells and organelles regulating the uptake/excretion of nutrients/waste.1,2 They are made primarily of a lipid bilayer with proteins anchored on it. Polar lipid head groups in these bilayers are exposed to internal/external environments with their acyl tails forming a non-polar core in between these environments. This creates a barrier that impedes the flow of polar molecules but allows the permeation of small non-polar molecules that are critical for life.3

Mechanisms by which molecules transport across cellular membranes can be classified into two main categories, active and passive.3 Active transport involves regulatory machinery which requires energy to transport target molecules in the opposite direction of their concentration gradient. Passive transport is an entropy-driven process of transporting molecules across cellular membranes. Cross-membrane transport of hydrophilic species is typically an activated process which requires assistance of specialized protein-based transporters or ion channels.4–6 Most small neutral molecules and drugs are transported across biological membranes passively.2,3,7,8 This includes oxygen (O2) that diffuses across cell membranes with a permeability coefficient of about 23 cm s−1. The passive diffusion of oxygen through the cell membrane is, however, strongly affected by lipid composition, which may help to rationalize the broad spectrum of oxygen levels observed across different cells/tissues.9 Here, we focus on providing insights into the atomic mechanisms accounting for oxygen permeation and the role played by cholesterol and unsaturated fatty acids in this process.

Cellular oxygenation is crucial for brain function,9 wound healing,10 anesthesia,11 adipose tissue dysfunction,12 reperfusion injury,13 and tumor radiation therapy.14 This process is so important that oxygen concentrations below normal levels have been related to pathologies such as cancer and heart diseases. In the newly emerging cold (room temperature) plasma therapy, a partially ionized gas composed of reactive oxygen and nitrogen species (RONS) is directly delivered to living tissues, and it is being explored for treatment of skin related disorders,15,16 cancer therapy,17–19 and wound healing.20,21 Among the RONS that cold plasma produces, is oxygen in its singlet delta state, O2(1Δg),22 which has been suggested to have a major role in cell apoptosis.23–25 Production of singlet delta oxygen (among other reactive oxygen species) is also the basis of photodynamic therapy in which a light source and a photosensitizer are used in combination with oxygen to kill malignant cells.26,27

Cellular membranes are composed of a diverse set of phospholipids. They are versatile molecules with different headgroups, charge, chain length, and saturation.1 Phosphatidylcholine (PC) is the main phospholipid in all mammalian cells (40–50%) and lipoprotein particles.28–30 The double bonds in the unsaturated acyl chains are prone to oxidative stress caused by reactive oxygen species.31–33 These reactions typically modify their chemical structure which consecutively disturbs the physiological properties of cellular membrane.34–36

In this work, we employ classical molecular dynamics (MD) simulations to study the mechanisms of permeability of O2 in six different phospholipid-based cell membrane models (details in Section 2). Although, singlet delta state of O2 is not described explicitly in our simulations, our study provides insights in the general mechanisms of its diffusion across cell membrane. We investigate how the O2 molecules interact with the double bonds of the lipid chains by exploring its location and distribution in the different lipid bilayers. As cholesterol has been suggested to protect lipid bilayers from oxidative stress,37,38 we investigate the effect(s) of cholesterol on the distribution and diffusional pathways of O2 in the lipid bilayers. In addition, biased umbrella sampling MD simulations were used to calculate the free energy profiles of O2 from the water phase to the lipid bilayers.

2 Methods

2.1 Unbiased MD simulations

To study the permeability of O2 in model cell membranes, we considered phospholipid bilayers at three different saturation levels, namely 1-palmitoyl-2-oleoyl-glycero-3-phosphocholine (POPC), 1-palmitoyl-2-linoleoyl-sn-glycero-3-phosphocholine (PLPC), and 1-palmitoyl-2-arachidonoyl-sn-glycero-3-phosphocholine (PAPC). The numbers of double bonds in POPC, PLPC, and PAPC lipids are 1, 2 and 4, respectively. In addition, binary mixtures of the above lipids and cholesterol (CHOL) at 50% CHOL concentration were studied to examine the effect(s) of cholesterol on the permeability of O2. The chemical structures of lipids and cholesterol are shown in Fig. 1 and the details with the numbers of molecules for each simulated system are listed in Table 1. All MD simulations were performed with the GROMACS software package version 5.1.2 (ref. 39) using the CHARMM36 all-atom force field.40 The TIP3P water model41,42 was used in all simulations.

Fig. 1 Chemical structures of the PC lipids (POPC, PLPC, and PAPC) and cholesterol. There are 1, 2 and 4 double bonds in the sn-2 chains of POPC, PLPC, and PAPC, respectively.

List of the simulated systems showing the total numbers of phospholipids, cholesterol (CHOL), water (SOL), and oxygen molecules (O2)

System	Lipid	Cholesterol	Water	O2	
POPC	128	0	6400	16	
PLPC	128	0	6400	16	
PAPC	128	0	6400	16	
POPC:CHOL	128	128	12 800	16	
PLPC:CHOL	128	128	12 800	16	
PAPC:CHOL	128	128	12 800	16	

First, solvated lipid bilayers were constructed using the CHARMM-GUI Membrane Builder.43 The pure bilayers had 128 phospholipid molecules and 6400 water molecules, and the cholesterol-containing ones had 128 lipids, 128 CHOL molecules and 12 800 water molecules. For all the systems, an equilibration simulation was run for 100 ns to optimize the bilayer structure.

After construction of the bilayers, O2 molecules were randomly inserted into the equilibrated systems. Energy minimization using the steepest descents algorithm was performed, followed by MD simulations under a constant particle number, pressure, and temperature conditions (the NPT ensemble). The pressure was set to 1 bar with semi-isotropic pressure coupling, time constant of 5.0 ps and compressibility of 4.5 × 10−5 bar−1 by using the Parrinello–Rahman method.44 The temperature was set to 310 K with a time constant of 1.0 ps using Nosé–Hoover thermostat technique.45–47 Periodic boundary conditions were applied in all directions. A cutoff distance of 1.2 nm was used for the van der Waals and the real space part of the electrostatic interactions. The Particle-Mesh Ewald (PME) method48,49 was applied to calculate the long-range part of the electrostatic interactions. The Lennard-Jones forces were smoothly switched to zero between 1.0 and 1.2 nm. All bond lengths were constrained by the parallel linear constraint solver (P-LINCS) algorithm.50 The integration time step was set to 2 fs, and the simulations were run up to 1 μs. The Visual Molecular Dynamics (VMD) software51 was used for molecular visualizations.

2.2 Biased umbrella MD simulations

We computed the free energy profiles for O2 transport in all the six bilayers. The potential of mean force (PMF) for moving an O2 molecule from the water phase to the bilayer center was computed using the biased umbrella sampling MD simulations52 with the Weighted Histogram Analysis Method.53,54

Initially, an O2 molecule was placed in the water phase at 3.4 nm away from the lipid bilayer center. During the umbrella sampling MD simulations, the distance between the center of mass of the O2 and the lipid bilayer was constrained in the z-direction (reaction coordinate) using a harmonic potential function with a force constant of 1000 kJ (mol nm2)−1. In the simulations, 35 windows with varying positions of O2 were performed to compute the PMF profile for each bilayer. The calculation was divided into 0.1 nm per window in the z-direction to pull O2 from the water phase (z = 3.4 nm) to the bilayer center (z = 0 nm). All simulations were performed in the NPT ensemble at a constant of temperature at 310 K and a constant of pressure at 1 bar. Each window was run for at least 100 ns, and the last 50 ns were used to compute the PMF profile. The statistical uncertainty of the PMF profile was estimated using bootstrap analysis.54

3 Results

The location and distribution of O2 in the different bilayers are considered to describe the potential of the O2 molecules to interact with the double bonds of the lipid chains. Under oxidative stress in the cell membrane, radiation, or photosensitized oxidation, an O2 molecules in its ground state gains energy and may jump to an excited singlet state (its lowest excited state).55 This may lead to lipid peroxidation in the cell membrane which has harmful effects on the biophysical properties and functions of the cell membrane. Previous studies37 have demonstrated that cholesterol help to protect lipid bilayers from oxidative attack by free radicals. Here, we examine the effect of cholesterol on the distribution and diffusional pathway of O2 in lipid bilayers. In addition, free energy calculations were performed using biased umbrella MD simulation techniques to determine the free energy of transfer of O2 from the water phase to the lipid bilayer. All the simulations showed that the O2 molecules passively translocate into the bilayers without causing pore formation or membrane rupture. Interestingly, the presence of cholesterol in the unsaturated lipid bilayers strongly decreases the permeability of O2.

3.1 Location and distribution of O2 in the lipid bilayer

Fig. 2 shows snapshots of the systems after 1.0 μs. O2 is a small and non-polar molecule, and it is capable to passively and easily diffuse through the lipid bilayer. This is consistent with numerous previous experimental and computational studies.56–58 We calculated the mass density profiles of the O2 and cholesterol molecules, and all the carbon double bonds in the unsaturated lipid chains as a function of the distance from the bilayer center (the z-direction). The results are shown in Fig. 3. The average distance of the phosphorus atoms (of the lipid head groups) from the bilayer center was used as the operational definition for the lipid–water interface in order to determine the bilayer thickness. Bilayer thicknesses decreased as the number of double bonds in the lipid chain increased. This is due to the increase in fluidity and disorder in the presence of double bonds. The data is shown in Table 2.

Fig. 2 Snapshots at the end of the 1 μs simulation of each of the systems. The lipid bilayers without cholesterol (A–C) consist 16 oxygen molecules (O2). The light blue region is water. The red spheres are O2, and yellow, black, green, and purple lines represent POPC, PLPC, PAPC, and cholesterol molecules, respectively. The green spheres are the phosphorus atoms of the lipid headgroups. For the cholesterol-containing lipid bilayers (D–F), the pink spheres are the oxygen atoms of the hydroxyl groups in the cholesterol headgroups. Note that all O2 molecules in the water phase can enter and exit the lipid bilayer.

Fig. 3 Mass density profiles of the carbon atoms in the double bonds of the lipids hydrocarbon chains both with and without cholesterol. (A) POPC, (B) PLPC, and (C) PAPC bilayer. The mass density profiles of O2 are also shown. The vertical solid and dashed lines refer to the average distances of the phosphorus atoms from bilayer center for the systems with and without cholesterol, respectively.

Bilayer thickness, average number of O2 contacts with the double bond carbons in the lipid chains (C <svg xmlns="http://www.w3.org/2000/svg" version="1.0" width="13.200000pt" height="16.000000pt" viewBox="0 0 13.200000 16.000000" preserveAspectRatio="xMidYMid meet"><metadata> Created by potrace 1.16, written by Peter Selinger 2001-2019 </metadata><g transform="translate(1.000000,15.000000) scale(0.017500,-0.017500)" fill="currentColor" stroke="none"><path d="M0 440 l0 -40 320 0 320 0 0 40 0 40 -320 0 -320 0 0 -40z M0 280 l0 -40 320 0 320 0 0 40 0 40 -320 0 -320 0 0 -40z"/></g></svg> C) and OH in cholesterol (OH-CHOL), and the lifetimes of O2 contacts with (CC) (ps) were calculated using eqn (2)

Systems	Thickness (nm)	Average area per lipid (nm2)	Average number of O2 in contact with CC	Average number of O2 in contact with OH-CHOL	Lifetime of O2 in contact with CC (ps)	
POPC	3.89 ± 0.06	0.646 ± 0.014	1.87 ± 1.74	—	9.74	
POPC:CHOL	4.59 ± 0.03	0.425 ± 0.004	0.60 ± 1.01	0.28 ± 0.64	7.85	
PLPC	3.79 ± 0.06	0.663 ± 0.013	4.13 ± 2.65	—	7.97	
PLPC:CHOL	4.49 ± 0.03	0.433 ± 0.004	1.90 ± 1.86	0.30 ± 0.65	7.92	
PAPC	3.65 ± 0.06	0.711 ± 0.015	7.35 ± 3.56	—	7.03	
PAPC:CHOL	4.49 ± 0.04	0.445 ± 0.005	3.31 ± 2.52	0.27 ± 0.63	6.99	

Previous experimental and computational studies have reported that cholesterol increases the membrane thickness while reducing the area per lipid.59 Our results are well in agreement with those studies for all lipid bilayer types studied here, see Table 2. Although the number of double bonds in the PLPC lipid chain (2 double bonds) is less than in the PAPC chain (4 double bonds), we noted no significant differences in the bilayer thicknesses between these two lipid types at 50% of cholesterol (Table 2).

In addition, we calculated the time evolution of all O2 molecules moving in the system with respect to the distance from the bilayer center in the z-direction (see Fig. S1–S6†). The results demonstrate that the O2 molecules passively enter and exit the lipid bilayer, the most favourable location of O2 is being at the bilayer center as shown by the mass density in Fig. 3. In the presence of cholesterol, a small peak in density was observed at the position of the cholesterol hydrocarbon rings, and the density of O2 decreased in the hydrocarbon lipid tail region especially in the double bond region. This result suggests that cholesterol decrease the possibility of O2 to be in contact with the double-bond regions of lipids.

3.2 Number of contacts and contact lifetimes

We calculated the number of contacts between the O2 molecules and the carbon atoms in the double bonds of the lipids. A contact was defined to exist when the distance was ≤0.35 nm; overlap distance criterion was taken from a carbon–oxygen hydrogen bonding.60 Additionally, the average contact lifetime (τcontact) was calculated from the autocorrelation functions C(t) of the lifetime distribution (P(τ)) of all O2 and CC contacts as1

2

As the results in Table 2 show, the number of contacts increases as the number of double bonds increases. For all the lipid bilayer types, the presence of cholesterol resulted in a decrease in the number of contacts; the number of contacts was roughly halved with 50% cholesterol present.

Table 2 also shows that the lifetimes for the O2 molecules in contact with the double bonds is <10 ps for all bilayers. The lifetimes were calculated using eqn (2). These results indicate that the presence of cholesterol helps to protect the lipid bilayer and to decrease the possibility of lipid oxidation by reducing molecular contacts between the O2 molecules and the double-bond regions of lipids.

3.3 Permeability of O2 in the lipid bilayer

O2 is a small and relatively hydrophobic molecule, see Table 1 in Endeward et al.,61 and it can diffuse through the lipid bilayers. We calculated the frequency of O2 diffusion through a lipid bilayer. The path of oxygen translocation in a lipid bilayer was defined with three lines that represent the boundaries and the center of the lipid bilayer, see Fig. 4. For the upper and lower leaflets, the average positions of the phosphorous atoms are assigned to the first and the third line in Fig. 4, respectively. The center of a bilayer is shown by the second line. Translocation of O2 across these three lines was counted as a successful permeation event. The durations of successful permeation events was also recorded. We monitored the systems with 16 O2 molecules. The number of events of diffusion of O2 molecules across the lipid bilayers during 1 μs is shown in Fig. 5.

Fig. 4 Definition of the permeation pathway for the oxygen molecules to move through a lipid bilayer. The lines show the bilayer center and the average positions of the phosphorous atoms at the head groups. The latter defines the interface between the bilayer and the water solution.

Fig. 5 The number of O2 molecules passing through the bilayer as a function of time.

The total number of O2 molecules permeations through the lipid bilayers per ns was determined by the slope of the total number of diffused O2 molecules through the lipid bilayers and are shown in Table 3.

The number of events that 16 O2 molecules in the system diffuse across the lipid bilayer

Systems	Total number of O2 passing the bilayer within 1 μs	The number of O2 permeations per ns (# per ns)	Average time for O2 passing through the lipid bilayer (ns)	
POPC	586	0.57 ± 0.17	12.70 ± 0.50	
POPC:CHOL	217	0.21 ± 0.10	24.90 ± 1.52	
PLPC	695	0.72 ± 0.19	9.54 ± 0.32	
PLPC:CHOL	284	0.30 ± 0.10	19.27 ± 1.03	
PAPC	1072	1.07 ± 0.11	6.73 ± 0.19	
PAPC:CHOL	365	0.35 ± 0.12	16.29 ± 0.82	

The results show that the more double bonds present in the lipid tails, the easier it is for the O2 molecules to permeate across the lipid bilayer. This was observed both in the absence and presence of CHOL, see Table 3.

The number of O2 permeation events through the bilayers is shown in Table 3. The results show that the presence of cholesterol significantly decreases the number of O2 permeation across the lipid bilayer. This is consistent with the decrease in the number of contacts between O2 and the double bond(s) due to the presence of cholesterol (Table 2).

We also calculated the average times that the O2 molecules reside inside the lipid bilayer, Table 3. The results show that an increase in the number of double bonds results in faster O2 passage through the lipid bilayer. Interestingly, the presence of cholesterol leads to roughly doubling the residence time in all systems. This increase in residence times caused by cholesterol, which traps oxygen in the middle of the lipid bilayer, prevents the O2 molecules from reaching the lipid double-bond region to exit the lipid bilayer, see Fig. 3. The reasons for this are elaborated in the next section where the free energy profiles are computed.

Interestingly, the concentration of about 50% is rather high, but there are tissues that have over 70% of cholesterol concentration. Of the lipids used in this study, POPC has the highest main phase transition temperature at around −2 °C (about 271 K). Mixtures of POPC and cholesterol have been studied quite extensively. For example, Dotson et al.62 performed MD simulations of POPC with different concentrations of cholesterol up to 100% cholesterol at 310 K (about 36.85 °C), the same temperature as in our simulations. They did not report any signs of a phase transition. Here, the systems were visualized, and we calculated the lipid tail's angle distributions as shown in Fig. 6. The emergence of an ordered phase was observed at 50% of cholesterol, this is also evident in the snapshots in Fig. 2. The effect of ordered phase on oxygen and water permeation to bilayer was studied by Ghysels et al.63 They showed that the bilayer with an ordered phase is less permeable for oxygen and water compared to the one in disordered phase. This implies that the presence of an order phase may play the role on O2 permeation.

Fig. 6 Average tail angle (sn1 and sn2) distributions of POPC, PLPC, and PAPC lipids were calculated between the vector from the second to the last carbon of each tail and the bilayer normal (z-direction).

3.4 Permeation scales with membrane thickness and area per lipid

The structure of a lipid membrane is strongly affected by the degree of unsaturation of its comprising lipids and the presence of cholesterol. In particular, these quantities affect the membrane thickness and the area per lipid (Table 2) that can be related to the distance travel by small molecules during a permeation event and the space available to them within the bilayer core, respectively. Accordingly, membranes with higher area per lipid can more easily accommodate molecules as they permeate the membrane. Conversely, the probability of a molecule to abort a permeation event increases with the membrane thickness. This suggests that the number of permeation events (PE) per unit of time, increases with the area per lipid and decrease with the membrane thickness: . The inset of Fig. 7 depicts this scaling separately for membranes without (black) and with (red) cholesterol.

Fig. 7 The number of permeation events (PE) scales with the area per lipid and membrane thickness. The fraction fCHOL of cholesterol molecules per lipids in the membrane is zero and one in simulations performed without (black) and with (red) cholesterol, respectively. The area per lipid and the membrane thickness of the six membranes studied here is given in Table 2. The number of permeation events per ns is given in Table 3. Best-fit lines to data points are shown using a dashed line.

To account for a scaling law that unifies membranes with and without cholesterol, one has to recognize that cholesterol occupies some of the space that is otherwise available to accommodate O2 in a membrane with equal thickness/area per lipid but depleted of this sterol. Accordingly, a factor needs to be added to the proposed empirical scaling law to account for a reduction in permeation due to cholesterol: . The constant γ enables superposition of best-fit-lines in the inset. For our two sets of membranes, γ is 0.53858.

The existence of a scaling law (Fig. 7) shows that despite the complex effects (e.g., Fig. 3 and 6) induced by unsaturated lipids and cholesterol on the membrane structure, the number of permeation events is dominated by the area per lipid and the membrane thickness.64

3.5 Free energy profiles

Fig. 8 shows the PMF profiles for O2 translocation into the bilayers. In all cases, the free energy barrier is at the lipid–water interface. After the maximum at the interfacial barrier, the free energy decreases toward the bilayer interior. The minimum value was found to be always at the bilayer center (at zero nm in Fig. 8). This represents the equilibrium position, or the most favorable location, of O2 within a lipid bilayer. This minimum position directly matches the mass density profiles of O2 measured from the unbiased MD simulations, Fig. 3, with the bilayer center having the highest density of O2.

Fig. 8 The free energy profiles for moving an O2 molecule from the water phase to the bilayer center. The vertical solid lines refer to the average positions of the phosphorus atoms in the systems with and without cholesterol. The vertical dot-dashed lines represent the hydroxyl groups of the cholesterols. The free energy was set to zero in the water phase (z = 3.4 nm). The bootstrap method54 was used to compute the error bars, shown as transparent shading.

The values for the free energy of O2 translocation across the lipid bilayer are reported in Table 4. The minimum free energy at the bilayer center (Gcenter) was −7.50 ± 0.26, −8.61 ± 0.25, and −6.97 ± 0.26 kJ mol−1 for POPC, PLPC, and PAPC lipid bilayer, respectively. In the presence of cholesterol, a deeper Gcenter was observed both in the POPC:CHOL bilayer (Gcenter = −10.00 ± 0.30) and the PLPC:CHOL bilayer (Gcenter = −9.34 ± 0.38) which may be related to the decrease in density at the bilayer center (see Fig. S7†). Except for the PAPC:CHOL bilayer, Gcenter was found to be increased to −5.76 ± 0.32 kJ mol−1. The decrease of free energy at the bilayer is related to the decrease of density at the bilayer center. It is worth noting that the local minima in the free energy profiles, located at approximately 1.75 nm, align with the positions of the cholesterol hydroxyl groups. The presence of cholesterol within the bilayer creates a dual energy barrier that significantly retards the diffusion of oxygen molecules through the central region of the bilayers. Note that the thermal energy in our simulations is 2.58 kJ mol−1, and thus thermal fluctuations play a role in aiding O2 to overcome the energy barriers in Fig. 8.

The free energy of O2 transfer from water to the lipid bilayer (kJ mol−1)

Systems	Free energy at bilayer center	Free energy barrier at lipid–water interface	Free energy at cholesterol	Total energy	
POPC	−7.50 ± 0.26	2.74 ± 0.15	—	2.74 ± 0.15	
POPC:CHOL	−10.00 ± 0.30	2.13 ± 0.11	1.50 ± 0.20	3.63 ± 0.23	
PLPC	−8.61 ± 0.25	2.4 ± 0.12	—	2.4 ± 0.12	
PLPC:CHOL	−9.34 ± 0.38	1.87 ± 0.16	1.07 ± 0.31	2.94 ± 0.35	
PAPC	−6.97 ± 0.26	2.05 ± 0.11	—	2.05 ± 0.11	
PAPC:CHOL	−5.76 ± 0.32	1.59 ± 0.12	1.43 ± 0.32	3.02 ± 0.35	

In a prior computational study, Dotson et al.62 studied mixtures of POPC and cholesterol (0–100% of cholesterol). Their free energy results for the case that can be compared, that is, POPC membrane with 50% of cholesterol are in excellent agreement. In another study, van der Paal et al.65 investigated systems with cholesterol and 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC). Similarly to our findings, they found the appearance of new free energy barriers in the presence of cholesterol.

In another study, Zuniga-Hertz and Patel66 concluded that as the number of unsaturations in the lipids increases, there is an increasing number of free volume pockets in the system. Such pockets were proposed to be able to accommodate O2 molecules. The current results indicate that the center of the bilayer is able to accommodate O2 molecules as it is the location with the most favourable free energy. However, additional free energy minima (Fig. 8) within the bilayer may also play an important role in trapping O2 molecules. The latter are more pronounced in membranes made from PLPC and PAPC lipids and in the presence of cholesterol.

4 Conclusions

In this work, we investigated the permeation of O2 through the polyunsaturated lipid bilayers, with and without cholesterol, using molecular dynamics (MD) simulations. The results show that O2 passively diffuses from the aqueous phase through the bilayer. O2 tends to favor the central region of bilayer, in agreement with the calculated free energy transfer from the aqueous phase to the lipid bilayer. The presence of cholesterol in the unsaturated lipid bilayers notably reduces the permeation of O2 through the membrane. This effect is accompanied by an increase in bilayer thickness and the emergence of a local energy minimum at the location of cholesterol hydroxyl group. Our results suggest that the presence of cholesterol serves as a protective shield for lipid bilayers by reducing the frequency (number of contacts) and the duration (lifetime) of interactions between the O2 molecules and the double-bond regions of the lipids. This new insight supports the antioxidant function of cholesterol in cellular membranes.

Although beyond the current study, we would like to point out that there the matter of permeation is even more complex and there are several factors that influence it, in particular the amount of cholesterol, the phase state of the membrane, and the area per lipid and thickness, both of which depend on the phase of the membrane.

Finally, we would also like to mention that no spontaneous flip-flops were observed in simulations. We have, however, studied flip-flops in oxidized PLPC systems in a prior study.38 Spontaneous flip-flops have been reported in ketosterone-containing systems,67 albeit in the absence of oxidized lipids. Molecular mechanisms of flip-flop have been discussed in detail, for example, by Gurtovenko and Vattulainen.68,69

Data availability

All data generated or analyzed during this study are included in this published article and ESI.†

Author contributions

All authors conceived the work, wrote the manuscript and provided critical feedback on the interpretation of data, analysis and discussion. PB carried out the simulations and data analysis.

Conflicts of interest

There are no conflicts to declare.

Supplementary Material

RA-014-D4RA04846F-s001

This work was financially supported by the National Research Council of Thailand (NRCT) through the Research Grants for Talented Mid-Career Researchers with grant no. N41A640080, the National Science Research and Innovation fund (NSRF) via the Program Management Unit for Human Resources and Institutional Development, Research and Innovation [grant number B42G670041] and Kasetsart University Research and Development Institute (KURDI) through Fundamental Fund (grant number FF(KU)33.67. MJ thanks the Development and Promotion of Science and Technology Talent Project (DPST). BB thanks the strategic alliance between TU/e, Utrecht University, and University Medical Center Utrecht for financial support. MK thanks the Natural Sciences and Engineering Research Council of Canada (NSERC) and the Canada Research Chairs Program for financial support, and the Digital Research Alliance of Canada for computational resources. CD is supported by the National Institute of General Medical Health under grant no. 1R15GM148982-01.
==== Refs
References

Van Meer G. Voelker D. R. Feigenson G. W. Nat. Rev. Molec. Cell. Biol. 2008 9 112 124 18216768
Subczynski W. K. Widomska J. Mainali L. Adv. Exp. Med. Biol. 2017 977 27 34 28685424
Shinoda W. Biochim. Biophys. Acta 2016 1858 2254 2265 27085977
Deamer D. W. Bramhall J. Chem. Phys. Lipids 1986 40 167 188 2427233
Gurtovenko A. A. Vattulainen I. Biophys. J. 2007 92 1878 1890 17208976
Wong-ekkabut J. Karttunen M. J. Biol. Phys. 2016 42 133 146 26264478
Nalakarn P. Boonnoy P. Nisoh N. Karttunen M. Wong-ekkabut J. Sci. Rep. 2019 9 1037 30705323
Nisoh N. Jarerattanachat V. Karttunen M. Wong-ekkabut J. Biomolecules 2022 12 639 35625567
Carreau A. Hafny-Rahbi B. E. Matejuk A. Grillon C. Kieda C. J. Cell. Mol. Med. 2011 15 1239 1253 21251211
Castilla D. M. Liu Z.-J. Velazquez O. C. Adv. Wound Care 2012 1 225 230
Dunn J.-O. Mythen M. Grocott M. BJA Educ. 2016 16 341 348
Netzer N. Gatterer H. Faulhaber M. Burtscher M. Pramsohler S. Pesta D. Biomolecules 2015 5 1143 1150 26061760
Granger D. N. Kvietys P. R. Redox Biol. 2015 6 524 551 26484802
Multhoff G. Radons J. Vaupel P. Cancers 2014 6 813 828 24717239
Wirtz M. Stoffels I. Dissemond J. Schadendorf D. Roesch A. J. Eur. Acad. Dermatol. Venereol. 2018 32 e37 e39 28695987
Heinlin J. Isbary G. Stolz W. Zeman F. Landthaler M. Morfill G. Shimizu T. Zimmermann J. Karrer S. J. Eur. Acad. Dermatol. Venereol. 2013 27 324 331 22188329
Ratovitski E. A. Cheng X. Yan D. Sherman J. H. Canady J. Trink B. Keidar M. Plasma Processes Polym. 2014 11 1128 1137
Rutkowski R. Schuster M. Unger J. Seebauer C. Metelmann H. Woedtke T. v. Weltmann K. Daeschlein G. Clin. Plasma Med. 2017 7 52 57
Metelmann H.-R. Seebauer C. Miller V. Fridman A. Bauer G. Graves D. B. Pouvesle J.-M. Rutkowski R. Schuster M. Bekeschus S. et al. Clin. Plasma Med. 2018 9 6 13
Ulrich C. Kluschke F. Patzelt A. Vandersee S. Czaika V. Richter H. Bob A. Hutten J. v. Painsi C. Hüge R. et al. J. Wound Care 2015 24 196 203 25970756
Bekeschus S. Schmidt A. Napp M. Kramer A. Kerner W. von Woedtke T. Wende K. Hasse S. Masur K. Exp. Dermatol. 2017 26 145 147 27673417
Sousa J. S. Niemi K. Cox L. Algwari Q. T. Gans T. O’connell D. J. Appl. Phys. 2011 109 123302
Riethmüller M. Burger N. Bauer G. Redox Biol. 2015 6 157 168 26225731
Bauer G. Sersenová D. Graves D. B. Machala Z. Sci. Rep. 2019 9 13931 31558835
Nasri Z. Memari S. Wenske S. Clemen R. Martens U. Delcea M. Bekeschus S. Weltmann K.-D. von Woedtke T. Wende K. Chem.–Eur. J. 2021 27 14702 14710 34375468
Li B. Lin L. Lin H. Wilson B. C. J. Biophotonics 2016 9 1314 1325 27136270
Agostinis P. Berg K. Cengel K. A. Foster T. H. Girotti A. W. Gollnick S. O. Hahn S. M. Hamblin M. R. Juzeniene A. Kessel D. et al. Ca-Cancer J. Clin. 2011 61 250 281 21617154
Fruhwirth G. O. Loidl A. Hermetter A. Biochim. Biophys. Acta, Mol. Basis Dis. 2007 1772 718 736
Makky A. Tanaka M. J. Phys. Chem. B 2015 119 5857 5863 25870900
van der Veen J. N. Kennelly J. P. Wan S. Vance J. E. Vance D. E. Jacobs R. L. Biochim. Biophys. Acta, Biomembr. 2017 1859 1558 1572 28411170
Frankel E. Prog. Lipid Res. 1984 23 197 221 6100997
Brodnitz M. H. Nawar W. W. Fagerson I. S. Lipids 1968 3 59 64 17805843
Brodnitz M. H. Nawar W. W. Fagerson I. S. Lipids 1968 3 65 71 17805844
Spiteller P. Kern W. Reiner J. Spiteller G. Biochim. Biophys. Acta, Mol. Cell Biol. Lipids 2001 1531 188 208
Pratt D. A. Mills J. H. Porter N. A. J. Am. Chem. Soc. 2003 125 5801 5810 12733921
Boonnoy P. Jarerattanachat V. Karttunen M. Wong-Ekkabut J. J. Phys. Chem. Lett. 2015 6 4884 4888 26673194
Parasassi T. Giusti A. M. Raimondi M. Ravagnan G. Sapora O. Gratton E. Free Radical Biol. Med. 1995 19 511 516 7590402
Boonnoy P. Jarerattanachat V. Karttunen M. Wong-Ekkabut J. Biophys. J. 2021 120 4525 4535 34478697
Abraham M. J. Murtola T. Schulz R. Páll S. Smith J. C. Hess B. Lindahl E. SoftwareX 2015 1 19 25
Lee S. Tran A. Allsopp M. Lim J. B. Hénin J. Klauda J. B. J. Phys. Chem. B 2014 118 547 556 24341749
Jorgensen W. L. Chandrasekhar J. Madura J. D. Impey R. W. Klein M. L. J. Chem. Phys. 1983 79 926 935
MacKerell A. D. J. Bashford D. Bellott M. Dunbrack R. L. J. Evanseck J. D. Field M. J. Fischer S. Gao J. Guo H. Ha S. Joseph-McCarthy D. Kuchnir L. Kuczera K. Lau F. T. K. Mattos C. Michnick S. Ngo T. Nguyen D. T. Prodhom B. Reiher W. E. Roux B. Schlenkrich M. Smith J. C. Stote R. Straub J. Watanabe M. Wiórkiewicz-Kuczera J. Yin D. Karplus M. J. Phys. Chem. B 1998 102 3586 3616 24889800
Jo S. Cheng X. Lee J. Kim S. Park S.-J. Patel D. S. Beaven A. H. Lee K. I. Rui H. Park S. Lee H. S. Roux B. MacKerell A. D. Klauda J. B. Qi Y. Im W. J. Comput. Chem. 2016 38 1114 1124 27862047
Parrinello M. Rahman A. J. Appl. Phys. 1981 52 7182 7190
Nosé S. J. Chem. Phys. 1984 81 511 519
Nosé S. Mol. Phys. 1984 52 255 268
Hoover W. G. Phys. Rev. A: At., Mol., Opt. Phys. 1985 31 1695 1697
Darden T. York D. Pedersen L. J. Chem. Phys. 1993 98 10089 10092
Essmann U. Perera L. Berkowitz M. L. Darden T. Lee H. Pedersen L. G. J. Chem. Phys. 1995 103 8577 8593
Hess B. J. Chem. Theory Comput. 2008 4 116 122 26619985
Humphrey W. Dalke A. Schulten K. J. Mol. Graphics 1996 14 33 38
Torrie G. M. Valleau J. P. J. Comput. Phys. 1977 23 187 199
Kumar S. Rosenberg J. M. Bouzida D. Swendsen R. H. Kollman P. A. J. Comput. Chem. 1992 13 1011 1021
Hub J. S. de Groot B. L. van der Spoel D. J. Chem. Theory Comput. 2010 6 3713 3720
Maharjan P. S. Bhattarai H. K. J. Oncol. 2022 2022 7211485 35794980
Subczynski W. K. Hyde J. S. Kusumi A. Proc. Natl. Acad. Sci. U.S.A. 1989 86 4474 4478 2543978
Al-Samir S. Itel F. Hegermann J. Gros G. Tsiavaliaris G. Endeward V. Cell. Mol. Life Sci. 2021 78 7649 7662 34694438
Subczynski W. K. Widomska J. Raguz M. Pasenkiewicz-Gierula M. Oxygen 2022 2 295 316 36852103
Róg T. Pasenkiewicz-Gierula M. Vattulainen I. Karttunen M. Biochim. Biophys. Acta, Biomembr. 2009 1788 97 121
Horowitz S. Trievel R. C. J. Biol. Chem. 2012 287 41576 41582 23048026
Endeward V. Al-Samir S. Itel F. Gros G. Front. Physiol. 2014 4 382 24409149
Dotson R. J. Smith C. R. Bueche K. Angles G. Pias S. C. Biophys. J. 2017 112 2336 2347 28591606
Ghysels A. Krämer A. Venable R. M. Teague W. E. Lyman E. Gawrisch K. Pastor R. W. Nat. Commun. 2019 10 5616 31819053
Andziak B. Buffenstein R. Aging Cell 2006 5 525 532 17129214
Van der Paal J. Verheyen C. Neyts E. C. Bogaerts A. Sci. Rep. 2017 7 39526 28059085
Zuniga-Hertz J. P. Patel H. H. Front. Physiol. 2019 10 1340 31736773
Róg T. Stimson L. M. Pasenkiewicz-Gierula M. Vattulainen I. Karttunen M. J. Phys. Chem. B 2008 112 1946 1952 18217743
Gurtovenko A. A. Vattulainen I. J. Phys. Chem. B 2007 111 13554 13559 17988118
Gurtovenko A. A. and Vattulainen I. , Biomembrane Frontiers: Nanostructures, Models, and the Design of Life, Humana Press, Totowa, NJ, 2009, pp. 121–139
