Quantifying the Heterogeneous Dynamics of a Simulated Dipalmitoylphosphatidylcholine (DPPC) Membrane

Nanoscale lipid domains determine the dynamic molecular portraits of mixed DOPC/DOPS bilayers in a fluid phase: A computational insight

Lateral heterogeneity, or mosaicity, is a fundamental property inherent to cell membranes that is crucial for their functioning. While microscopic inhomogeneities (e.g. rafts) are easily detected experimentally, lipid domains with nanoscale dimensions (nanoclusters of nanodomains, NDs) resist reliable characterization by instrumental methods. In such a case, important insight can be gained via computer modeling. Here, NDs composed of lipid's head groups in the mixed zwitterionic dioleoylphosphatidylcholine (DOPC) and negatively charged dioleoylphosphatidylserine (DOPS) bilayers were studied by molecular dynamics. A new algorithm has been developed to identify NDs. Unlike most similar methods, it implicitly considers the heterogeneous distribution of lipid head atomic density and does not require subjectively chosen parameters. In DOPS-rich membranes, lipids form more compact and stable NDs due to strong interlipid interactions. In DOPC-rich systems, NDs arise due to the "packing" effect of weakly bound lipid heads. The clustering picture is related to the physical properties of the bilayer surface: DOPS-rich systems show more pronounced surface heterogeneity of hydrophilic/hydrophobic regions compared to DOPC-rich ones. The results obtained are important for the effective quantitative characterization of the "dynamic molecular portrait" of a membrane surface - its "fingerprint" characterizing dynamical distribution of its physicochemical properties.

Relaxation time scales of interfacial water upon fluid to ripple to gel phase transitions of bilayers

The slow relaxation of interface water (IW) across three primary phases of membranes is relevant to understand the influence of IW on membrane functions at supercooled conditions. To this objective, a total of ∼16.26μs all-atom molecular dynamics simulations of 1,2-dimyristoyl-sn-glycerol-3-phosphocholine lipid membranes are carried out. A supercooling-driven drastic slow-down in heterogeneity time scales of the IW is found at the fluid to the ripple to the gel phase transitions of the membranes. At both fluid-to-ripple-to-gel phase transitions, the IW undergoes two dynamic crossovers in Arrhenius behavior with the highest activation energy at the gel phase due to the highest number of hydrogen bonds. Interestingly, the Stokes-Einstein (SE) relation is conserved for the IW near all three phases of the membranes for the time scales derived from the diffusion exponents and the non-Gaussian parameters. However, the SE relation breaks for the time scale obtained from the self-intermediate scattering functions. The behavioral difference in different time scales is universal and found to be an intrinsic property of glass. The first dynamical transition in the α relaxation time of the IW is associated with an increase in the Gibbs energy of activation of hydrogen bond breaking with locally distorted tetrahedral structures, unlike the bulk water. Thus, our analyses unveil the nature of the relaxation time scales of the IW across membrane phase transitions in comparison with the bulk water. The results will be useful to understand the activities and survival of complex biomembranes under supercooled conditions in the future.

Quantifying dynamical heterogeneity length scales of interface water across model membrane phase transitions

All-atom molecular dynamics simulations of 1,2-dimyristoyl-sn-glycero-3-phosphocholine lipid membranes reveal a drastic growth in the heterogeneity length scales of interface water (IW) across fluid to ripple to gel phase transitions. It acts as an alternate probe to capture the ripple size of the membrane and follows an activated dynamical scaling with the relaxation time scale solely within the gel phase. The results quantify the mostly unknown correlations between the spatiotemporal scales of the IW and membranes at various phases under physiological and supercooled conditions.

Dehydration induced dynamical heterogeneity and ordering mechanism of lipid bilayers

Understanding the influence of dehydration on the membrane structure is crucial to control membrane functionality related to domain formation and cell fusion under anhydrobiosis conditions. To this end, we perform all-atom molecular dynamic simulations of 1,2-dimyristoyl-sn-glycero-3-phosphocholine dimyristoylphosphatidylcholine lipid membranes at different hydration levels at 308 K. As dehydration increases, the lipid area per head group decreases with an increase in bilayer thickness and lipid order parameters indicating bilayer ordering. Concurrently, translational and rotational dynamics of interfacial water (IW) molecules near membranes slow down. On the onset of bilayer ordering, the IW molecules exhibit prominent features of dynamical heterogeneity evident from non-Gaussian parameters and one-dimensional van Hove correlation functions. At a fully hydrated state, diffusion constants (D) of the IW follow a scaling relation, D∼τ_α ^-1, where the α relaxation time (τ_α) is obtained from self-intermediate scattering functions. However, upon dehydration, the relation breaks and the D of the IW follows a power law behavior as D∼τ_α ^-0.57, showing the signature of glass dynamics. τ_α and hydrogen bond lifetime calculated from intermittent hydrogen bond auto-correlation functions undergo a similar crossover in association with bilayer ordering on dehydration. The bilayer ordering is accompanied with an increase in fraction of caged lipids spanned over the bilayer surface and a decrease in fraction of mobile lipids due to the non-diffusive dynamics. Our analyses reveal that the microscopic mechanism of lipid ordering by dehydration is governed by dynamical heterogeneity. The fundamental understanding from this study can be applied to complex bio-membranes to trap functionally relevant gel-like domains at room temperature.

Quantitative Evidence of Mobile Ion Hopping in Polymerized Ionic Liquids

Atomistic molecular dynamics simulations were performed, and an extensive set of analyses were undertaken to understand the ion transport mechanism in the polymerized ionic liquid poly(C₂VIm)Tf₂N. The ion hopping events were investigated at different time scales. Ion hopping was examined by monitoring the instantaneous cation-anion association and dissociation. Ion diffusion was subsequently evaluated with correlation functions and the calculation of relaxation times at different time scales. Dynamical heterogeneity in the mobility of the ions was observed with only a small portion of the anions classified as fast mobile ions. The mobile ions were characterized as the ones traveling farther than a certain distance during a characteristic period, which was much longer than the time scale of the instant ion pair dissociation. Effective hopping of the mobile ions contributed to the diffusivity which was dominated by interchain hopping and generally facilitated with five associating cations from two different polymer chains. Mobile anions had relatively fewer associating cations from more associating chains than immobile anions. The stringlike cooperative motion was observed in the mobile anions. The string length was determined to decrease with increasing temperature. These findings provided an in-depth understanding of the ion transport in polymerized ionic liquids and important information for the rational design of novel materials.

Dynamic coupling of a hydration layer to a fluid phospholipid membrane: intermittency and multiple time-scale relaxations

Understanding the coupling of a hydration layer and a lipid membrane is crucial to gaining access to membrane dynamics and understanding its functionality towards various biological processes.

Quantification of spatio-temporal scales of dynamical heterogeneity of water near lipid membranes above supercooling

A hydrated 1,2-dimyristoyl-sn-glycero-3-phosphorylcholine (DMPC) lipid membrane is investigated using an all atom molecular dynamics simulation at 308 K to determine the physical sources of universal slow relaxations of hydration layers and length-scale of the spatially heterogeneous dynamics. Continuously residing interface water (IW) molecules hydrogen bonded to different moieties of lipid heads in the membrane are identified. The non-Gaussian parameters of all classes of IW molecules show a cross-over from cage vibration to translational diffusion. A significant non-Gaussianity is observed for the IW molecules exhibiting large length correlations in translational van Hove functions. Two time-scales for the ballistic motions and hopping transitions are obtained from the self intermediate scattering functions of the IW molecules with an additional long relaxation, which disappears for bulk water. The long relaxation time-scales for the IW molecules obtained from the self intermediate scattering functions are in good accordance with the hydrogen bond relaxation time-scales irrespective of the nature of the chemical confinement and the confinement lifetime. Employing a block analysis approach, the length-scale of dynamical heterogeneities is captured from a transition from non-Gaussianity to Gaussianity in van Hove correlation functions of the IW molecules. The heterogeneity length-scale is comparable to the wave-length of the small and weak undulations of the membrane calculated by Fourier transforms of lipid tilts. This opens up a new avenue towards a possible correlation between heterogeneity length-scale and membrane curvature more significant for rippled membranes. Thus, our analyses provide a measure towards the spatio-temporal scale of dynamical heterogeneity of confined water near membranes.

Universal nature of dynamic heterogeneity in glass-forming liquids: A comparative study of metallic and polymeric glass-forming liquids

Glass-formation is a ubiquitous phenomenon that is often observed in a broad class of materials ranging from biological matter to commonly encountered synthetic polymer, as well as metallic and inorganic glass-forming (GF) materials. Despite the many regularities in the dynamical properties of GF materials, the structural origin of the universal dynamical properties of these materials has not yet been identified. Recent simulations of coarse-grained polymeric GF liquids have indicated the coexistence of clusters of mobile and immobile particles that appear to be directly linked, respectively, to the rate of molecular diffusion and structural relaxation. The present work examines the extent to which these distinct types of "dynamic heterogeneity" (DH) arise in metallic GF liquids (Cu-Zr, Ni-Nb, and Pd-Si alloys) having a vastly different molecular structure and chemistry. We first identified mobile and immobile particles and their transient clusters and found the DH in the metallic alloys to be remarkably similar in form to polymeric GF liquids, confirming the "universality" of the DH phenomenon. Furthermore, the lifetime of the mobile particle clusters was found to be directly related to the rate of diffusion in these materials, while the lifetime of immobile particles was found to be proportional to the structural relaxation time, providing some insight into the origin of decoupling in GF liquids. An examination of particles having a locally preferred atomic packing, and clusters of such particles, suggests that there is no one-to-one relation between these populations of particles so that an understanding of the origin of DH in terms of static fluid structure remains elusive.

Complex dynamics at the nanoscale in simple biomembranes

Nature is known to engineer complex compositional and dynamical platforms in biological membranes. Understanding this complex landscape requires techniques to simultaneously detect membrane re-organization and dynamics at the nanoscale. Using super-resolution stimulated emission depletion (STED) microscopy coupled with fluorescence correlation spectroscopy (FCS), we reveal direct experimental evidence of dynamic heterogeneity at the nanoscale in binary phospholipid-cholesterol bilayers. Domain formation on the length scale of ~200–600 nm due to local cholesterol compositional heterogeneity is found to be more prominent at high cholesterol content giving rise to distinct intra-domain lipid dynamics. STED-FCS reveals unique dynamical crossover phenomena at length scales of ~100–150 nm within each of these macroscopic regions. The extent of dynamic heterogeneity due to intra-domain hindered lipid diffusion as reflected from the crossover length scale, is driven by cholesterol packing and organization, uniquely influenced by phospholipid type. These results on simple binary model bilayer systems provide novel insights into pathways leading to the emergence of complex nanodomain substructures with implications for a wide variety of membrane mediated cellular events.

Emergence of compositionally tunable nanoscale dynamical heterogeneity in model binary lipid biomembranes

While the existence of nanoscale dynamical heterogeneity in biological membranes has been suggested to act as an active functional platform for enabling various cellular processes like signal transduction and viral or bacterial entry, it has been extremely difficult to detect the existence of such domains. Model lipid bilayer membranes have been widely used to detect such dynamical heterogeneity in order to avoid complications arising from the compositional heterogeneity of cellular membranes. However, even in model biological membranes the issue of nanoscale lipid dynamics has remained controversial and unresolved due to the difficulty of detecting the existence of such dynamical heterogeneity on the scale of 10-300 nm. Here we report direct evidence of nanoscale lipid dynamical heterogeneity in model binary lipid bilayer membranes using a combination of super-resolution stimulated emission depletion (STED) microscopy and fluorescence correlation spectroscopy (FCS). We control the phase behavior of the lipid bilayers by varying their composition and discuss how this leads to the emergence of dynamical lipid domains on the scale of 80-150 nm, which is also dependent on the lipid phase in which such dynamics are observed. Notably, our work shows that the presence of cholesterol is not required for the existence of such domains even in fluid like bilayers, as has been widely believed, and specifies the minimal conditions required for the emergence of such dynamical heterogeneity in cellular membranes. Our work will thus not only be of great significance towards understanding the nanoscale dynamic organizing principles of cellular membranes but could also be useful in understanding the dynamics of related soft matter systems and nanoparticle-cell membrane interactions.

Comparative Study of the Collective Dynamics of Proteins and Inorganic Nanoparticles

Molecular dynamics simulations of ubiquitin in water/glycerol solutions are used to test the suggestion by Karplus and coworkers that proteins in their biologically active state should exhibit a dynamics similar to 'surface-melted' inorganic nanoparticles (NPs). Motivated by recent studies indicating that surface-melted inorganic NPs are in a 'glassy' state that is an intermediate dynamical state between a solid and liquid, we probe the validity and significance of this proposed analogy. In particular, atomistic simulations of ubiquitin in solution based on CHARMM36 force field and pre-melted Ni NPs (Voter-Chen Embedded Atom Method potential) indicate a common dynamic heterogeneity, along with other features of glass-forming (GF) liquids such as collective atomic motion in the form of string-like atomic displacements, potential energy fluctuations and particle displacements with long range correlations ('colored' or 'pink' noise), and particle displacement events having a power law scaling in magnitude, as found in earthquakes. On the other hand, we find the dynamics of ubiquitin to be even more like a polycrystalline material in which the α-helix and β-sheet regions of the protein are similar to crystal grains so that the string-like collective atomic motion is concentrated in regions between the α-helix and β-sheet domains.

Temperature dependent heterogeneous rotational correlation in lipids

Lipid structures exhibit complex and highly dynamic lateral structure; and changes in lipid density and fluidity are believed to play an essential role in membrane targeting and function. The dynamic structure of liquids on the molecular scale can exhibit complex transient density fluctuations. Here the lateral heterogeneity of lipid dynamics is explored in free standing lipid monolayers. As the temperature is lowered the probes exhibit increasingly broad and heterogeneous rotational correlation. This increase in heterogeneity appears to exhibit a critical onset, similar to those observed for glass forming fluids. We explore heterogeneous relaxation in in a single constituent lipid monolayer of 1, 2-dimyristoyl-sn-glycero-3-phosphocholine  by measuring the rotational diffusion of a fluorescent probe (1-palmitoyl-2-[1]-sn-glycero-3-phosphocholine), which is embedded in the lipid monolayer at low labeling density. Dynamic distributions are measured using wide-field time-resolved fluorescence anisotropy. The observed relaxation exhibits a narrow, liquid-like distribution at high temperatures (τ ∼ 2.4 ns), consistent with previous experimental measures (Dadashvand et al 2014 Struct. Dyn. 1 054701, Loura and Ramalho 2007 Biochim. Biophys. Acta 1768 467-478). However, as the temperature is quenched, the distribution broadens, and we observe the appearance of a long relaxation population (τ ∼ 16.5 ns). This supports the heterogeneity observed for lipids at high packing densities, and demonstrates that the nanoscale diffusion and reorganization in lipid structures can be significantly complex, even in the simplest amorphous architectures. Dynamical heterogeneity of this form can have a significant impact on the organization, permeability and energetics of lipid membrane structures.