Impact of Oxygen on DNA Damage Distribution in 3D Genome and Its Correlation to Oxygen Enhancement Ratio under High LET Irradiation

Basic Information

• Paper ID: 2503.21837
• Title: Impact of Oxygen on DNA Damage Distribution in 3D Genome and Its Correlation to Oxygen Enhancement Ratio under High LET Irradiation
• Authors: Ankang Hu, Wanyi Zhou, Xiyu Luo, Rui Qiu, Junli Li
• Affiliated Institutions: Department of Engineering Physics, Tsinghua University; Institute of Fluid Physics, China Academy of Engineering Physics
• Classification: physics.bio-ph physics.med-ph
• Paper Link: https://arxiv.org/abs/2503.21837

Abstract

The variation of oxygen enhancement ratio (OER) under different linear energy transfer (LET) values currently lacks comprehensive mechanistic explanation and theoretical models. Previous research by the authors revealed significant correlations between double-strand break (DSB) distribution within 3D genomes and radiation-induced cell death, providing valuable insights for understanding oxygen effects. This study proposes a model representing oxygen response as the probability of inducing DNA strand breaks and integrates it into track structure Monte Carlo simulations to investigate oxygen's impact on DSB distribution within 3D genomes. Results demonstrate that the ratio of occurrence rates for clustered DSBs within individual topologically associating domains (TADs) (Case 2) and DSBs in frequently interacting TADs (Case 3) under both aerobic and hypoxic conditions is consistent with the trend of cellular survival OER across different LET values.

Research Background and Motivation

Core Problems

1. Missing OER Mechanism: The variation of oxygen enhancement ratio (OER) with linear energy transfer (LET) lacks comprehensive mechanistic explanation, particularly the mechanism of OER reduction under high LET irradiation remains unclear
2. Limitations of Traditional Models: Existing models primarily rely on empirical parameter fitting, lacking mechanistic explanations from the perspective of DNA damage distribution
3. Importance of 3D Genome Structure: Traditional research focuses mainly on total DSB numbers, overlooking the importance of DSB spatial distribution within 3D genome architecture

Research Significance

• Radiotherapy Optimization: Understanding OER mechanisms is crucial for optimizing radiotherapy strategies, particularly in hypoxic tumor treatment
• Radiosensitizer Development: Provides theoretical foundation for developing radiosensitizers based on oxygen effect principles
• FLASH Effect Mechanism: Supports understanding of oxygen effects in the current research hotspot of FLASH radiotherapy

Limitations of Existing Methods

1. Empirical Models: Most OER models rely on empirical parameters fitted directly to cellular survival data
2. Simplified Chemical Reactions: Although some models attempt to incorporate DNA radical-oxygen reaction corrections, these adjustments remain empirical
3. Scale Mismatch: DSB yield curves cannot match relative biological effectiveness (RBE) curves, suggesting the need for analysis at larger genomic structure scales

Core Contributions

1. Innovative Mechanistic Model: Proposes an OER mechanistic explanation model based on 3D genome DSB distribution, modeling oxygen response as the probability of inducing DNA strand breaks
2. Track Structure Simulation Integration: Integrates oxygen effects into track structure Monte Carlo simulations to study oxygen's impact on DSB distribution within 3D genomes
3. TAD-Scale Damage Classification: Establishes a classification system for three DSB distribution patterns based on topologically associating domains (TADs)
4. Experimental Validation: Calculated OER results show good correspondence with experimental data, validating model effectiveness
5. High LET Mechanism Explanation: Provides potential mechanistic explanation for OER reduction under high LET irradiation

Methodology Details

Task Definition

Input: Radiation particles of different LET (electrons, protons, He-3 ions, C-12 ions) and their energies Output: DSB distribution patterns within 3D genome, occurrence rates of various damage scenarios, OER values Constraints: Based on IMR90 cell nucleus model, considering TAD structure and interactions

Model Architecture

1. 3D Genome-Based Cell Survival Model

DSB distribution is classified into three cases based on TAD structure:

• Case 1: Isolated DSBs - only one DSB within a single TAD, with no DSBs in its frequently interacting TADs
• Case 2: Clustered DSBs within a single TAD - two or more DSBs within a single TAD
• Case 3: DSBs in frequently interacting TADs - multiple frequently interacting TADs each containing one or more DSBs

Survival fraction calculation formula:

-ln S = p₁n₁(D) + p₂n₂(D) + p₃n₃(D)

where p₃ > p₂ > p₁, indicating that Cases 3 and 2 have higher cell lethality probabilities.

2. Oxygen Effect Mechanistic Modeling

Oxygen's impact on radiation damage is primarily achieved through reactions with DNA radicals, involving three competing processes:

1. Oxygen Fixation: Reaction with oxygen to form peroxy radicals
2. Natural Fixation: Conversion to lesions requiring biological repair
3. Chemical Repair: Repair by antioxidants such as glutathione

Probability of DNA radical conversion to lesions:

p_lesion([O₂]) = (k₁·[O₂] + k₂)/(k₁·[O₂] + k₂ + k₃·[antioxidants])

3. Track Structure Monte Carlo Simulation

• Physical Simulation: Energy deposition events simulated using Geant4-DNA
• Cell Nucleus Model: IMR90 cell nucleus model based on Hi-C data, containing 15,282 TADs
• Damage Analysis: DSB distribution and clustering analyzed using DBSCAN algorithm

Technical Innovations

1. Multi-Scale Integration: Combines molecular-scale oxygen chemistry reactions with cell nucleus-scale 3D genome structure
2. Probabilistic Treatment: Simplifies complex oxygen chemistry reactions into probability p_DNA of energy deposition events inducing strand breaks
3. TAD-Based Analysis: Classifies damage based on biologically meaningful TAD structure rather than simple spatial distance
4. Statistical Methods: Employs Poisson distribution assumptions to handle statistical variability in simulation results

Experimental Setup

Simulation Parameters

• Cell Nucleus Model: IMR90 cells, radius 5.0 μm, volume fraction 0.15
• TAD Parameters: 15,282 TADs, average radius 89.12 nm, average DNA content 397 kbp
• Radiation Types: 0.9 MeV electrons (low LET), 0.4-66.46 MeV protons, 0.8-23.0 MeV/u He-3 ions, 2.98-209.46 MeV/u C-12 ions
• Oxygen Effect Parameters: p_DNA range 0.04-0.22, aerobic condition p_DNA=0.14, hypoxic condition p_DNA=0.09

Evaluation Metrics

• OER_α: Oxygen enhancement ratio based on α parameter
• OER₁₀: Oxygen enhancement ratio at 10% survival fraction
• Damage Yield Ratio: Ratio of SSB, DSB, and various cases under aerobic/hypoxic conditions

Simulation Workflow

1. Physical simulation using Geant4 to generate energy deposition events
2. Determine potential lesions based on energy threshold and p_DNA probability
3. Randomly assign to two DNA strands
4. Define DSBs: two or more strand breaks on opposite strands within 3.2 nm
5. Analyze DSB distribution using DBSCAN algorithm
6. Calculate occurrence rates of three cases and compute survival fractions

Experimental Results

Main Results

1. Impact of p_DNA on Damage Distribution

• Low LET: SSB yield shows linear relationship with p_DNA; DSB and complex damage show nonlinear relationships
• High LET: All damage types exhibit more complex nonlinear relationships
• Classical Alper Formula: Valid only for simple strand breaks, deviates from predictions for complex damage scenarios

2. Oxygen Impact on Damage Distribution at Different LET

With increasing LET:

• Yield ratio of SSB and DSB decreases
• Occurrence rate ratio of Cases 2 and 3 decreases significantly
• Trend of Cases 2 and 3 ratio curve highly consistent with cellular survival OER trend

3. OER Calculation Results

• OER_α: Remains relatively constant in LET range of tens of keV/μm
• OER₁₀: Shows slight increase in same LET range, eventually approaching OER_α
• Experimental Comparison: Calculated OER curves show good correspondence with experimental data

Key Findings

1. Mechanistic Explanation: High LET radiation produces dense ionization, causing excessive lesions; even with partial chemical repair, remaining lesions are sufficient to induce Cases 2 and 3
2. Scale Importance: DSB distribution within 3D genome better predicts cell death than total DSB numbers
3. TAD Correlation: TAD-based damage classification highly correlates with OER variation trends, validating the importance of 3D genome structure

Related Work

Traditional Oxygen Effect Research

1. Oxygen Fixation Hypothesis: Oxygen reacts with DNA radicals to form irreparable peroxy radicals
2. Track Structure Simulation: Incorporates oxygen effects into Monte Carlo simulations, but primarily focuses on total DSB numbers
3. High LET Special Mechanisms: Hypothesis of multiple ionization producing oxygen, but with limited effects

3D Genome Damage Research

• Cell nucleus modeling based on Hi-C data
• Impact of TAD structure on DNA repair
• Relationship between chromatin 3D structure and radiation sensitivity

OER Modeling Methods

• Empirical formula fitting
• Chemical reaction kinetics models
• Microdosimetry approaches

Conclusions and Discussion

Main Conclusions

1. New Mechanistic Explanation: Provides potential mechanistic explanation for OER variation with LET based on 3D genome DSB distribution
2. TAD-Scale Importance: Confirms the importance of TAD-scale DNA damage distribution in predicting cell death
3. High LET Characteristics: High LET radiation reduces the relative importance of oxygen effects by producing excessive lesions

Limitations

1. Model Simplification: Simplifies complex oxygen chemistry reactions into probability models, potentially overlooking certain details
2. Single Cell Type: Uses only IMR90 cell model; applicability to other cell types requires verification
3. Scattered Experimental Data: High LET experimental OER data are highly dispersed, limiting quantitative comparison
4. Chemical Reaction Details: Overlooks secondary factors such as oxygen production from multiple ionization under high LET

Future Directions

1. Multi-Cell Type Verification: Extend to different cell lines to verify model universality
2. Refined Chemical Reactions: Incorporate more detailed chemical reaction mechanisms
3. FLASH Effect Application: Apply to understanding oxygen effects in FLASH radiotherapy
4. Clinical Translation: Explore potential applications in clinical radiotherapy planning

In-Depth Evaluation

Strengths

1. Strong Innovation: First to combine 3D genome structure with OER mechanism, providing new research perspective
2. Rigorous Methodology: Multi-disciplinary approach combining track structure simulation, statistical analysis, and biological modeling
3. Experimental Validation: Good correspondence between calculated and experimental data enhances model credibility
4. Mechanistic Explanation: Provides reasonable explanation for long-unexplained OER changes under high LET

Weaknesses

1. Model Complexity: Involves multiple scales and disciplines, increasing complexity of model verification and application
2. Parameter Dependence: Model performance depends on multiple empirical parameters requiring more experimental data for calibration
3. Computational Cost: Monte Carlo simulation and 3D genome analysis have high computational costs
4. Biological Assumptions: TAD-based damage classification requires more biological experimental verification

Impact

1. Academic Contribution: Provides new theoretical framework and research direction for radiation biology
2. Application Potential: Can optimize radiotherapy planning, particularly for heavy ion therapy and hypoxic tumor treatment
3. Interdisciplinary Value: Promotes cross-disciplinary research integration of physics, biology, and medicine
4. Methodological Value: Provides modeling paradigm for integrating multi-scale biophysical processes

Applicable Scenarios

1. Radiotherapy Planning: Optimization considering tumor oxygenation status
2. Radiation Protection: Assessment of biological effects of different radiation types
3. Drug Development: Mechanistic research and efficacy prediction of radiosensitizers
4. Basic Research: Deeper understanding of DNA damage repair mechanisms

References

The paper cites 38 important references covering:

• Classical oxygen effect theory (Alper formula, oxygen fixation hypothesis)
• Track structure Monte Carlo simulation methods
• 3D genome modeling techniques
• Experimental OER data
• DNA damage repair mechanism research

Overall Assessment: This is an innovative interdisciplinary research paper that successfully introduces 3D genome structure into OER mechanism research, providing reasonable mechanistic explanation for the long-puzzling OER changes under high LET irradiation. Despite some model simplifications and verification requirements, its innovative methodology and important findings open new research directions for radiation biology and radiotherapy fields.