📄 synthetic biology

Complex-phase stochastic modeling of mitochondrial heteroplasmy

This paper introduces a novel complex-phase stochastic model that represents mitochondrial heteroplasmy as a complex number to simultaneously track copy numbers and mutant load, thereby enabling a quantitative, personalized assessment of disease risk by integrating selection, genetic drift, and tissue migration dynamics.

Original authors: Nurbaev, S., Pocheshkhova, E.

Published 2026-06-09

📖 3 min read☕ Coffee break read

CC BY 4.0

Original authors: Nurbaev, S., Pocheshkhova, E.

Original paper licensed under CC BY 4.0 (https://creativecommons.org/licenses/by/4.0/). ⚕️ This is an AI-generated explanation of a preprint that has not been peer-reviewed. It is not medical advice. Do not make health decisions based on this content. Read full disclaimer

Imagine your body is a massive city made up of billions of tiny neighborhoods (cells). Inside each neighborhood, there are thousands of tiny power plants called mitochondria. These power plants have their own instruction manuals, known as mitochondrial DNA (mtDNA).

Usually, all these manuals are perfect copies of the "wild-type" (healthy) version. But sometimes, some manuals get typos or damage, becoming "mutant" versions. When a single cell contains a mix of both healthy and damaged manuals, it's called heteroplasmy. Think of it like a library where some books are pristine and others are torn up.

The Problem with Old Methods
Scientists used to measure this mix by just counting the percentage of bad books. If 30% of the books were damaged, they'd say the risk is 30%. But the paper argues this is too simple. It's like trying to predict a traffic jam just by counting cars, without knowing how fast they are moving, if they are changing lanes, or if the road is getting narrower. The old method misses the "stochastic" (random) nature of how these manuals shuffle around and the fact that different tissues (like the brain vs. the skin) react differently to the same amount of damage.

The New "Complex-Phase" Idea
The authors propose a clever new way to look at this using a bit of math magic. Instead of just a single number, they describe the state of a cell's power plants using a complex number (a number with two parts: a real part and an imaginary part).

The Analogy: Imagine a clock face.
- The length of the hand pointing from the center represents the total number of power plant manuals (both healthy and damaged).
- The angle (or "phase") of the hand represents the ratio of damaged manuals to healthy ones.
- If the hand points straight up, you have a mix. If it points far to the right, you have mostly healthy manuals. If it swings far to the left, you have a dangerous load of damaged ones.

This "phase" angle is the key. It captures both how many manuals there are and how bad the mix is in one single, elegant measurement.

The Simulation: A Random Walk
The researchers built a computer model (running 1,000 simulations) to watch how this "clock hand" moves over time. They included three main forces:

Selection: Some manuals are better at copying themselves than others.
Genetic Drift: Random chance causes the mix to fluctuate, like shuffling a deck of cards.
Migration: Manuals can move between different neighborhoods (tissues) in the body.

What They Found
When they ran the simulations, they discovered something surprising about neuronal tissues (the brain and nerves). Even if the body's overall average of damaged manuals is low, the brain is like a volatile neighborhood where the mix can swing wildly by chance. The "clock hand" can randomly swing into the "danger zone" (the pathological threshold) much more easily in the brain than in other tissues.

The Conclusion
The paper concludes that getting sick isn't a simple "yes or no" switch based on a fixed percentage. Instead, it's a time-to-event process, similar to waiting for a bus that might arrive at any random moment. By using this new "complex-phase" clock method, doctors can better calculate the specific, individual risk of a person developing a mitochondrial disease, acknowledging that the brain is particularly sensitive to these random swings.

Technical Summary: Complex-phase stochastic modeling of mitochondrial heteroplasmy

Problem Statement
Mitochondrial heteroplasmy, defined as the coexistence of wild-type and mutant mitochondrial DNA (mtDNA) within a single cell, is a critical determinant in the pathogenesis of mitochondrial diseases. Current classical approaches rely exclusively on the scalar fraction of mutant DNA to model these dynamics. However, the paper identifies significant limitations in these traditional methods: they fail to adequately account for threshold effects, the inherent stochastic nature of heteroplasmy dynamics, and the specific variations observed across different tissues. Consequently, existing models are insufficient for a precise quantitative assessment of disease risk.

Methodology
To address these gaps, the authors propose a novel complex-phase formalism. In this framework, the state of a cell's mitochondrial genome is represented by a complex number $Z = a + ib$, where:

$a$ represents the absolute number of normal mtDNA copies.
$b$ represents the absolute number of mutant mtDNA copies.

This mathematical representation allows for the simultaneous integration of total copy number and heteroplasmy level. The argument of this complex number, $\phi = \arctan(b/a)$ , is interpreted as a "phase" that characterizes the mutant load.

Building on this formalism, the authors developed a stochastic model of tissue dynamics that explicitly incorporates three biological processes:

Selection: The differential replication or survival of mtDNA variants.
Genetic Drift: The random fluctuations in mtDNA copy numbers.
Intertissue Migration: The movement of mitochondrial genomes between different tissues.

The model was evaluated using Monte Carlo simulations (specifically 1,000 simulations) to analyze the probabilistic behavior of the system. Additionally, Kaplan-Meier survival analysis was employed to treat the onset of pathology as a time-to-event process.

Key Results
The simulation results yielded several specific findings regarding tissue-specific dynamics:

Neuronal Tissue Variability: Neuronal tissues were found to exhibit high heteroplasmy variability.
Threshold Crossing: There is a significant probability that neuronal tissues will reach a pathological threshold even when the systemic mutant load is relatively low.
Probabilistic Pathology: The development of pathology is not deterministic but probabilistic. The Kaplan-Meier analysis confirms that disease onset can be effectively described as a time-to-event process, dependent on the stochastic accumulation of mutant loads over time.

Significance and Claims
The paper claims that this complex-phase approach offers a more comprehensive framework than scalar models by naturally combining copy number and heteroplasmy level information. The primary significance of this work lies in its ability to enable a quantitative assessment of the individual risk of developing mitochondrial diseases. By modeling the stochastic nature of heteroplasmy and tissue-specific migration, the approach opens the door to personalized prognosis, moving beyond general population averages to individualized risk prediction. The authors present this as a foundational step toward understanding the complex mechanisms driving mitochondrial disease pathogenesis.

Technical Summary: Complex-phase stochastic modeling of mitochondrial heteroplasmy

More like this