Development of an improved inhibitor of Lats kinases to promote regeneration of mammalian organs.

The Hippo signaling pathway acts as a brake on regeneration in many tissues. This cascade of kinases culminates in the phosphorylation of the transcriptional cofactors Yap and Taz, whose concentration in the nucleus consequently remains low. Various types of cellular signals can reduce phosphorylation, however, resulting in the accumulation of Yap and Taz in the nucleus and subsequently in mitosis. We earlier identified a small molecule, TRULI, that blocks the final kinases in the pathway, Lats1 and Lats2, and thus elicits proliferation of several cell types that are ordinarily postmitotic and aids regeneration in mammals. In the present study, we present the results of chemical modification of the original compound and demonstrate that a derivative, TDI-011536, is an effective blocker of Lats kinases in vitro at nanomolar concentrations. The compound fosters extensive proliferation in retinal organoids derived from human induced pluripotent stem cells. Intraperitoneal administration of the substance to mice suppresses Yap phosphorylation for several hours and induces transcriptional activation of Yap target genes in the heart, liver, and skin. Moreover, the compound initiates the proliferation of cardiomyocytes in adult mice following cardiac cryolesions. After further chemical refinement, related compounds might prove useful in protective and regenerative therapies.

Development of a Highly Selective Plasmodium falciparum Proteasome Inhibitor with Anti-malaria Activity in Humanized Mice.

Plasmodium falciparum proteasome (Pf20S) inhibitors are active against Plasmodium at multiple stages-erythrocytic, gametocyte, liver, and gamete activation stages-indicating that selective Pf20S inhibitors possess the potential to be therapeutic, prophylactic, and transmission-blocking antimalarials. Starting from a reported compound, we developed a noncovalent, macrocyclic peptide inhibitor of the malarial proteasome with high species selectivity and improved pharmacokinetic properties. The compound demonstrates specific, time-dependent inhibition of the ß5 subunit of the Pf20S, kills artemisinin-sensitive and artemisinin-resistant P. falciparum isolates in vitro and reduces parasitemia in humanized, P. falciparum-infected mice.

Dosimetric validation of a magnetic resonance image gated radiotherapy system using a motion phantom and radiochromic film.

PURPOSE: Magnetic resonance image (MRI) guided radiotherapy enables gating directly on the target position. We present an evaluation of an MRI-guided radiotherapy system's gating performance using an MRI-compatible respiratory motion phantom and radiochromic film. Our evaluation is geared toward validation of our institution's clinical gating protocol which involves planning to a target volume formed by expanding 5 mm about the gross tumor volume (GTV) and gating based on a 3 mm window about the GTV. METHODS: The motion phantom consisted of a target rod containing high-contrast target inserts which moved in the superior-inferior direction inside a body structure containing background contrast material. The target rod was equipped with a radiochromic film insert. Treatment plans were generated for a 3 cm diameter spherical planning target volume, and delivered to the phantom at rest and in motion with and without gating. Both sinusoidal trajectories and tumor trajectories measured during MRI-guided treatments were used. Similarity of the gated dose distribution to the planned, motion-frozen, distribution was quantified using the gamma technique. RESULTS: Without gating, gamma pass rates using 4%/3 mm criteria were 22-59% depending on motion trajectory. Using our clinical standard of repeated breath holds and a gating window of 3 mm with 10% target allowed outside the gating boundary, the gamma pass rate was 97.8% with 3%/3 mm gamma criteria. Using a 3 mm window and 10% allowed excursion, all of the patient tumor motion trajectories at actual speed resulting in at least 95% gamma pass rate at 4%/3 mm. CONCLUSIONS: Our results suggest that the device can be used to compensate respiratory motion using a 3 mm gating margin and 10% allowed excursion results in conjunction with repeated breath holds. Full clinical validation requires a comprehensive evaluation of tracking performance in actual patient images, outside the scope of this study.

A Dose Accumulation Assessment of Alignment Errors During Spatially Fractionated Radiation Therapy.

PURPOSE: Spatially fractionated radiation therapy (SFRT) techniques produce high-dose peaks and low-dose valleys within a tumor. Lattice stereotactic body radiation therapy (SBRT) is a form a SFRT delivered across 5 fractions. Because of the high spatial dose gradients associated with SFRT, it is critical for fractionated SFRT patients to be aligned correctly for treatment. Here we investigate the dosimetric effect of daily alignment uncertainty through a dose accumulation study. METHODS AND MATERIALS: Dose accumulation was retrospectively performed for 10 patients enrolled on a phase 1 trial. Lattice stereotactic body radiation therapy was completed in 5 fractions with 20 Gy prescribed to the entire tumor and a simultaneous integrated boost of 66.7 Gy prescribed to a set of regularly spaced high-dose spheres. Daily alignment error was quantified through manually selected landmarks in both the planning computed tomography scan and daily cone beam computed tomography. The dosimetric effect of alignment errors was quantified by translating the isocenter in the treatment planning system by the daily average alignment error. Large errors were simulated by translating isocenter 5 and 10 mm for 1 and 2 fractions, independently assessing errors in the superior-inferior and axial directions. The reduction of dose gradients was quantified using the dose ratio (DR) of the mean dose in the high-dose and low-dose spheres. RESULTS: The average alignment error was 1.8 mm across the patient population resulting in minor smoothing of the high- and low-dose distributions in the dose accumulation. Quantitatively, the DR decreased from 3.42 to 3.32 (P = .093) in the dose accumulation study. The simulated worst case was an inferior-superior shift of 10 mm for 2 fractions where the average DR decreased to 2.72 (P = .0001). CONCLUSIONS: The dose accumulation study revealed on average DR only decreased from 3.42 to 3.32. However, setup errors >5 mm resulted in larger dosimetric degradation, reflecting a larger effect for individual high-dose spheres within regions exhibiting larger displacements.

A clinical and time savings evaluation of a deep learning automatic contouring algorithm.

Automatic contouring algorithms may streamline clinical workflows by reducing normal organ-at-risk (OAR) contouring time. Here we report the first comprehensive quantitative and qualitative evaluation, along with time savings assessment for a prototype deep learning segmentation algorithm from Siemens Healthineers. The accuracy of contours generated by the prototype were evaluated quantitatively using the Sorensen-Dice coefficient (Dice), Jaccard index (JC), and Hausdorff distance (Haus). Normal pelvic and head and neck OAR contours were evaluated retrospectively comparing the automatic and manual clinical contours in 100 patient cases. Contouring performance outliers were investigated. To quantify the time savings, a certified medical dosimetrist manually contoured de novo and, separately, edited the generated OARs for 10 head and neck and 10 pelvic patients. The automatic, edited, and manually generated contours were visually evaluated and scored by a practicing radiation oncologist on a scale of 1-4, where a higher score indicated better performance. The quantitative comparison revealed high (> 0.8) Dice and JC performance for relatively large organs such as the lungs, brain, femurs, and kidneys. Smaller elongated structures that had relatively low Dice and JC values tended to have low Hausdorff distances. Poor performing outlier cases revealed common anatomical inconsistencies including overestimation of the bladder and incorrect superior-inferior truncation of the spinal cord and femur contours. In all cases, editing contours was faster than manual contouring with an average time saving of 43.4% or 11.8 minutes per patient. The physician scored 240 structures with > 95% of structures receiving a score of 3 or 4. Of the structures reviewed, only 11 structures needed major revision or to be redone entirely. Our results indicate the evaluated auto-contouring solution has the potential to reduce clinical contouring time. The algorithm's performance is promising, but human review and some editing is required prior to clinical use.

Knowledge-based adaptive planning quality assurance using dosimetric indicators for stereotactic adaptive radiotherapy for pancreatic cancer.

INTRODUCTION: We aimed to develop knowledge-based tools for robust adaptive radiotherapy (ART) planning to determine on-table adaptive DVH metric variations or planning process errors for stereotactic pancreatic ART. We developed volume-based dosimetric identifiers to identify deviations of ART plans from simulation plans. MATERIALS AND METHODS: Two patient cohorts who were treated on MR-Linac for pancreas cancer were included in this retrospective study; a training cohort and a validation cohort. All patients received 50 Gy in 5 fractions. PTV-OPT was generated by subtracting the critical organs plus a 5 mm-margin from PTV. Several metrics that potentially can identify failure-modes were calculated including PTV & PTV_OPT V95% and PTV & PTV_OPT D95%/D5%. The difference between each DVH metric in each adaptive plan with the DVH metric in simulation plan was calculated. The 95% confidence interval (CI) of the variations in each DVH metric was calculated for the patient training cohort. Variations in DVH metrics that exceeded the 95% CI for all fractions in training and validation cohort were flagged for retrospective investigation for root-cause analysis to determine their predictive power for identifying failure-modes. RESULTS: The CIs for the PTV & PTV_OPT V95% and PTV & PTV_OPT D95%/D5% were ± 13%, ± 5%, ± 0.1, ± 0.03, respectively. We estimated the positive predictive value and negative predictive value of our method to be 77% and 89%, respectively, for the training cohort, and 80% for both in the validation cohort. DISCUSSION: We developed dosimetric indicators for ART planning QA to identify population-based deviations or planning errors during online adaptive process for stereotactic pancreatic ART. This technology may be useful as an ART clinical trial QA tool and improve overall ART quality at an institution.

High-throughput screening identifies small molecule inhibitors of thioesterase superfamily member 1: Implications for the management of non-alcoholic fatty liver disease.

OBJECTIVE: Thioesterase superfamily member 1 (Them1) is a long chain acyl-CoA thioesterase comprising two N-terminal HotDog fold enzymatic domains linked to a C-terminal lipid-sensing steroidogenic acute regulatory transfer-related (START) domain, which allosterically modulates enzymatic activity. Them1 is highly expressed in thermogenic adipose tissue, where it functions to suppress energy expenditure by limiting rates of fatty acid oxidation, and is induced markedly in liver in response to high fat feeding, where it suppresses fatty acid oxidation and promotes glucose production. Them1-/- mice are protected against non-alcoholic fatty liver disease (NAFLD), suggesting Them1 as a therapeutic target. METHODS: A high-throughput small molecule screen was performed to identify promising inhibitors targeting the fatty acyl-CoA thioesterase activity of purified recombinant Them1.Counter screening was used to determine specificity for Them1 relative to other acyl-CoA thioesterase isoforms. Inhibitor binding and enzyme inhibition were quantified by biophysical and biochemical approaches, respectively. Following structure-based optimization, lead compounds were tested in cell culture. RESULTS: Two lead allosteric inhibitors were identified that selectively inhibited Them1 by binding the START domain. In mouse brown adipocytes, these inhibitors promoted fatty acid oxidation, as evidenced by increased oxygen consumption rates. In mouse hepatocytes, they promoted fatty acid oxidation, but also reduced glucose production. CONCLUSION: Them1 inhibitors could prove attractive for the pharmacologic management of NAFLD.

Shape-Based Virtual Screening of a Billion-Compound Library Identifies Mycobacterial Lipoamide Dehydrogenase Inhibitors.

Lpd (lipoamide dehydrogenase) in Mycobacterium tuberculosis (Mtb) is required for virulence and is a genetically validated tuberculosis (TB) target. Numerous screens have been performed over the last decade, yet only two inhibitor series have been identified. Recent advances in large-scale virtual screening methods combined with make-on-demand compound libraries have shown the potential for finding novel hits. In this study, the Enamine REAL library consisting of ∼1.12 billion compounds was efficiently screened using the GPU Shape screen method against Mtb Lpd to find additional chemical matter that would expand on the known sulfonamide inhibitor series. We identified six new inhibitors with IC50 in the range of 5-100 µM. While these compounds remained chemically close to the already known sulfonamide series inhibitors, some diversity was found in the cores of the hits. The two most potent hits were further validated by one-step potency optimization to submicromolar levels. The co-crystal structure of optimized analogue TDI-13537 provided new insights into the potency determinants of the series.

Structure-Activity Relationship Studies of Antimalarial Plasmodium Proteasome Inhibitors─Part II.

With increasing reports of resistance to artemisinins and artemisinin-combination therapies, targeting the Plasmodium proteasome is a promising strategy for antimalarial development. We recently reported a highly selective Plasmodium falciparum proteasome inhibitor with anti-malarial activity in the humanized mouse model. To balance the permeability of the series of macrocycles with other drug-like properties, we conducted further structure-activity relationship studies on a biphenyl ether-tethered macrocyclic scaffold. Extensive SAR studies around the P1, P3, and P5 groups and peptide backbone identified compound TDI-8414. TDI-8414 showed nanomolar antiparasitic activity, no toxicity to HepG2 cells, high selectivity against the Plasmodium proteasome over the human constitutive proteasome and immunoproteasome, improved solubility and PAMPA permeability, and enhanced metabolic stability in microsomes and plasma of both humans and mice.

Design, Synthesis, and Optimization of Macrocyclic Peptides as Species-Selective Antimalaria Proteasome Inhibitors.

With over 200 million cases and close to half a million deaths each year, malaria is a threat to global health, particularly in developing countries. Plasmodium falciparum, the parasite that causes the most severe form of the disease, has developed resistance to all antimalarial drugs. Resistance to the first-line antimalarial artemisinin and to artemisinin combination therapies is widespread in Southeast Asia and is emerging in sub-Saharan Africa. The P. falciparum proteasome is an attractive antimalarial target because its inhibition kills the parasite at multiple stages of its life cycle and restores artemisinin sensitivity in parasites that have become resistant through mutation in Kelch K13. Here, we detail our efforts to develop noncovalent, macrocyclic peptide malaria proteasome inhibitors, guided by structural analysis and pharmacokinetic properties, leading to a potent, species-selective, metabolically stable inhibitor.

A Chemical Strategy toward Novel Brain-Penetrant EZH2 Inhibitors.

Aberrant gene-silencing through dysregulation of polycomb protein activity has emerged as an important oncogenic mechanism in cancer, implicating polycomb proteins as important therapeutic targets. Recently, an inhibitor targeting EZH2, the methyltransferase component of PRC2, received U.S. Food and Drug Administration approval following promising clinical responses in cancer patients. However, the current array of EZH2 inhibitors have poor brain penetrance, limiting their use in patients with central nervous system malignancies, a number of which have been shown to be sensitive to EZH2 inhibition. To address this need, we have identified a chemical strategy, based on computational modeling of pyridone-containing EZH2 inhibitor scaffolds, to minimize P-glycoprotein activity, and here we report the first brain-penetrant EZH2 inhibitor, TDI-6118 (compound 5). Additionally, in the course of our attempts to optimize this compound, we discovered TDI-11904 (compound 21), a novel, highly potent, and peripherally active EZH2 inhibitor based on a 7 member ring structure.

Implementing stereotactic accelerated partial breast irradiation using magnetic resonance guided radiation therapy.

INTRODUCTION: Accelerated partial breast irradiation (APBI) seeks to reduce irradiated volumes and radiation exposure for patients while maintaining acceptable clinical outcomes. Magnetic resonance image-guided radiotherapy (MRgRT) provides excellent soft-tissue contrast for treatment localization, which can reduce setup uncertainty, thus reducing margins in the external beam setting. Additionally, stereotactic body radiotherapy (SBRT)-style regimens with high gradients can also be executed. This MR-guided stereotactic APBI (MRgS-APBI) approach can be utilized for a lower number of fractions and spare a greater volume of healthy tissues compared to conventional 3D external beam APBI. METHODS: Our MRgS-APBI program was developed for two prospective non-randomized phase I/II clinical trials (20Gyx1 and 8.5Gyx3). Both breast SBRT treatment planning and MRgRT delivery techniques were described in this study. Simulation included both CT and MRI with specialized immobilization to accommodate MR-guided setup and cine-MRI treatment gating. Dosimetry data from 48 single-fraction and 19 three-fraction patients were collected and evaluated. This included planning objectives and SBRT-specific indices. During treatment, setup errors were calculated to evaluate setup reproducibility and duty cycle was calculated using cine-MRI data during gated delivery. RESULTS: In both the single- and three- fraction trials combined, 88.5% of the possible dosimetric objectives across all patients were met during planning. The majority of the planning objectives were easily achievable indicating the potential for stricter objectives for subsequent S-APBI treatments. The average magnitude of setup uncertainties was 1.0 cm ±â€¯0.6 cm across all treatments. In the three-fraction trial, the average beam-on duty-cycle for the MRI-gated delivery was 83.0 ±â€¯13.0%. There were no technical MRgS-APBI related issues that resulted in discontinuation of treatment across all patients. CONCLUSION: SBRT-style dosimetry and delivery for APBI is feasible using MR-guidance. The program development and dosimetric outcomes reported here can serve as a guide for other institutions considering the clinical implementation of MR-guided stereotactic APBI.

Whole Cell Active Inhibitors of Mycobacterial Lipoamide Dehydrogenase Afford Selectivity over the Human Enzyme through Tight Binding Interactions.

Tuberculosis remains a leading cause of death from a single bacterial infection worldwide. Efforts to develop new treatment options call for expansion into an unexplored target space to expand the drug pipeline and bypass resistance to current antibiotics. Lipoamide dehydrogenase is a metabolic and antioxidant enzyme critical for mycobacterial growth and survival in mice. Sulfonamide analogs were previously identified as potent and selective inhibitors of mycobacterial lipoamide dehydrogenase in vitro but lacked activity against whole mycobacteria. Here we present the development of analogs with improved permeability, potency, and selectivity, which inhibit the growth of Mycobacterium tuberculosis in axenic culture on carbohydrates and within mouse primary macrophages. They increase intrabacterial pyruvate levels, supporting their on-target activity within mycobacteria. Distinct modalities of binding between the mycobacterial and human enzymes contribute to improved potency and hence selectivity through induced-fit tight binding interactions within the mycobacterial but not human enzyme, as indicated by kinetic analysis and crystallography.

Macrocyclic Peptides that Selectively Inhibit the Mycobacterium tuberculosis Proteasome.

Treatment of tuberculosis (TB) currently takes at least 6 months. Latent Mycobacterium tuberculosis (Mtb) is phenotypically tolerant to most anti-TB drugs. A key hypothesis is that drugs that kill nonreplicating (NR) Mtb may shorten treatment when used in combination with conventional drugs. The Mtb proteasome (Mtb20S) could be such a target because its pharmacological inhibition kills NR Mtb and its genetic deletion renders Mtb unable to persist in mice. Here, we report a series of macrocyclic peptides that potently and selectively target the Mtb20S over human proteasomes, including macrocycle 6. The cocrystal structure of macrocycle 6 with Mtb20S revealed structural bases for the species selectivity. Inhibition of 20S within Mtb by 6 dose dependently led to the accumulation of Pup-tagged GFP that is degradable but resistant to depupylation and death of nonreplicating Mtb under nitrosative stress. These results suggest that compounds of this class have the potential to develop as anti-TB therapeutics.

Substituted 2H-isoquinolin-1-ones as potent Rho-kinase inhibitors: part 3, aryl substituted pyrrolidines.

The discovery and SAR of a series of beta-aryl substituted pyrrolidine 2H-isoquinolin-1-one inhibitors of Rho-kinase (ROCK) derived from 2 is herein described. SAR studies have shown that aryl groups in the beta-position are optimal for potency. Our efforts focused on improving the ROCK potency of this isoquinolone class of inhibitors which led to the identification of pyrrolidine 32 which demonstrated a 10-fold improvement in aortic ring (AR) potency over 2.

Substituted 2H-isoquinolin-1-ones as potent Rho-kinase inhibitors: part 2, optimization for blood pressure reduction in spontaneously hypertensive rats.

Phenylglycine substituted isoquinolones 1 and 2 have previously been described as potent dual ROCK1/ROCK2 inhibitors. Here we describe the further SAR of this series to improve metabolic stability and rat oral exposure. Piperidine analog 20 which demonstrates sustained blood pressure normalization in an SHR blood pressure reduction model was identified through this effort.

A motion prediction confidence estimation framework for prediction-based radiotherapy gating.

PURPOSE: Motion prediction can compensate for latency in image-guided radiotherapy and has been an active area of research. However, motion predictions are subject to error and variations. We have developed and evaluated a novel motion prediction confidence estimation framework to improve the efficacy and robustness of prediction-based radiotherapy gating decision-making. The specific scenario of adaptive gating in magnetic resonance imaging (MRI)-guided radiotherapy is studied as an example, but the method generalizes to other modalities and motion management setups. METHODS: The proposed prediction confidence estimator is based on a generic training/testing paradigm and consists of a weighted combination of three components: the prediction model's goodness of fit, variation in the prediction using a leave-one-out process and the velocity of the tracked target. Roughly, these terms quantify respectively the consistency between prediction and the training data, the robustness of model inference, and the stability due to target speed. The weight parameters and the action level in triggering beam-off decision are optimized. The method is assessed and validated in 8 healthy volunteer and 13 patient studies using a 0.35T MRI-guided radiotherapy system predicting 0.25-0.33 s ahead. The effect of the action level on the predicted gating decision accuracy, beam-on positive predictive value (PPV) and median distance between the predicted and ground-truth target centroids were evaluated. Statistical significance was evaluated using a paired t-test. The tradeoff between these performance metrics and gating duty cycle was assessed. RESULTS: Use of the confidence estimator threshold increased gating accuracy by up to 2.42%, increased PPV by up to 3.00%, and reduced the median centroid distance up to 0.28 mm. The confidence estimator threshold on average increased gating accuracy to 96.5% (P = 2.08 × 10-4 ), increased PPV to 96.7% (P = 1.46 × 10-5 ), reduced the median centroid distance to 0.54 mm (P = 1.71 × 10-5 ) at the cost of reducing the gating duty cycle by 14.3% to 48.5%. Hyperparameter tuning revealed that contrary to intuition, the velocity term offered only minimal performance improvement in some cases but also introduced potential stability issues. The combination of goodness of fit and leave-one-out prediction variation provided the most effective confidence estimator, yielding universally better performance in gating decisions. CONCLUSION: Confidence estimation utilizing prediction model fitness criterion and validation principles can complement prediction methods to guide MRI-guided radiotherapy gating. Results from both volunteer and patient studies showed improved gating quality.

An image regression motion prediction technique for MRI-guided radiotherapy evaluated in single-plane cine imaging.

PURPOSE: To develop and evaluate a novel motion prediction method for magnetic resonance image (MRI)-guided radiotherapy applications. This method, which we deem "image regression," predicts future tissue motion based on a weighted combination of previously observed motion states. Motion predictions are derived from a sliding window of recent motion states which are defined by a temporal sequence of images. A key advantage of this method compared to other motion prediction methods is that its computational complexity scales weekly with the number of spatial points predicted. Applications of gating latency reduction and improvement in deformable registration-based target tracking are demonstrated. METHODS: The image regression (IR) motion prediction method was developed and evaluated using 26.9 h of real-time imaging acquired from eight healthy volunteers and 13 patients using a 0.35 T MRI-guided radiotherapy system. Motion predictions were performed 0.25-0.33 s into the future using a weighted sum of previously observed motion states with image similarity-derived weights. The set of previously observed motion states were continuously updated to incorporate the changes in breathing patterns. The accuracy of the predicted radiotherapy gating decision, beam-on positive predictive value (PPV), and predicted vs ground-truth target centroid position errors are reported. The IR technique was compared against no prediction, linear extrapolation, and an established autoregressive linear prediction algorithm. The usage of IR to initialize the deformable registration and enhance the target tracking was demonstrated in the healthy volunteer studies. Deformable registration with IR initialization was compared to the initialization performed by current clinical software: no initialization, previous image registration initialization and linear motion extrapolation initialization. RESULTS: The average IR-predicted radiation beam gating decision accuracy was 95.8%, with a PPV of 95.7%, and median and 95th percentile centroid position errors of 0.63 and 2.08 mm, respectively. Compared to the autoregressive linear prediction method, gating accuracy was 1.15% greater, PPV was 1.61% greater, and median and 95th percentile centroid distances were 0.21 and 0.23 mm smaller. The IR-initialized registration on average converged within 0.50 mm of the ground-truth position in fewer than 10 iterations whereas the next best initialization method required more than 25 iterations. CONCLUSIONS: Image regression motion prediction has the potential to reduce the gating latencies and improve the speed and accuracy of deformable registration-based target tracking in MRI-guided radiotherapy.

The discovery of thienopyridine analogues as potent IkappaB kinase beta inhibitors. Part II.

An SAR study that identified a series of thienopyridine-based potent IkappaB Kinase beta (IKKbeta) inhibitors is described. With focuses on the structural optimization at C4 and C6 of structure 1 (Fig. 1), the study reveals that small alkyl and certain aromatic groups are preferred at C4, whereas polar groups with proper orientation at C6 efficiently enhance compound potency. The most potent analogues inhibit IKKbeta with IC50s as low as 40 nM, suppress LPS-induced TNF-alpha production in vitro and in vivo, display good kinase selectivity profiles, and are active in a HeLa cell NF-kappaB reporter gene assay, demonstrating that they directly interfere with the NF-kappaB signaling pathway.

Multislice motion modeling for MRI-guided radiotherapy gating.

PURPOSE: On-board magnetic resonance imaging (MRI) greatly enhances real-time target tracking capability during radiotherapy treatments. However, multislice and volumetric MRI techniques are frame rate limited and introduce unacceptable latency between the target moving out of position and the beam being turned off. We present a technique to estimate continuous volumetric tissue motion using motion models built from a repeated acquisition of a stack of MR slices. Applications including multislice target visualization and out-of-slice motion estimation during MRI-guided radiotherapy are demonstrated. METHODS: Eight healthy volunteer studies were performed using a 0.35 T MRI-guided radiotherapy system. Images were acquired at three frames per second in an interleaved fashion across ten adjacent sagittal slice positions covering 4.5 cm using a balanced steady-state-free precession sequence. A previously published five-dimensional (5D) linear motion model used for MRI-guided radiotherapy gating was extended to include multiple slices. This model utilizes an external respiratory bellows signal recorded during imaging to simultaneously estimate motion across all imaged slices. For comparison to an image-based approach, the manifold learning technique local linear embedding (LLE) was used to derive a respiratory surrogate for motion modeling. Manifolds for every slice were aligned during LLE in a group-wise fashion, enabling motion estimation outside the current imaged slice using a motion model, a process which we denote as mSGA. Additionally, a method is developed to evaluate out-of-slice motion estimates. The multislice motion model was evaluated in a single slice with each newly acquired image using a leave-one-out approach. Model-generated gating decision accuracy and beam-on positive predictive value (PPV) are reported along with the median and 95th percentile distance between model and ground truth target centroids. RESULTS: The average model gating decision accuracy and PPV across all volunteer studies was 93.7% and 92.8% using the 5D model, and 96.8% and 96.1% using the mSGA model, respectively. The median and 95th percentile distance between model and ground truth target centroids was 0.91 and 2.90 mm, respectively, using the 5D model and 0.58 and 1.49 mm using the mSGA model, averaged over all eight subjects. The mSGA motion model provided a statistically significant improvement across all evaluation metrics compared to the external surrogate-based 5D model. CONCLUSION: The proposed techniques for out-of-slice target motion estimation demonstrated accuracy likely sufficient for clinical use. Results indicate the mSGA model may provide higher accuracy, however, the external surrogate-based model allows for unbiased in vivo accuracy evaluation.