PubMed Health⌕ Search

PubMed · 11260661

A model for optimizing normal tissue complication probability in the spinal cord using a generalized incomplete repair scheme.

Abstract

The purpose of this study was to determine the treatment protocol, in terms of dose fractions and interfraction intervals, which minimizes normal tissue complication probability in the spinal cord for a given total treatment dose and treatment time. We generalize the concept of incomplete repair in the linear-quadratic model, allowing for arbitrary dose fractions and interfraction intervals. This is incorporated into a previously presented model of normal tissue complication probability for the spinal cord. Equations are derived for both mono-exponential and bi-exponential repair schemes, regarding each dose fraction and interfraction interval as an independent parameter, subject to the constraints of fixed total treatment dose and treatment time. When the interfraction intervals are fixed and equal, an exact analytical solution is found. The general problem is nonlinear and is solved numerically using simulated annealing. For constant interfraction intervals and varying dose fractions, we find that optimal normal tissue complication probability is obtained by two large and equal doses at the start and conclusion of the treatment, with the rest of the doses equal to one another and smaller than the two dose spikes. A similar result is obtained for bi-exponential repair. For the general case where the interfraction intervals are discrete and also vary, the pattern of two large dose spikes is maintained, while the interfraction intervals oscillate between the smallest two values. As the minimum interfraction interval is reduced, the normal tissue complication probability decreases, indicating that the global minimum is achieved in the continuum limit, where the dose delivered by the "middle" fractions is given continuously at a low dose rate. Furthermore, for bi-exponential repair, it is seen that as the slow component of repair becomes increasingly dominant as the magnitude of the dose spikes decreases. Continuous low-dose-rate irradiation with dose spikes at the start and end of treatment yields the lowest normal tissue complication probability in the spinal cord, given a fixed total dose and total treatment time, for both mono-exponential and bi-exponential repair. The magnitudes of the dose spikes can be calculated analytically, and are in close agreement with the numerical results.

Explore related subjects

Keep this discovery

Explore connections, maps & timelines

BibTeXRIS

D Levin-Plotnik, R J Hamilton, A Niemierko, S Akselrod. 2001. A model for optimizing normal tissue complication probability in the spinal cord using a generalized incomplete repair scheme.. https://doi.org/10.1667/0033-7587(2001)155%5B0593%3Aamfont%5D2.0.co%3B2

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related citations

A note on a generalized single step theory for any number of hierarchical genomic matrices.

BACKGROUND: The Single Step algorithm allows combining information from genotyped and un-genotyped individuals, provided they are connected by a pedigree. However, current single step theory is limited to a single list of markers. RESULTS: We present a generalized single step (GSS) method that can accommodate any number of hierarchical molecular datasets (e.g. sequence, high and low density arrays) and pedigree, avoiding imputation. We prove that a similar efficient inversion algorithm exists. The method is recursive, starting with the highest marker density scenario. We illustrate the method with simulation and show that GSS can increase predictive accuracy compared to standard single step. R code is provided so that custom scenarios can be easily compared, either with simulated or real data. CONCLUSION: The method developed generalizes extant single step theory to any number of hierarchical molecular relationship matrices, broadening the scenarios where single step can be applied. A topic of particular interest can be ecology field data or human populations where pedigree is not available, but where samples sequenced and genotyped at different densities can exist. GSS can also be a useful tool to optimize allocation of genotyping and / or sequencing resources.

Algorithms↗

cgDist: Nucleotide-level distance calculation from cgMLST allelic profiles.

Bacterial genomic surveillance requires balancing computational efficiency with genetic resolution for effective cluster investigation. cgMLST distance calculations treat all allelic differences as equivalent units, obscuring nucleotide-level variation. Furthermore, single nucleotide polymorphism-based pipelines provide finer resolution at substantially higher computational cost, which limits their routine deployment in surveillance laboratories. We present cgDist, an algorithm that calculates nucleotide-level distances directly from cgMLST allelic profiles, providing finer resolution than allele-count distances by leveraging within-allele nucleotide variation. The cache architecture stores alignment statistics, enabling distance calculation modes without computation and supporting both dataset-specific and schema-complete cache generation. This design enables incremental surveillance analysis, with performance benefits as laboratories accumulate alignment data. cgDist functions as a precision 'zoom lens' for the investigation of clusters identified through initial cgMLST screening. Rather than restructuring population relationships, this targeted approach concentrates enhanced resolution where it is most informative. The algorithm ensures that cgDist distances are greater than or equal to corresponding cgMLST distances, preserving epidemiological interpretability while adding genetic discrimination. By increasing resolution within identified clusters, cgDist may also support outbreak investigation, a potential application that remains to be evaluated on outbreak-derived data.

Algorithms↗

Theseus: fast and optimal affine-gap sequence-to-graph alignment.

MOTIVATION: Sequence-to-graph alignment is a central problem in bioinformatics, with applications in multiple sequence alignment (MSA) and pangenome analysis, among others. However, current algorithms for optimal affine-gap alignment impose high memory and computational requirements, limiting their scalability to aligning long sequences to complex graphs. Practical solutions partially address this problem using heuristic strategies that ultimately trade off optimality for speed. RESULTS: This work presents Theseus, a novel, fast, and optimal affine-gap sequence-to-graph alignment algorithm. Theseus leverages similarities between genomic sequences to accelerate the alignment computation and reduces the overall memory requirements without compromising optimality. To that end, Theseus processes only a subset of the dynamic programming cells, using a sparse-data strategy that enables efficient sequence-to-graph alignment. Moreover, our algorithm supports optimal affine-gap alignment on arbitrary directed graphs, including those with cycles. We evaluate Theseus on two key problems: MSA and pangenome read mapping. For MSA, we compare it against SPOA, abPOA, and POASTA. Theseus is 1.6× to 17.6× faster than POASTA, and 7.3× faster, on average, than SPOA, both optimal aligners. Compared with abPOA, Theseus ensures optimality and scales to the largest problems. For pangenome read mapping, we benchmark Theseus against the alignment stage of the mapping tool vg map, along with the alignment kernels of SPOA, abPOA, and POASTA. Theseus outperforms the other methods, showing a 1.9× to 16.9× speedup on short reads. Moreover, Theseus is 1.5× to 36.3× faster than vg when aligning against synthetic cyclic graphs. AVAILABILITY AND IMPLEMENTATION: Theseus code and documentation are publicly available at https://github.com/albertjimenezbl/theseus-lib.

Algorithms↗