A theoretical stochastic control framework for adapting radiotherapy to hypoxia

Fatemeh Saberian, Archis Ghate, Minsun Kim

Research output: Contribution to journalArticlepeer-review

13 Scopus citations

Abstract

Hypoxia, that is, insufficient oxygen partial pressure, is a known cause of reduced radiosensitivity in solid tumors, and especially in head-and-neck tumors. It is thus believed to adversely affect the outcome of fractionated radiotherapy. Oxygen partial pressure varies spatially and temporally over the treatment course and exhibits inter-patient and intra-tumor variation. Emerging advances in non-invasive functional imaging offer the future possibility of adapting radiotherapy plans to this uncertain spatiotemporal evolution of hypoxia over the treatment course. We study the potential benefits of such adaptive planning via a theoretical stochastic control framework using computer-simulated evolution of hypoxia on computer-generated test cases in head-and-neck cancer. The exact solution of the resulting control problem is computationally intractable. We develop an approximation algorithm, called certainty equivalent control, that calls for the solution of a sequence of convex programs over the treatment course; dose-volume constraints are handled using a simple constraint generation method. These convex programs are solved using an interior point algorithm with a logarithmic barrier via Newton's method and backtracking line search. Convexity of various formulations in this paper is guaranteed by a sufficient condition on radiobiological tumor-response parameters. This condition is expected to hold for head-and-neck tumors and for other similarly responding tumors where the linear dose-response parameter is larger than the quadratic dose-response parameter. We perform numerical experiments on four test cases by using a first-order vector autoregressive process with exponential and rational-quadratic covariance functions from the spatiotemporal statistics literature to simulate the evolution of hypoxia. Our results suggest that dynamic planning could lead to a considerable improvement in the number of tumor cells remaining at the end of the treatment course. Through these simulations, we also gain insights into when and why dynamic planning is likely to yield the largest benefits.

Original languageEnglish (US)
Pages (from-to)7136-7161
Number of pages26
JournalPhysics in Medicine and Biology
Volume61
Issue number19
DOIs
StatePublished - Sep 19 2016
Externally publishedYes

Bibliographical note

Publisher Copyright:
© 2016 Institute of Physics and Engineering in Medicine.

Keywords

  • convex optimization
  • functional imaging
  • intensity modulated radiation therapy

Fingerprint

Dive into the research topics of 'A theoretical stochastic control framework for adapting radiotherapy to hypoxia'. Together they form a unique fingerprint.

Cite this