Understanding the Continual Reassessment Method for Dose Finding Studies: An Overview for Non-Statisticians

The Continual Reassessment Method (CRM) has gained popularity since its proposal by O’Quigley et al. [1]. Many variations have been published and discussed in the statistical literature, but there has been little attention to making the design considerations accessible to non-statisticians. As a result, some clinicians or reviewers of clinical trials tend to be wary of the CRM due to safety concerns. This paper presents the CRM in a non-technical way, describing the original CRM with some of its modified versions. It also describes the specifications that define a CRM design, along with two simulated examples of CRMs for illustration.

Up-regulation of somatostatin receptor density on rat CA20948 tumors escaped from low dose [(177)Lu-DOTA(0),Tyr(3)]octreotate therapy

AIM: Peptide receptor radionuclide therapy using the somatostatin analogue [(177)Lu-DOTA(0),Tyr(3)]octreotate is a convincing treatment modality for metastasized neuroendocrine tumors. Therapeutic doses are administered in 4 cycles with 6-10 week intervals. A high somatostatin receptor density on tumor cells is a prerequisite at every administration to enable effective therapy. In this study, the density of the somatostatin receptor subtype 2 (sst2) was investigated in the rat CA20948 pancreatic tumor model after low dose [(177)Lu-DOTA(0), Tyr(3)]octreotate administration resulting in approximately 20 Gy tumor radiation absorbed dose, whereas 60 Gy is needed to induce complete tumor regression in these and the majority of tumors. METHODS: Sixteen days after inoculation of the CA20948 tumor, male Lewis rats were injected with 185 MBq [(177)Lu-DOTA(0),Tyr(3)]octreotate to initiate a decline in tumor size. Approximately 40 days after injection, tumors re-grew progressively after initial response. Quantification of sst2 expression was performed using in vitro autoradiography on frozen sections of three groups: control (not-treated) tumors, tumors in regression and tumors in re-growth. Histology and proliferation were determined using HE- and anti-Ki-67-staining. RESULTS: The sst2 expression on CA20948 tumor cells decreased significantly after therapy to 5% of control level. However, tumors escaping from therapy showed an up-regulated sst2 level of 2-5 times higher sst2 density compared to control tumors. CONCLUSION: After a suboptimal therapeutic dose of [(177)Lu-DOTA(0),Tyr(3)]octreotate, escape of tumors is likely to occur. Since these cells show an up-regulated sst2 receptor density, a next therapeutic administration of radiolabelled sst2 analogue can be expected to be highly effective.

Copper and human health: biochemistry, genetics, and strategies for modeling dose-response relationships

Copper (Cu) and its alloys are used extensively in domestic and industrial applications. Cu is also an essential element in mammalian nutrition. Since both copper deficiency and copper excess produce adverse health effects, the dose-response curve is U-shaped, although the precise form has not yet been well characterized. Many animal and human studies were conducted on copper to provide a rich database from which data suitable for modeling the dose-response relationship for copper may be extracted. Possible dose-response modeling strategies are considered in this review, including those based on the benchmark dose and categorical regression. The usefulness of biologically based dose-response modeling techniques in understanding copper toxicity was difficult to assess at this time since the mechanisms underlying copper-induced toxicity have yet to be fully elucidated. A dose-response modeling strategy for copper toxicity was proposed associated with both deficiency and excess. This modeling strategy was applied to multiple studies of copper-induced toxicity, standardized with respect to severity of adverse health outcomes and selected on the basis of criteria reflecting the quality and relevance of individual studies. The use of a comprehensive database on copper-induced toxicity is essential for dose-response modeling since there is insufficient information in any single study to adequately characterize copper dose-response relationships. The dose-response modeling strategy envisioned here is designed to determine whether the existing toxicity data for copper excess or deficiency may be effectively utilized in defining the limits of the homeostatic range in humans and other species. By considering alternative techniques for determining a point of departure and low-dose extrapolation (including categorical regression, the benchmark dose, and identification of observed no-effect levels) this strategy will identify which techniques are most suitable for this purpose. This analysis also serves to identify areas in which additional data are needed to better define the characteristics of dose-response relationships for copper-induced toxicity in relation to excess or deficiency.

Effect of patient size on radiation dose for abdominal MDCT with automatic tube current modulation: phantom study

OBJECTIVE: The purpose of this study was to evaluate in a phantom study the effect of patient size on radiation dose for abdominal MDCT with automatic tube current modulation. MATERIALS AND METHODS: One or two 4-cm-thick circumferential layers of fat-equivalent material were added to the abdomen of an anthropomorphic phantom to simulate patients of three sizes: small (cross-sectional dimensions, 18 x 22 cm), average size (26 x 30 cm), and oversize (34 x 38 cm). Imaging was performed with a 64-MDCT scanner with combined z-axis and xy-axis tube current modulation according to two protocols: protocol A had a noise index of 12.5 H, and protocol B, 15.0 H. Radiation doses to three abdominal organs and the skin were assessed. Image noise also was measured. RESULTS: Despite increasing patient size, the image noise measured was similar for protocol A (range, 11.7-12.2 H) and protocol B (range, 13.9-14.8 H) (p > 0.05). With the two protocols, in comparison with the dose of the small patient, the abdominal organ doses of the average-sized patient and the oversized patient increased 161.5-190.6%and 426.9-528.1%, respectively (p < 0.001). The skin dose increased as much as 268.6% for the average-sized patient and 816.3% for the oversized patient compared with the small patient (p < 0.001). CONCLUSION: Oversized patients undergoing abdominal MDCT with tube current modulation receive significantly higher doses than do small patients. The noise index needs to be adjusted to the body habitus to ensure dose efficiency.

Time- and dose-related regional fluxes of tissue-type plasminogen activator in anesthetized endotoxemic pigs

BACKGROUND: Acute endotoxinemia elicits an early fibrinolytic response. This study analyzes the effects of the dose and duration of endotoxin infusion on arterial levels of tissue-type plasminogen activator (tPA) and pulmonary, mesenteric and hepatic plasma tPA fluxes. METHODS: Pigs were randomized to receive an acute, high-dose (for 6 h, n=13, high ETX) or a prolonged, low-dose (for 18 h, n=18, low ETX) infusion of endotoxin or saline vehicle alone (for 18 h, n=14, control). All animals were fluid resuscitated to maintain a normodynamic circulation. Systemic and regional blood flows were measured and arterial, pulmonary arterial, portal and hepatic venous blood samples were analyzed to calculate regional net fluxes of tPA. Plasma tumor necrosis factor (TNF-alpha) levels were analyzed. RESULTS: Mesenteric tPA release and hepatic uptake increased maximally at 1.5 h in ETX groups related to dose. Maximal mesenteric tPA release [high ETX 612 (138-1185) microg/min/kg, low ETX 72 (32-94) microg/min/kg, median+/-interquartile range] and hepatic tPA uptake [high ETX -1549 (-1134 to -2194) microg/min/kg, low ETX -153 (-105 to -307) microg/min/kg] correlated to TNF-alpha levels. Regional tPA fluxes returned to baseline levels at 6 h in both ETX groups and also remained low during sustained low ETX. No changes were observed in control animals. CONCLUSIONS: Endotoxemia induces an early increase in mesenteric tPA release and hepatic tPA uptake related to the severity of endotoxemia. The time patterns of changes in mesenteric and hepatic tPA fluxes are similar in acute high-dose endotoxemia and sustained low-dose endotoxemia.

Multilead quantitative electroencephalogram profile and cognitive evoked potentials (P300) in healthy subjects after a single dose of olanzapine

RATIONALE: Olanzapine is an atypical antipsychotic drug with a more favourable safety profile than typical antipsychotics with a hitherto unknown topographic quantitative electroencephalogram (QEEG) profile. OBJECTIVES: We investigated electrical brain activity (QEEG and cognitive event related potentials, ERPs) in healthy subjects who received olanzapine. METHODS: Vigilance-controlled, 19-channel EEG and ERP in an auditory odd-ball paradigm were recorded before and 3 h, 6 h and 9 h after administration of either a single dose of placebo or olanzapine (2.5 mg and 5 mg) in ten healthy subjects. QEEG was analysed by spectral analysis and evaluated in nine frequency bands. For the P300 component in the odd-ball ERP, the amplitude and latency was analysed. Statistical effects were tested using a repeated-measurement analysis of variance. RESULTS: For the interaction between time and treatment, significant effects were observed for theta, alpha-2, beta-2 and beta-4 frequency bands. The amplitude of the activity in the theta band increased most significantly 6 h after the 5-mg administration of olanzapine. A pronounced decrease of the alpha-2 activity especially 9 h after 5 mg olanzapine administration could be observed. In most beta frequency bands, and most significantly in the beta-4 band, a dose-dependent decrease of the activity beginning 6 h after drug administration was demonstrated. Topographic effects could be observed for the beta-2 band (occipital decrease) and a tendency for the alpha-2 band (frontal increase and occipital decrease), both indicating a frontal shift of brain electrical activity. There were no significant changes in P300 amplitude or latency after drug administration. Conclusion: QEEG alterations after olanzapine administration were similar to EEG effects gained by other atypical antipsychotic drugs, such as clozapine. The increase of theta activity is comparable to the frequency distribution observed for thymoleptics or antipsychotics for which treatment-emergent somnolence is commonly observed, whereas the decrease of beta activity observed after olanzapine administration is not characteristic for these drugs. There were no clear signs for an increased cerebral excitability after a single-dose administration of 2.5 mg and 5 mg olanzapine in healthy controls.

High-dose chemotherapy followed by autologous stem cell transplantation for relapsed or refractory Hodgkin lymphoma: prognostic features and outcomes

Between January 1990 and April 2001, 115 patients received high-dose chemotherapy (HDT) followed by autologous stem cell transplantation (ASCT) for relapsed or refractory Hodgkin lymphoma (HL). With a median follow-up of 58 months (range, 1 - 175 months), 5-year progression-free survival (PFS) and overall survival (OS) were 46% and 58%, respectively. Twelve patients with primary refractory disease had a 5-year PFS of 41% and OS of 58%, not significantly different from those of the remaining cohort. Early and overall regimen related mortality were 7% and 16%, respectively. Male gender (P = 0.04) and a time to relapse (TTR) < 12 months (P = 0.03) were associated with decreased OS by univariate analysis. In multivariate analysis, TTR < 12 months remained statistically significant (P = 0.04). We have confirmed that HDT and ASCT result in long-term survival for a proportion of patients with relapsed or refractory HL. All patients, including those with primary refractory disease, benefited from HDT and ASCT.

Collapsed cone convolution and analytical anisotropic algorithm dose calculations compared to VMC++ Monte Carlo simulations in clinical cases

The purpose of this work was to study and quantify the differences in dose distributions computed with some of the newest dose calculation algorithms available in commercial planning systems. The study was done for clinical cases originally calculated with pencil beam convolution (PBC) where large density inhomogeneities were present. Three other dose algorithms were used: a pencil beam like algorithm, the anisotropic analytic algorithm (AAA), a convolution superposition algorithm, collapsed cone convolution (CCC), and a Monte Carlo program, voxel Monte Carlo (VMC++). The dose calculation algorithms were compared under static field irradiations at 6 MV and 15 MV using multileaf collimators and hard wedges where necessary. Five clinical cases were studied: three lung and two breast cases. We found that, in terms of accuracy, the CCC algorithm performed better overall than AAA compared to VMC++, but AAA remains an attractive option for routine use in the clinic due to its short computation times. Dose differences between the different algorithms and VMC++ for the median value of the planning target volume (PTV) were typically 0.4% (range: 0.0 to 1.4%) in the lung and -1.3% (range: -2.1 to -0.6%) in the breast for the few cases we analysed. As expected, PTV coverage and dose homogeneity turned out to be more critical in the lung than in the breast cases with respect to the accuracy of the dose calculation. This was observed in the dose volume histograms obtained from the Monte Carlo simulations.

Conversion of CT numbers into tissue parameters for Monte Carlo dose calculations: a multi-centre study

The conversion of computed tomography (CT) numbers into material composition and mass density data influences the accuracy of patient dose calculations in Monte Carlo treatment planning (MCTP). The aim of our work was to develop a CT conversion scheme by performing a stoichiometric CT calibration. Fourteen dosimetrically equivalent tissue subsets (bins), of which ten bone bins, were created. After validating the proposed CT conversion scheme on phantoms, it was compared to a conventional five bin scheme with only one bone bin. This resulted in dose distributions D(14) and D(5) for nine clinical patient cases in a European multi-centre study. The observed local relative differences in dose to medium were mostly smaller than 5%. The dose-volume histograms of both targets and organs at risk were comparable, although within bony structures D(14) was found to be slightly but systematically higher than D(5). Converting dose to medium to dose to water (D(14) to D(14wat) and D(5) to D(5wat)) resulted in larger local differences as D(5wat) became up to 10% higher than D(14wat). In conclusion, multiple bone bins need to be introduced when Monte Carlo (MC) calculations of patient dose distributions are converted to dose to water.

A multiple source model for 6 MV photon beam dose calculations using Monte Carlo

A multiple source model (MSM) for the 6 MV beam of a Varian Clinac 2300 C/D was developed by simulating radiation transport through the accelerator head for a set of square fields using the GEANT Monte Carlo (MC) code. The corresponding phase space (PS) data enabled the characterization of 12 sources representing the main components of the beam defining system. By parametrizing the source characteristics and by evaluating the dependence of the parameters on field size, it was possible to extend the validity of the model to arbitrary rectangular fields which include the central 3 x 3 cm2 field without additional precalculated PS data. Finally, a sampling procedure was developed in order to reproduce the PS data. To validate the MSM, the fluence, energy fluence and mean energy distributions determined from the original and the reproduced PS data were compared and showed very good agreement. In addition, the MC calculated primary energy spectrum was verified by an energy spectrum derived from transmission measurements. Comparisons of MC calculated depth dose curves and profiles, using original and PS data reproduced by the MSM, agree within 1% and 1 mm. Deviations from measured dose distributions are within 1.5% and 1 mm. However, the real beam leads to some larger deviations outside the geometrical beam area for large fields. Calculated output factors in 10 cm water depth agree within 1.5% with experimentally determined data. In conclusion, the MSM produces accurate PS data for MC photon dose calculations for the rectangular fields specified.

Simple beam models for Monte Carlo photon beam dose calculations in radiotherapy

Monte Carlo (code GEANT) produced 6 and 15 MV phase space (PS) data were used to define several simple photon beam models. For creating the PS data the energy of starting electrons hitting the target was tuned to get correct depth dose data compared to measurements. The modeling process used the full PS information within the geometrical boundaries of the beam including all scattered radiation of the accelerator head. Scattered radiation outside the boundaries was neglected. Photons and electrons were assumed to be radiated from point sources. Four different models were investigated which involved different ways to determine the energies and locations of beam particles in the output plane. Depth dose curves, profiles, and relative output factors were calculated with these models for six field sizes from 5x5 to 40x40cm2 and compared to measurements. Model 1 uses a photon energy spectrum independent of location in the PS plane and a constant photon fluence in this plane. Model 2 takes into account the spatial particle fluence distribution in the PS plane. A constant fluence is used again in model 3, but the photon energy spectrum depends upon the off axis position. Model 4, finally uses the spatial particle fluence distribution and off axis dependent photon energy spectra in the PS plane. Depth dose curves and profiles for field sizes up to 10x10cm2 were not model sensitive. Good agreement between measured and calculated depth dose curves and profiles for all field sizes was reached for model 4. However, increasing deviations were found for increasing field sizes for models 1-3. Large deviations resulted for the profiles of models 2 and 3. This is due to the fact that these models overestimate and underestimate the energy fluence at large off axis distances. Relative output factors consistent with measurements resulted only for model 4.