Alpha-hemolysin is an exotoxin secreted by the bacterium Staphylococcus aureus. In its active form, the protein is built from 293 water-soluble amino acids that assemble into a heptameric pore, a seven-part barrel that punches through cell membranes. That structure, first resolved in crystallographic detail in the 1990s, has made alpha-hemolysin the workhorse of stochastic sensing: the first and still most widely used biological nanopore for detecting individual molecules. When a single pore is embedded in an artificial membrane bathed in electrolyte and a voltage is applied, a steady ionic current flows through the opening. When an analyte molecule wanders into the pore, it partially obstructs that current, producing a characteristic downward blip. Each blip corresponds to one molecule, and the pattern of blips over time carries information about the analyte’s identity and concentration.

In the new experiments, the team used monodisperse polyethylene glycol of molecular weight 1294, abbreviated PEG 1294, as the analyte. The experimental platform was a solvent-free planar lipid bilayer assembled by the classic Montal and Mueller technique: two monomolecular lipid films were apposed across a hole roughly 100 micrometers in diameter in a Teflon partition, forming a barrier that separated two identical compartments of a Teflon chamber. The synthetic lipid 1,2-diphytanoyl-sn-glycero-3-phosphocholine was used at 2 percent by weight in hexane, and the hole was pretreated with a dilute hexadecane solution to promote stable bilayer formation. Both compartments held a concentrated electrolyte of 4 molar potassium chloride buffered with 5 millimolar Tris-citric acid at pH 7.5. Alpha-hemolysin was added to one side, the cis side, at a concentration tuned so that exactly one pore inserted into the membrane.

The researchers then added PEG 1294 to the opposite, trans side at two concentrations, 400 and 1000 micromolar, and applied voltages ranging from +20 to +100 millivolts in 20 millivolt steps. Ionic currents were recorded with an Axopatch 200B amplifier in voltage clamp mode, using silver/silver chloride electrodes connected through agar salt bridges. The signals were low-pass filtered at 15 kilohertz and digitized at a sampling rate of 250 kilohertz, all at a controlled temperature of 25 degrees Celsius. Under these conditions, the current trace showed the textbook signature of stochastic sensing: stepwise transitions between a fully open level and a partially blocked level, each step corresponding to the entry or exit of a single PEG molecule from the pore lumen.

The novelty of the study lay not in the electrophysiology but in the analysis. Rather than focusing on conventional metrics such as event duration or blocking depth, the team turned to Lempel-Ziv complexity, a measure introduced in 1976 by Abraham Lempel and Jacob Ziv to quantify the randomness of finite sequences. The method works by converting a time series into a binary string: each data point is compared with the average of the entire series, receiving a 1 if it exceeds the average and a 0 if it falls below. The resulting string of ones and zeros is then scanned from left to right, and a counter increments every time a new, previously unseen substring appears. Finally, the counter is normalized by a theoretical upper bound, n divided by the base-2 logarithm of n, yielding a complexity value that ranges from 0 to 1 and is independent of sequence length.

The interpretation of this number is intuitive. A value close to 1 indicates a highly random, unpredictable signal with little internal self-similarity, while a value close to 0 indicates a repetitive, ordered series. Lempel-Ziv complexity has previously found applications in biomedical signal analysis, including studies of electroencephalogram background activity in Alzheimer’s disease patients and the detection of ventricular tachycardia and fibrillation in cardiac recordings. Applying it to nanopore data, however, allowed the researchers to ask a question that conventional event statistics cannot easily answer: how complex, in an information-theoretic sense, is the pattern of molecular blocking events?

The answer depended strongly on how much analyte was present. At the lower PEG concentration of 400 micromolar, complexity values were consistently lower than at 1000 micromolar, across every voltage tested. The full range of measured values stretched from a minimum of 0.45, recorded at 400 micromolar and +40 millivolts, to a maximum of 0.97, recorded at 1000 micromolar and +80 millivolts. The physical explanation is straightforward: higher analyte concentrations produce more frequent blocking events, filling the time series with rapid alternations between open and blocked states and pushing the complexity measure toward its upper limit. At lower concentrations, longer stretches of open-pore current dominate, and the signal becomes more ordered and less random.

To probe whether the apparent order at low concentration reflected genuine structure or simply sparse sampling, the researchers performed a clever control. They segmented the original time series recorded at 400 micromolar, shuffled the segments into a random order, and recalculated the complexity of the resulting randomized series. For every applied potential, the shuffled data yielded higher Lempel-Ziv values than the original recordings. This result is diagnostic: if the original series had already been maximally random, shuffling it would have changed nothing. The fact that randomization increased complexity means the unshuffled sequences contained residual structure, and that the biosensor’s output at low analyte concentration falls short of full stochastic behavior.

The team also compared series built exclusively from blocked periods with series built exclusively from unblocked periods, at both concentrations and all voltages. Statistical analysis using Student’s t-test, with a significance threshold of p less than 0.05, found no meaningful differences between the complexity of the blocked-only and unblocked-only series, nor between those segmented series and the original combined recordings. In other words, the complexity of the signal does not reside preferentially in either the moments when PEG occupies the pore or the intervals between occupations; the stochastic character is distributed across the entire record.

The broader implications reach into a long-standing debate in biology about the role of randomness. Molecular biology has traditionally leaned on deterministic principles, yet stochasticity is now recognized as a genuine feature of living systems rather than mere noise. Ion channel kinetics are widely treated as stochastic events, and some theorists argue that intrinsic randomness in the nervous system may even enable flexible decision-making. As the authors note, biological studies often lack the tools to distinguish chaotic, noisy, deterministic, and probabilistic behavior. Lempel-Ziv complexity offers one such tool. Their conclusion is nuanced: the alpha-hemolysin biosensor behaves in a purely stochastic fashion only at higher analyte concentrations, while at lower concentrations its stochasticity diminishes, though it can be mathematically restored by randomizing the event sequence. For engineers designing nanopore sensors, that finding matters, because the reliability of single-molecule detection depends on understanding when the digital current signature is truly random and when it carries hidden structure that smarter algorithms might exploit.

Subject of Research: Stochastic analysis of alpha-hemolysin nanopore ionic current blocking events using Lempel-Ziv complexity

Article Title: Analyzing alpha-hemolysin nanopore behavior using Lempel-Ziv complexity

Article References: Neves, G. F., Machado, D. C., Consoni, L. H. A., Costa, E. V. L., Carneiro, C. M. M., Rodrigues, C. G., & Nogueira, R. A. (2026). Analyzing alpha-hemolysin nanopore behavior using Lempel-Ziv complexity. Discover Chemistry, 3(1), Article 571. https://doi.org/10.1007/s44371-026-01022-8

Image Credits: AI Generated

DOI: 10.1007/s44371-026-01022-8

Keywords: alpha-hemolysin, nanopore, Lempel-Ziv complexity, stochastic sensing, polyethylene glycol, ion channel, planar lipid bilayer, biosensor, information theory, single-molecule detection, biophysics, stochasticity

The work, published in Nature by Chenyu Tang, Mayank Prakash Pandey and colleagues at the Université de Lorraine and CNRS in Nancy, France, with collaborators in Madrid, Chicago and Urbana-Champaign, addresses a stubborn bottleneck in computational biology. Standard molecular dynamics simulations propagate atoms step by step according to the laws of classical physics, and they do so faithfully. But the transitions that matter most, such as protein folding, allosteric signaling and membrane transport, are rare events. A simulation might churn through millions of time steps without ever witnessing the molecule cross from one stable state to another, because the probability of stumbling into the transition region by chance is vanishingly small.

Computational scientists have long responded with enhanced-sampling strategies, which apply artificial forces or cleverly designed biases to push molecules over energetic barriers. These methods have produced genuine insights, but they come with serious caveats. They are computationally demanding, and they typically require the researcher to choose, in advance, a set of collective variables, essentially a guess about which combination of atomic coordinates describes the reaction. If that guess is wrong or incomplete, the resulting free-energy landscape and apparent mechanism can be distorted. A recent study cited by the team put the problem bluntly in its title: convergence is not correctness, and enhanced-sampling methods can perform inconsistently as biological complexity increases.

Gen-COMPAS takes a different route. The framework couples a denoising diffusion probabilistic model, the same class of generative architecture that powers modern image synthesis, with a rigorous statistical concept from transition state theory known as the committor. The diffusion model is trained to produce structurally plausible intermediate conformations that sit between the two known end points of a transition. These generated structures are not accepted on faith. Instead, they are filtered using committor-based criteria, which identify the configurations from which a molecule is equally likely to fall forward to the product state or backward to the reactant state, the hallmark of a true transition state.

The committor is a powerful but demanding quantity. For any given molecular configuration, it answers a deceptively simple question: if the system were released from this point with random thermal velocities, what is the probability that it would commit to the product rather than the reactant? Configurations with a committor of one half define the transition state ensemble. Estimating this probability traditionally requires launching many short trajectories from each candidate configuration, which is exactly the kind of brute-force calculation that makes rare-event statistics so expensive. Gen-COMPAS turns this into an iterative loop: the generative model proposes candidate intermediates, short unbiased simulations from those intermediates test whether they genuinely belong to the transition region, and the results refine both the model and the ensemble.

The payoff is a dramatic compression of the required sampling. According to the authors, short unbiased simulations launched from the generated intermediates yield transition-region ensembles at aggregate sampling scales ranging from nanoseconds to submicroseconds, whereas conventional approaches would require orders of magnitude more simulation time to assemble equivalent statistics. In practical terms, the method replaces the hopeless task of waiting for a rare event to happen spontaneously with the tractable task of verifying that generated candidate structures are genuine members of the transition ensemble.

To demonstrate the framework, the team applied it to three systems of increasing size and biological complexity. The first was Trp-cage, a miniprotein of roughly twenty residues that has become a standard benchmark for folding studies because it folds quickly enough to be simulated directly. Decades of work, including landmark long-trajectory studies of fast-folding proteins, have mapped its folding landscape in detail, giving the team an independent yardstick against which to judge their results. Gen-COMPAS recovered committors, transition states and free-energy landscapes consistent with this established knowledge, starting only from the folded and unfolded end-point structures.

The second test system was the mitochondrial ADP/ATP carrier, a membrane transporter that shuttles the energy currency ATP across the inner mitochondrial membrane in exchange for ADP. This protein operates through a switching mechanism between conformational states that has been studied with structural biology and a variety of simulation techniques, but its full transition pathway remains challenging to sample because the carrier is embedded in a lipid membrane and undergoes large-scale domain movements. Gen-COMPAS reconstructed the transition pathway between the carrier’s known states without any predefined reaction coordinate, offering a mechanistic picture that previously would have demanded extensive biased simulations and careful, potentially contentious, choices of collective variables.

The third and most ambitious application was the nicotinic acetylcholine receptor, a pentameric ligand-gated ion channel that mediates synaptic transmission. Channel activation involves coordinated conformational changes across five protein subunits, and recent work has suggested that these subunit transitions occur asynchronously, priming the receptor for activation in a stepwise fashion. Simulating the full gating transition of such a large membrane protein is far beyond the reach of straightforward molecular dynamics. Gen-COMPAS nevertheless recovered the transition pathway and free-energy landscape from the known resting and activated structures alone, with no prior mechanistic knowledge encoded into the calculation.

The significance of the method lies not only in its speed but in what it does not require. No collective variables need to be chosen by hand, no reaction coordinates need to be guessed, and no prior hypothesis about the mechanism needs to be supplied. The generative model and the committor filter together let the data define the pathway. This removes one of the most persistent sources of bias in enhanced sampling, the researcher’s own assumptions about how the transition ought to proceed, and it opens the door to systems where no expert intuition exists to guide variable selection.

The method also arrives amid a broader convergence of machine learning and molecular simulation. Generative models have been used to emulate equilibrium protein ensembles, to sample molecular dynamics trajectories directly, and to guide transition path sampling in systems ranging from lipid flip-flop to membrane nanoporation. Gen-COMPAS distinguishes itself by combining generative proposal with a physically rigorous acceptance criterion rooted in transition rate theory, ensuring that the structures it produces are not merely plausible but statistically meaningful members of the true transition ensemble. The authors acknowledge prior machine-guided path-sampling work and committor-learning approaches as foundations for this synthesis.

Accessibility may prove as important as performance. Gen-COMPAS is released as an installable open-source package on GitHub, complete with a unified command-line interface, ready-to-use configuration files, molecular-simulation inputs for representative systems, and a graphical user interface that guides users through the available options. Input files for all systems analyzed in the paper are publicly available, and source data accompany the publication. The team also gratefully acknowledges reference molecular dynamics trajectories provided by D. E. Shaw Research, whose specialized hardware has long set the benchmark for long unbiased simulations that methods like this one can now be validated against.

For the field, the implications extend across drug discovery, basic biophysics and beyond. Ion channels and membrane transporters are among the most heavily targeted protein families in pharmacology, and understanding the exact pathway by which they switch states can reveal where small molecules might bind to stabilize one state or another. Allosteric transitions in signaling proteins, another class of rare events, could similarly be mapped without the guesswork that has historically limited such studies. If the framework generalizes as its demonstration systems suggest, the computational microscope that biophysicists have been building for half a century may finally be able to zoom in on the fleeting moments that matter most, the instants in which a molecule decides what it will become next.

Subject of Research: Generative committor-guided path sampling of biomolecular conformational transitions

Article Title: Breaking timescales with generative sampling of conformational transitions

Article References: Breaking timescales with generative sampling of conformational transitions. (n.d.). https://doi.org/10.1038/s41586-026-11025-1

Image Credits: AI Generated

DOI: 10.1038/s41586-026-11025-1

Keywords: Gen-COMPAS, molecular dynamics, transition path sampling, denoising diffusion model, committor, protein folding, enhanced sampling, ion channel, membrane transporter, free-energy landscape, machine learning, computational biophysics