🔬 condensed matter

Effects of Divalent Cations on Diffusion Dynamics of Biological Water Confined between Lipid Membranes

Using molecular dynamics simulations, this study reveals that calcium and magnesium ions exert contrasting, concentration-dependent effects on the diffusion dynamics of water confined between lipid membranes, driven by their distinct hydration radii and impacts on interfacial water structure.

Original authors: Minho Lee, Jinwon Park, Ji-Hyun Kim, Minhaeng Cho, Jaeyoung Sung

Published 2026-03-23

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Minho Lee, Jinwon Park, Ji-Hyun Kim, Minhaeng Cho, Jaeyoung Sung

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine a cell membrane not as a solid wall, but as a bustling, crowded hallway between two rows of people (the lipid molecules). In the middle of this hallway flows "biological water"—the lifeblood of the cell. Usually, we think of water as just a fluid that flows smoothly, but inside a cell, it's more like a chaotic dance floor where the water molecules are constantly bumping into ions (charged particles) and the membrane walls.

This paper investigates what happens to this water dance when you add two specific types of "dancers" to the floor: Calcium (Ca²⁺) and Magnesium (Mg²⁺). Even though both are divalent cations (they have a double positive charge), they behave like two completely different characters in a play, causing the water to move in opposite ways.

Here is the breakdown of their effects using simple analogies:

1. The Two Characters: Calcium vs. Magnesium

Think of the ions as guests at a party who interact with the water and the membrane walls differently based on their "hydration radius" (how big their personal space bubble is).

Calcium (Ca²⁺) is the "Tight-Grip Dancer."
- The Analogy: Calcium has a small personal bubble. It sheds some of its water coat to get very close to the membrane walls. It grabs onto the wall's "headgroups" (the lipid heads) tightly, almost like a hug.
- The Effect: Because it hugs the walls so tightly, it pulls the walls closer together. Imagine two rows of people in a hallway suddenly hugging each other; the hallway gets narrower. This squeezes the water out of the tight corners and forces it into the open, spacious middle of the hallway.
- The Result: Since the water is pushed into the open, spacious middle (the "bulk" region), it has more room to run. As you add more Calcium, the water flows faster and faster.
Magnesium (Mg²⁺) is the "Bulky Dancer."
- The Analogy: Magnesium has a huge personal bubble (a large hydration shell). It holds onto its water coat tightly and refuses to get close to the membrane walls. It stays in the middle of the hallway, floating like a big, heavy balloon.
- The Effect: Because it stays in the middle, it doesn't squeeze the walls together. Instead, it crowds the open space. It acts like a giant obstacle in the middle of a running track.
- The Result: At first, adding a little Magnesium might clear a path, but as you add more, the hallway becomes a traffic jam. The water molecules get stuck behind these bulky Magnesium balloons. As you add more Magnesium, the water flow slows down and gets chaotic.

2. The "Traffic Jam" vs. The "Highway"

The researchers measured how fast the water molecules moved (diffusion) and how predictably they moved.

With Calcium: It's like opening a highway. As you add more Calcium, the "squeeze" pushes all the water into the fast lane. The water moves faster and more predictably.
With Magnesium: It's like a traffic jam that gets worse the more cars you add. The water molecules are constantly bumping into the bulky Magnesium ions. The movement becomes erratic and slower.

3. The "Random Walk" and the "Memory"

Water molecules don't just move in a straight line; they wiggle and bounce. The paper looked at how long it takes for this wiggling to settle into a smooth, predictable pattern (called a "Gaussian distribution").

The Finding: In both cases, adding salt makes the water take longer to "calm down" and move smoothly. It's like adding more people to a dance floor; it takes longer for everyone to find a rhythm.
The Difference: With Magnesium, the water gets stuck in a "super-Gaussian" state (very erratic movement) for a long time because the bulky ions create a very uneven environment. With Calcium, the environment becomes more uniform (all water is in the middle), so the chaos settles down differently.

4. The "Hidden Variable" (The Z-Axis)

The paper explains that the main reason for these weird behaviors is that the water isn't the same everywhere.

Near the walls: The water is slow and sticky.
In the middle: The water is fast and free.
Calcium's Trick: It forces the water to move from the "slow sticky zone" (near the walls) to the "fast free zone" (the middle).
Magnesium's Trick: It stays in the "fast free zone" and clogs it up, making the whole system slower.

Summary: Why Does This Matter?

Cells rely on water moving quickly to send signals, transport nutrients, and keep the cell alive.

If a cell needs to speed up a process, it might use Calcium to squeeze the membranes and push the water into the fast lane.
If a cell needs to slow down or stabilize a structure, Magnesium might be the tool, acting as a heavy anchor that slows the flow.

In a nutshell: This paper reveals that not all "salt" is created equal. Calcium acts like a traffic cop directing cars to a highway, making water flow faster. Magnesium acts like a giant roadblock, making water flow slower and more chaotically. Understanding this helps us understand how cells control their internal environment and how diseases related to ion imbalances might disrupt these delicate flows.

1. Problem Statement

Biological water is an ionic solution containing both monovalent and divalent ions, which play critical roles in cellular functions such as membrane fusion, signal transmission, and enzyme activity. While the structural effects of divalent cations (specifically Calcium, $Ca^{2+}$ , and Magnesium, $Mg^{2+}$ ) on lipid membranes are well-documented, their influence on the transport dynamics of water molecules nanoconfined between these membranes remains poorly understood.

Previous studies have shown that $Ca^{2+}$ and $Mg^{2+}$ induce distinct structural changes in phosphatidylcholine (PC) membranes (e.g., $Ca^{2+}$ condenses membranes while $Mg^{2+}$ does not), but the microscopic mechanisms governing how these ions alter water diffusion and the non-Gaussian nature of water displacement in confined biological environments were elusive.

2. Methodology

The authors employed a multi-faceted approach combining All-Atom Molecular Dynamics (MD) simulations with a Generalized Transport Equation framework.

Simulation System:
- Membrane Model: Two Dimyristoylphosphatidylcholine (DMPC) lipid bilayers.
- Solvent: Water nanoconfined between the bilayers with a water-to-lipid ratio of 37:1.
- Ions: Systems were simulated with varying concentrations (0 M, 0.1 M, 0.4 M, 0.8 M) of $CaCl_2$ and $MgCl_2$ .
- Force Fields: AMBER Lipid21 for lipids, SPC/E for water, and Li-Merz 12-6-4 parameters for ions.
- Protocol: Simulations were performed at 318 K using the AMBER 21 package. Production runs included short-time (10 fs sampling) and long-time (1 ps sampling) trajectories up to 1 $\mu$ s.
Theoretical Framework:
- The study utilized a generalized transport equation for biological water to analyze the lateral displacement distribution $P(r_\parallel, t)$ .
- Key metrics extracted included the Mean Square Displacement (MSD), the Non-Gaussian Parameter (NGP), and the Diffusion Kernel Correlation (DKC).
- The authors decomposed the system into seven discrete layers along the $z$ -axis (perpendicular to the membrane) to analyze layer-dependent diffusion coefficients ( $D(n)$ ) and population fractions.
- Trapped Water Identification: Water molecules exhibiting long-time subdiffusive behavior (exponent $\alpha < 0.75$ ) were identified and excluded from the main ensemble average to focus on the bulk-like transport dynamics.

3. Key Contributions

Quantification of Ion-Specific Dynamics: The study provides the first comprehensive atomistic explanation for the contrasting effects of $Ca^{2+}$ and $Mg^{2+}$ on the lateral diffusion of confined water.
Decoupling Structure and Dynamics: It links the macroscopic diffusion behavior to microscopic structural changes, specifically the "squeeze-out" effect of water from interfacial regions and the redistribution of ions.
Theoretical Extraction of Fluctuations: By applying the generalized transport equation, the authors successfully extracted the time-correlation function of diffusion coefficient fluctuations and the relative variance of these fluctuations, linking them to the longitudinal thermal motion of water molecules.

4. Key Results

A. Divergent Diffusion Trends

$Ca^{2+}$ (Calcium): The mean lateral diffusion coefficient ( $D_\parallel$ ) monotonically increases with increasing $CaCl_2$ concentration.
$Mg^{2+}$ (Magnesium): The mean lateral diffusion coefficient exhibits a non-monotonic dependence, initially increasing slightly at low concentrations (0.1 M) but then decreasing significantly as concentration rises to 0.4 M and 0.8 M.

B. Non-Gaussian Behavior

The deviation of water displacement from a Gaussian distribution (measured by the peak height of the NGP) behaves oppositely for the two ions:
- $Mg^{2+}$ : The NGP peak height increases with concentration, indicating growing dynamic heterogeneity.
- $Ca^{2+}$ : The NGP peak height decreases with concentration, suggesting a more homogeneous dynamic environment.
The relaxation time of the displacement distribution toward a Gaussian increases for both ions as salt concentration rises, but this slowing is more pronounced for $Mg^{2+}$ .

C. Mechanistic Origins (The "Why")

The contrasting behaviors are attributed to the hydration radii and binding characteristics of the ions:

$Ca^{2+}$ Mechanism (Squeeze-Out Effect):
- $Ca^{2+}$ has a small hydration radius and binds tightly to lipid headgroups, bridging them and condensing the membrane (reducing Area Per Lipid).
- This condensation effectively "squeezes out" water molecules from the interfacial region (Layers 1-3) into the central bulk region (Layer 4).
- Since the bulk region has a higher diffusion coefficient than the interfacial region, the shift in water population toward the bulk causes the overall mean diffusion to increase.
- The relative variance of diffusion fluctuations decreases because the system becomes more dominated by the high-diffusion bulk water.
$Mg^{2+}$ Mechanism (Hydration Shell Retention):
- $Mg^{2+}$ has a large hydration radius and retains a strong hydration shell, preventing it from penetrating the lipid headgroup region as deeply as $Ca^{2+}$ .
- Consequently, $Mg^{2+}$ ions remain concentrated in the bulk region rather than the interface.
- The high concentration of ions in the bulk region slows down the water molecules there (due to ion-water friction).
- Since the bulk region is where most water resides, the overall mean diffusion decreases as $Mg^{2+}$ concentration rises.
- The relative variance increases because the heterogeneity between the slow bulk water (high ion density) and the interfacial water increases.

D. Source of Fluctuations

The study identifies that the primary source of lateral diffusion coefficient fluctuations is the thermal motion of water molecules in the longitudinal ( $z$ ) direction. As water molecules move between layers with different microscopic environments (lipid functional groups, ion concentrations, hydrogen bonding), their local diffusion coefficients fluctuate. The relaxation of these fluctuations is governed by the time it takes for water to traverse the intermembrane space.

5. Significance

Biological Relevance: The findings explain how cells might regulate water transport and membrane dynamics by modulating specific divalent ion concentrations. The distinct behaviors of $Ca^{2+}$ and $Mg^{2+}$ suggest that biological systems can fine-tune membrane fluidity and water permeability through ion selection.
Methodological Advancement: The successful application of the generalized transport equation to extract diffusion kernel correlations from MD simulations provides a robust framework for analyzing complex, heterogeneous fluids where simple Fickian diffusion models fail.
Resolution of Experimental Discrepancies: The results align with and explain previous Quasi-Elastic Neutron Scattering (QENS) experiments that observed opposite effects of $Ca^{2+}$ and $Mg^{2+}$ on hydration water populations, offering a microscopic mechanism (ion-specific redistribution of water between bulk and interfacial layers) for these observations.