2026-01-22 | | Total: 7
The relationship between the physical characteristics of the radiation field and biological damage is central to both radiotherapy and radioprotection, yet the link between spatial scales of energy deposition and biological effects remains not entirely understood. To address this, we developed an interpretable deep learning model that predicts cell survival after proton and carbon ion irradiation, leveraging sequential attention to highlight relevant features and provide insight into the contribution of different energy deposition scales. Trained and tested on the PIDE dataset, our model incorporates, beside LET, nanodosimetric and microdosimetric quantities simulated with MC-Startrack and Open-TOPAS, enabling multi-scale characterization. While achieving high predictive accuracy, our approach also emphasizes transparency in decision-making. We demonstrate high accuracy in predicting RBE for in vitro experiments. Multiple scales are utilized concurrently, with no single spatial scale being predominant. Quantities defined at smaller spatial domains generally have a greater influence, whereas the LET plays a lesser role.
The Hodgkin-Huxley equations have described neuronal excitability for seventy years, yet their mathematical structure-gating exponents m3h and n4, exponential voltage dependencies, and bounded activation variables, has remained empirically justified rather than theoretically derived. Hodgkin and Huxley introduced voltage-dependent conductances controlled by gating variables. While these equations reproduce experimental observations, they were derived through curve-fitting without theoretical justification. Modern theoretical physics derives governing equations from symmetry principles through Lie group theory. We prove that the complete Hodgkin-Huxley equations necessarily follow from three fundamental symmetries: (1) compact conformational state spaces, (2) multiplicative conductance scaling, and (3) temporal translation invariance. These symmetries uniquely determine a Lie group structure isomorphic to SO(2) semidirect product with R2. From representation theory, we derive: boundedness from SO(2) compactness, exponential Boltzmann factors from scale invariance, specific integer exponents m3h and n4 from irreducible representations, and first-order kinetics from Lie algebra flows. This demonstrates that the HH equations are not empirical curve-fits but the unique mathematical structure mandated by fundamental symmetries. We reveal why gating variables must be bounded, voltage dependencies must be exponential, sodium requires three activation gates and one inactivation gate, potassium requires four activation gates, and kinetics must be first-order. This establishes that neural electrophysiology obeys the same theoretical framework as modern physics, where symmetries determine dynamics, providing a foundation for understanding channel mutations and network dynamics through group theory.
Developing epithelial tissues coordinate cell proliferation and mechanical forces to achieve proper size and shape. As epithelial cells tightly adhere together to form the confluent tissue, the distribution of cell areas significantly influences possible patterns of cellular packing and thereby also the mechanics of the epithelium. Therefore, it is important to understand the origin of cell area heterogeneity in developing tissues and, if possible, how to control it. Previous models of cell growth and division have been successful in accounting for experimentally observed area distributions in cultured cells and bacterial colonies, but developing tissues present additional complexity due to self-organized patterns of mechanical stresses that guide morphogenesis. Here, we address this challenge focusing on the D. melanogaster wing disc epithelium. We consider a simple model that couples cell cycle dynamics to tissue mechanics. From time-lapse imaging of the cellular network, we extract all model parameters - cell growth rates, division rates, and mechanical fluctuations - revealing that they all depend on cell size. With these independently measured parameters, our model quantitatively reproduces the observed cell area distribution without any fitting parameters and further predicts tissue pressure gradients, in quantitative agreement with previously published data. Importantly, we find that cell proliferation accounts for 85% of cell area variance, establishing it as the dominant source of packing disorder that influences tissue mechanics and organization.
Proton minibeam (pMB) radiotherapy, delivers highly heterogeneous dose distributions alternating high-dose peaks and low-dose valleys. This aims to widen the therapeutic window by improving normal tissue sparing while maintaining the same or even better tumour control. The performance of pMB strongly depends on the collimator design and physical parameters. To better understand the physical and radiobiological drivers of this enhanced therapeutic window, we perform a detailed microdosimetric characterization of proton minibeams and assess their impact. We characterize radiation quality with microdosimetry through Monte Carlo simulations. Then we extend the Generalized Stochastic Microdosimetric Model to predict the normal tissue complication probability (NTCP) at different depths in water, 1cm, 2cm, and 4cm, for 100MeV proton minibeams realized with varying configurations of collimator. Results are compared with conventional homogeneous field (HF) irradiation after dose normalization to the tumor. The developed model is applied by considering tissues as divided into several functional subunits, connected by a seriality parameter. Microdosimetric characterization of proton minibeam irradiation shows differences between peak and valley regions in shaping lineal energy spectra, especially at low depth, while radiation quality uniforms progressively getting closer to the tumor. NTCP calculations results suggest an increased sparing effect for pMB over conventional HF. A strong dependence is found on the peak-to-valley dose ratio (PVDR), and on the seriality parameter. Predictions indicate substantial sparing from pMB, especially for PVDR > 15, including relatively serial organs with seriality around 0.7. This integrated dose-microdosimetric-biological framework elucidates how spatial fractionation, radiation quality, and organ architecture collectively shape tissue sparing in pMB.
LiDRoSIS is an automated MATLAB-Python software suite for the segmentation and quantification of lipid droplets (LDs) and reactive oxygen species (ROS) in fluorescence microscopy images of irradiated A549 and MCF7 cells exposed to gold-based nanoparticles. It combines classical image processing algorithms with statistical post-analysis through a companion Python tool, StatLysis. The platform enables reproducible, high-throughput analysis of morphological and spectral parameters linked to nanoparticle-enhanced radiobiological responses. By bridging imaging and quantitative analytics, LiDRoSIS provides a robust framework for nanomedicine and radiation biology research.
Cadaveric studies have provided important insights into the mechanics of the human foot arch and plantar fascia. However, repeatedly probing posture-dependent viscoelastic responses immediately after landing impact is difficult in biological specimens, leaving the contribution of skeletal architecture to landing dynamics incompletely understood. In this study, we developed an anthropomimetic foot joint structure aimed at replicating the skeletal geometry of the human foot. Using a vertical drop apparatus that simulates landing and a viscoelastic system-identification model, we investigated how skeletal structure and posture modulate the apparent post-impact viscoelastic response. The results show that the multi-jointed anthropomimetic structure exhibited a higher damping ratio than simplified flat and rigid feet. Moreover, ankle dorsiflexion and toe extension systematically shifted the identified parameters, reducing the damping ratio under the tested conditions. Taken together, these findings indicate that an arch-like, multi-jointed skeletal architecture can enhance impact attenuation in an anthropomimetic mechanical foot, and that morphology and passive posture alone can tune the trade-off between attenuation and rebound. The observed posture-dependent trends are qualitatively consistent with reported differences in human landing strategies, suggesting that skeletal architecture may partly account for the modulation. Furthermore, these results highlight the engineering advantage of anatomically informed skeletal replication for achieving human-like apparent viscoelastic behavior through postural adjustment during landing.
In contrast to inert materials, living cells use molecular motors to generate forces independently of elastic strain. How do local active forces translate into large-scale shape, and what are the mechanical properties of the resulting "living matter"? Here, we address these questions within the active tension network (ATN) model for the mechanics of 2d epithelia. We represent the configuration of active forces geometrically by a tension triangulation dual to the cell tessellation. The Voronoi dual of the tension triangulation is shown to be macroscopically stress-free -- the tensions hence define an emergent reference state for the tissue. Two soft modes, curl-free and conformal deformations, map this reference to the internally stressed, but force-balanced, physical configuration of the tissue. Conformal deformations parametrize pressure gradients, leading to a generalization of von Neumann's law for the pressure in a foam. Via finite-element-like interpolation and a notion of "discrete" conformal maps, we both construct the mechanically balanced cell tessellations exactly on the microscopic level, and pass to the continuum limit. We systematically incorporate cell rearrangement into our theory by representing cell-scale topology in terms of circle packings. The results of this bottom-up coarse-graining study match a top-down continuum analysis presented in a companion paper. Thus, the geometry of force balance bridges between micro- and macro-scales, elucidating how cell-level active forces program shape. The present formalism may be useful in the study of foams, granular matter, and metamaterials, as well as in numerical simulations.