Signatures of glassy dynamics in highly ordered lipid bilayers with emergence of soft dynamic channels

Harini Sureshkumar,Sahithya Iyer,Atreyee Banerjee,Pratyush Poduval,Edward Lyman,Anand Srivastava

DOI: https://doi.org/10.1101/2024.11.22.624930

2024-11-25

Abstract:Over the last few decades, extensive investigations on spatial and dynamic heterogeneity have been performed on carefully reconstituted biological lipid membranes. Characterising the molecular features in heterogeneous membranes is extremely challenging due to experimentally inaccessible time- and length-scales of these emergent systems. In this context, simulations can provide important insights into molecular-level interactions leading to membrane heterogeneity and associated functions. To that end, we use the non-affine displacement (NAD) framework (a concept borrowed from Physics of granular material) to faithfully capture molecular-scale local membrane order in simulated heterogeneous bilayers. In our latest application of NAD, we investigate the temperature-dependent spatial and temporal organisation on microseconds trajectories of liquid-ordered bilayer systems at all-atom resolution (DPPC/DOPC/CHOL: 0.55:0.15:0.30; 40 nm x 40 nm with a total of 5600 lipids and 2 million atoms). Lateral organisation in these large bilayer patches show noticeable dynamic heterogeneity despite their liquid-ordered nature. Moreover, our NAD analyses reveal soft fluid channels within the tightly packed membrane reminiscent of the classical two-component Kob-Andersen glass-forming binary mixture. Hence, we characterised these systems using classical glass physics markers for dynamic heterogeneities such as Overlap, Four-point Susceptibility, van Hove, and Intermediate Scattering functions to quantify the multiple time scales underlying the lipid dynamics. Our analyses reveal that highly ordered membrane systems can have glass-like dynamics with distinct soft fluid channels inside them. Biologically, these dynamic channels could act as conduits for facilitating molecular encounters for biological functions even in highly ordered phases such as lipid nanodomains and rafts.

Biophysics

What problem does this paper attempt to address?

The problem that this paper attempts to solve is to explore the characteristics of glassy dynamics and the emergence of soft dynamic channels in highly ordered lipid bilayer membranes. Specifically, the researchers are concerned with the significant dynamic heterogeneity in these highly ordered membrane systems, despite their liquid - ordered nature. In addition, they also found that soft fluid channels similar to those in classical binary mixtures emerged within these closely - packed membranes. To gain a deep understanding of these phenomena, the researchers used classical glass physics indicators, such as the Overlap function, Four - point Susceptibility, van Hove function and Intermediate Scattering Function, to quantify multiple time scales behind lipid dynamics.

### Research Background
Cell membranes are quasi - two - dimensional fluids composed of complex mixtures of lipids and proteins. In mammalian cell membranes, approximately 30 - 40 mole percent of the lipid mixture is cholesterol, an extremely important lipid whose chemical properties and interactions with other lipids determine a complex phase diagram. In model membranes containing cholesterol and two other lipids (usually saturated and unsaturated lipids), a liquid - liquid coexistence region, known as the liquid - ordered (L_o) and liquid - disordered (L_d) phases, is observed. This behavior is also observed in vesicles directly extracted from living cell membranes, which contain all the chemical complexity of cell membranes, including hundreds of different lipids and proteins.

### Research Methods
The researchers used the Non - Affine Displacement (NAD) framework to capture molecular - scale local membrane order in simulated heterogeneous bilayer membranes. NAD is a concept borrowed from the physics of amorphous materials and is used to identify viscoplastic deformation in amorphous materials. The researchers characterized the dynamic heterogeneity of these systems through the following methods:

1. **Mean - Square Displacement (MSD)**:
\[
\langle \Delta r(t)^2 \rangle=\left\langle\frac{1}{NP} \sum_{i = 1}^{N}|r_i(t)-r_i(0)|^2\right\rangle
\]
where \(N\) is the number of particles and \(r(t)\) represents the position at time \(t\). The diffusion index \(\alpha(t)\) is defined as:
\[
\alpha(t)=\frac{d(\ln(\langle \Delta r(t)^2 \rangle))}{d(\ln(\Delta t))}
\]

2. **Self - Intermediate Scattering Function**:
\[
F_s(q, t)=\frac{1}{N} \sum_{i = 1}^{N}\left\langle\exp[-i\mathbf{q} \cdot(\mathbf{r}_i(t)-\mathbf{r}_i(0))]\right\rangle
\]
where the wave vector \(q = \frac{2\pi}{\lambda}\).

3. **Overlap Function**:
\[
Q_s(t)=\sum_{i = 1}^{N} w(|\mathbf{r}_i(t)-\mathbf{r}_i(0)|)
\]
where \(w(|\mathbf{r}_1-\mathbf{r}_2|)\) is 1 when the particle spacing is less than or equal to the distance \(a\), and 0 otherwise.

4. **Dynamic Susceptibility**:
\[
\chi_4(t)=\beta V[\langle Q(t)^2 \rangle-\langle Q(t) \rangle^2]