Next Article in Journal
Assessment of Oxidative Stress-Induced Oral Epithelial Toxicity
Next Article in Special Issue
A New Theory about Interfacial Proton Diffusion Revisited: The Commonly Accepted Laws of Electrostatics and Diffusion Prevail
Previous Article in Journal
Effects of Tropomodulin 2 on Dendritic Spine Reorganization and Dynamics
Previous Article in Special Issue
Fatty Acid-Activated Proton Transport by Bisaryl Anion Transporters Depolarises Mitochondria and Reduces the Viability of MDA-MB-231 Breast Cancer Cells
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Biomolecules
Volume 13
Issue 8
10.3390/biom13081238
biomolecules-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Graph-Based Analyses of Dynamic Water-Mediated Hydrogen-Bond Networks in Phosphatidylserine: Cholesterol Membranes

by
Honey Jain
1,2,
Konstantina Karathanou
2 and
Ana-Nicoleta Bondar
1,3,*
1
Faculty of Physics, University of Bucharest, Atomiștilor 405, 077125 Măgurele, Romania
2
Department of Physics, Freie Universität Berlin, Arnimallee 14, 14195 Berlin, Germany
3
IAS-5/INM-9, Forschungszentrum Jülich, Institute of Computational Biomedicine, Wilhelm-Johnen Straße, 52428 Jülich, Germany
*
Author to whom correspondence should be addressed.
Biomolecules 2023, 13(8), 1238; https://doi.org/10.3390/biom13081238
Submission received: 20 June 2023 / Revised: 3 August 2023 / Accepted: 7 August 2023 / Published: 11 August 2023
(This article belongs to the Special Issue Proton and Proton-Coupled Transport)
Browse Figures
Review Reports Versions Notes

Abstract

:
Phosphatidylserine lipids are anionic molecules present in eukaryotic plasma membranes, where they have essential physiological roles. The altered distribution of phosphatidylserine in cells such as apoptotic cancer cells, which, unlike healthy cells, expose phosphatidylserine, is of direct interest for the development of biomarkers. We present here applications of a recently implemented Depth-First-Search graph algorithm to dissect the dynamics of transient water-mediated lipid clusters at the interface of a model bilayer composed of 1-palmytoyl-2-oleoyl-sn-glycero-2-phosphatidylserine (POPS) and cholesterol. Relative to a reference POPS bilayer without cholesterol, in the POPS:cholesterol bilayer there is a somewhat less frequent sampling of relatively complex and extended water-mediated hydrogen-bond networks of POPS headgroups. The analysis protocol used here is more generally applicable to other lipid:cholesterol bilayers.

1. Introduction

Hydrogen (H)-bonds and H-bond networks shape the structure and dynamics of lipid–lipid interactions [1] and the fluidity of the lipid bilayer [2,3,4]. H-bonding properties depend on the nature of the lipid headgroup. The anionic phosphatidylserine (PS) has more inter-lipid H-bonds than, e.g., the zwitterionic phosphatidylethanolamine (PE) and phosphatidylcholine (PC), or than the anionic phosphatidylglycerol (PG) and phosphatidic acid (PA) [5,6,7]. This propensity for intermolecular H-bonding of PS may be associated with somewhat larger bilayer thicknesses [7].
Lipid clusters can be defined as groups of lipids H-bonded via direct H-bonds, water molecules, or ion interactions. The clustering of PS appears to be associated with the formation of nanometer-sized domains [8,9,10] that could provide, e.g., a platform for peripheral protein recruitment [11], viral assembly [12], and the binding to the membrane of the influenza A virus matrix protein [9]. Altered distribution of PS has been associated with a number of human diseases [13], including cancer [14,15,16]. The formation of lipid clusters, particularly in the case of PS [17,18,19], PG [17], and phosphatidylinositol 4,5-bisphosphate (PIP2) [20], might be promoted by the presence of cations. Experiments using, e.g., time-resolved fluorescence spectroscopy [18], fluorescence correlation spectroscopy [20,21], steady-state probe-partitioning fluorescence resonance energy transfer (SP-FRET) [20], and molecular dynamics (MD) simulations [18,19,22] suggested that the presence of sodium [9,18,21,23,24,25,26], potassium [18,22], calcium [19,22,25,26,27,28], or magnesium ions [28] may be associated with the clustering of lipids. Details of ion-mediated lipid cluster formation depend both on the cation and on the lipid headgroup. For example, compared to potassium, the presence of sodium ions tends to be associated with larger clusters in PIP2 bilayers [22] and shows a higher affinity to a membrane composed of 1-palmytoyl-2-oleoyl-sn-glycero-2-phosphatidylcholine (POPC) and 1-palmytoyl-2-oleoyl-sn-glycero-2-phosphatidylserine (POPS) [18]. At the interface of a POPS/POPC bilayer, sodium ions have reduced mobility relative to POPC [18]. Likewise, sodium and potassium ions have weaker binding at a PC than a PE lipid membrane interface [29].
Cholesterol favors interactions with anionic lipids (PS and PG) as compared to the zwitterionic PC and PE [30,31,32,33], and it can H-bond to lipid headgroups [34]. Cholesterol prefers H-bond to single phospholipids [35], typically with the lipid carbonyl or phosphate oxygen atoms [35,36,37]. Intercalation of cholesterol molecules could influence the packing of anionic lipids [30], looser at the lipid headgroups and tighter at the tails [38] and be associated with reduced sodium ion binding at the interface [39].
MD studies have shown that lipids can form H-bond networks mediated by direct or water-mediated H-bonds between lipid headgroups [5,17]. The water dynamics at the lipid interface is slower than the bulk [1,40,41]. The slower relaxation rates of interfacial waters H-bonded to lipids suggest that H-bond networks in the hydration layer facilitate lipid–lipid interactions [1]. The number of H-bonds between water molecules at the interface of DMPC [1] and DOPS [42] bilayers decreases with an increasing membrane depth.
To identify and characterize dynamic water-mediated clusters of lipid headgroups, we have recently implemented a Depth-First-Search (DFS)-based graph algorithm that uses Connected Component searches to explore the nodes of the H-bond graphs and identify four main types of lipid H-bond clusters [17]. Clusters are then characterized according to properties such as cluster size, which is given by the number of nodes (lipid headgroups) that constitute the cluster, and cluster length, which is given by the number of H-bonds within the longest H-bond path within the cluster [17]. Using the DFS-based algorithm we implemented, we dissected the water-mediated H-bond clusters at the interface of hydrated lipid bilayers composed of POPE, POPG, 3:1 POPE:POPG, 5:1 POPE:POPG, POPS, and an Escherichia coli lipid membrane model, all without cholesterol [17]. We found that, regardless of the lipid membrane composition, short linear H-bond arrangements, typically two lipid headgroups that H-bond directly or via one water molecule, were sampled frequently during atomistic MD simulations [17]; relative to the other lipids studied, POPS could engage more frequently in circular arrangements of three lipid headgroups [17].
Here, we apply the DFS graph algorithm to characterize dynamic water-mediated interactions in a POPS:cholesterol membrane with 10% cholesterol, which is in the range of cholesterol concentrations relevant to the plasma membrane of a healthy eukaryotic cell [43]. As a reference, we use a hydrated POPS bilayer without cholesterol. We compute separately H-bond clusters contributed by phosphate, serine vs. ester groups; likewise, we compute separately H-bond clusters that involve direct POPS:POPS H-bonds, vs. one-water bridges or ion-mediated bridges between POPS headgroups. We find that cholesterol disfavors the transient formation of H-bonded POPS clusters.

2. Methods

2.1. MD Simulations of Hydrated Lipid Membranes

CHARMM-GUI [44,45,46] was used to generate coordinates for two distinct hydrated POPS bilayers (Table 1): (i) without cholesterol; (ii) with 10% cholesterol. All simulation systems have 0.15 M neutralizing KCl salt and a starting size of approximately 81 Å × 81 Å × 85 Å.
We used the CHARMM36 force field for lipids and ions [47,48,49,50,51,52] and the TIP3P water model [53]. All simulations were performed with NAMD 2.13 [54,55]. The standard CHARMM-GUI scheme for the initial equilibration was used. For the production runs, we used a Langevin dynamics scheme [56,57] with an oscillation period of 200 fs, a damping time scale of 100 fs, and a damping coefficient of 5 ps−1 in the NPT ensemble (constant number of particles N, constant pressure P = 1 bar, and constant temperature T = 310.15 K). We used the smooth particle mesh Ewald summation to compute Coulomb interactions [58,59] and a 10–12 Å switch function for real-space interactions; covalent bonds to H-atoms were fixed. An integration step of 1 fs was used throughout equilibration and the first 1 ns of the production runs; for the production runs, we used a multiple time-step integration scheme [60,61] with 1 fs for bonded interactions and 2 fs for short-range nonbonded forces and long-range electrostatics. Coordinates were saved every 10 ps for data analyses. The total sampling time of the two simulations was 1 μs.

2.2. H-bond Criteria, H-bond Graphs, H-bond Paths, and H-bond Occupancies

We identified H-bonds using a distance criterion whereby two groups were considered as H-bonded if the distance between the H-atom and acceptor hetero-atom was within 2.5 Å. This single distance criterion gives results largely equivalent to a combined distance and angle criterion of H-bond distances ≤3.5 Å between the H-bond donor and acceptor hetero-atoms, and H-bond angles within 60° [62]. For potassium ion-mediated bridges between two POPS phosphate groups, we used a distance criterion of 4 Å.
A H-bond graph consists of nodes, which here are the POPS or cholesterol headgroups, and the edges, which are direct, 1 water-mediated, or ion-mediated H-bond connections between the nodes. A H-bond path between two lipid headgroups connects these two headgroups via the shortest distance H-bond path, i.e., via the smallest number of intermediate H-bonds. The occupancy of a H-bond is the percentage of coordinate sets utilized in analyses in which the H-bond criterion is satisfied [24,63]. We used Bridge [63,64] to compute H-bond occupancies.

2.3. H-bond Clusters, Topology, Cluster Size, Path Length, and Occupancy of Lipid Clusters

A H-bond cluster consists of nodes and edges that interconnect with each other. The topology of a H-bond lipid cluster is defined as the geometric arrangement of the nodes and edges in that cluster. To identify lipid clusters, we used the recently developed DFS algorithm [17] that analyses the H-bond graph to identify transient lipid clusters and the topologies of these clusters. Briefly, the DFS algorithm treats the nodes and edges of the H-bond graph as, respectively, data points and relationships between these data points. The algorithm explores all graph nodes, monitors H-bond distances in the coordinate dataset (MD trajectory), and constructs adjacency or connection matrices that represent the H-bond connections between pairs of nodes. The nodes of the graph are clustered, according to the H-bonding relationships between them, into four main topologies—linear, star and linear, circular, and complex combinations of these three arrangements. A lipid cluster is then characterized by its size, which is given by the number of lipid headgroups that constitute the cluster, and by the path length L, which is given by the number of edges (H-bonds) in the longest path connecting two nodes of the H-bond cluster.
The occupancy of a lipid cluster is the percentage of coordinate sets utilized in analyses in which the cluster is present. The length of the water wire that connects two lipid headgroups is given by the number of bridging H-bonded water molecules [63]; for simplicity, and because water bridges between lipids tend to have very small occupancies, here we considered only one-water bridges. All topology calculations were performed using the VMD [65] and MATLAB data analysis scripts of the DFS algorithm [17] as deposited in the Mendeley repository [66]. For data analyses, the complete trajectories of the production runs were used.

2.4. Structure Factors and the Number of H-bonds per POPS

We used the MEMBPLUGIN [67] in VMD [65] to evaluate the average thickness of each bilayer, area per lipid, and order parameters. The H-bonds per lipid are reported as the average number of H-bonds sampled by the specific lipid group in each coordinate set.
All molecular graphics were prepared using VMD 1.9.4 [65]. Boxplots of H-bond occupancy values extracted with Bridge [63,64] were prepared using Python. Files pertaining to the simulations reported here are released as a public Mendeley repository (see the Data Availability Statement).

3. Results and Discussion

The simulations reported here have average membrane thickness, area per lipid, and order parameter values close to previous experiments and simulations. For the POPS membrane without cholesterol, the average membrane thickness we calculate is 41.2 ± 0.5 Å, which is close to values reported from previous MD simulations, e.g., 43.2 ± 0.5 Å [68], 40.6 Å [69], 42.4 ± 0.2 Å [70], 42.1 ± 0.5 Å [17], and to the 42.2 Å thickness measured by neutron and X-ray scattering in the presence of NaCl [69]. In the presence of 10% cholesterol, the POPS membrane is ~1.5 Å thicker (Table 1, Figure 1A,C), which is compatible with the ~2.2 Å increase, relative to pure DMPC, in the thickness of DMPC/cholesterol membranes with 10% cholesterol [71]. The calculated area per lipid for the POPS membrane is 59.8 ± 1.9 Å2, which is relatively close to previously reported values from MD simulations (58.4 Å2 [72], 57.5 ± 1.2 Å2 [70], and 62.0 Å2 [69]). A decrease of ~5 Å2 in the area per lipid and an increase in the order parameters found here in the presence of cholesterol (Figure 1B,D) are in reasonable qualitative agreement with the decrease by ~6–7 Å2 [73,74] in the area per lipid of a DOPC/cholesterol membrane with 10% cholesterol. Likewise, the order parameters for the POPS membrane (Figure 1D) are compatible with previous studies [70,72].

3.1. Direct and One-Water-Mediated H-bonding of POPS Headgroups

On average, each POPS headgroup has 1.5 ± 0.1 H-bonds with another POPS (Figure 2A), which is close to the 1.1–1.2 H-bonds/POPS headgroup reported previously from CHARMM36 MD simulations [7]. Each POPS headgroup also has, on average, 2.3 ± 0.2 H-bonds from one-water bridges (Figure 2B) and 1.5 ± 0.1 potassium ion interactions (Figure 2C). Each phosphate and serine group has about 1.3–1.4 H-bonds to a bridging water molecule (Figure 2), slightly higher than the ~1.1 one-water-mediated H-bonds per ester group (Figure 2B). When the ester groups are included in the H-bond computation, the average number of one-water-mediated H-bonds increases from 1.9 ± 0.1 (when only the POPS phosphate and serine groups are included in the H-bond computation) to 2.3 ± 0.2 (Figure 2), but these ester H-bonds have very low occupancies (Figure 3). The total average number of H-bonds per POPS is similar in membranes with and without cholesterol (Figure 2).
Overall, both the direct and one-water-mediated bridges between POPS have small occupancies: 83–93% of the direct H-bonds have occupancies within 10–15% (Figure 3A,C), and 69–100% of the one-water bridges computed for the different POPS H-bonding groups have occupancies within about 5% (Figure 3B,D), suggesting these interactions are quite dynamic. In the distribution of the occupancies of direct H-bonds and one-water-mediated bridges, there are also outliers: for each of the H-bonds and water-mediated bridges presented in Figure 3, about 4–8% of the data points are outliers, i.e., these are H-bonds and water-mediated bridges with occupancy values significantly larger than the average (Figure 3). The highest-occupancy H-bonds are, in the POPS membrane, a pair of POPS headgroups that remain within H-bond distance for 306.2ns (i.e., 61.2% POPS simulation) and, in the POPS:cholesterol simulation, a pair of POPS headgroups within H-bond distance for 360.2 ns (72% of the simulation).

3.2. Water-Mediated H-bond Clusters in the POPS:Cholesterol Membrane

We used the DFS algorithm as illustrated in Scheme 1 to characterize POPS clusters in membranes with cholesterol, using as a reference the results obtained for the POPS membrane without cholesterol [17,66].
We first considered all unique H-bonds sampled during the simulations. At any given time, about 15–18 clusters are likely to be sampled in the POPS:cholesterol membrane, that is, because each leaflet has 109 lipids (Table 1), at any moment of time about 81% of the lipids participate in dynamic, low-occupancy H-bond paths of at least two lipids. About half of these H-bond paths (8–9) are linear (Scheme 1, Figure 4A). Linear, relatively short H-bond paths are favoured by POPS irrespective of whether cholesterol is present or not (Figure 4A and Figure 5) and are exclusively sampled in cholesterol clusters (Figure 6). Circular and star and linear clusters are rare and, when sampled, have minimal sizes (3 and 4–5 lipids, respectively, Figure 4A,B). Star and linear and circular arrangements may also be sampled transiently, such that, overall, most of the H-bond paths sampled are either linear or a more complex topology that includes a linear path segment (Figure 4).
As reported previously for a hydrated POPS bilayer [17], linear H-bond paths (Scheme 1) mediated by one water molecule are sampled throughout the entire POPS simulation: on average, at least one linear path is sampled at any moment in time (Figure 5E). The one-water-mediated linear H-bond paths tend to be relatively short, composed of 2–3 lipids, that is, these are singular H-bonds, or two H-bonds with one common POPS headgroup (Figure 5E). These results for the reference POPS simulation are in qualitatively good agreement with our previous computations on hydrated POPS [17] and 4:1 POPC:POPG membranes [23]. In the presence of cholesterol, direct singular H-bonds (H-bond paths with length L = 1) between POPS headgroups are slightly less likely than in the POPS membrane (Figure 5A), though such paths are sampled throughout the entire simulation (Figure 5F). Taken together with the H-bond occupancy analysis summarized in Figure 3, the time series of the H-bond paths of different lengths (Figure 5) suggest that, overall, most H-bond clusters at the POPS headgroup interface are short and transient. It should also be noted that the exact number of H-bond paths, or clusters, sampled throughout the simulation trajectories (Figure 5), depends on the size of the lipid bilayer, and the time series presented in Figure 5 are meant as an illustration of the H-bond dynamics at the lipid headgroup interface of the hydrated membrane patches used here (Table 1).
As noted before for a hydrated POPS bilayer [17], linear paths with more than 2–3 lipids are sampled but only very rarely (Figure 4 and Figure 5A,B). Most (~60%) of the linear H-bond clusters sampled by POPS headgroups have a path length L = 1, and a significant percentage (25%) of the linear clusters have a path length L = 2. That is, when sampled, linear H-bond paths are typically short, only 2–3 lipids; longer linear paths with more lipids may be sampled but very rarely (Figure 5A,B). Cholesterol molecules may, infrequently, sample small clusters (L = 1, i.e., two cholesterol molecules); longer linear H-bond clusters with 3–4 cholesterol molecules H-bonded via one-water bridges are very rarely sampled (Figure 6). The infrequent sampling of cholesterol clusters can be attributed to its preference for H-bonding to single phospholipids [35]. Even in 1:1 concentrations, cholesterol–lipid is more likely than cholesterol–cholesterol interactions [35,36,37]. We also note that the sampling frequency of specific clusters can depend on the length of the simulation.

4. Conclusions

Anionic PS lipid headgroups have a high propensity to H-bond [5,6,7], and topology analyses have revealed preferred arrangements of POPS bridged by water molecules [17]. Because the H-bond dynamics of PS can be affected by cholesterol [30,31,32,34], here we studied the dynamics of H-bond clusters sampled in atomistic simulations of a POPS membrane with cholesterol and compared the results with those for a cholesterol-free POPS bilayer.
The concentration of cholesterol used here, 10%, approximates that of a healthy eukaryotic plasma membrane [43]. We stress, however, that the POPS:cholesterol membrane studied here is but a model membrane, because the eukaryotic plasma membrane is a complex mixture of no fewer than 11 different types of lipids [43]. As anticipated based on previous results reported in the literature [75], the POPS membrane is about 1.5 Å thicker when cholesterol is present (Figure 1A,C).
Regardless of the presence of cholesterol, POPS samples dynamic (low-occupancy) water-mediated H-bond networks (Figure 3). As reported before, linear H-bond clusters are sampled more frequently than more complex topologies (Scheme 1, Figure 4), yet more complex arrangements may be sampled, particularly in the absence of cholesterol (Figure 4).
Regardless of the presence of cholesterol, POPS tends to form small linear lipid clusters with 2–3 lipids directly H-bonded or one-water-bridged (Figure 4 and Figure 5). Cholesterol is associated with fewer POPS being part of complex H-bond clusters (Figure 4).
POPS-cholesterol clusters are sampled less frequently than POPS-POPS clusters, and cholesterol–cholesterol clusters are rare (Figure 6). When sampled, cholesterol-only clusters are small, typically only two cholesterol molecules (Figure 6); this finding is compatible with the previous study [76] reporting a preference of cholesterol molecules to form dimers. POPS clusters mediated by potassium ions prefer a linear arrangement of 3–4 lipids, whereas the more complex topologies are much less frequent (Figure 4). This is compatible with previous work [17] indicating that sodium-mediated POPS lipid clusters tend to prefer linear topologies with two lipid clusters; we note, however, that the number of cation-mediated clusters sampled during the simulations, and the size of the clusters, might depend on the length of the simulations and on the size of the lipid membrane patch.
We anticipate that the protocol used here to evaluate the H-bonding of a simple cholesterol-containing model bilayer could become useful to dissect dynamic H-bond clusters sampled in simulations of more complex lipid mixtures, including for water-mediated lipid headgroup H-bond paths of potential interest for lateral proton transfer along lipid membrane interfaces, and to evaluate how different cations influence the sampling of transient H-bond clusters, and the preferred topologies of these clusters, in hydrated membranes with distinct lipid composition.

Author Contributions

Conceptualization, A.-N.B.; methodology, H.J., K.K. and A.-N.B.; software, K.K.; validation, H.J., K.K. and A.-N.B.; formal analysis, H.J.; investigation, H.J.; resources, A.-N.B.; data curation, H.J. and K.K.; writing—original draft preparation, H.J.; writing—review and editing, H.J. and A.-N.B.; visualization, H.J. and A.-N.B.; supervision, A.-N.B.; project administration, A.-N.B.; funding acquisition, A.-N.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the European Union’s Horizon 2020 Research and Innovation Program under the Marie Sklodowska-Curie grant agreement No H2020 MSCA-ITN-2019-860592, Innovative Training Network ‘Proton and proton-coupled transport’ (to A.-N.B.).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Input coordinate and protein structure files, and trajectory files with coordinate sets from the simulations reported here, are openly available at Jain, Honey (2023), “Graph-based analyses of dynamic water-mediated hydrogen-bond networks in phosphatidylserine:cholesterol membranes”, Mendeley Data, V1, doi: 10.17632/h8y6zjckcw.1 as accessed on 11 August 2023.

Acknowledgments

Research was supported by funds from the European Union’s Horizon 2020 Research and Innovation Program under the Marie Sklodowska-Curie grant agreement No. H2020 MSCA-ITN-2019-860592, Innovative Training Network ‘Proton and proton-coupled transport’ (to A-NB). MD simulations were performed using computing resources of the Physics Department at the Freie Universität Berlin, and data analyses were performed using computer resources of the Faculty of Physics of the University of Bucharest. The authors are grateful to Jens Dreger (Freie Universität Berlin) for technical support.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Srivastava, A.; Debnath, A. Hydration Dynamics of a Lipid Membrane: Hydrogen Bond Networks and Lipid-Lipid Associations. J. Chem. Phys. 2018, 148, 094901. [Google Scholar] [CrossRef]
  2. Boggs, J.M. Lipid Intermolecular Hydrogen Bonding: Influence on Structural Organization and Membrane Function. Biochim. Biophys. Acta 1987, 906, 353–404. [Google Scholar] [CrossRef]
  3. Seu, K.J.; Cambrea, L.R.; Everly, R.M.; Hovis, J.S. Influence of Lipid Chemistry on Membrane Fluidity: Tail and Headgroup Interactions. Biophys. J. 2006, 91, 3727–3735. [Google Scholar] [CrossRef] [Green Version]
  4. Pasenkiewicz-Gierula, M.; Takaoka, Y.; Miyagawa, H.; Kitamura, K.; Kusumi, A. Hydrogen Bonding of Water to Phosphatidylcholine in the Membrane As Studied by a Molecular Dynamics Simulation: Location, Geometry, and Lipid-Lipid Bridging via Hydrogen-Bonded Water. J. Phys. Chem. A 1997, 101, 3677–3691. [Google Scholar] [CrossRef]
  5. Pyrkova, D.V.; Tarasova, N.V.; Krylov, N.A.; Nolde, D.E.; Pentkovsky, V.M.; Efremov, R.G. Dynamic Clustering of Lipids in Hydrated Two-Component Membranes: Results of Computer Modeling and Putative Biological Impact. J. Biomol. Struct. Dyn. 2013, 31, 87–95. [Google Scholar] [CrossRef]
  6. Lewis, R.N.A.H.; McElhaney, R.N. Calorimetric and Spectroscopic Studies of the Thermotropic Phase Behavior of Lipid Bilayer Model Membranes Composed of a Homologous Series of Linear Saturated Phosphatidylserines. Biophys. J. 2000, 79, 2043–2055. [Google Scholar] [CrossRef] [Green Version]
  7. Jiang, L.; Dávila-Contreras, E.M.; Beaven, A.H.; Im, W.; Klauda, J.B. An Extensive Simulation Study of Lipid Bilayer Properties with Different Head Groups, Acyl Chain Lengths, and Chain Saturations. Biochim. Biophys. Acta Biomembr. 2016, 1858, 3093–3104. [Google Scholar] [CrossRef]
  8. Fairn, G.D.; Schieber, N.L.; Ariotti, N.; Murphy, S.; Kuerschner, L.; Webb, R.I.; Grinstein, S.; Parton, R.G. High-Resolution Mapping Reveals Topologically Distinct Cellular Pools of Phosphatidylserine. J. Cell Biol. 2011, 194, 257–275. [Google Scholar] [CrossRef]
  9. Bobone, S.; Hilsch, M.; Storm, J.; Dunsing, V.; Herrmann, A.; Chiantia, S. Phosphatidylserine Lateral Organization Influences the Interaction of Influenza Virus Matrix Protein 1 with Lipid Membranes. J. Virol. 2017, 91, 10–1128. [Google Scholar] [CrossRef] [Green Version]
  10. Hirama, T.; Das, R.; Yang, Y.; Ferguson, C.; Won, A.; Yip, C.M.; Kay, J.G.; Grinstein, S.; Parton, R.G.; Fairn, G.D. Phosphatidylserine Dictates the Assembly and Dynamics of Caveolae in the Plasma Membrane. J. Biol. Chem. 2017, 292, 14292–14307. [Google Scholar] [CrossRef] [Green Version]
  11. Lucas, N.; Cho, W. Phosphatidylserine Binding Is Essential for Plasma Membrane Recruitment and Signaling Function of 3-Phosphoinositide-Dependent Kinase-1. J. Biol. Chem. 2011, 286, 41265–41272. [Google Scholar] [CrossRef] [Green Version]
  12. Husby, M.L.; Amiar, S.; Prugar, L.I.; David, E.A.; Plescia, C.B.; Huie, K.E.; Brannan, J.M.; Dye, J.M.; Pienaar, E.; Stahelin, R.V. Phosphatidylserine Clustering by the Ebola Virus Matrix Protein Is a Critical Step in Viral Budding. EMBO Rep. 2022, 23, e51709. [Google Scholar] [CrossRef]
  13. Zwaal, R.F.A.; Comfurius, P.; Bevers, E.M. Surface Exposure of Phosphatidylserine in Pathological Cells. Cell. Mol. Life Sci. 2005, 62, 971–988. [Google Scholar] [CrossRef]
  14. Ran, S.; Thorpe, P.E. Phosphatidylserine Is a Marker of Tumor Vasculature and a Potential Target for Cancer Imaging and Therapy. Int. J. Radiat. Oncol. Biol. Phys. 2002, 54, 1479–1484. [Google Scholar] [CrossRef]
  15. Riedl, S.; Rinner, B.; Asslaber, M.; Schaider, H.; Walzer, S.M.; Novak, A.; Lohner, K.; Zweytick, D. In Search of a Novel Target—Phosphatidylserine Exposed by Non-Apoptotic Tumor Cells and Metastases of Malignancies with Poor Treatment Efficacy. Biochim. Biophys. Acta Biomembr. 2011, 1808, 2638–2645. [Google Scholar] [CrossRef] [Green Version]
  16. Tan, L.T.H.; Chan, K.-G.; Pusparajah, P.; Lee, W.L.; Chuah, L.H.; Khan, T.M.; Lee, L.H.; Goh, B.H. Targeting Membrane Lipid a Potential Cancer Cure? Front. Pharmacol. 2017, 8, 12. [Google Scholar] [CrossRef] [Green Version]
  17. Karathanou, K.; Bondar, A.-N. Algorithm to Catalogue Topologies of Dynamic Lipid Hydrogen-Bond Networks. Biochim. Biophys. Acta Biomembr. 2022, 1864, 183859. [Google Scholar] [CrossRef]
  18. Jurkiewicz, P.; Cwiklik, L.; Vojtiskova, A.; Jungwirth, P.; Hof, M. Structure, Dynamics, and Hydration of POPC/POPS Bilayers Suspended in NaCl, KCl, and CsCl Solutions. Biochim. Biophys. Acta Biomembr. 2012, 1818, 609–616. [Google Scholar] [CrossRef] [Green Version]
  19. Boettcher, J.M.; Davis-Harrison, R.L.; Clay, M.E.; Nieuwkoop, A.J.; Ohkubo, Y.Z.; Tajkhorshid, E.; Morrissey, J.H.; Rienstra, C.M. Atomic View of Calcium-Induced Clustering of Phosphatidylserine in Mixed Lipid Bilayers. Biochemistry 2011, 50, 2264–2273. [Google Scholar] [CrossRef] [Green Version]
  20. Wang, Y.H.; Slochower, D.R.; Janmey, P.A. Counterion-Mediated Cluster Formation by Polyphosphoinositides. Chem. Phys. Lipids 2014, 182, 38–51. [Google Scholar] [CrossRef] [Green Version]
  21. Mukhopadhyay, P.; Monticelli, L.; Tieleman, D.P. Molecular Dynamics Simulation of a Palmitoyl-Oleoyl Phosphatidylserine Bilayer with Na+ Counterions and NaCl. Biophys. J. 2004, 86, 1601–1609. [Google Scholar] [CrossRef] [Green Version]
  22. Böckmann, R.A.; Hac, A.E.; Heimburg, T.; Grubmüller, H. Effect of Sodium Chloride on a Lipid Bilayer. Biophys. J. 2003, 85, 1647–1655. [Google Scholar] [CrossRef] [Green Version]
  23. Han, K.; Kim, S.H.; Venable, R.M.; Pastor, R.W. Design Principles of PI(4,5)P2 Clustering under Protein-Free Conditions: Specific Cation Effects and Calcium-Potassium Synergy. Proc. Natl. Acad. Sci. USA 2022, 119, e2202647119. [Google Scholar] [CrossRef]
  24. Karathanou, K.; Bondar, A.-N. Dynamic Water Hydrogen-Bond Networks at the Interface of a Lipid Membrane Containing Palmitoyl-Oleoyl Phosphatidylglycerol. J. Membr. Biol. 2018, 251, 461–473. [Google Scholar] [CrossRef]
  25. Filippov, A.V.; Orädd, G.; Lindblom, G. Effect of NaCl and CaCl2 on the Lateral Diffusion of Zwitterionic and Anionic Lipids in Bilayers. Chem. Phys. Lipids 2009, 159, 81–87. [Google Scholar] [CrossRef]
  26. Porasso, R.D.; Cascales, J.J.L. Study of the Effect of Na+ and Ca2+ Ion Concentration on the Structure of an Asymmetric DPPC/DPPC+DPPS Lipid Bilayer by Molecular Dynamics Simulation. Colloids Surf. B Biointerfaces 2009, 73, 42–50. [Google Scholar] [CrossRef]
  27. Mao, L.; Yang, L.; Zhang, Q.; Jiang, H.; Yang, H. Effects of Ion Interactions with a Cholesterol-Rich Bilayer. Biochem. Biophys. Res. Commun. 2015, 468, 125–129. [Google Scholar] [CrossRef]
  28. Bradley, R.; Slochower, D.R.; Janmey, P.A.; Radhakrishnan, R. Divalent Cations Bind to Phosphoinositides to Induce Ion and Isomer Specific Propensities for Nano-Cluster Initiation in Bilayer Membranes. R. Soc. Open Sci. 2020, 7, 192208. [Google Scholar] [CrossRef]
  29. Gurtovenko, A.A.; Vattulainen, I. Effect of NaCl and KCl on Phosphatidylcholine and Phosphatidylethanolamine Lipid Membranes: Insight from Atomic-Scale Simulations for Understanding Salt-Induced Effects in the Plasma Membrane. J. Phys. Chem. B 2008, 112, 1953–1962. [Google Scholar] [CrossRef] [Green Version]
  30. Yesylevskyy, S.O.; Demchenko, A.P. Cholesterol Behavior in Asymmetric Lipid Bilayers: Insights from Molecular Dynamics Simulations. In Methods in Molecular Biology; Springer Science + Business Media: New York, NY, USA, 2014; pp. 291–306. [Google Scholar]
  31. Silvius, J.R. Role of Cholesterol in Lipid Raft Formation: Lessons from Lipid Model Systems. Biochim. Biophys. Acta Biomembr. 2003, 1610, 174–183. [Google Scholar] [CrossRef] [Green Version]
  32. Simons, K.; Vaz, W.L.C. Model Systems, Lipid Rafts, and Cell Membranes. Annu. Rev. Biophys. Biomol. Struct. 2004, 33, 269–295. [Google Scholar] [CrossRef]
  33. Van Dijck, P.W.M. Negatively Charged Phospholipids and Their Position in the Cholesterol Affinity Sequence. Biochim. Biophys. Acta Biomembr. 1979, 555, 89–101. [Google Scholar] [CrossRef]
  34. Rowlands, L.J.; Marks, A.; Sanderson, J.E.; Law, R.V. 17O NMR Spectroscopy as a Tool to Study Hydrogen Bonding of Cholesterol in Lipid Bilayers. Chem. Commun. 2020, 56, 14499–14502. [Google Scholar] [CrossRef]
  35. Bhide, S.Y.; Zhang, Z.; Berkowitz, M.L. Molecular Dynamics Simulations of SOPS and Sphingomyelin Bilayers Containing Cholesterol. Biophys. J. 2007, 92, 1284–1295. [Google Scholar] [CrossRef] [Green Version]
  36. Hénin, J.; Chipot, C. Hydrogen-Bonding Patterns of Cholesterol in Lipid Membranes. Chem. Phys. Lett. 2006, 425, 329–335. [Google Scholar] [CrossRef]
  37. Chiu, S.C.; Jakobsson, E.; Mashl, R.J.; Scott, H.L. Cholesterol-Induced Modifications in Lipid Bilayers: A Simulation Study. Biophys. J. 2002, 83, 1842–1853. [Google Scholar] [CrossRef] [Green Version]
  38. Zocher, F.; Van Der Spoel, D.; Pohl, P.; Hub, J.S. Local Partition Coefficients Govern Solute Permeability of Cholesterol-Containing Membranes. Biophys. J. 2013, 105, 2760–2770. [Google Scholar] [CrossRef] [Green Version]
  39. Magarkar, A.; Dhawan, V.; Kallinteri, P.; Viitala, T.; Elmowafy, M.; Róg, T.; Bunker, A. Cholesterol Level Affects Surface Charge of Lipid Membranes in Saline Solution. Sci. Rep. 2014, 4, 5005. [Google Scholar] [CrossRef] [Green Version]
  40. Flanagan, J.C.; Valentine, M.L.; Baiz, C.R. Ultrafast Dynamics at Lipid–Water Interfaces. Acc. Chem. Res. 2020, 53, 1860–1868. [Google Scholar] [CrossRef]
  41. Higuchi, Y.; Asano, Y.; Kuwahara, T.; Hishida, M. Rotational Dynamics of Water at the Phospholipid Bilayer Depending on the Head Groups Studied by Molecular Dynamics Simulations. Langmuir 2021, 37, 5329–5338. [Google Scholar] [CrossRef]
  42. Efremov, R.G. Dielectric-Dependent Strength of Interlipid Hydrogen Bonding in Biomembranes: Model Case Study. J. Chem. Inf. Model. 2019, 59, 2765–2775. [Google Scholar] [CrossRef]
  43. Symons, J.; Cho, K.-J.; Chang, J.T.; Du, G.; Waxham, M.N.; Hancock, J.F.; Levental, I.; Levental, K.R. Lipidomic Atlas of Mammalian Cell Membranes Reveals Hierarchical Variation Induced by Culture Conditions, Subcellular Membranes, and Cell Lineages. Soft Matter 2021, 17, 288–297. [Google Scholar] [CrossRef]
  44. Jo, S.; Kim, T.H.; Iyer, V.; Im, W. CHARMM-GUI: A Web-Based Graphical User Interface for CHARMM. J. Comput. Chem. 2008, 29, 1859–1865. [Google Scholar] [CrossRef]
  45. Lee, J.; Cheng, X.; Swails, J.M.; Yeom, M.S.; Eastman, P.; Lemkul, J.A.; Wei, S.; Buckner, J.; Jeong, J.J.; Zhang, J.Z.H.; et al. CHARMM-GUI Input Generator for NAMD, GROMACS, AMBER, OpenMM, and CHARMM/OpenMM Simulations Using the CHARMM36 Additive Force Field. J. Chem. Theory Comput. 2015, 12, 405–413. [Google Scholar] [CrossRef]
  46. Wu, E.L.; Cheng, X.; Jo, S.; Rui, H.; Song, K.W.; Dávila-Contreras, E.M.; Zhang, J.Z.H.; Lee, J.; Monje-Galvan, V.; Venable, R.M.; et al. CHARMM-GUI Membrane Builder toward Realistic Biological Membrane Simulations. J. Comput. Chem. 2014, 35, 1997–2004. [Google Scholar] [CrossRef] [Green Version]
  47. Brooks, B.R.; Bruccoleri, R.E.; Olafson, B.D.; States, D.J.; Swaminathan, S.; Karplus, M. CHARMM: A Program for Macromolecular Energy, Minimization, and Dynamics Calculations. J. Comput. Chem. 1983, 4, 187–217. [Google Scholar] [CrossRef]
  48. Feller, S.E.; MacKerell, A.D. An Improved Empirical Potential Energy Function for Molecular Simulations of Phospholipids. J. Phys. Chem. B 2000, 104, 7510–7515. [Google Scholar] [CrossRef] [Green Version]
  49. Klauda, J.B.; Venable, R.M.; Freites, J.A.; O’Connor, J.K.; Tobias, D.J.; Mondragon-Ramirez, C.; Vorobyov, I.; MacKerell, A.D.; Pastor, R.W. Update of the CHARMM All-Atom Additive Force Field for Lipids: Validation on Six Lipid Types. J. Phys. Chem. B 2010, 114, 7830–7843. [Google Scholar] [CrossRef] [Green Version]
  50. Beglov, D.; Roux, B. Finite Representation of an Infinite Bulk System: Solvent Boundary Potential for Computer Simulations. J. Chem. Phys. 1994, 100, 9050–9063. [Google Scholar] [CrossRef] [Green Version]
  51. MacKerell, A.D.; Feig, M.; Brooks, C.L. Extending the Treatment of Backbone Energetics in Protein Force Fields: Limitations of Gas-Phase Quantum Mechanics in Reproducing Protein Conformational Distributions in Molecular Dynamics Simulations. J. Comput. Chem. 2004, 25, 1400–1415. [Google Scholar] [CrossRef]
  52. MacKerell, A.D.; Bashford, D.; Bellott, M.; Dunbrack, R.L.; Evanseck, J.D.; Field, M.J.; Fischer, S.M.; Gao, J.R.; Guo, H.; Ha, S.; et al. All-Atom Empirical Potential for Molecular Modeling and Dynamics Studies of Proteins. J. Phys. Chem. B 1998, 102, 3586–3616. [Google Scholar] [CrossRef]
  53. Jorgensen, W.L.; Chandrasekhar, J.; Madura, J.D.; Impey, R.; Klein, M.L. Comparison of Simple Potential Functions for Simulating Liquid Water. J. Chem. Phys. 1983, 79, 926–935. [Google Scholar] [CrossRef]
  54. Kale, L.V.; Skeel, R.D.; Bhandarkar, M.; Brunner, R.J.; Gursoy, A.; Krawetz, N.; Phillips, J.G.; Shinozaki, A.; Varadarajan, K.; Schulten, K. NAMD2: Greater Scalability for Parallel Molecular Dynamics. J. Comput. Phys. 1999, 151, 283–312. [Google Scholar] [CrossRef]
  55. Phillips, J.G.; Braun, R.; Wang, W.; Gumbart, J.C.; Tajkhorshid, E.; Villa, E.; Chipot, C.; Skeel, R.D.; Kale, L.V.; Schulten, K. Scalable Molecular Dynamics with NAMD. J. Comput. Chem. 2005, 26, 1781–1802. [Google Scholar] [CrossRef] [Green Version]
  56. Martyna, G.J.; Tobias, D.J.; Klein, M.L. Constant Pressure Molecular Dynamics Algorithms. J. Chem. Phys. 1994, 101, 4177–4189. [Google Scholar] [CrossRef] [Green Version]
  57. Feller, S.E.; Zhang, Y.; Pastor, R.W.; Brooks, B.R. Constant Pressure Molecular Dynamics Simulation: The Langevin Piston Method. J. Chem. Phys. 1995, 103, 4613–4621. [Google Scholar] [CrossRef]
  58. Darden, T.; York, D.M.; Pedersen, L.G. Particle Mesh Ewald: An(N)⋅log(N) Method for Ewald Sums in Large Systems. J. Chem. Phys. 1993, 98, 10089–10092. [Google Scholar] [CrossRef] [Green Version]
  59. Essmann, U.; Perera, L.; Berkowitz, M.L.; Darden, T.; Lee, H.F.; Pedersen, L.G. A Smooth Particle Mesh Ewald Method. J. Chem. Phys. 1995, 103, 8577–8593. [Google Scholar] [CrossRef]
  60. Grubmüller, H.; Heller, H.; Windemuth, A.; Schulten, K. Generalized Verlet Algorithm for Efficient Molecular Dynamics Simulations with Long-Range Interactions. Mol. Simul. 1991, 6, 121–142. [Google Scholar] [CrossRef]
  61. Tuckerman, M.E.; Berne, B.J.; Martyna, G.J. Reversible Multiple Time Scale Molecular Dynamics. J. Chem. Phys. 1992, 97, 1990–2001. [Google Scholar] [CrossRef] [Green Version]
  62. Karathanou, K.; Lazaratos, M.; Bertalan, É.; Siemers, M.; Buzar, K.; Schertler, G.F.X.; Del Val, C.; Bondar, A.-N. A Graph-Based Approach Identifies Dynamic H-Bond Communication Networks in Spike Protein S of SARS-CoV-2. J. Struct. Biol. 2020, 212, 107617. [Google Scholar] [CrossRef]
  63. Siemers, M.; Lazaratos, M.; Karathanou, K.; Guerra, F.; Brown, L.S.; Bondar, A.-N. Bridge: A Graph-Based Algorithm to Analyze Dynamic H-Bond Networks in Membrane Proteins. J. Chem. Theory Comput. 2019, 15, 6781–6798. [Google Scholar] [CrossRef]
  64. Siemers, M.; Bondar, A.-N. Interactive Interface for Graph-Based Analyses of Dynamic H-Bond Networks: Application to Spike Protein S. J. Chem. Inf. Model. 2021, 61, 2998–3014. [Google Scholar] [CrossRef]
  65. Humphrey, W.; Dalke, A.; Schulten, K. VMD: Visual Molecular Dynamics. J. Mol. Graph. 1996, 14, 33–38. [Google Scholar] [CrossRef]
  66. Karathanou, K. Graph-Based Algorithm for Common Topologies of Dynamic Lipid Clusters. Mendeley Data 2022, V2. [Google Scholar] [CrossRef]
  67. Guixà-González, R.; Giorgino, T.; Ramírez-Anguita, J.M.; Carrió-Gaspar, P.; Martinez-Seara, H.; Hansen, T.V.; Selent, J. MEMBPLUGIN: Studying Membrane Complexity in VMD. Comput. Appl. Biosci. 2014, 30, 1478–1480. [Google Scholar] [CrossRef]
  68. Shahane, G.; Ding, W.; Palaiokostas, M.; Orsi, M. Physical Properties of Model Biological Lipid Bilayers: Insights from All-Atom Molecular Dynamics Simulations. J. Mol. Model. 2019, 25, 76. [Google Scholar] [CrossRef] [Green Version]
  69. Pan, J.; Cheng, X.; Monticelli, L.; Katsaras, J.; Kučerka, N.; Tieleman, D.P.; Katsaras, J. The Molecular Structure of a Phosphatidylserine Bilayer Determined by Scattering and Molecular Dynamics Simulations. Soft Matter 2014, 10, 3716. [Google Scholar] [CrossRef]
  70. Skjevik, Å.A.; Madej, B.D.; Dickson, C.J.; Lin, C.P.; Teigen, K.; Walker, R.C.; Gould, I.R. Simulation of Lipid Bilayer Self-Assembly Using All-Atom Lipid Force Fields. Phys. Chem. Chem. Phys. 2016, 18, 10573–10584. [Google Scholar] [CrossRef] [Green Version]
  71. Róg, T.; Pasenkiewicz-Gierula, M. Cholesterol Effects on the Phospholipid Condensation and Packing in the Bilayer: A Molecular Simulation Study. FEBS Lett. 2001, 502, 68–71. [Google Scholar] [CrossRef] [Green Version]
  72. Venable, R.M.; Luo, Y.; Gawrisch, K.; Roux, B.; Pastor, R.W. Simulations of Anionic Lipid Membranes: Development of Interaction-Specific Ion Parameters and Validation Using NMR Data. J. Phys. Chem. B 2013, 117, 10183–10192. [Google Scholar] [CrossRef] [Green Version]
  73. Hung, W.-C.; Lee, M.-T.; Chen, F.-Y.; Huang, H.W. The Condensing Effect of Cholesterol in Lipid Bilayers. Biophys. J. 2007, 92, 3960–3967. [Google Scholar] [CrossRef] [Green Version]
  74. Alwarawrah, M.; Dai, J.S.; Huang, J. A Molecular View of the Cholesterol Condensing Effect in DOPC Lipid Bilayers. J. Phys. Chem. B 2010, 114, 7516–7523. [Google Scholar] [CrossRef] [Green Version]
  75. Doktorova, M.; Katsaras, J.; Kingston, R.L.; Khelashvili, G.; Cuendet, M.A.; Wen, Y.; Feigenson, G.W.; Vogt, V.M.; Dick, R.P. Cholesterol Promotes Protein Binding by Affecting Membrane Electrostatics and Solvation Properties. Biophys. J. 2017, 113, 2004–2015. [Google Scholar] [CrossRef] [Green Version]
  76. Bandara, A.; Panahi, A.; Pantelopulos, G.A.; Straub, J.E. Exploring the Structure and Stability of Cholesterol Dimer Formation in Multicomponent Lipid Bilayers. J. Comput. Chem. 2016, 38, 1479–1488. [Google Scholar] [CrossRef]
Biomolecules 13 01238 g001
Figure 1. Basic features of the hydrated POPS and POPS:chol membranes. (A) Time series and average values of the lipid bilayer thickness computed from POPS membranes (cyan profile) vs. POPS:chol (orange profile). (B) Time series of the area per lipid computed from POPS and POPS:chol membranes. (C) Molecular graphics based on coordinate snapshots from simulations of POPS (left) vs. POPS:chol (right). POPS are shown as yellow bonds; when present, cholesterol is shown as orange bonds. (D) SCD order parameter profiles. (E) Molecular graphics illustrating POPS with selected atoms labeled.
Figure 1. Basic features of the hydrated POPS and POPS:chol membranes. (A) Time series and average values of the lipid bilayer thickness computed from POPS membranes (cyan profile) vs. POPS:chol (orange profile). (B) Time series of the area per lipid computed from POPS and POPS:chol membranes. (C) Molecular graphics based on coordinate snapshots from simulations of POPS (left) vs. POPS:chol (right). POPS are shown as yellow bonds; when present, cholesterol is shown as orange bonds. (D) SCD order parameter profiles. (E) Molecular graphics illustrating POPS with selected atoms labeled.
Biomolecules 13 01238 g001
Biomolecules 13 01238 g002
Figure 2. Number of H-bonds per POPS. (A–C) Time series of interactions sampled in simulations of the POPS membrane without cholesterol. We show the number of direct lipid H-bonds (panel A) and 1-water-mediated (panel B) and ion-mediated bridges between POPS (panel C). The total number of H-bonds per POPS is shown by the blue profile. For groups of atoms of POPS, we use the following color code: serine, green; phosphate, red; ester, purple; phosphate and serine, orange; all, blue. For clarity, coordinate sets were read each 1 ns. (D–F) Interactions sampled in POPS:chol membrane simulations. We show the number of direct lipids H-bonds (panel C), 1-water-mediated (panel D) vs. ion-mediated bridges between POPS (panel E). (G–I) POPS-cholesterol interactions. We monitor the number of direct POPS-cholesterol H-bonds (panel G), 1-water-mediated (panel H) vs. ion-mediated bridges (panel I). (J) Molecular graphics illustrating, from left to right, direct POPS H-bonds, 1-water bridging, and potassium ion-mediated bridges.
Figure 2. Number of H-bonds per POPS. (A–C) Time series of interactions sampled in simulations of the POPS membrane without cholesterol. We show the number of direct lipid H-bonds (panel A) and 1-water-mediated (panel B) and ion-mediated bridges between POPS (panel C). The total number of H-bonds per POPS is shown by the blue profile. For groups of atoms of POPS, we use the following color code: serine, green; phosphate, red; ester, purple; phosphate and serine, orange; all, blue. For clarity, coordinate sets were read each 1 ns. (D–F) Interactions sampled in POPS:chol membrane simulations. We show the number of direct lipids H-bonds (panel C), 1-water-mediated (panel D) vs. ion-mediated bridges between POPS (panel E). (G–I) POPS-cholesterol interactions. We monitor the number of direct POPS-cholesterol H-bonds (panel G), 1-water-mediated (panel H) vs. ion-mediated bridges (panel I). (J) Molecular graphics illustrating, from left to right, direct POPS H-bonds, 1-water bridging, and potassium ion-mediated bridges.
Biomolecules 13 01238 g002
Biomolecules 13 01238 g003
Figure 3. H-bonds and water-mediated bridges have low occupancies. We present boxplots of H-bond occupancy values for interactions between phosphate, serine, and ester groups, and for all three groups (see Scheme 1A). (A–D) Occupancies (%) of direct H-bonds (panels A, C) and for one-water bridges (panels B, D). In panels A and C, note that direct H-bonds between phosphate or ester groups are not found. (E,F) Occupancies (%) of direct POPS-cholesterol H-bonds (panel E) and one-water bridges between POPS and cholesterol (panel F). Overall, the same pairs of headgroups tend to give the high-occupancy H-bonding and water bridging presented here for the distinct moieties.
Figure 3. H-bonds and water-mediated bridges have low occupancies. We present boxplots of H-bond occupancy values for interactions between phosphate, serine, and ester groups, and for all three groups (see Scheme 1A). (A–D) Occupancies (%) of direct H-bonds (panels A, C) and for one-water bridges (panels B, D). In panels A and C, note that direct H-bonds between phosphate or ester groups are not found. (E,F) Occupancies (%) of direct POPS-cholesterol H-bonds (panel E) and one-water bridges between POPS and cholesterol (panel F). Overall, the same pairs of headgroups tend to give the high-occupancy H-bonding and water bridging presented here for the distinct moieties.
Biomolecules 13 01238 g003
Biomolecules 13 01238 g004
Figure 4. POPS clusters sampled during simulations without and with cholesterol. (A,B) Bar graphs of the average number of lipids clusters (panel A) and of the average cluster size (panel B) computed separately for H-bond clusters mediated by direct POPS and POPS-cholesterol H-bonds, one-water bridges, and potassium ion bridges. (C) Illustration of lipid clusters mediated by direct POPS H-bonds in coordinate snapshots from simulations of POPS bilayers without (left) and with cholesterol (right). POPS headgroups are shown with yellow bonds, and cholesterol oxygen atoms are shown as purple spheres.
Figure 4. POPS clusters sampled during simulations without and with cholesterol. (A,B) Bar graphs of the average number of lipids clusters (panel A) and of the average cluster size (panel B) computed separately for H-bond clusters mediated by direct POPS and POPS-cholesterol H-bonds, one-water bridges, and potassium ion bridges. (C) Illustration of lipid clusters mediated by direct POPS H-bonds in coordinate snapshots from simulations of POPS bilayers without (left) and with cholesterol (right). POPS headgroups are shown with yellow bonds, and cholesterol oxygen atoms are shown as purple spheres.
Biomolecules 13 01238 g004
Biomolecules 13 01238 g005
Figure 5. Path lengths for linear clusters in the POPS and POPS:chol bilayer. (A,B) Percentage of number of linear clusters with different path lengths mediated by direct H-bonds (panel A) vs. one-water bridges (panel B) in POPS (blue) vs. POPS:chol (orange) bilayers. Note that most of these paths are transient, as most H-bonds and water-mediated bridges have small occupancies (Figure 3). (C) Water-mediated linear clusters with path length L = 4; the four connecting H-bonds are shown as thick green lines. (D,E) Time series of the number of direct (panel D) and one-water bridges (panel E) linear clusters with path lengths between 1 and 6, computed from simulations of POPS without cholesterol. (F,G) Time series of the number of direct (panel F) and one-water bridges (panel G) linear clusters with path lengths between 1 and 6, computed from simulations of POPS with cholesterol.
Figure 5. Path lengths for linear clusters in the POPS and POPS:chol bilayer. (A,B) Percentage of number of linear clusters with different path lengths mediated by direct H-bonds (panel A) vs. one-water bridges (panel B) in POPS (blue) vs. POPS:chol (orange) bilayers. Note that most of these paths are transient, as most H-bonds and water-mediated bridges have small occupancies (Figure 3). (C) Water-mediated linear clusters with path length L = 4; the four connecting H-bonds are shown as thick green lines. (D,E) Time series of the number of direct (panel D) and one-water bridges (panel E) linear clusters with path lengths between 1 and 6, computed from simulations of POPS without cholesterol. (F,G) Time series of the number of direct (panel F) and one-water bridges (panel G) linear clusters with path lengths between 1 and 6, computed from simulations of POPS with cholesterol.
Biomolecules 13 01238 g005
Biomolecules 13 01238 g006
Figure 6. Cholesterol H-bond clusters. (A) Time series of the number of linear clusters mediated by direct cholesterol–cholesterol H-bonds, one-water bridges, and ion bridges between cholesterol molecules with path lengths 1–2. (B) Illustration of linear paths mediated via one-water bridging between 2 and 3 cholesterol molecules. (C) Occupancy (%) of cholesterol–cholesterol H-bonds mediated by 2, 3, 4, or 5 water molecules. (D) Cluster size. (E) Illustration of linear cholesterol paths formed via water wires of length 2, 3, 4, or 5 water molecules.
Figure 6. Cholesterol H-bond clusters. (A) Time series of the number of linear clusters mediated by direct cholesterol–cholesterol H-bonds, one-water bridges, and ion bridges between cholesterol molecules with path lengths 1–2. (B) Illustration of linear paths mediated via one-water bridging between 2 and 3 cholesterol molecules. (C) Occupancy (%) of cholesterol–cholesterol H-bonds mediated by 2, 3, 4, or 5 water molecules. (D) Cluster size. (E) Illustration of linear cholesterol paths formed via water wires of length 2, 3, 4, or 5 water molecules.
Biomolecules 13 01238 g006
Biomolecules 13 01238 sch001
Scheme 1. H-bonds and lipid cluster topologies sampled in hydrated POPS bilayers. (A) H-bonding groups used for topology analysis and H-bond calculations of POPS. (B–D) Molecular graphics illustrating 1-water mediated H-bond interaction between POPS headgroups (Panel B), POPS-cholesterol (Panel C), and between two cholesterol molecules (Panel D). (E) Illustrations of the topologies studied here with the DFS algorithm [17]. H-bond clusters with path lengths L = 3 and L = 4 are shown.
Scheme 1. H-bonds and lipid cluster topologies sampled in hydrated POPS bilayers. (A) H-bonding groups used for topology analysis and H-bond calculations of POPS. (B–D) Molecular graphics illustrating 1-water mediated H-bond interaction between POPS headgroups (Panel B), POPS-cholesterol (Panel C), and between two cholesterol molecules (Panel D). (E) Illustrations of the topologies studied here with the DFS algorithm [17]. H-bond clusters with path lengths L = 3 and L = 4 are shown.
Biomolecules 13 01238 sch001
Table 1. MD simulations performed. We report the length of the production runs performed without any constraints. The total number of atoms we report includes atoms of the lipids and water molecules, and the ions.
Table 1. MD simulations performed. We report the length of the production runs performed without any constraints. The total number of atoms we report includes atoms of the lipids and water molecules, and the ions.
Sim#POPS/Chol Molecules per Leaflet#AtomsLength (μs)
POPSPOPS (109)53,6800.5
POPS:chol (10%)POPS/Chol (98/11)51,8590.5
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Jain, H.; Karathanou, K.; Bondar, A.-N. Graph-Based Analyses of Dynamic Water-Mediated Hydrogen-Bond Networks in Phosphatidylserine: Cholesterol Membranes. Biomolecules 2023, 13, 1238. https://doi.org/10.3390/biom13081238

AMA Style

Jain H, Karathanou K, Bondar A-N. Graph-Based Analyses of Dynamic Water-Mediated Hydrogen-Bond Networks in Phosphatidylserine: Cholesterol Membranes. Biomolecules. 2023; 13(8):1238. https://doi.org/10.3390/biom13081238

Chicago/Turabian Style

Jain, Honey, Konstantina Karathanou, and Ana-Nicoleta Bondar. 2023. "Graph-Based Analyses of Dynamic Water-Mediated Hydrogen-Bond Networks in Phosphatidylserine: Cholesterol Membranes" Biomolecules 13, no. 8: 1238. https://doi.org/10.3390/biom13081238

APA Style

Jain, H., Karathanou, K., & Bondar, A. -N. (2023). Graph-Based Analyses of Dynamic Water-Mediated Hydrogen-Bond Networks in Phosphatidylserine: Cholesterol Membranes. Biomolecules, 13(8), 1238. https://doi.org/10.3390/biom13081238

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop