Ken Muldrew

A model for the time dependent three-dimensional thermal distribution within iceballs surrounding multiple cryoprobes

A time dependent three-dimensional finite difference model of iceball formation about multiple cryoprobes has been developed and compared to experimental data. Realistic three-dimensional probe geometry is specified and the number of cryoprobes, the cryoprobe cooling rates, and the locations of the probes are arbitrary inputs by the user. The simulation accounts for observed longitudinal thermal gradients along the cryoprobe tips. Thermal histories for several points around commercially available cryoprobes have been predicted within experimental error for one, three, and five probe configurations. The simulation can be used to generate isotherms within the iceball at arbitrary times. Volumes enclosed by the iceball and any isotherms may also be computed to give the ablative ratio, a measure of the iceballs killing efficiency. This ratio was calculated as the volume enclosed by a critical isotherm divided by the total volume of the iceball for assumed critical temperatures of -20 and -40 degrees C. The ablative ratio for a single probe is a continuously decreasing function of time but when multiple probe configurations are used the ablative ratio increases to a maximum and then essentially plateaus. Maximum values of 0.44 and 0.55 were observed for three and five probe configurations, respectively, with an assumed critical temperature of -20 degrees C. Assuming a critical temperature of -40 degrees C, maximum ablative ratios of 0.21 and 0.3 for three and five probe configurations, respectively, were observed.

Freeze-thaw treatment effects on the dynamic mechanical properties of articular cartilage

BackgroundAs a relatively non-regenerative tissue, articular cartilage has been targeted for cryopreservation as a method of mitigating a lack of donor tissue availability for transplant surgeries. In addition, subzero storage of articular cartilage has long been used in biomedical studies using various storage temperatures. The current investigation studies the potential for freeze-thaw to affect the mechanical properties of articular cartilage through direct comparison of various subzero storage temperatures.MethodsBoth subzero storage temperature as well as freezing rate were compared using control samples (4°C) and samples stored at either -20°C or -80°C as well as samples first snap frozen in liquid nitrogen (-196°C) prior to storage at -80°C. All samples were thawed at 37.5°C to testing temperature (22°C). Complex stiffness and hysteresis characterized load resistance and damping properties using a non-destructive, low force magnitude, dynamic indentation protocol spanning a broad loading rate range to identify the dynamic viscoelastic properties of cartilage.ResultsStiffness levels remained unchanged with exposure to the various subzero temperatures. Hysteresis increased in samples snap frozen at -196°C and stored at -80°C, though remained unchanged with exposure to the other storage temperatures.ConclusionsMechanical changes shown are likely due to ice lens creation, where frost heave effects may have caused collagen damage. That storage to -20°C and -80°C did not alter the mechanical properties of articular cartilage shows that when combined with a rapid thawing protocol to 37.5°C, the tissue may successfully be stored at subzero temperatures.

A semi-empirical treatment planning model for optimization of multiprobe cryosurgery

A model is presented for treatment planning of multiprobe cryosurgery. In this model a thermal simulation algorithm is used to generate temperature distribution from cryoprobes, visualize isotherms in the anatomical region of interest (ROI) and provide tools to assist estimation of the amount of freezing damage to the target and surrounding normal structures. Calculations may be performed for any given freezing time for the selected set of operation parameters. The thermal simulation is based on solving the transient heat conduction equation using finite element methods for a multiprobe geometry. As an example, a semi-empirical optimization of 2D placement of six cryoprobes and their thermal protocol for the first freeze cycle is presented. The effectiveness of the optimized treatment protocol was estimated by generating temperature-volume histograms and calculating the objective function for the anatomy of interest. Two phantom experiments were performed to verify isotherm locations predicted by calculations. A comparison of the predicted 0 °C isotherm with the actual iceball boundary imaged by x-ray CT demonstrated a spatial agreement within ±2 mm.

Cryopreservation of intact human articular cartilage

Damaged articular cartilage (AC) impairs joint function and many treatment techniques are being investigated to determine their long term results. Successful cryopreservation of AC can provide a reliable source of intact matrix with viable chondrocytes to maintain the cartilage over long periods of time. This study investigated the application of an established cryopreservation protocol to determine the recovery of intact chondrocytes from human AC. Ten millimeter diameter osteochondral dowels were harvested from two human donors. The cryopreservation protocol was performed and the samples were rapidly warmed from varying experimental holding temperatures (−10, −20, −30, −40°C), with and without plunging into liquid nitrogen, using 1 M dimethyl sulfoxide as cryoprotectant. The cartilage was stained with membrane integrity dyes and viewed under fluorescence microscopy. The percent of intact chondrocytes was compared to fresh controls. Low recovery of intact chondrocytes was recorded from all temperature levels with and without cryoprotectant. The results of this experiment demonstrated that the cryopreservation procedure used to achieve moderate success with intact sheep AC was not successful with intact human AC and further investigation is required.

Simultaneous optimization of cryoprobe placement and thermal protocol for cryosurgery.

We demonstrate that it is possible to simultaneously optimize multiple cryoprobe placements and their thermal protocol for one freeze-thaw cycle. A numerical optimization algorithm is used and three different forms of objective function are examined in terms of algorithm convergence rate, minimum value of the chosen objective function, temperature-volume histograms and isotherm distributions. The optimization results depend on the initial values of the variables, the form of the objective function, optimization goals and the mathematical method adopted for gradient calculation. The proposed optimization model offers significant advantages over the previously reported semi-empirical approach to conformal cryotherapy, such as the ability to handle an unlimited number of variables and eliminating the need for the user input between iterations, thereby reducing, if not removing, the subjectivity of cryosurgery treatment planning.

Suppression of high-density artefacts in x-ray CT images using temporal digital subtraction with application to cryotherapy

Image guidance in cryotherapy is usually performed using ultrasound. Although not currently in routine clinical use, x-ray CT imaging is an alternative means of guidance that can display the full 3D structure of the iceball, including frozen and unfrozen regions. However, the quality of x-ray CT images is compromised by the presence of high-density streak artefacts. To suppress these artefacts we applied temporal digital subtraction (TDS). This TDS method has the added advantage of improving the grey-scale contrast between frozen and unfrozen tissue in the CT images. Two sets of CT images were taken of a phantom material, cryoprobes and a urethral warmer (UW) before and during the cryoprobe freeze cycle. The high-density artefacts persisted in both image sets. TDS was performed on these two image sets using the corresponding mask image of unfrozen material and the same geometrical configuration of the cryoprobes and the UW. The resultant difference image had a significantly reduced artefact content. Thus TDS can be used to significantly suppress or eliminate high-density CT streak artefacts without reducing the metallic content of the cryoprobes. In vivo study needs to be conducted to establish the utility of this TDS procedure for CT assisted prostate or liver cryotherapy. Applying TDS in x-ray CT guided cryotherapy will facilitate estimation of the number and location of all frozen and unfrozen regions, potentially making cryotherapy safer and less operator dependent.

The possible influence of osmotic poration on cell membrane water permeability

It has been hypothesized that pores in the plasma membrane form under conditions of rapid water efflux, allowing extracellular ice to grow into the cytoplasm under conditions of rapid freezing. When cells with intracellular ice are thawed slowly, the transmembrane ice crystal expands through recrystallization causing the cell to lyse. One of the implications of this hypothesis is that osmotic pores will provide an alternative route for water movement under conditions of osmotically induced flow. We show that the plasma membrane water permeability of a fibroblast cell changes as a function of the osmotic pressure gradient that is used to drive water movement. It is further shown that cell volume is more important than the magnitude of water flux in causing this departure from a uniform water permeability. We suggest that these data provide evidence of a transient route for water movement across cell membranes.

Osteoarthritis as an inevitable consequence of the structure of articular cartilage

The hypothesis that osteoarthritis is an inevitable consequence of the structure of articular cartilage is proposed. Cartilage structure is viewed as an evolutionary solution to the optimization of articular cartilage function. A simplified model of the structural elements and constraints will be used to illustrate how the hypothesis leads to a mathematical solution and then a more realistic situation will be described, but not optimized rigorously. The goal for finding a solution will be to optimize the function of articular cartilage subject to the constraints of a biological tissue. In comparison, the structure of genuine articular cartilage will be discussed in the context of an optimal solution to the functional problem and, particularly, in relation to the disease of osteoarthritis.

Mechanisms of intracellular ice formation

The phenomenon of intracellular freezing in cells was investigated by designing experiments with cultured mouse fibroblasts on a cryomicroscope to critically assess the current hypotheses describing the genesis of intracellular ice: (a) intracellular freezing is a result of critical undercooling; (b) the cytoplasm is nucleated through aqueous pores in the plasma membrane; and (c) intracellular freezing is a result of membrane damage caused by electrical transients at the ice interface. The experimental data did not support any of these theories, but was consistent with the hypothesis that the plasma membrane is damaged at a critical gradient in osmotic pressure across the membrane, and intracellular freezing occurs as a result of this damage. An implication of this hypothesis is that mathematical models can be used to design protocols to avoid damaging gradients in osmotic pressure, allowing new approaches to the preservation of cells, tissues, and organs by rapid cooling.

Suppression of high-density artifacts in x-ray CT images using temporal digital subtraction with application to cryotherapy

Image guidance in cryotherapy is usually performed using ultrasound. Although not currently in routine clinical use, x-ray CT imaging is an alternative means of guidance that can display the full 3D structure of the iceball, including frozen and unfrozen regions. However, the quality of x-ray CT images is compromised by the presence of high-density streak artefacts. To suppress these artefacts we applied temporal digital subtraction (TDS). This TDS method has the added advantage of improving the grey scale contrast between frozen and unfrozen tissue in the CT images. Two sets of CT images were taken of a phantom material, cryoprobes and a urethral warmer (UW) before and during the cryoprobe freeze cycle. The high density artefacts persisted in both image sets. TDS was performed on these two image sets using the corresponding mask image of unfrozen material and the same geometrical configuration of the cryoprobes and the UW. The resultant difference image had a significantly reduced artefact content. Thus TDS can be used to significantly suppress or eliminate high-density CT streak artefacts without reducing the metallic content of the cryoprobes. In vivo study needs to be conducted to establish the utility of this TDS procedure for CT assisted prostate or liver cryotherapy. Applying TDS in x-ray CT guided cryotherapy will facilitate estimation of the number and location of all frozen and unfrozen regions, potentially making cryotherapy safer and less operator dependent.