Equilibrium state dynamics‐based modeling of temporal dose delivery dependencies of FLASH skin sparing

TL;DR

This study introduces a new model to understand how FLASH radiotherapy spares normal skin tissue based on how radiation is delivered over time.

Contribution

The novel perturbed equilibrium state dynamics (ESD) framework quantitatively models temporal dose delivery effects on FLASH skin sparing.

Findings

An ESD model with three parameters accurately reproduced 52 averaged FMF values with low error.

The model captures experimentally observed features like increased sparing with dose and saturation at high doses.

Three ESD variants showed strong support, while others were less favored based on BIC metrics.

Abstract

FLASH radiotherapy shows promise in sparing normal tissues while maintaining tumor control, with the magnitude of its sparing effect being strongly influenced by the temporal dose delivery (TDD) structure of ultra‐high dose rate (UHDR) irradiation. Quantitatively describing these dependencies is critical for optimizing UHDR treatments and guiding preclinical and clinical applications.

This study introduces perturbed equilibrium state dynamics (ESD) as a minimal phenomenological framework to describe TDD dependencies of FLASH normal tissue sparing.

ESD‐based modeling assumes a transient perturbation of a generic equilibrium state during irradiation, modulating the biological system's instantaneous radiosensitivity. Isoeffective dose ratios (FLASH‐modifying factors, FMF) were derived from a collection of acute skin toxicity data, consisting of both previously published and newly acquired data, from 721 mice across 60 experimental groups irradiated with pulsed electron beams and scanned proton beams spanning various TDD patterns. Eight ESD model variants were evaluated based on their ability to reproduce experimental FMF values, optimizing their three or four free parameters and assessing goodness‐of‐fit metrics.

The best‐ranked model based on the Bayesian Information Criterion (BIC) used three parameters to reproduce accurately 52 averaged FMF values (mean absolute error: 0.034, mean absolute percentage error: 4.3%). It consistently reproduces experimentally observed features of FLASH sparing in terms of FMF, including an increased sparing with dose followed by saturation at high doses. Two additional model variants received similar support by BIC, while the remaining five variants were less favored.

ESD‐based modeling offers an effective framework to describe TDD dependencies of FLASH murine skin sparing while limiting complexity and number of free parameters to avoid overfitting. Through parameter tuning, the ESD modeling aligns with key experimental data, supporting its potential as a predictive, parsimonious framework for experimental and translational research. By highlighting the potential role of a transient equilibrium perturbation modulating instantaneous radiosensitivity and quantifying respective measures, ESD‐based models may also guide the exploration and modeling of underlying mechanisms. While the present application is tissue‐specific, the approach is broadly adaptable, with clinical translation requiring refitting and validation in other contexts.

Figures

(a) Schematic representation: experimental dose ratios for isoeffective normal tissue complication probability (NTCP), referred to as experimental FLASH modifying factors (FMFexp), were computed utilizing the experimental data point for the test temporal dose delivery pattern (𝐷test) and the reference dose (𝐷ref) obtained from a logistic fit to the dose‐response curve for the respective reference temporal dose delivery pattern. (b) NTCP for skin ulcers (>3 mm) after conventional dose rate (CDR) and ultra‐high dose rate (UHDR) electron irradiation (i.e., datasets e1 and e2). A logistic NTCP curve (line) was fitted to the experimental single mouse data for CDR data.

(a–c) Changes in normal tissue sparing for a fixed single fraction dose ((a) 10 Gy, 26 (b) 10 Gy, 25 (c) 11.2 Gy (9–10 weeks)/12.5 Gy (30–31 weeks) 4 ) as a function of exposure time for different in vivo assays. Data points were digitized from the original publications. 4 , 25 , 26 Dashed grey lines (CTRL) correspond to reference level for control groups. (d) Schematic graph illustrating perturbed equilibrium state dynamics of a generic quantity 𝑄(𝑡) that may be pushed away from its equilibrium state 𝑄0 by an irradiation at instantaneous dose rate D˙(t) inducing a change Δ(𝑄). 𝑄(𝑡) returns toward equilibrium 𝑄0 at a given recovery rate R˙(Q). (e) Evaluated dependencies for the instantaneous relative radiosensitivity Srel(𝑄) include a Heaviside function, a power function, and a Michaelis–Menten (M–M) type behavior. For details, see text and Table 2.

Simulated and experimental isoeffective dose ratios (FLASH‐modifying factors, FMF) and corresponding residuals for model variant ESD1_RRS1 as a function of dose (a) and as a function of time‐averaged dose rate (TADR) (b). For the residuals, dashed horizontal lines indicate ± 0.05.

Four representative temporal irradiation profiles from the data collection are shown (left) alongside the simulated dynamics for 𝑄(𝑡) and 𝑆rel(𝑡) for the ESD1_RRS1 model variant (right). The corresponding accumulated RRS‐weighted dose, 𝐷RRS(𝑡), is also displayed (left).

Tables

TABLE 1. Summary of acute skin toxicity data used for model assessment.

Reference	Dataset ID	Label	Skin toxicity	Evaluation time	Mice groups	Mice total	Delivery type	Nominal energy (MeV)	Source	Total dose (Gy)	TADR (Gy/s)	Exposure time a) (s)
This work	e1	CDR	Ulcers (>3 mm)	Max over 1–50 d	6	60	Pulsed broad electron beam (collimated)	5	Oriatron eRT6, PMB	18–31.9	0.11	165–292
This work	e2	UHDR	Ulcers (>3 mm)	Max over 1–50 d	7	70	Pulsed broad electron beam (collimated)	6	Oriatron eRT6, PMB	23.1–37	109–111	0.21–0.34
Sør22 19	p1	CDR	Grade 1.5–3.5	Max over 11–25 d	11	154	Scanned proton pencil beam	244	ProBeam, Varian	23.2–39.2	0.35–0.40	61–107
Sør22 19	p2	UHDR	Grade 1.5–3.5	Max over 11–25 d	14	147	Scanned proton pencil beam	250	ProBeam, Varian	31.4–52.8	65–92	0.35–0.73
Sør24 2	p3	Varying DR no repaint	Grade 1.5–3.5	Max over 11–25 d	7	109	Scanned proton pencil beam	250	ProBeam, Varian	39.3	0.7–80	0.5–60.5
Sør24 2	p4	Varying DR with repaint	Grade 1.5–3.5	Max over 11–25 d	8	86	Scanned proton pencil beam	250	ProBeam, Varian	39.3	0.37–80	0.44–107.8
Sør24 2	p5	Total dose split into 1–6 equal deliveries with 2 min pauses	Grade 1.5–3.5	Max over 11–25 d	7	125	Scanned proton pencil beam	244–250	ProBeam, Varian	39.3	0.065–60	0.68–607.3

TABLE 2. Evaluated functional forms of (a) ordinary differential equations (ODE) (Equation 2) and (b) Srel(Q).

(a)
ODE ID	Δ(Q)	R.(Q)
ESD1	−gΘ(Q(t))	rΘ(Q0−Q(t))
ESD2	−gQ(t)	r(Q0−Q(t))
ESD3	−gΘ(Q(t))	r(Q0−Q(t))

TABLE 3. Goodness of fit evaluations for all eight evaluated combinations of equilibrium state dynamics (ESD) and relative radiosensitivity (RRS) Srel.

Model ID a)	R2	MAE	MAPE	RMSE	BIC	ΔBIC	Optimized parameters
ESD1_RRS1	0.911	0.034	4.3%	0.043	−314.7	0.0	3
ESD1_RRS2	0.910	0.033	4.2%	0.043	−310.7	4.0	4
ESD1_RRS3	0.909	0.034	4.3%	0.044	−310.1	4.6	4
ESD2_RRS2	0.917	0.032	4.1%	0.042	−314.4	0.3	4
ESD2_RRS3	0.911	0.034	4.3%	0.043	−311.3	3.4	4
ESD3_RRS1	0.908	0.034	4.3%	0.044	−313.3	1.4	3
ESD3_RRS2	0.907	0.034	4.3%	0.044	−308.9	5.8	4
ESD3_RRS3	0.908	0.034	4.3%	0.044	−309.1	5.6	4

Full text

INTRODUCTION

1

Experimental evidence suggests that the normal tissue‐sparing effect of FLASH radiotherapy (FLASH‐RT) may be influenced by temporal dose delivery parameters, including average dose rate, exposure duration, beam pulse structure, and delivery pauses.1, 2, 3, 4, 5, 6 For preclinical research and clinical translation, it is therefore essential to understand these dependencies to optimize the design of ultra‐high dose rate (UHDR) devices and treatment approaches, as well as to guide future experimental and translational research.

The modeling of the FLASH effect is currently hampered by a limited, elusive understanding of its underlying mechanisms 5, 7, 8, 9 and, from a modeling perspective, by its scarce experimental characterization. Mechanistic, process‐oriented models of the FLASH effect risk being overly complex and tend to obscure the essential temporal role of the irradiations, reducing their practical applicability and requiring extensive parameter tuning. Oxygen enhancement ratio (OER)‐based modeling has been shown to capture temporal dose delivery dependencies,10, 11 but is difficult to reconcile with oxygen depletion measurements and modeling.12, 13, 14

In this study, we propose a phenomenological modeling approach based on perturbed equilibrium state dynamics (ESD), where a biological system's instantaneous radiosensitivity is modulated by a transient perturbation of a generic equilibrium state. This modeling approach focuses on capturing experimentally observed effect behaviors in terms of isoeffective dose ratios as a function of the temporal dose delivery profiles and integrates well into current dose delivery frameworks without invoking complex biochemical pathways. It should be understood as a generalization of existing dose‐modifying approaches, rather than as an entirely new modeling approach. Its mathematical structure is inspired by dynamic equilibrium models commonly used in biology and ecology, where a latent quantity is transiently saturated or depleted and recovers toward equilibrium.11, 15, 16 We validate the model against an extensive dataset of acute murine skin reactions, drawing on both published studies and new in vivo measurements. This data‐driven methodology provides a simplified yet flexible foundation for quantifying temporal dose delivery dependencies of FLASH sparing, which are crucial for the clinical implementation of the FLASH effect and its further preclinical investigation.

METHODS

2

In vivo experimentation

2.1

Experiments were conducted in compliance with the guidelines established by the Animal Ethics Committee of Vaud, Switzerland, and were approved by license nos. VD3934 and VD3962 (DGAV Canton Vaud). Mice used within this research were acclimatized for 1 week prior to experimentation. Mice were maintained and handled in pathogen‐free conditions in cages with a maximum of five mice per cage, under a controlled 12 h light/dark cycle, relative humidity of 55%, and controlled temperature of 21°C, with food and sterile water provided daily, all in accordance with the guidelines provided by the institution. The experiments were conducted on female BALB/c mice (Charles River; France) at 8 weeks of age. The dorsal skin of 8‐week‐female old BALB/c mice was placed under ketamine anesthesia (75 µg/g body weight), supplemented with medetomidine (0.5 µg/g body weight) using intraperitoneal (IP) injection and followed by atipamezole (1 µg/g body weight) IP injection post‐irradiation.

Electron irradiations were conducted at our institute using the 5–6 MeV pulsed electron beam from the prototype linear accelerator Oriatron eRT6.17 The dorsal skin was gently lifted and secured to the back of a 3‐cm‐thick carbon applicator, featuring a square aperture measuring 15 mm × 15 mm, covered by a 0.5 mm thick PMMA layer. Only the portion of skin extending up to 10 mm mark behind the aperture was exposed to either UHDR or CDR treatments, while the remainder of the animal was shielded by the applicator. Considering that the extended skin was folded, the total irradiated area was 20 mm × 15 mm. A summary of key setup and irradiation parameters for electron skin irradiations for the two radiation modalities is provided in Table S1, following recommendations.18 Absolute dose measurements were performed using EBT‐XD GafChromic films, positioned identically to the extended mouse skin. Given that the folded mouse skin remained thinner than 1 mm, the dose measured by the film served as an accurate proxy for the dose absorbed by the skin layers. Skin reactions were monitored for 50 days post‐irradiation unless euthanasia was required earlier due to mice reaching severity endpoint of >3 mm skin ulceration, occurring 16–30 days post CDR‐ and UHDR‐irradiation. Normal tissue complication probability (NTCP) is expressed as the fraction of animals in a group that developed skin ulcers exceeding 3 mm in any direction. This criterion was established as the humane endpoint in compliance with the approved animal license.

Mouse skin toxicity data for model assessment

2.2

A collection of acute mouse skin toxicity data was used to assess the models’ capability of describing FLASH normal tissue sparing for different temporal dose delivery patterns. The data collection evaluated skin toxicity scores of 721 mice divided into 60 experimental groups irradiated either with a pulsed electron beam (see previous section) or with the entrance dose plateau of a scanned proton beam of 244–250 MeV and included dose‐response curves for conventional dose rate (CDR) and ultra‐high dose rate (UHDR) irradiations for both electron and proton beams. The proton experiments have been described in detail elsewhere2, 19 and currently represent the largest consistently acquired data series evaluating acute murine skin reactions after CDR and UHDR irradiation. A summary of the datasets used for model assessment is given in Table 1. They comprise diverse temporal dose delivery structures spanning a dose range of 18–53 Gy, time‐averaged dose rates (TADR) between 0.11 and 111 Gy/s, and exposure times of 0.21 s to 10 min, including dose rate scans and split‐beam deliveries.

Derivation of experimental isoeffective dose ratios

2.3

The FLASH‐modifying factor (FMF) is the ratio of doses required to be delivered at a reference temporal dose delivery pattern (Dref, i.e. a CDR irradiation) and at a test temporal dose delivery pattern (Dtest, i.e. any temporal irradiation pattern) to achieve an isoeffect for a given biological endpoint1

$$[eqn]$$

This follows a definition analogous to the one of relative biological effectiveness for different radiation qualities.20 Experimental data provided group‐wise NTCP. Dataset e1 (TADR of 0.11 Gy/s) and dataset p1 (TADR of 0.35‐0.40 Gy/s), see Table 1, were used to establish reference dose‐response curves for electron and proton irradiations, respectively. For this purpose, a logistic NTCP curve was fitted to the experimental single mouse data for each evaluated skin toxicity grade. The experimental isoeffective dose ratios (FMF_exp_) were then calculated, as illustrated in Figure 1a, using the experimental data points to determine Dtest and obtaining the corresponding Dref for an isoeffective NTCP from the logistic fit to the respective reference dataset. Since isoeffective dose ratios are undefined for cases of complete absence or saturation of effect, FMF_exp_ was only determined for data points with an NTCP between 1% and 99%. For proton data, multiple grades of skin reactions were evaluated and resulted often in multiple FMF_exp_ values pertaining to different toxicity grades for the same mouse group. These FMF_exp_ values were averaged for model evaluation.

(a) Schematic representation: experimental dose ratios for isoeffective normal tissue complication probability (NTCP), referred to as experimental FLASH modifying factors (FMFexp), were computed utilizing the experimental data point for the test temporal dose delivery pattern (𝐷test) and the reference dose (𝐷ref) obtained from a logistic fit to the dose‐response curve for the respective reference temporal dose delivery pattern. (b) NTCP for skin ulcers (>3 mm) after conventional dose rate (CDR) and ultra‐high dose rate (UHDR) electron irradiation (i.e., datasets e1 and e2). A logistic NTCP curve (line) was fitted to the experimental single mouse data for CDR data.

Equilibrium state dynamics‐based modeling and simulation of experimental isoeffective dose ratios

2.4

A recent meta‐analysis examining in vivo data on NT sparing by FLASH as a function of dose found that available experimental data can be described by a piecewise linear function of two parts when expressed as FMF‐weighted dose versus dose, also referred to as sudden effect transition (SET) function.1 This implies for UHDR irradiations a decreasing sparing effect in terms of FMF for lower doses and a saturation for very high doses. These features have been subsequently confirmed for data acquired over a large dose range,21, 22, 23 see Figure S1. Such a behavior is compatible with an instantaneous relative radiosensitivity Srel that is initially one (compared to a low dose rate irradiation) and that may reach a minimum value Smin during UHDR irradiation. Furthermore, experiments have shown that the magnitude of normal tissue sparing by FLASH is typically modulated by irradiations on time scales of ≲1 s (see Figure 2a–c).4, 24, 25 The observed key dependencies and relevant time scales for the modulation of the FLASH normal tissue sparing magnitude suggest that the radiosensitivity of the biological system behaves as if it is transiently pushed away from its conventional state by swift dose delivery. Such behavior can be effectively modeled using the framework of perturbed equilibrium state dynamics (ESD). For an effect description based on ESD, we introduce a generic latent quantity Q(t), which represents a pool of chemical or biological ‘resources’ that modulates the instantaneous relative radiosensitivity Srel(Q) of the biological system during irradiation relative to a low dose rate irradiation. Srel(Q) is a dimensionless scaling factor representing the relative weight of incremental dose contributions under perturbed equilibrium versus steady state Q0. Q(t) may be pushed away from its equilibrium state Q0 by irradiation with an instantaneous dose rate D.(t) inducing a change Δ(Q) and recovers toward equilibrium at rate R.(Q), reflecting ‘replenishment’ processes. Such a behavior is illustrated in Figure 2d. Not considering any spatial components, the dynamics of such a system can be described by a 1st‐order ordinary differential equation (ODE)

$$[eqn]$$

with the boundary condition Q(t=0)=Q0 and a system state between equilibrium and complete saturation/depletion, i.e. Q0≥Q(t)≥0. While not at equilibrium, Q(t) may alter the instantaneous relative radiosensitivity Srel(Q) of the biological system between 1 and a minimum Smin, with Srel(t=0)=1 and Smin>0. The dose D delivered during the irradiation session of time T is given by

$$[eqn]$$

(a–c) Changes in normal tissue sparing for a fixed single fraction dose ((a) 10 Gy, 26 (b) 10 Gy, 25 (c) 11.2 Gy (9–10 weeks)/12.5 Gy (30–31 weeks) 4 ) as a function of exposure time for different in vivo assays. Data points were digitized from the original publications. 4 , 25 , 26 Dashed grey lines (CTRL) correspond to reference level for control groups. (d) Schematic graph illustrating perturbed equilibrium state dynamics of a generic quantity 𝑄(𝑡) that may be pushed away from its equilibrium state 𝑄0 by an irradiation at instantaneous dose rate D˙(t) inducing a change Δ(𝑄). 𝑄(𝑡) returns toward equilibrium 𝑄0 at a given recovery rate R˙(Q). (e) Evaluated dependencies for the instantaneous relative radiosensitivity Srel(𝑄) include a Heaviside function, a power function, and a Michaelis–Menten (M–M) type behavior. For details, see text and Table 2.

Assuming that radiobiological damage is altered by the status of Q(t) during dose delivery, the relative radiosensitivity (RRS)‐weighted dose DRRS provides a measure for the effective radiobiological damage to the system compared to an irradiation at equilibrium

$$[eqn]$$ $$[eqn]$$

In this study, we evaluated multiple ESD models with different dependencies for Δ(Q), R.(Q), and Srel(Q) based on their ability to describe experimental FMF_exp_ values. For this purpose, combinations of simple behaviors for Δ(Q), R.(Q), and Srel(Q) were evaluated as detailed in Table 2. Those behaviors are compatible with the previously discussed observations of FLASH normal tissue sparing and boundary conditions. For Srel(Q), this included a step‐like behavior (RRS1), a power function (RRS2) that accommodates proportional as well as under‐ and over‐proportional behaviors, and a Michaelis–Menten‐type function (RRS3) that results in an intermediate behavior between RRS1 and a proportional behavior (i.e., under‐proportional), see Figure 2e. In the following, Q0 was fixed to a generic value of 1. Because ODE ESD2 cannot reach Q=0, the combination ESD2_RRS1 was not evaluated, as it would leave RRS unchanged. Consequently, we evaluated in total eight combinations of ESD and RRS (see Tables 2 and 3) that had three or four free parameters. The tested functional forms were deliberately restricted to the simplest plausible dependencies (step‐like, linear, under/over‐proportional saturation/depletion responses), chosen to represent principal experimental trends while minimizing unnecessary complexity and reducing the risk of overfitting.

Isoeffective dose ratios predicted by ESD models were calculated as

$$[eqn]$$

as follows from Equation 1 and isoeffectiveness (i.e., RRStestDtest=RRSrefDref). RRS‐weighted dose DRRS, RRS and FMFpred (Equations (2), (3), (4), (5), (6)) were computed by numerically solving the corresponding ODE and integrals for the respective experimental temporal dose delivery profiles D.(t). For this purpose, experimental irradiation profiles D.(t) were calculated at the dosimetric reference points. Broad beam electron irradiations were spatiotemporally constant18 and temporal structure of UHDR dose delivery was derived from the temporal pulse beam structure parameters and the total dose while assuming rectangular pulses, see Table S1. For scanned proton irradiations, D.(t) at the field center was calculated based on simulated intensities and delivery times of each spot, its distance to the field center and the spot scan order, while taking into account the non‐Gaussian spot dose profile previously measured in vivo during the experiments.27 D.(t) of reference dose rate irradiations were approximated by a continuous irradiation at constant dose rates of 0.11 and 0.37 Gy/s for datasets e1 and p1 (Table 1), respectively.

Model fitting, evaluation, and ranking

2.5

The eight ESD model variants (see Tables 2 and 3) were evaluated based on their capability to fit the collection of experimental FMFexp. The mean squared error objective function

$$[eqn]$$

was optimized with the Nelder–Mead simplex algorithm28 (scipy.optimize.minimize v1.6.2). Model ranking was performed using the Bayesian Information Criterion (BIC) with ΔBIC being the BIC difference between a given model and the best‐performing model.29, 30 ΔBIC values of [0–2], [2–6], and [6–10] are referred to hereafter as weak, positive, and strong support for model rejection, respectively.30 The coefficient of determination (R2), the mean absolute error (MAE), the mean absolute percentage error (MAPE), and the root mean squared error (RMSE) were evaluated as additional goodness‐of‐fit indicators. The Shapiro–Wilk test was used to assess normal distribution of the residuals with a significance level of 0.05.

RESULTS

3

Dose‐response curves for electron irradiations (e1 and e2) are shown in Figure 1b. A total of 123 experimental FMF_exp_ values were derived from the experimental data collection (Table 1) and resulted in 52 FMF_exp_ values after averaging FMF_exp_ values obtained from the same mouse group for different toxicity levels (Figure 3). Model variant ESD1_RRS1 was ranked best based on BIC, see Table 3, with g=0.22Gy−1,r=0.29s−1, and Smin=0.65. ΔBIC indicated weak support for model rejection for ESD2_RRS2 and ESD3_RRS1, whereas positive support was indicated for the remaining five model variants. All eight candidate ESD models were able to reproduce experimental FMF_exp_ values with good accuracy (R2 within 0.907–0.917, MAE within 0.032–0.034, MAPE within 4.1–4.3%, and RMSE within 0.042–0.044) using three or four free parameters, see Table 3, and resulted in similar values for Smin (range: 0.55–0.65). Simulated FMF_pred_ and experimental FMF_exp_ values and corresponding residuals for the best ranked model are shown in Figure 3a,b as a function of dose and TADR, respectively. The Shapiro–Wilk test statistic for residuals was 0.98 with p = 0.4, indicating no significant deviation from normality. Variation of residuals were similar for the reference data sub‐group (e1, p1: MAE=0.030, RMSE=0.046) and the non‐reference data sub‐group (e2, p2, p3, p4, p5: MAE=0.035, RMSE=0.042), suggesting that residuals are driven by uncertainties in the experimental data, rather than by inabilities of the model to capture key data trends. Figure 4 displays examples of temporal irradiation profiles together with the simulated dynamics for Q(t), Srel(t), and the corresponding accumulated RRS‐weighted dose DRRS(t) for the selected model variant ESD1_RRS1.

Simulated and experimental isoeffective dose ratios (FLASH‐modifying factors, FMF) and corresponding residuals for model variant ESD1_RRS1 as a function of dose (a) and as a function of time‐averaged dose rate (TADR) (b). For the residuals, dashed horizontal lines indicate ± 0.05.

Four representative temporal irradiation profiles from the data collection are shown (left) alongside the simulated dynamics for 𝑄(𝑡) and 𝑆rel(𝑡) for the ESD1_RRS1 model variant (right). The corresponding accumulated RRS‐weighted dose, 𝐷RRS(𝑡), is also displayed (left).

DISCUSSION

4

Equilibrium state dynamics‐based modeling offers a parsimonious, data‐driven, and mechanistically agnostic framework that generalizes established dose‐modifying concepts (e.g., OER‐weighted dose) to describe FLASH normal tissue sparing quantitatively, emphasizing the role of transient equilibrium perturbations in modulating radiosensitivity. Unlike our previous work focusing only on single‐fraction dose dependence,1 the present study extends to capture key experimental dependencies on temporal structure using only three or four free parameters in a simple first‐order ODE. The limited parameterization inherent to ESD combined with its simplicity addresses common challenges in radiobiological modeling, such as overfitting of sparse data, and provides a framework for parameter extraction from experimental studies. ESD does not assume a priori that the temporal dose delivery structure on a certain time scale determines FLASH sparing. Instead, the modeling framework allows such dependencies to emerge naturally from the fitted parameters if supported by data. ESD models were shown in this work to be able to describe consistently an extensive collection of experimental isoeffect dose ratios (FMF_exp_), spanning diverse temporal dose delivery patterns, as well as pulsed electron and scanned proton beams, covering thereby a large domain of relevance for the clinical transfer of FLASH‐RT. The best‐performing model ESD1_RRS1 could reproduce experimental FMF with an MAE of 0.034 using three free parameters.

If the extent of FLASH normal tissue sparing is modulated according to ESD, this holds important implications for its clinical transfer. It underscores the critical importance of delivering hypo‐ or even extremely hypofractionated FLASH‐RT in order to maximize normal tissue sparing by FLASH, assuming equitoxicity for tumor tissues.31 However, in many clinical scenarios, such fractionation schemes may compromise the sparing of late‐reacting tissues.32, 33, 34 Furthermore, it suggests that any FLASH‐RT fraction should be delivered on time scales well below recovery times of critical normal tissues to maximize their sparing. A linear recovery of r=0.29s−1 was obtained for model variant ESD1_RRS1 for acute murine skin reactions, though recovery times likely vary by tissue and endpoint.

Out of the eight evaluated ESD‐based models, all reproduced the experimental isoeffect dose ratios (FMF_exp_) with good accuracy using three or four free parameters (Table 3). Models with four parameters achieved similar RMSE and MAE compared to three‐parameter variants but generally had higher BIC values, indicating mild overfitting—apart from ESD2_RRS2, which performed comparably. The fact that several variants yielded similarly low BIC values highlights that the dataset lacks sufficient discriminative power to decisively favor one functional form, reflecting both the limited precision and domain of the data, as well as similarities between the tested functional dependencies (e.g., Michaelis–Menten vs. factor‐type forms). Within this context, ESD1_RRS1 can be regarded as a straightforward and practical three‐parameter model for capturing the primary temporal dose‐delivery dependencies of acute murine FLASH skin sparing. However, a constant change Δ(Q) and recovery R.(Q) combined with a non‐continuous, step‐like Srel(Q) may be overly simplistic, and future investigations should refine these behaviors to provide a more accurate representation and investigate predictive performances for unseen data.

ESD1_RRS1 (Equations (3), (4), (5), (6) and Table 2) reduces for irradiations at a constant dose rate D.c with gD.c>r to

$$[eqn]$$

and, similarly, for irradiations for which exposure times T are short so that recovery r is negligible, it follows

$$[eqn]$$

Pulsed irradiations with many pulses can be approximated by the constant dose rate case and representations similar to Equations 8 and 9 can also be obtained for other ESD model variants for high doses and in some cases also for lower doses. Assuming a RRSref (Equations (4), (5), (6)) of 1, Equations 8 and 9 are both equivalent to the SET function, which was shown to describe the single fraction dose‐dependency of FLASH normal tissue sparing for multiple experimental datasets in a recent meta‐analysis.1 For RRSref=1, the following association with parameters of the SET function can be made: FMFmin=Smin and DT=D.cQ0(gD.c−r). For a negligible r, this reduces to DT=Q0g. A FMFmin of 0.67 and 0.51 was found for all pooled data and for skin data, respectively.1 This is in reasonable agreement with Smin=0.65, as determined in this study for model ESD1_RRS1. In contrast, the DT values of 9.6 Gy (all pooled data) and 16.6 Gy (skin data) found by the meta‐analysis are considerably higher than a DT of 4.3 Gy obtained by ESD1_RRS1 for a negligible r and RRSref=1 (Equation 9). This might be attributed to limitations when extrapolating the ESD model to lower dose regions for which it was not fitted.

In this study, ESD was fitted and evaluated exclusively on acute murine skin reactions obtained with pulsed electron beams and scanned proton beams. Although the framework is general and broadly adaptable, its parameters are context‐specific and not assumed to transfer across tissues, endpoints, or species. Re‐fitting will be required for new applications. Consequently, claims of generalizability are limited, and any potential clinical translation remains speculative pending validation with independent datasets. ESD models are bound to the data domain for which they have been validated, and as a phenomenological modeling approach, they do not rule out that additional, complementary processes may further influence normal tissue sparing by FLASH. For instance, in contrast to the electron and proton skin toxicity data utilized in this work, recent studies identified considerable differences in intestinal responses to UHDR proton and electron irradiation35 and suggest a strong biological component in the mechanisms of the FLASH sparing effect for different systems and endpoints.21 The absence of a sparing effect for tumoral tissues is also not directly deducible from ESD models. However, while these aspects and differences cannot be explained by the ESD‐based modeling approach presented here, they can in principle be incorporated by endpoint‐ and modality‐specific model parameters. Furthermore, we have omitted explicit spatial recovery dynamics (e.g., diffusion) due to a lack of suitable data, even though such spatial factors are likely to affect FLASH sparing when irradiation and system recovery times occur on comparable time scales (∼1s).

ESD‐based models obtained in this study are compatible with hypotheses that sub‐second physical and chemical processes may contribute to the mechanism of the FLASH effect, such as radical recombination, reactive oxygen species scavenging, lipid peroxidation, and Fenton chemistry, as well as oxygen depletion and recovery.7, 8, 9, 36 Taken together, the ESD‐based data‐driven description of FLASH sparing aligns with a saturation or depletion phenomenon that modifies the instantaneous radiosensitivity and that has for acute murine skin reactions a characteristic recovery rate on a sub‐second time scale. In fact, as a phenomenological modeling approach, ESD may help narrow down the search for eligible mechanisms of action for the FLASH effect and may serve as a forerunner and benchmark for more mechanistic/process‐oriented descriptions. Although ESD models reproduced the acute murine skin dataset with high fidelity, the inability to decisively favor one variant over another limits the extent to which biological inferences can be drawn. This underscores that current datasets are better suited for phenomenological description than for mechanistic discrimination. Future experiments should therefore be designed to probe parameter‐sensitive regimes, for example via split‐dose studies with controlled gap times, different temporal beam structure parameters, such as a varying pulse repetition frequency, or intermediate dose ranges where predictions diverge. Incorporating high‐precision readouts and cross‐tissue validations may further increase discriminative power.

It should be underscored that ESD‐based modeling is mathematically similar to a description of FLASH normal tissue sparing based on the conventional radiolytic oxygen depletion hypothesis10, 11 that is difficult to reconcile with oxygen depletion measurements and modeling,12, 13, 14 as initially mentioned. In fact, the model variant ESD2_RRS3 in the current study is mathematically identical to the OER‐weighted dose model, which has already previously been shown to describe the presented proton dataset successfully,10 however, with a different naming of the involved parameters (O(t), g, λ, O env, OER/OER_env_, 1/OER_env_, K/O env, and m in Ref. 10 correspond to Q(t), g, r, Q 0, S rel, S min, K, and 1/S min(1+K)‐K in ESD2_RRS3). Despite yielding similar residuals, ESD2_RRS3 was penalized by the BIC because of its larger number of free parameters and was therefore less favored than ESD1_RRS1. The novelty of the current study lies in extending the experimental dataset with pulsed electron beam measurements and demonstrating that a broader variety of ESD‐based models can reproduce the experimental data consistently. Furthermore, the ESD‐based description and corresponding model parameters as formulated in this work are not attached to specific mechanistic processes, such as oxygen depletion.

The functional forms used to link the latent quantity Q to radiosensitivity were deliberately selected to represent the simplest plausible dependencies that are consistent with experimental observations: a step function to mimic a threshold‐like response and power‐law and Michaelis–Menten‐type forms to capture proportional or supra‐/sub‐proportional effects and to represent saturation/depletion behavior. These forms are widely used across biological modeling because they balance descriptive power with parsimony, and they allow the data to identify whether responses behave more like thresholds, proportional shifts, or saturating processes. By starting with these minimal formulations, the ESD framework avoids introducing unnecessary mechanistic assumptions while still providing flexibility to capture the principal experimental trends. In this sense, the purpose of the model is not to prescribe a single mechanistic explanation but to offer a parsimonious, phenomenological abstraction that can unify different temporal‐dose dependencies within one consistent mathematical framework. Such abstractions may support hypothesis generation, enable quantitative comparison across datasets, and help design experiments that can further test where different functional forms diverge.

CONCLUSIONS

5

ESD‐based modeling can capture key temporal dose‐delivery dependencies of FLASH acute skin sparing using a parsimonious, data‐driven approach. Further investigations using independent datasets are needed to assess the predictive performances of ESD and whether modeled behaviors extend to other tissues, endpoints, modalities, and delivery patterns. If these behaviors generalize, ESD‐based modeling could serve as a practical framework for evaluating and optimizing UHDR irradiations in preclinical and translational settings.

CONFLICT OF INTEREST STATEMENT

The authors declare no conflict of interest.

HUMAN/ANIMAL REGULATORY APPROVAL

For new murine electron data, published first within this study: the Animal Ethics Committee of Vaud, Switzerland (Approval No. VD3934 and VD3962). For previously published research data that are used by this study, the respective publications apply.

Supporting information

Supporting Information