Thermostability of Biological Systems: Fundamentals, Challenges, and Quantification

Xiaoming He*
Multiscale Biothermostability Engineering Laboratory, Department of Mechanical Engineering and Biomedical Engineering Program, University of South Carolina, 300 Main Street, Columbia, SC 29208, USA

Article Metrics

CrossRef Citations:
Total Statistics:

Full-Text HTML Views: 2060
Abstract HTML Views: 356
PDF Downloads: 209
Total Views/Downloads: 2625
Unique Statistics:

Full-Text HTML Views: 444
Abstract HTML Views: 199
PDF Downloads: 148
Total Views/Downloads: 791

© Xiaoming He; Licensee Bentham Open.

open-access license: This is an open access article distributed under the terms of the Creative Commons Attribution Non-Commercial License ( which permits unrestricted, non-commercial use, distribution and reproduction in any medium, provided the work is properly cited.

* Address correspondence to this author at the Multiscale Biothermostability Engineering Laboratory, Department of Mechanical Engineering and Biomedical Engineering Program, University of South Carolina, 300 Main Street, Columbia, SC 29208, USA; Tel: 803-777-0070; Fax: 803-777-0106; E-mail:


This review examines the fundamentals and challenges in engineering/understanding the thermostability of biological systems over a wide temperature range (from the cryogenic to hyperthermic regimen). Applications of the bio-thermostability engineering to either destroy unwanted or stabilize useful biologicals for the treatment of diseases in modern medicine are first introduced. Studies on the biological responses to cryogenic and hyperthermic temperatures for the various applications are reviewed to understand the mechanism of thermal (both cryo and hyperthermic) injury and its quantification at the molecular, cellular and tissue/organ levels. Methods for quantifying the thermophysical processes of the various applications are then summarized accounting for the effect of blood perfusion, metabolism, water transport across cell plasma membrane, and phase transition (both equilibrium and non-equilibrium such as ice formation and glass transition) of water. The review concludes with a summary of the status quo and future perspectives in engineering the thermostability of biological systems.

Key Words: Thermotherapy, thermoablation, cryotherapy, cryoablation, cryosurgery, hyperthermia, cryopreservation, vitrification, desiccation, lyopreservation, freeze-drying, membrane, vascular injury, stasis, cell injury, protein denaturation, phase transition, intracellular ice formation (IIF), water transport, glass transition, blood perfusion, kinetics, nanotechnology..


Efforts on engineering the thermostability of biological systems could be traced back to thousands of years ago according to the records of early civilizations. For example, high temperatures have been used to treat a variety of diseases including back pain in the traditional Chinese medicine [1-3]; the ancient Hindus heated metal bars to stop bleeding [4]; in 500 B.C., the Greek physician Parmenides believed that if he could create fever, he could cure all illness [5]; dried mummies have been found in Egypt and other places to be well-preserved at ambient temperatures for thousands of years [6]; and the use of cold was recommended by Hippocrates to reduce bleeding and swelling [7]. Nothing speaks for the unlimited potential of temperature in medicine better than an aphorism of Hippocrates [8,9]: “Those diseases which medicines do not cure, iron (author’s note: iron means scalpel here) cures; those which iron cannot cure, fire (author’s note: fire means temperature here) cures; and those which fire cannot cure, are to be reckoned wholly incurable.”

Indeed, temperature (either heat or cold) is still one of the most important tools for the treatment of diseases in modern medicine. Recent research on engineering the (thermo) stability of biologicals at various temperatures for biomedical applications has been focused on either destroying unwanted/diseased cells/tissues or stabilizing cells/tissues for future use in the treatment of diseases as summarized in Fig. (1). The use of abnormally low and high temperatures to locally destroy unwanted tissues is called cryosurgery/cryoablation/cryotherapy and thermal surgery/thermal ablation/thermal therapy/thermotherapy (including hyperthermia for which the temperature is usually between 39 and 45 oC), respectively. For thermal stabilization (or biopreservation) applications, different terminologies including cryopreservation, hypothermic preservation and dry or lyopreservation are used for temperatures below -80 oC, at ~ 4 oC, and at ambient temperatures, respectively. There is a significant overlap in terms of temperature range between hypothermic and lyopreservation. However, they are quite different since hypothermic preservation is designed for preservation of fully hydrated tissues/organs in a short period usually from days to weeks while lyopreservation is aimed at stabilizing biologicals in the long term (i.e., up to years) in a dehydrated state. More detailed discussions of the different applications are given below.

Fig. (1).

A summary of typical biomedical applications for which the goal is to engineer the thermostability (for either destruction or stabilization) of biological systems from cryogenic to hyperthermic temperatures.

Fig. (2).

Sketches of a typical cryoprobe (A) designed based on the Joule-Thomson effect for destroying tissue by freezing and a typical thermal probe (B) using energy carried either in an electromagnetic (e.g., radio frequency, microwave, or laser) or acoustic wave field (HIFU, high intensity focused ultrasound) to destroy tissue by heating.

Fig. (3).

An illustration on the extended phase diagram of the various approaches to achieve thermal stabilization on biologicals at either cryogenic or ambient temperature from an initial liquid state (A) to a final glassy phase (e.g., I, II, III, or IV): The phase diagram is divided into four thermodynamic regimens by the liquidus, extended liquidus, solidus, and the glass transition curve; the four regimens are the liquid, subcooled liquid, supersaturated liquid, and the glassy phase; cells (in a liquid solution initially) must enter the glassy phase for long-term preservation; CPA represents cryoprotectant, and θg and θm represent the glass transition and melting temperature, respectively. Of note, the diagram is not to scale (for example, the melting temperature of pure CPAs is usually below 20 oC). Figure reprinted from reference [53] with permission from Xiaoming He (the author).

Fig. (4).

Tardigrade (A and B) and resurrection plant (C and D) can survive extreme drought in nature upon rehydration - a phenomenon called anhydrobiosis or life without water. Figure reprinted and redrawn from references [85,86] with permission from John H. Crowe at U of California-Davis (for panels A and B) and Wayne P. Armstrong at Palomar College (for panels C and D).

Fig. (5).

An illustration (not to scale) of some important structures/confirmations of a typical biological tissue (A), mammalian cell (B), cell plasma membrane (C), and protein (D): A biological tissue typically consists of cells (1), extracellular matrix (2) and blood vessels (3) through which oxygen and nutrients are transported to the cells; a mammalian cell consists of many important sub-cellular organelles such as the rough endoplasmic reticulum (4), ribosome (5), smooth endoplasmic reticulum (6), mitochondria (7), plasma membrane (8), cytosol (9), lysosome (10), centriole (11), Golgi apparatus (12), and nucleus that are demonstrated as nuclear membrane (13), nucleolus (14), and nuclear membrane pore (15); the cell plasma membrane contains a phospholipid bilayer with hydrophilic heads (16) and hydrophobic tails (17), and membrane proteins (18) that works as transport channels, receptors, and structural connections between intracellular skeleton and extracellular matrix; and the most important secondary structures of a protein are α helix (19) and β sheet (20).

Fig. (6).

Fourier transform infrared (FTIR) and differential scanning calorimetry (DSC) signatures of thermally induced changes in cellular proteins in situ: (A), FTIR spectra at different temperatures showing the increase of extended β sheet structure and a simultaneous decrease of the α helix structure in cellular protein during heating mammalian cells slowly (2 °C/min); (B), DSC endotherm demonstrating heat absorption as a result of protein denaturation during heating mammalian cells slowly also at 2 °C/min; (C), FTIR spectra showing extensive change of protein secondary structure from α helix to extended β sheet in thawed cells after freezing slowly to -80 °C, but not -20 °C; and (D), the DSC and FTIR measurements match well, suggesting protein denaturation is one of the major events that result in cell injury during heating between 40-70 °C. Figure reprinted and redrawn from references [158,195] with permission from Elsevier (for panels A and C) and Springer (for panels B and D).

Fig. (7).

Viability and function of mammalian cells post thermal treatment evaluated by exemplary assays: (A), immediate viability of cells assayed using calcein AM (green fluorescence indicating cellular metabolic activity and intact membrane of viable cells) and ethidium homodimer (red fluorescence indicating compromised cell membrane of injured and potentially dead cells); (B) and (C), expression of green fluorescent protein (GFP) under the control of Oct-4 gene (B) and a membrane surface glycoprotein (SSEA-1, C) indicating the undifferentiated properties of R1 murine embryonic stem cells post thermal treatment; and (D) histological difference between viable and damaged renal cells in normal porcine kidney tissue after 2 day’s culture in media post thermal treatment showing the transition (arrow) from the lower left zone of significant necrotic tissue to the upper right zone with mainly intact tissue: Damaged cells lack nuclei. Figure reprinted and redrawn from references [38,64] with permission from Elsevier (for panels B and C) and Informa (for panel D).

Fig. (8).

Both histology and blood perfusion defects have been used to evaluate tissue injury in vivo: (A), a micrograph showing the patterns of in vivo hyperthermic injury in porcine renal tissue 7 days after thermal treatment and (B) in vivo blood perfusion defects due to cryothermic injury in prostate tissue at different times (i.e. 0, 3 and 7 days) measured in situ by fluorescence contrast. Although wound healing was observed in the necrotic and partially ablated zone, a central lesion named thermal fixation knowing for its resistance to wound healing response were observed in kidney tissue exposed to hyperthermic temperatures (A). Wound healing was observed throughout the cryogenic lesion leading to its shrinkage from day 3 to day 7 (B). Figure reprinted and redrawn from references [38,247] with permission from Informa (for panel A) and ASME (for panel B).

Fig. (9).

Cell survival as a function of cooling rate due to intracellular ice formation (IIF) only (➀), freeze concentration induced cell dehydration (solute effect) only (➁), and the combination of the two factors (➂ and ➃): During slow-freezing, the solute effect and IIF determines cell survival at low and high cooling rates, respectively; the combined effect of IIF and freeze concentration (i.e., solute effect) results in the classical inverted U curve of cell survival vs. cooling rate and an optimal cooling rate (CRSF) is observable for slow-freezing; the optimal cooling rate is cell type and also cryoprotectant concentration dependent; and when the cooling rate is higher than a critical cooling rate (CRV), cells are vitrified without freezing (or ice formation) and high cell survival ensues. Figure reprinted from reference [53] with permission from Xiaoming He (the author).

Fig. (10).

An illustration of energy-state relationship of the 1st order kinetic process of thermally induced protein denaturation: The native protein is activated after absorbing energy followed by denaturation and potentially aggregation and agglomeration: Aggregation and agglomeration make the process irreversible. Figure reprinted and redrawn from reference [143] with permission from John Wiley & Sons.

Fig. (11).

An illustration of the procedures for determining the activation energy and frequency factor using Equations 4 and 5: (A), determining the injury rate at various temperatures by fitting Equation 5 to experimental data on cell survival vs. exposure time and (B) determining activation energy and frequency factor by fitting Equation 4 to the natural logarithm of injury rate vs. the inverse of temperature obtained from (A). Figure reprinted and redrawn from reference [187] with permission from ASME.

Fig. (12).

A summary of the kinetic parameters in the first order kinetic model for various proteins, cells, and tissues: A simple, phenomenological, linear relationship exists between activation energy and the natural logarithm of frequency factor in the first order kinetic model for hyperthermic injury. Figure reprinted and redrawn from reference [142] with permission from Begell House.

Fig. (13).

Dependence on cooling rate of intracellular ice formation (IIF) when freezing hepatocytes from -4 (after seeding extracellular ice) to -40 °C showing the darkening of the cellular interior as a result of significant IIF under a bright field of light microscope: When the cooling rate is very low (10 °C/min), cell dehydration dominates (A); when the cooling rate is high (400 °C/min), IIF dominates (C); and when the cooling rate is intermediate (100 °C/min), both IIF and cell dehydration can occur (B). Figure reprinted and redrawn from reference [272] with permission from John Wiley & Sons.

Fig. (14).

A typical SAR pattern of a dipole microwave probe showing the penetration of microwave energy into the tissue space (in r direction). The frequency of the microwave used is 921 MHz. Figure reprinted and redrawn from reference [38] with permission from Informa.

Fig. (15).

Change of normalized blood perfusion (w/w37C) in normal porcine kidney tissue during heating as a function of vascular stasis predicted using the 1st order kinetic model: Blood perfusion increases initially when the predicted vascular stasis is less than ~ 10% followed by a gradual decrease and eventually complete stasis, indicating the first order model is insufficient in predicting the change of blood perfusion in response to heating in normal tissue. It is reported that the initial increase in blood perfusion from the baseline value is absent in heating tumor tissue in vivo. Figure reprinted and redrawn from reference [38] with permission from Informa.

1.1. Thermal Destruction

On the thermal destruction side, whole body hyperthermia (usually with brain cooling to avoid thermal damage to the extremely heat susceptible cerebral tissue [10,11]) with temperatures usually between 39-45 oC has been used for many years as an adjuvant treatment to enhance the effect of chemo and radio therapy of malignant diseases (i.e., cancer) [12-16]. More recently, thermal therapy and cryosurgery have been developed as single treatment modalities for various benign and malignant diseases [17-22]. With the advancement of modern surgical techniques such as laparoscopy, both thermal therapy and cryosurgery can be performed in a minimally invasive manner to reduce the patient suffering and hospital stay during and post surgical operation [23-25]. Consequently, they are becoming increasingly popular and have been studied in essentially all surgical sub-specialties as a minimally invasive alternative to the conventional radical surgical intervention [26-28].

1.1.1. Thermal Destruction at Cryogenic Temperatures (Cryotherapy): Cooling Approaches

In cryosurgery, unwanted or diseased tissue is destroyed by freezing it to subzero temperatures (usually below -20 oC) typically using one or multiple cryoprobes (up to a few millimeters in diameter, Fig. 2A). The tissue is frozen through direct contact with the cryoprobe which is cooled to a cryogenic temperature by circulating cryogenic media at a temperature usually below -100 oC inside the cryoprobe lumen [29,30]. Liquid nitrogen (-196 oC) is one of the commonly used cryogenic media, particularly in cryoprobes designed before the 1990s. Gases such as argon, carbon oxide, and nitrous oxide at cryogenic temperatures obtained by throttling (i.e., from high to low pressure) as a result of the Joule-Thomson effect [31] are becoming increasingly popular. Liquid nitrogen has been replaced with these gases in more recent designs of cryoprobes for better control of the probe temperature. An additional advantage of cryoprobes designed based on the Joule-Thomson effect is that the circulating gas for freezing in the cryoprobe can be switched to a different gas such as helium, whose temperature actually increases after throttling, due to a negative Joule-Thomson coefficient (i.e., µJT defined in Fig. 2A) of the gas at around room temperature. This coefficient is positive for argon at room temperature or below. This feature is utilized to achieve fast, active (instead of slow, passive) thawing of the diseased tissue after freezing, which is important for accurate control of the frozen tissue volume. Other gas mixtures (mainly Freon) at cryogenic temperatures (usually higher than -130 oC) obtained utilizing a refrigeration cycle system essentially like that of a household refrigerator/freezer, are also used in some cryoprobe designs.

1.1.2. Thermal Destruction at Hyperthermic Temperatures (Thermotherapy): Heating Approaches

In thermal therapy, unwanted or diseased tissue is thermally ablated or destroyed by locally heating the tissue to above 50 oC for minutes or even seconds typically using one (for removing a small unwanted tissue volume) or an array (for removing a large unwanted tissue volume) of thermal probes (also up to a few millimeters in diameter, Fig. 2B). Unlike the cryoprobes for cryosurgery that must have a cold surface to freeze tissue, a thermal probe for thermal therapy is not necessary to be hot and even not necessary to be in direct contact with the tissue, particularly when laser or high intensity focused ultrasound (HIFU) is used as the energy source [32-37]. This is because it is the electromagnetic or acoustic wave field applied on the tissue through the thermal probe that heats up and thus destroys the tissue. Moreover, a cooling mechanism might be employed to protect the thermal probe from being overheated when the probe does have direct contact with the tissue to be thermally destroyed [38,39]. This is because the efficacy of the probe for emitting electromagnetic or acoustic wave might be significantly compromised at the high temperature necessary for destroying unwanted tissue.

The electromagnetic wave over a wide range of frequency (or wavelength) including radiofrequency (RF), microwave, and laser, has been used for heating. The terms of radiofrequency (RF) and microwave have often been used interchangeably in the hyperthermic literature for electromagnetic waves with frequency over a broad range from several hundred kilohertz (kHz) to several hundred gigahertz (GHz). Traditionally, RF is referred as electromagnetic wave with frequency from 300 kHz to 1GHz and microwave is from 300 MHz to 300 GHz. The corresponding wavelength of RF and microwave in air is ~ 1000 – 0.3 m and ~ 1000 – 1 mm, respectively. Considering the speed of light in water (the major building component of most soft tissue) is approximately half of that in air (3×108 m/s), the wavelength of RF/microwave in most soft tissue should decrease by half as well. Typical RF/microwave frequencies used for heating biological tissue in the hyperthermic literature are around 2.45 GHz, 915 MHz, 40.68 MHz, 27.12 MHz, and 13.56 MHz, which are the frequencies allocated by the Federal Communications Commission (FCC) for medical use in USA. The two highest frequencies (2.45 GHz and 915 MHz) are usually (although not always) referred as microwave while the rest are called RF in the thermal therapy literature.

Both conductive (more important at low frequencies) and dielectric (more important at high frequencies) heating mechanisms contribute to the absorption of electromagnetic energy by biological tissue [40-42]. The former is due to the electrical resistance of biological tissue and is the major heat generation mechanism at low frequencies (e.g., < 1 MHz). The latter is due to friction between molecules in tissue as a result of the alternating movement of polar molecules (particularly the water dipoles in a biological tissue) trying to align themselves with the alternating RF/microwave electric field. Therefore, dielectric heating is frequency dependent and is the major heating mechanism at high frequencies (e.g., > 100 MHz), particularly in the nonconductive biological tissue. A more comprehensive review of the biological effect of RF/microwave can be found elsewhere [40].

Lasers (acronym for light amplification by stimulated emission of radiation) of infrared, visible light, and ultraviolet with a wave length much less than 1 mm (typically ~ 150 nm – 100 µm in the infrared domain) have also been used for thermal destruction. Correspondingly, the wave frequency is much higher than 300 GHz but usually less than 1016 Hz. In this case, the tissue is heated by absorbing photons carried in the laser wave field [32]. Due to the small wave length, laser wave could be scattered significantly along its propagation path in biological tissue which is generally inhomogeneous. Light scattering can significantly enhance tissue absorption of the laser energy on one hand and limit the penetration depth of the laser wave in tissue on the other. Generally, the heating depth of a single interstitial laser probe is only approximately one tenth of that of an interstitial RF/microwave probe under comparable working conditions. Therefore, multiple probes or multiple laser shots at different locations in the targeted tissue are necessary to destroy a large tissue volume. It is interesting and important to note that lasers with wavelength from 800 to 1,200 nm can be absorbed much more strongly by metals than water [33]. Therefore, lasers in this wavelength range can be utilized to specifically destroy unwanted tissue by heating metallic (e.g. gold) nanoparticles delivered into the tissue [43,44].

The acoustic (or mechanical) wave used for heating and destroying unwanted tissue is generally referred as high intensity focused ultrasound (HIFU). The frequency of medical acoustic waves (a mechanical wave) is typically from 1 to 3 MHz. As a result, the wavelength of the HIFU mechanical wave is typically less than ~ 1.5 mm since the speed of sound in water is ~ 1500 m/s between 0-100 oC. Therefore, tissue scattering of the mechanical wave is significant, which can enhance its absorption by tissue and at the same time limit its penetration depth in tissue. The latter is due to undesired heating of the surface tissue such as the dermis. Consequently, multiple HIFU shots at different locations in the targeted tissue are necessary to destroy a large diseased tissue volume. At a low intensity, tissue absorption of acoustic energy can be attributed to interstitial mechanical friction as a result of the alternating movement in tissue driven by the alternating acoustic (mechanical) wave while at a very high intensity, cavitation will occur as a result of microbubble formation/collapse and its interaction with the acoustic wave [35]. A significant enhancement of heat generation has been observed when cavitation occurs [36,37]. In addition, the generation of shock wave associated with microbubble collapse during cavitation can destroy tissue mechanically [35].

1.2. Thermal Stabilization

The goal of thermal stabilization (or biopreservation) is to preserve important biologicals such as biomolecules, cells, tissues, and organs for future use (e.g., transplantation) in modern medicine. For short term stabilization, biologicals are usually cooled to around ice temperature (hypothermic preservation) without ice formation or freezing. The purpose is to slow down metabolism and degradation of the biologicals to extend their shelf life to days and up to weeks after explantation [45-47]. Studies on short term hypothermic preservation are mainly focused on organs, for which long-term stabilization has not been practically successful. For long-term stabilization, however, biologicals (e.g., proteins, cells and small tissue slices) are usually processed and stabilized in an amorphous or glassy phase, which can be best demonstrated in an extended phase diagram shown in Fig. (3). The glassy phase is a thermodynamically metastable state with an extremely high viscosity and low molecular mobility and activity. Therefore, any processes that require the diffusion of participating molecules including both degradative and life supporting processes in a biological system are essentially arrested/suspended in a desired glassy phase. Practically, five different approaches including the conventional slow-freezing, conventional vitrification using a high concentration of cryoprotectant (or cryoprotective agent, or CPA for short), low-CPA vitrification, freeze-drying, and evaporative/convective drying (or desiccation) have been studied to bring biologicals from an aqueous liquid state at room temperature to a glassy phase at either a cryogenic or ambient temperature as shown in Fig. (3). The former is usually called cryopreservation while the latter is usually called dry or lyopreservation. A more detailed discussion of each of the approaches for long-term cell/tissue preservation is given in the following sub-sections.

1.2.1. Thermal Stabilization at Cryogenic Temperature (Cryopreservation)

Thermal stabilization of biologicals at a cryogenic temperature can be achieved with two major approaches: the conventional slow-freezing with inevitable ice formation and the ice-free vitrification [48-53]. For the conventional slow-freezing approach, the following steps are typically adopted (A→B→C→D→E→I in Fig. 3): (1), biologicals such as cells in an aqueous solution (state A) are first loaded with CPA(s) such as DMSO (dimethylsulfoxide), glycerol, propylene glycol (or 1,2-propanediol), and ethylene glycol at a concentration usually up to 15 wt% (or up to 2 M, state B); (2), the samples are sub-cooled usually to between -2 and -7 oC (phase C: note, phase rather than state is used to indicate that it is not in equilibrium) to seed ice in the extracellular space by touching the sample contained usually in a cryovial with a deeply cooled (e.g., in liquid nitrogen) object (C→D); (3), the samples are further cooled slowly (typically, < 10 oC/min) along the liquidus to between -40 and -100 oC, a process called freeze concentration (D→E); and (4), the samples are transferred into liquid nitrogen for long-term storage (E→I). In this approach, the formation of extracellular ice leads to freeze concentration of the unfrozen solutions by ejecting solutes and cells from the frozen to unfrozen phases. As a result, dehydration of cells in the unfrozen phase ensues, which minimizes intracellular water available for ice formation inside the cells so that the cells can enter the glassy phase (I) easily when transferring into liquid nitrogen. This approach typically requires a specialized machine usually called controlled rate freezer (CRF) to achieve freezing in a controllable manner. The time required for the slow-freezing process is typically hours.

Vitrification by definition is ice free. In other words, no (or negligible) ice formation or freezing will occur in the sample during cooling [48,51,54]. Conventional vitrification (A→F→II in Fig. 3) has also been studied for both cells and tissue. In this approach, biological samples (state A) are first loaded with a very high concentration of CPA(s) (up to ~ 7 M, state F) [48,51,54]. The samples are then cooled directly from ambient temperature to a cryogenic temperature usually in liquid nitrogen (state II) and stored there for future use. Although the conventional vitrification approach can be used to eliminate the detrimental effect of ice formation altogether, the unusually high concentration of CPA required by the approach is toxic to most mammalian cells even in a short period of exposure (ranging from seconds to minutes dependent on the specific cells and tissues) [55-60]. Therefore, the samples should be cooled as soon as possible after loading with CPAs. A mixture of multiple CPAs is often used to reduce the cytotoxicity of the high CPA concentration required [61]. In addition, large, membrane impermeable molecules such as sugars (typically sucrose and trehalose) have been used to minimize ice formation and protect cell membrane from injury during cooling [62-64]. Vitrification can be done without a specialized machine and the time required is generally much shorter than that for slow-freezing.

Low-CPA vitrification (A→B→III in Fig. 3) is a further advancement of the conventional vitrification with the goal to reduce the CPA concentration (e.g., at state F vs. B) required for vitrification to a low, nontoxic level (similar to that used for slow-freezing). This can be done by creating an ultrafast cooling rate to cool the cells for cryopreservation. This is because the higher the cooling rate, the less the amount of cryoprotectants is required for achieving vitrification [64-69]. For example, even pure water can be vitrified without any ice formation when the cooling rate is about one million degree Celsius per second [69-71].

Various devices have been utilized to achieve fast cooling rates (~ 20,000 oC/min) such as the traditional French type straw, open pulled straw, electron grid, and cryoloops [56,72-74]. As a result of the fast cooling rate, the amount of cryoprotectants required for vitrification can be reduced to around 4 M. To achieve an ultrafast cooling rate, two recent studies reported the use of a micro-fabricated oscillating heat pipe (OHP) device [75,76]. Although their theoretical analysis shows that an ultrafast cooling rate of ~ 106oC/min could be achieved, testing of the device for low-CPA vitrification using living cells has not been reported to date. Another recent study reported that a cooling rate as high as 200,000 oC/min can be achieved by plunging an ultra-thin walled (10 µm) quartz microcapillary (QMC, 180 µm inner diameter, slightly bigger than the diameter of a human oocyte) into liquid nitrogen [64]. As a result, the CPA concentration required for vitrification of mouse embryonic stem cells and mouse oocytes was found to be as low as 2.5 M altogether [64,77], which is close to the upper boundary of CPA concentrations used for slow-freezing. More recently, it was found that the required CPA concentration could be further reduced to ~ 1.5 M when mouse mesenchymal stem cells were encapsulated in ~ 100 µm alginate microcapsules for vitrifying using the thin-walled (380 µm inner diameter) QMC [78]. Therefore, the QMC-based low-CPA vitrification technique is promising to provide a solution to the major challenge facing the conventional vitrification approach associated with cytotoxicity of high CPA concentration. Because low-CPA vitrification avoids the shortcomings of both the conventional slow-freezing and vitrification approaches while combining their advantages, it is becoming increasingly popular, particularly for the cryopreservation of important and osmotically sensitive mammalian cells such as oocytes, sperm, stem cells, and many primary cells important for tissue engineering, cell-based therapy, and assisted reproduction.

1.2.2. Thermal Stabilization at Ambient Temperatures (Lyopreservation)

With the development of modern cell-based medicine, the demand on wide distribution of biopreserved products to end users in medical units at both urban and remote locations is increasingly high. The difficulty to widely distribute cryopreserved biologicals in liquid nitrogen has motivated cryobiologist and engineers to develop approaches for long-term stabilization of biologicals at ambient temperature. The idea of dry or lyopreservation at ambient temperature is actually not new since many lower organisms (e.g., tardigrade shown in Fig. 4A and B), resurrection plants (e.g., Selaginella Lepidophylla shown in Fig. 4C and D), and seeds can survive extreme drought in nature upon rehydration, a phenomenon called anhydrobiosis or life without water [79-86]. A high concentration of sugars (typically sucrose for plants and trehalose for lower organisms) has been found in these organisms and plants when they are in the anhydrobiotic state. Learning from nature, both sucrose and trehalose have been used as the protective agent (also termed lyoprotectant) in protocols for lyopreservation [87-91]. In practice, two approaches have been studied to bring biologicals in aqueous samples to a dry state: freeze-drying (or lyophilization) and evaporative or convective drying (or desiccation) as shown in Fig. (3), as well.

A typical freeze-drying protocol is as follows (A→B→C→D→G→H→IV in Fig. 3): (1), a biological sample such as proteins and cells in an aqueous solution (state A) is first supplemented with lyoprotectants (e.g. sucrose and trehalose) at a concentration of up to ~ 15 wt% (state B); (2), the sample is sub-cooled to usually between -3 and -7 oC (phase C) to seed ice in the solution by touching the samples with a deeply cooled object (C→D); (3), the sample is further cooled to between -30 and -50 oC slowly at a cooling rate usually less than 10 oC/min (D→G); (4), the ice formed in the sample during freezing is then sublimated by exposing the sample to a vacuum usually less than 10 Pa at phase G (primary drying); and (5) a secondary drying process is then done by heating the sample in vacuum slowly to ambient temperature to further dehydrate the sample for additional hours to days (G→H→IV). The samples are then sealed and preserved in the dry phase (IV) at ambient temperature for future use. Freeze-drying has been used successfully in achieving lyopreservation of many biomacromolecules such as proteins and lipids and many pharmaceutical drugs. It has also been used for achieving lyopreservation of prokaryotic cells such as bacterium, red blood cells, and platelets, but not eukaryotic cells at this time. For freeze-drying, it is crucial to keep the temperature low enough (below the so-called collapse temperature [92-106]) during primary drying (at phase G). Otherwise, the sample may collapse (can’t maintain the morphology of the frozen sample) resulting in incomplete drying and heterogeneity in the freeze-dried product, which could significantly decrease the biostability of the freeze-dried biologicals [107-109]. More importantly, a recent study reported that the collapse temperature of cell culture medium-based trehalose solutions important for freeze-drying mammalian cells can be much lower than that of a simple binary trehalose-water solution and trehalose solutions used for freeze-drying pharmaceuticals and prokaryotes [110].

Unlike freeze-drying, during desiccation by evaporative/convective drying (A→B→IV in Fig. 3) water in an aqueous sample is removed by exposing the sample to a dry environment (e.g., dry air, inert gas such as nitrogen, and vacuum) without freezing (or ice formation) after loading with up to 15 wt% lyoprotectants. Desiccation by evaporative/convective drying has been used to achieve lyopreservation of both biomacromolecules such as proteins and lipids, pharmaceutical drugs, and prokaryotic cells (such as bacterium, red blood cells, and platelets), but not eukaryotic mammalian cells. A major engineering difficulty to dry the glass-forming trehalose based solution for cell lyopreservation by convective drying is that a thin glassy skin can easily form on the interface between the solution and the dry environment, resulting in incomplete drying and heterogeneity in the dried product. This problem might be minimized by breaking down the solution into micron-sized droplets or thin-films [108].

Beside the engineering challenge to effectively dry the trehalose solutions, effective delivery of the small hydrophilic lyoprotectants (e.g., trehalose and sucrose) into mammalian cells has been challenging as the first step toward cell preservation at ambient temperature. This is because lyoprotectant such as trehalose must be present both intra and extracellularly to provide the maximum protection during drying, but mammalian cells lack a mechanism to synthesize trehalose endogenously and their plasma membrane is impermeable to the sugars [111-113]. Over the past decades, a number of approaches have been explored to introduce trehalose into living cells for preservation purpose. The most straightforward approach is to deliver exogenous trehalose into the cytosol of living cells by direct microinjection. This approach has been successfully used for intracellular delivery of trehalose to cryopreserve mammalian oocytes that have a large size (~ 100 µm in diameter) and are generally in a small quantity (less than a few hundred) [113-116].

However, the microinjection approach is difficult (if not impossible) to apply for most living cells that are generally much smaller (< ~ 20 µm) than mammalian oocytes and usually present in a large quantity (millions). Small living cells have been genetically engineered to synthesize trehalose endogenously. This approach requires the constant production of adenoviral vectors that exhibit significant cytotoxicity, particularly at high multiplicities of infection [117-119]. Trehalose has also been introduced into mammalian cells or their organelles through engineered or native transmembrane pores [91,112,120-122], electroporation [123,124], fluid-phase endocytosis [125-127], and lipid phase transition [127,128].

In spite of the various approaches being explored, a consistent report of cell preservation using trehalose for small eukaryotic living cells is still absent [88,111,129,130]. This could be due to the inability to deliver a sufficient amount of intracellular trehalose (≥ ~ 0.1 M) for cell preservation using some of the approaches (e.g., fluid phase endocytosis). In addition, cells could be too severely compromised during the delivery step to withstand further freezing/dehydration stresses during preservation, considering the highly invasive nature of some of the approaches (e.g., electroporation). Recently, research has been sought to use liposomal and polymeric nanoparticles as the intracellular delivery vehicles of small hydrophilic molecules including the lyoprotectants with promising outcomes [131-135].

Besides the non-reducing disaccharide (trehalose and sucrose), small stress proteins particularly, the late embryogenesis abundant (LEA) proteins have been suggested to be an important part of the molecular repertoire that renders desiccation tolerance in anhydrobiotic organisms and are attracting more and more research attention [79-84,136-141].


The biological responses to either abnormally low or high temperatures have been reported at the molecular, cellular and tissue (both in vitro and in vivo) levels [142-149]. Some important structures/conformations of important biologicals (e.g., tissue, cell, cell membrane, and protein) that have been reported to be the target of both cryo and hyperthermic injury are shown in Fig. (5). The driving forces leading to the biological responses at abnormally high temperature are different from those at abnormally low temperatures due to the different biophysical events involved. The biological responses at hyperthermic temperature appear to be dominantly thermotropic. At abnormally low temperatures, however, additional biophysical events such as ice formation and freeze concentration (see D→E→G in Fig. 3) as a result of ice formation at subzero temperatures can lead to additional biological alterations at the molecular, cell, and tissue levels.

2.1. Hyperthermic Injury

A bevy of studies have been performed to investigate the alterations in biologicals exposing to hyperthermic temperatures. At the molecular level, lipid, protein, DNA, and RNA are the major structural and functional macromolecules of living cells and tissues. DNA and some structured RNAs do not undergo conformational change until above ~ 85–90 °C [150-153]. Therefore, they are unlikely to play an important role in cell killing at temperatures ranging from 37 to 85 °C for most thermal therapy protocols. Other RNAs that can perform enzymatic functions and form protein/RNA complexes such as ribosomes [151,154] could be potential targets in thermal therapy applications [153]. Both protein denaturation and lipid alteration (phase transition) have been commonly observed after thermal treatment [143]. Although many studies on protein denaturation were done using purified proteins, the attention of recent studies have gradually shifted to study such changes in situ in living cells using various analytical techniques including differential scanning calorimetry (DSC), Fourier transform infrared (FTIR) spectroscopy, and circular dichroism (CD) spectroscopy [155-159]. These studies have shown that alterations to lipid mainly occur below ~ 45 oC while major protein denaturation starts from 40 to 45 oC (Fig. 6A and B) dependent on the heating rate and continue to occur at more than 100 oC. In view of the fact that major cell injury occurs at temperatures above 45 oC, protein denaturation has been proposed to play a more important role in defining cell injury than the change in lipids. This argument may be further strengthened by identifying the denaturation of specific proteins or protein groups that is the rate limiting step for thermal injury to cells [159]. However, lipid alterations are still considered to be important in defining the rate of hyperthermic injury because many cellular proteins (especially membrane proteins) are embedded in lipids and their stability is strongly affected by the lipids surrounding them [160-165].

Of note, with much attention being focused on studying the effect of hyperthermic temperatures on macromolecules, recent studies have found that a variety of molecules can be used to augment the performance of high temperatures for thermal destruction [166]. For example, the small peptide of tissue necrosis factor-alpha (TNF-α, 12.6 kD) has been shown to sensitize tumor for thermal therapy when delivered to the tumor site using functionalized gold nanoparticles [167]. Anticancer drugs encapsulated in liposomes have been shown to significantly increase the effect of hyperthermia in inhibiting tumor growth in vivo [168-170]. Using therapeutic adjuvants such as TNF-α and anticancer drugs to augment thermal therapy and the targeted delivery of these adjuvants specifically into tumor using various nanoparticles (gold, magnetic, liposomal, and polymeric nanoparticles and carbon nanotubes for heating, as well) are attracting more and more attention in the field of thermal therapy [43,170-176].

At the cellular level, thermal damage to essentially all the subcellular organelles (e.g., the plasma membrane, mitochondria, Golgi apparatus, and endoplasmic reticulum in Fig. 5B) have been reported in the hyperthermic literature, presumably due to the damage to their two important building blocks, the cellular proteins and lipids as discussed above [142]. The plasma membrane that provides cells with the first protection from its extracellular environment has been the major focus of many studies on hyperthermic injury. Hyperthermic injury to the plasma membrane usually is manifested as bleb formation and hyper-permeability to large molecules that are impermeable to the plasma membrane of intact cells (Fig. 7A). Damage to the plasma membrane and other sub-cellular organelles can further cause cells to lose their capability of attaching to a substrate and the subsequent proliferation and therefore, cell death ensues. Of note, when cells are exposed to a hyperthermic temperature below 43 oC, they can adapt themselves to the sub-lethal temperature and become more thermally resistant to subsequent heat treatments, a phenomenon called thermotolerance [14,177]. Cells obtain this capability by increasing the concentration of a group of constitutive stress proteins called the heat shock proteins inside the cells [178-184]. Heat shock proteins are known to protect critical cellular proteins for normal cellular function under stress conditions including that as a result of abnormal hyperthermic temperatures [180,181].

At the in vitro tissue level, living cells are embedded in their native extracellular matrix (mainly collagen, Fig. 5A), which is different from the substrate or culture medium for attached and suspended cells, respectively. The effect of different extracellular matrices on the thermal sensitivity of mammalian cells has been reported in a few studies [185-188]. The results however, are still inconsistent. Some studies show that cells are more thermally sensitive in their native matrix while others show no significant difference or even the opposite [185-188]. Further studies in a more controlled fashion are required to clarify the discrepancy in the existing literature in this regard.

At the in vivo tissue level, blood perfusion (Fig. 5A) and inflammatory wounding healing responses have been shown to play an important role in determining the extent of thermal injury: The thermal threshold required to destroy cells embedded in in vivo tissue is lower than that in vitro [38,167,189,190]. This observation has been attributed to vascular stasis as a result of thermal damage to tissue vasculature during and post thermal therapy, which leads to the deprivation of nutrition and oxygen, a secondary mechanism that kills the cells in addition to direct cell injury. Inflammatory wound healing response has been reported to be limited to the peripheral (i.e., the partial ablation and necrotic zone in Fig. 8A) of the thermal lesion in vivo even after more than one week [38,191]. The lack of wound healing response in the central zone (i.e., the thermal fixation zone in Fig. 8A) of a hyperthermic lesion in vivo is presumably due to the lack of blood perfusion as a result of vascular stasis post thermal therapy in the observed time period. Sufficient blood perfusion is necessary for recruiting inflammatory cells such as neutrophils and macrophages for the wound healing response to occur [192]. As a result of the lack of wound healing response in the central lesion, it is reported that the host actually takes the central lesion as a foreign body and tries to wall out the central lesion from its surrounding tissue by dystrophic calcification in the peripheral zone [38,191,193]. More studies are necessary to clarify the long-term (i.e., in months and years) fate of the central zone in a hyperthermic lesion in vivo.

2.2. Cryothermic Injury

The response of biologicals to cryogenic temperature may not be simply thermotropic in that two additional biophysical events including ice formation and freeze concentration/dehydration may occur and can result in significant damage to cells [50,68,194]. As with hyperthermic injury, cryothermic injury has been investigated at the molecular, cellular and tissue (both in vitro and in vivo) levels [142-144]. At the molecular level, alterations in both lipid and proteins in response to cryogenic temperatures have been observed [195-197]. Dehydration and excursion in temperature are the thermodynamic driving forces for lipid phase separation and phase change among various states including liquid crystal lamellar (hydrated), gel (less hydrated), and even inverted micellar (Hex II) in bulk lipid membranes (composed of many types of lipids) [144]. Both phase separation and the Hex II phase change in lipids have been reported to be highly correlated with cell membrane damage which ultimately destroys cells upon thawing in various plant systems [198,199]. Presumably, lipid phase change and separation are also important in defining cryothermic injury to the plasma membranes of mammalian cells, although it has not been as clearly demonstrated as that in plant systems.

Protein denaturation during freezing is also inevitable because a freezing process is generally associated with change in pH and solute/electrolyte concentration, both of which can result in conformational alterations of proteins [143]. It has been reported that denaturation of cellular protein was minimal in thawed cells after freezing them to -20 oC (Fig. 6C) [195]. However, significant denaturation of cellular proteins was observed when the cells were frozen to -80 oC and thawed (Fig. 6C) [195]. In the latter case, an apparent shift of the secondary structure of cellular protein from α-helix to extended β-sheet (Fig. 6C) is observable in the FTIR spectra, indicating denaturation and aggregation of cellular proteins. Interestingly, cell viability also correlates well with the denaturation of cellular proteins at -80 oC according to the study [195]. Finally, as with hyperthermic injury, various molecules including TNF-α [200-205], antifreeze proteins [206,207], excessive electrolytes [208,209], immunopotentiator [210-213], and various anticancer drugs [214-218] have been shown to augment the efficacy of low temperatures for thermal destruction. For example, TNF-α has been found to increase the thermal threshold of cryoinjury to prostate cancer cells from ~ -20 to 0 oC [200], which is significant because it indicates that all cancer cells could be killed in a frozen tumor iceball and it is much more convenient to monitor the size of the frozen tumor iceball (e.g., using ultrasonography) than the subzero temperature inside the iceball. Therefore, using therapeutic adjuvants such as TNF-α to augment cryoinjury and the controlled/targeted delivery of the adjuvants into tumor and tumor cells using nanotechnology have attracted more and more attention in the field of thermal destruction at cryogenic temperatures [201,219,220].

At the cellular level, cryothermic injury can be attributed to ice formation induced mechanical damage and freeze concentration induced physiochemical deviation from the physiological state. The latter can be connected to the alteration to cellular proteins and lipids at the molecular level as discussed previously. Ice formation can occur both intracellularly and extracellularly. Although significant intracellular ice formation (IIF) is generally considered to be a lethal event, extracellular ice is not as damaging [50,194]. Significant IIF could destroy the cellular plasma membrane and other sub-cellular organelles mechanically as a result of volumetric expansion associated with the phase change or solidification of cellular water to ice, and the stochastic event of crack formation in the frozen cellular interior during freezing and thawing. The growing dendritic front of ice during IIF could also dismantle the native structure of cellular proteins, lipids, and subcellular organelles. In addition, significant IIF could lead to freeze concentration of the cellular interior which is detrimental to cellular proteins and lipids. Therefore, multiple mechanisms contribute to the lethality of significant IIF to cells.

For applications with a low concentration of CPA (e.g., cryopreservation by the conventional slow-freezing) or without CPA (e.g., thermal destruction or cryosurgery), significant IIF and freeze concentration induced excessive dehydration are well-established to be the two primary biophysical events that cause cryoinjury. Since freeze concentration requires sufficient time for dehydration of cellular water across the cell plasma membrane, its cell damaging effect is apparent only when cooling cells at slow cooling rates. For IIF, however, a higher cooling rate (tens to hundreds of degree Celsius per minute) is required so that intracellular water has no time to diffuse out of the cell and is available for IIF. This differential dependence on cooling rate of the damaging effect of IIF and freeze concentration results in the classical inverted U curve (dark solid line ➂ in Fig. 9) of cell viability vs. cooling rate for the conventional slow-freezing (for cryopreservation) and cryosurgery applications [50]. Of note, this inverted U curve does not apply to thermal stabilization by vitrification. For conventional vitrification, the cell injury mechanism is related to the high concentration of CPA used, which leads to cell damage by applying a significant osmotic stress and disturbing the normal metabolic pathways at super-zero temperatures before cooling and after warming the cells [55-60,221]. For low-CPA vitrification, cells loaded with a low concentration of CPA are brought to the vitrified phase without any ice formation or freeze concentration. Therefore, the traditional inverted U curve for cell survival vs. cooling rate can be modified and extended into the ultrafast cooling rate domain where cell viability increases with increasing cooling rate (dark dashed line ➃ in Fig. 9), which typically occurs when the cooling rate is as high as thousands to millions of degree Celsius per second when the CPA concentration is lower than 2-4 M.

At the in vitro tissue level, the cell injury mechanism discussed at the molecular and cellular level should still apply. In general, cells isolated from a tissue tend to transport water more quickly than cells within the intact tissue [144]. Therefore, cells form intracellular ice more readily in tissues due to enhanced water trapping as well as cell–ECM (extracellular matrix) and cell–cell contact issues. For example, IIF in primary hepatocytes, fibroblast, smooth muscle cells were observed to be enhanced after placing them in a matrix (of either collagen or fibrin) over suspension [222-224]. IIF has been suggested to propagate through the gap junction between adjacent cells [225-227]. Mechanically, abnormal thermal stress higher than the yield stress has been shown to accumulate in tissue and micro and macro crack formation in the tissue is commonly observed, which can lead to damage to the cells and the integrity of the extracellular matrix in both cryosurgery and cryopreservation applications [228-237]. Freezing has been shown to result in extensive damage to the vasculature in renal tissue, which presumably is responsible for the side effect of significant bleeding post cryosurgical operation of the kidney [238-243]. Blood vessel is the target of freezing induced mechanical damage, presumably because it resembles the defect in an otherwise more homogeneous tissue [229]. Freezing has also been shown to dramatically affect the biomechanical properties such as the porosity of the extracellular matrix of soft tissues [244] and the Young’s modulus and ultimate tensile strength (UTS) of blood vessels [236]. Considering the tight connection between cells and extracellular matrix through focal adhesion, mechanical damage to the extracellular matrix may further disrupt the integrity of cellular membrane. Of note, the presence of CPA can greatly reduce the mechanical damage, particularly when the CPA concentration is high enough to result in vitrification [245,246].

At the in vivo tissue level, several studies have shown that the subzero temperature required for causing tissue necrosis is higher than that in vitro [200,239,247]. Again, this observation has been attributed to the vascular stasis as a result of freezing induced damage to tissue vasculature during freezing and thawing, which results in the deprivation of nutrition and oxygen to cells, a secondary mechanism that kills cells in addition to direct cell injury. The wound healing response in an in vivo cryolesion however, is quite different from that in an in vivo hyperthermic lesion in that wound healing has been observed throughout the whole cryolesion [191]. For example, the cryolesion shrinks from day 3 to day 7 (Fig. 8B) indicating wound healing response that assimilates the necrotic tissue and allow the regeneration of normal tissue in the cryolesion [144,247]. This is probably because the frozen tissue has generally been observed to be reperfused immediately after thawing and vascular stasis does not occur usually until a few hours after thawing in cryosurgical applications [239,247]. Presumably, the short period of temporary reperfusion of the frozen tissue after thawing allows the recruitment of inflammatory cells for the subsequent wound healing response in the whole damaged tissue domain.

2.3. Measurement of Thermal (Cryo and Hyperthermic) Injury

Various methods have been utilized to measure thermal injury to biologicals. At the molecular level, x-ray crystallography, nuclear magnetic resonance (NMR) spectroscopy, Fourier transform infrared (FTIR) spectroscopy, circular dichroism (CD), electron microscopy, and differential scanning calorimetry (DSC), have been used to investigate thermally induced alterations to many proteins and lipids and sometimes DNA/RNAs [143,187]. Both x-ray crystallography and NMR spectroscopy can determine definitively the 3D structure of pure proteins. However, x-ray crystallography requires the protein sample to be in the crystal state and NMR spectroscopy can only be used to study small proteins or protein domains (≤ 20 kD) [151]. Both FTIR (Fig. 6A and C) and CD have been widely used to study the structural change of proteins (secondary structures such as α-helix and β-sheet in Fig. 5D) and lipids in response to a thermal challenge [143,187]. DSC is the only technique that can be used to investigate the calorimetric effect associated with the structural change in proteins (denaturation and agglomeration, Fig. 6B), lipids (phase transition and separation), and DNA/RNAs (denaturation, breaking up the intra-molecular hydrogen bond). Although isolated proteins and lipids have been used in many studies of thermal injury, there is a shift of research interests to investigate the in situ protein denaturation and lipid change in living cells recently. This is important because protein and lipid stability is strongly affected by the biochemical microenvironment in which they reside [143]. This is particularly important for determining the correlation between cell injury and thermally induced alterations in proteins and lipids. A number of studies using FTIR and DSC in this regard have reported strong correlations between cell injury and in situ thermal denaturation of protein in a number of cell types including prostate cancer cells (Fig. 6D), primary liver cells, and red blood cells [156,158,159,195].

Immediate cell viability (or injury) post thermal treatment is usually judged by the integrity of cell plasma membrane. The ability of an intact plasma membrane to exclude some membrane impermeable dye (e.g., ethidium homodimer showing red fluorescence in Fig. (7A), propidium iodide, and trypan blue) has been the most convenient way to visualize thermally induced damage to mammalian cells using microscopy (fluorescence capability is required for using ethidium homodimer and propidium iodide) [187,248]. The hydrolysis of calcein-AM (a non-fluorescent hydrophobic compound that easily permeates intact cells) by intracellular esterases (important enzymes for intracellular hydrolytic reactions) to produce calcein (a hydrophilic compound with strong green fluorescence that is impermeable to the plasma membrane of intact cells) has also been widely used to judge immediate cell viability using fluorescence microscopy (Fig. 7A). However, it has been shown that many cells with compromised membrane (judged by impermeable dyes) also show strong green fluorescence of calcein. Therefore, caution should be taken when calcein AM is used solely to judge immediate cell viability [187].

Since thermal damage to the other subcellular organelles can also determine the eventual fate of cells exposed to abnormally high and low temperatures, a better assay to evaluate the true (or long-term) cell viability can be done by checking the capability of the cells to attach to a substrate and subsequently proliferate [187]. The former is usually termed attachment efficiency and is evaluated by calculating the percentage of attached cells with regard to total cells seeded 3-12 hrs after thermal treatment and seeding. Proliferation is obtained by monitoring the change in cell number in a sample in several consecutive days after treatment. Clonogenics, another assay for determining the long-term cell survival, is calculated as the ratio (or percentage) of the number of cell colonies formed for thermally treated cells to that of control (i.e., no thermal treatment) samples usually after 7-10 days of culture after seeding [249]. This assay avoids counting cells in multiple days for the proliferation assay. For thermal stabilization applications, it is also important to verify that the cells retain important cell specific functions. For example, it is important to make sure hepatocytes (or liver cells) retain the normal function of urea production and albumin synthesis after biopreservation [222]. For stem cells, it is important to make sure that they retain pluripotency (e.g., the undifferentiated properties indicated by the expression of some stem cells specific markers shown in Fig. 7B and C for murine embryonic stem cells), multipotency, and/or the capability to differentiate into more specialized adult cells post biopreservation [250].

Histology (hemotoxylin&eosin or H&E stain) is commonly used to judge cell (Fig. 7D) and vascular injury in tissue [192]. Although being qualitative initially, histology has been used to quantify thermal injury to cells in tissue in several recent studies (Fig. 8A) [38,186,251]. Membrane impermeable fluorescence dyes have also been used to quantify cell injury in tissue after thermal treatment [186,251]. For some tissues abundant of collagen such as rat tail tendon, joint capsule, and myocardium, thermal injury is reflected by collagen denaturation which can be measured by a variety of techniques such as DSC, loss of birefringence observable using polarized light microscopy, loss of transparency, and shrinkage as reviewed elsewhere [143,148,252]. For injury assessment in situ in the host, perfusion defect (or absence of blood perfusion) in the damaged tissue is commonly used, which can be measured by MR contrast agent (e.g., gadolinium) exclusion or by fluorescence contrast in situ (Fig. 8B) [200,253-255]. Recent studies have shown that tissue injury assessed by perfusion defect corresponds quantitatively with tissue necrosis determined by histology using biopsy [256,257].

Among the various methods, clonogenics (or proliferation) and histology are the gold standard for judging cell and tissue viability, respectively. Besides the methods used to determine it, the cell viability (0-100%) is also dependent on many other factors including the specific thermal history applied, the attachment state (i.e., suspended, attached on dish, attached in 3D matrices) of the cells, and the cell culture microenvironment (e.g., pH) [187].


Besides injury measurement, there is a need of injury models to predict the outcome of thermal treatment for better understanding of the injury mechanism in research and pretreatment planning in the clinic. Injury models need to be verified or parameters in the model need to be extracted using experimental data a priori. With the verified model, injury to biologicals under any thermal history can be predicted. Furthermore, models may help us understand the injury mechanism by comparing the model parameters with those of known biological processes such as protein denaturation, lipid alteration, and DNA breakage [142,146,252]. Modeling hyperthermic injury at the molecular, cellular and tissue level has been reported in the literature for a wealth of biologicals [142,146,252]. However, modeling of cryothermic injury is not as well-developed due to the complication of two additional biophysical events: freeze concentration (leading to cell dehydration) and IIF in addition to the thermotropic effect of cryogenic temperature per se.

3.1. Modeling of Hyperthermic Injury

Several models have been developed to predict hyperthermic injury to biologicals including the statistical model, enzyme denaturation model, thermal isoeffective dose (TID) model, and various kinetic models [142,146,258-261]. Among the various models, the TID and kinetic models are the most widely used.

3.1.1. Kinetic Models

The first order irreversible kinetic model has been the most commonly used for predicting hyperthermic injury. In this model, biologicals including biomacromolecules, cells, and tissue are considered to be either in a native (N) or injured (I) state with an injury rate (or rate constant in the chemistry literature) k as follows:

(1)  NkI

An energy-state illustration of the first order kinetic model for describing the denaturation of proteins is given in Fig. (10), where the denatured state of proteins is equivalent to the injured state of cells or tissue in Equation 1. Thermally induced protein denaturation is often irreversible as a result of aggregation/coagulation and agglomeration following denaturation, which eventually results in irreversible cell death.

Assuming all biologicals being in the native state initially, the fraction of biologicals in the native state after being exposed to a thermal history T(t) can be calculated as follows :

(2)  FN=exp0tkdt

where F represents fraction, t is time, and the subscript ‘N’ represents native state. The injury rate (or rate constant) can be calculated using the Arrhenius model as follows [262-264]:

(3)  k= AexpE/RgT

where Rg is the universal gas constant (8.314 J mole-1 K-1), T (in K) is the thermal history (time dependent), and ΔE and A are the two kinetic model parameters (constants) usually called activation energy and frequency factor, respectively. This equation can be rewritten in the following format by taking the natural logarithm on both sides of the equation:

(4)  lnk = lnAE/RgT

Equation 4 shows that the injury rate (k) is constant for a given temperature. Therefore, Equation 2 can be rewritten as follows under isothermal heating conditions:

(5)  FN=expkt

Equation 5 has been frequently used to fit to data of fractional cell survival vs. isothermal heating time under various constant temperatures, through which the cell injury rate at various constant temperatures can be determined. With the known injury rate at various constant temperatures, the two model parameters (i.e., ΔE and A) can be determined by fitting to the data of k vs. T-1 at various constant temperatures using the linear relationship between the natural logarithm of k and T-1 given by Equation 4. The frequency factor and activation energy can be determined from the intercept and the slope of the linear fit, respectively. The procedures for determining the activation energy and frequency factor are demonstrated in Fig. (11).

In the hyperthermic literature, a terminology called injury and often represented by Ω is defined as follows [262-264]:

(6)   = 0tkdt

According to Equation 6, k is the first derivative of injury (Ω) with respect to time, which explains why k is often called injury rate in the hyperthermic literature. The relationship between Ω and the fraction of biologicals in the native state for a thermal treatment can be determined by combining Equations 6 and 2 as follows:

(7)  FN = exp

From this equation, it is clear that an accumulated injury of 1 (i.e., Ω=1) implies that the fractional survival is e-1 or ~ 36.7%.

The reaction rate k can also be calculated using the absolute rate theory model as follows [265,266]:

(8)  k = k'/h'Texpg"/RgT

(9)  g" = h"Ts"

where k' is the Boltzmann constant (1.38×10-23J K-1), h' is the Planck’s constant (6.625×10-34 J s), and Δh″, Δs″, and Δg″ are enthalpy, entropy, and Gibb’s free energy of activation for the thermal injury process, respectively. The two kinetic model parameters for the absolute rate theory model are the enthalpy and entropy of activation. The relationship between the parameters of the absolute rate theory model and those of the Arrhenius model for a first order reaction such as that described by Equation 1 is as follows [148]:

(10)  E = h"+RgT

(11)  A = k'/h'Texps"/Rg1

According to Equation 10, the difference (i.e. RgT) between ΔE and Δh″ is less than 3 kJ/mole when temperature is below 80°C. This small difference is negligible for most situations due to the relatively high value of the activation energy for cell/tissue injury and protein denaturation (typically, more than 50 kJ/mole) [142]. The frequency factor is linearly dependent on temperature and exponentially dependent on the entropy of activation. Since the unit for temperature is Kelvin, it is usually more than 293.15 K (assuming room temperature is 20 oC) for thermal therapy applications. Therefore, the dependence of the frequency factor on temperature is weak considering that the temperature for cell killing is usually between 50 and 70 oC in most thermal therapy applications.

Many studies have been performed to determine the parameters in the first order kinetic model for many different types of biologicals at the molecular (mainly proteins), cell, and tissue levels. These kinetic parameters for a wealth of biologicals have been summarized/tabulated and the following conclusion can be drawn by a careful examination of the kinetic data (Fig. 12) [142]:

  • The literature values of activation energy varies over a broad range from 50 to 1600 kJ/mole dependent on the type of biologicals (i.e., lipids, proteins, cells and/or tissues), temperature range, pH, and assay used. The activation energy for lipid change is typically lower than that for protein denaturation and cell/tissue injury. The activation energy for protein denaturation and cell/tissue injury however, is on the same order of magnitude, suggesting protein denaturation should be an important mechanism of thermally induced cell/tissue injury.
  • Clonogenic assay for single cells usually results in very high activation energy (~ 350-1600 kJ/mole) whereas it is generally much lower for histology and fluorescence viability dye assays (~ 100-800 kJ/mole).
  • An approximate, phenomenological, linear relationship between the natural logarithm of frequency factor and activation energy exists (see Fig. 12). This approximate linear relationship is attributed to the small variation of Gibbs free energy of activation given by Equation 9 for a thermal injury process despite the large variation of activation enthalpy (related to activation energy by Equation 10) and entropy (related to the frequency factor by Equation 11). The Gibb’s free energy of activation (Δg″) typically varies between 100 and 115 kJ/mole.
  • Temperature break points below and above which the activation energy changes significantly are mainly located at ~ 43°C and ~ 50°C if they do exist [142,187]. The break point at 43°C has been suggested to be a result of thermotolerance. The reason for a break point at 50°C is unclear yet.
  • With few exceptions, the literature values of frequency factor are greater than 2.0 × 1013 s–1 which is larger than the frequency factor predicted using Equation 11 assuming 0 activation entropy. In other words, the activation entropy for the kinetic process of thermally induced alterations in protein, cells, and tissue are typically positive.

Besides the first order kinetic model, higher order kinetic models have also been developed for modeling thermal injury [259,260]. However, higher order kinetic models have not been widely used probably due to the difficulty to definitively determine the multiple (> 2) parameters in those models using experimental data.

3.1.2. Thermal Isoeffective Dose (TID) Model

The thermal isoeffective dose (TID) model was proposed nearly 30 years ago [267] and has been widely used in both the traditional and many recent hyperthermic literatures [149]. In the TID model, a thermal history (i.e., time-temperature data) is converted to a cumulative equivalent minutes (isoeffect) at 43oC. As a result, the comparison of thermal doses and treatment planning of thermal therapy become straightforward, which may explain why the TID model has been widely used [268,269]. The time required to achieve the isoeffect of 43 oC at a constant temperature in terms of thermal injury is calculated as follows [149]:

(12)  t = CEM43CS/R43

where θ is temperature in oC, CEM43CS is the cumulative equivalent minutes required to decrease cell survival from the initial (100%) to a given criterion (S can be calculated as FN × 100%) at 43 oC, and R is the number of minutes required to compensate for a 1 oC temperature change to achieve an isoeffect of cell injury. The parameter R (not dimensionless) is a function of temperature and the activation energy of a thermal injury process [267]:

(13)  R = expE/RgTT1

If assuming 36.7% (e-1) cell survival as the criterion, the threshold time CEM43C can be calculated from the Arrhenius model as follows:

(14)  CEM 43C36.7%=1/k43C = A1expE/Rg43273.15

Note that CEM43C36.7% is the same as the inverse of the injury rate at 43 oC (k43C). To calculate CEM43C assuming a different cell survival S (CEM43CS), the following equation can be used [187]:

(15)  CEM43CS = CEM43C36.7%lnS1

where the subscript S and 36.7% represent the percentage of survival that is assumed as the criterion for calculating the cumulative equivalent time at 43 oC.

In many traditional literature of hyperthermia, a simplified version of the TID model was used that assumes a universal CEM43C of either 120 or 240 min for complete injury of all cells at 43 oC. The parameter R in the simplified TID model is taken as 0.25 and 0.5 for temperatures below and above 43 oC, respectively. According to Equations 13-15, both the CEM43C and R are dependent on the activation energy and frequency factor (R is also temperature dependent), both of which are strongly cell type dependent. Therefore, caution should be taken when using the simplified TID model for predicting thermal injury in different types of cells. The CEM43C for a variety of cells and tissue of urologic origin together with the temperature and cell type-dependent R values for thermal therapy applications have been tabulated elsewhere [187].

3.2. Modeling of Cryothermic Injury

Modeling of cryothermic injury is more challenging since it is complicated by two additional but important biophysical events: freeze concentration and ice formation. To the author’s knowledge, quantitative and mechanistic modeling of the effect of low temperature, freeze concentration, and ice formation on the stability of biomacromolecules such as lipids and proteins has not been reported. At the cell and tissue levels, formation of extracellular ice leads to freeze concentration along the liquidus in the phase diagram (see Fig. 3) which further results in cell dehydration at a slow cooling rate (e.g., < 10 oC/min). On the other hand, IIF will occur when the cooling rate is fast so that the intracellular water has no time to leave the cells through exosmosis, but not fast enough (< CRV shown in Fig. 9) to supersede the kinetics of ice formation. As with the molecular level injury, quantitative and mechanistic modeling of cryoinjury at the cellular and tissue levels has not been well studied. So far, much of the effort in modeling cryoinjury has been focused on modeling intracellular ice formation (IIF) and cell dehydration during freezing. IIF has been studied using both phenomenological and mechanistic models [194,270]. The mechanistic model has been widely used and delineates the ice formation as two consecutive events: Nucleation of ice nuclei and the subsequent growth of the nuclei [194,271]. Nucleation of intracellular ice can be catalyzed by either a surface (surface catalyzed nucleation, SCN) such as the cell plasma membrane or a volume of subcooled solution (volume catalyzed nucleation, VCN) such as the cytoplasm [194,271]. The rate of ice nucleation (I) due to either VCN or SCN can be modeled as follows [68,272]:

(16)  I = 0NN00TTf00.5exp0Tf/Tf04TTf2T3

where Tf is the equilibrium freezing point of the intracellular solution, N is the number of water molecules either in the cells (for VCN) or in contact with the cell plasma membrane (for SCN), η is viscosity, and Ω and κ are two model parameters that are usually called the kinetic and thermodynamic coefficient, respectively, and the subscript 0 represents the isotonic solution state before freezing (ice formation). The two model parameters under isotonic solution state (Ω0 and κ0) need to be determined by experimental studies and have been reported in the literature for a number of cell types as summarized elsewhere [142,194,273]. The cumulative probability of intracellular ice formation (PIIF) can then be calculated as follows [68,194]:

(17)  PIIFVCN = 1exp0tVIVCNdt

(18)  PIIFSCN = 1exp0tAISCNdt


where V and A are the cell volume and surface area available for catalyzing the nucleation of intracellular ice, respectively. Significant IIF is usually manifested as darkening of the cell cytoplasm when observed under a bright field of light microscopy as demonstrated in Fig. (13), which has been used widely to quantify the kinetics of intracellular ice formation [274,275].

Of note, the above IIF model is valid only when ice nucleation (the first step of IIF) is the rate-limiting step of IIF which is often true when freezing cells/tissues at not very high cooling rates (e.g., less than a few hundred Celsius per minute) and in the absence of a high concentration of CPA (e.g., less than 10 wt%) [194]. Under these conditions, ice will immediately propagate (or grow) throughout the entire intracellular space once the ice nuclei form, which results in the darkening event. When cooling cells/tissue with much higher cooling rates and/or a high concentration of intracellular CPA (e.g., during vitrification and after significant cell dehydration during slow-freezing), the rate limiting-step of IIF is the growth of the ice nuclei (the second step of IIF). The IIF under these conditions is said to be diffusion-limited for which more complicated models are needed to predict the amount of ice in cells [67,69,276-278]. To accurately predict the diffusion-limited ice nucleation and growth, an advanced model such as the free volume model that can be used to predict the effect of glass transition on solution viscosity and diffusion coefficient might be necessary. Such free volume model for a couple of cryo and lyoprotectants has been reported in [107].

In order to predict the probability of IIF using the above model during slow-freezing where freeze concentration induced cell dehydration is significant, information on the cell volume, V (or cell surface area, A, related to the diameter of the cells when the cells assume a spherical geometry), during freezing is required. The following model has been commonly used to predict the cell volume change during freezing [67,279]:

(20)  dVdt = LpARgTwlnVVbVsVVbVssnsvwhfRg1Tref1T

where νw is the partial molar volume of water, n represents amount (in mole), φs is dissociation constant of solutes (e.g., 2 for NaCl), Vs is volume of solutes, Vb is osmotically inactive volume in cells, Δhf is latent heat of fusion of water, Tref is a reference temperature usually taken as the equilibrium melting point of the intracellular solution, T is thermal history, the subscripts s and w represent solute (including CPAs) and water, respectively, and Lp is the cell plasma membrane permeability to water that can be calculated as follows [280]:

(21)  Lp = Lpg expELpRg1T1Tref

where Lpg is the permeability of the cell membrane to water at the reference temperature (Tref) and ELp is the activation energy for the water transport across the cell plasma membrane. In the equation, Lpg and ELp are two model parameters that need to be determined a priori using experimental data. Cell dehydration during freezing can be measured using either a specialized cryostage mounted on a light microscope [275] or differential scanning calorimetry [281-284]. Many studies have been performed to determine the two model parameters for various cells, which were reviewed and tabulated elsewhere [142,273,284,285].

Although cell dehydration and intracellular ice formation (IIF) can be predicted using the above models, a quantitative understanding of the mechanistic link between the two biophysical events and cell injury has not been well established despite some early efforts in this respect as reviewed elsewhere [142]. The incidence of significant IIF (defined as darkening of cell cytoplasm) correlates strongly with cell death in many cell types (i.e., 50% of IIF in many cell populations yields 50% of dead cells) [194]. However, the exact amount or percentage of intracellular ice that is significant enough to result in irreversible cell death is still unclear. Some studies even suggest that a small amount of intracellular ice might be beneficial to cell survival [286,287]. Therefore, further studies to establish mechanistic models capable of accounting for the effect of all the freezing induced biophysical events including intracellular ice formation and freeze concentration (i.e., the so-called solute effect), and low temperatures per se is important to further our understanding of low temperature biology and its biomedical application such as cryosurgery and cryopreservation.


Heat transfer in living systems is unique because of the effect on heat transfer of blood perfusion (and sometimes metabolism) (Fig. 5A) [21]. Generally, the thermal response in a bioheat transfer process (either heating or cooling) can be modeled using the following energy balance equation together with appropriate initial and boundary conditions:

(22)   cTLt = kTTQEEQBPQM

where ρ is density, c is specific heat, L is latent heat (per unit mass), ∧ is the latent heat release pattern that could be temperature and/or time dependent [142], kT is thermal conductivity, Q represents heat source/sink term, and the subscripts EE, BP, and M represent the source/sink term due to the applied external energy, blood perfusion, and metabolism, respectively. A detailed discussion about the different terms is given below. Note that the effect on heat transfer due to convection is omitted in the above energy equation since convection in the bulk (excluding the blood vessel that is accounted for by the heat source/sink term due to blood perfusion (QBP)) soft tissue is generally considered to be negligible.

4.1. Metabolic Heat Generation

In many studies of thermal therapy and cryosurgery, metabolic heat generation is taken as zero because it is generally much smaller than that due to blood perfusion and the externally applied energy source. In some studies particularly for the modeling of heating or cooling without freezing (no ice formation), the metabolic source term is taken to be temperature dependent as follows [288,289]:

(23)  QM = QM,37C37/10

where QM,37C is the baseline metabolic heat generation at 37 oC that is tissue type dependent and is summarized in [146] for a number of tissues/organs, ϕ is a model parameter that should be tissue type dependent as well although it is usually taken as 3 in the literature [288,289].

4.2. External Energy Source Term

The external energy source term (QEE) is related to the so-called specific absorption rate (SAR) in the literature, particularly when RF/microwave is used as the energy source for thermally destroying unwanted tissue or warming cryopreserved biologicals. The unit of QEE is power per volume while SAR is in power per weight. Therefore, QEE is the product of SAR and density. The direct approach for determining this external source term is to solve either the electromagnetic (for microwave, RF, and laser) or mechanical (for high intensity focused ultrasound or HIFU) wave equation to determine the electrical or pressure field (i.e., E and p), respectively [35,290,291]. For RF/microwave applications, QEE can be calculated from the electric field as follows [40]:

(24)  QEE = "E2

where σ is electrical conductivity, ω is angular frequency, ε" is the imaginary part of electrical permittivity, and ||E|| represents the magnitude of electric field. For laser and HIFU heating, the governing equations for wave propagation become very nonlinear as a result of wave scattering and cavitation (for HIFU). Therefore, a simplified model derived based on the well-known Beer-Lambert law of absorption has been utilized to predict the QEE of a plane wave as follows [292]:

(25)  QEE = I0ez fr

where α is absorption coefficient that accounts for all the mechanisms of heat generation, I0 is the incident wave intensity, z is the coordinate in the direction of laser or HIFU wave propagation, and f(r) is a function of the intensity distribution in the radial direction which is usually taken in the form of Gaussian distribution as follows [292]:

(26)  fr = exp2r/r02

where r is the radial coordinate and ro is known as the ‘e-2 radius’ of the wave beam since at r = ro, f(r)=e-2.

Besides theoretical prediction, experimental measurement of SAR (or QEE) has also been widely adopted for complicated probe designs, particularly when RF/microwave is used as the energy source. Since the experimental method can only determine SAR values at discrete, representative locations, analytical expressions are usually used to fit to the experimental data to predict SAR in a continuous domain [38,39]. For example, Fig. (14) shows the measured SAR pattern together with its fit for a typical axisymmetric dipolar helical microwave probes using the following equations [38]:

(27)  SAR = QE/=Ae2rrn1azbz2cz3dz4ez/z02

where r is radial coordinate originated from the probe centerline, z is the axial coordinate originated from the middle of the microwave thermal probe, β (= 0.413 cm-1 for kidney tissue) is attenuation constant, |z| represents the absolute of z, and the constant model parameters A, n, a, b, c, d, and z0 are 1.24×103, 0.9523, -0.3065, 2.783, -0.941, 0.2295, and 1.0515, respectively. The two SAR peaks at r = 0 are the result of the dipole design of the microwave thermal probe. The SAR distribution along the axial direction of the probe (except near the two ends) is more homogeneous than that along the radial direction. SAR generally decreases exponentially in the radial direction due to tissue absorption of the electromagnetic energy.

4.3. Thermal Effect of Blood Perfusion

The effect of blood perfusion on thermal modeling in unfrozen tissue particularly during heating has been the focus of many studies in the literature. In frozen tissue however, blood perfusion is stopped and has no thermal effect. The classical and probably the most frequently used model of the blood perfusion term is given as follows (note, Equation 22 with the blood perfusion term given below is often called the classical Pennes bioheat equation) [293]:

(28)  QBP = wBcBTaT

where w is blood perfusion rate and the subscripts B and a represent blood and artery, respectively. This model assumes that blood enters tissue at the arterial temperature and leaves at an equilibrated temperature which is the same as the local tissue temperature, regardless of the tissue vasculature. This assumption has been questioned by many studies in the literature. It has been shown that thermal equilibration between blood and tissue occurs before the blood reaches the capillary bed, and heat exchange between blood and tissue in aorta and large arteries (D > 3 mm) is negligible in terms of affecting the blood temperature [294,295]. Further theoretical studies have shown that thermal equilibration between blood and tissue actually occurs in precapillary arteries (or arterioles) that are approximately 50 μm in diameter [296]. On the venous side, blood temperature begins to deviate from its surrounding tissue in approximately the terminal vein class (or venules also ~ 50 µm in diameter) [296]. Blood vessels with a diameter between 50 µm and 3 mm are therefore termed thermally significant vessels. These results show that the heat transfer process between blood and tissue is very different from the mass transfer process. The latter occurs dominantly in the capillary bed. Therefore, it is important to account for the effect on heat transfer of different classes (or diameters) of blood vessels differently for accurate thermal modeling. The model for blood perfusion in the classical Pennes bioheat equation seems inadequate for accurate thermal modeling.

Further investigation have shown that arterioles and venules that lead to and from the capillary bed lie in closely matched pairs with a countercurrent flow pattern, and models are therefore developed to take into account of this observation [296-305]. In this case, up to three separate but coupled governing equations are required to model the coupled heat transfer processes in arteries, veins, and the tissue region as reviewed elsewhere [146,306,307]. Although models accounting for the effects of vessel size and countercurrent flow are more theoretically rigorous, their application is very limited [146]. This is probably because the more mechanistic models require knowledge of the complex anatomy of the vasculature in a specific tissue, which unfortunately, is often unknown [306]. On the other hand, many recent studies have verified the applicability of the Pennes bioheat equation in a variety of tissues/organs by experimental data, especially in tissues/organs (e.g., kidney cortex) where small blood vessels (D < 0.3 mm) dominate [308-311]. Blood vessels of 0.3 – 3 mm in diameter could have a significant local cooling or heating effect to the adjacent tissue, which is not considered in the Pennes bioheat equation and should be accounted for in thermal modeling with special attention.

With much attention being focused on modeling the thermal effect of blood perfusion using more realistic vasculature data, the magnitude of the blood perfusion (i.e., w) has been often taken as constant. However, the blood perfusion during heating has been observed to be strongly dependent on both heating time and temperature [312-314], which has been shown to significantly affect the thermal and injury modeling [38,315,316]. For example, the change of blood perfusion has been modeled as follows [317]:

(29)  w = w37Ce

where w37C is the baseline blood perfusion rate at 37 oC that can be found in the literature for many types of tissues [146] and Ω is injury to the vasculature which can be calculated using Equation 6. Therefore, this model basically takes the blood perfusion change in response to heat as a first order kinetic process. This model may be sufficient to predict the change of blood perfusion in tumor, which has been observed to decrease monotonically with the increase of thermal dose (i.e., the accumulation of heat with time reflected by injury, Ω) [312-314]. However, it is inadequate to model blood perfusion response to heating in normal tissue, which has been observed to increase initially followed by a gradual decrease when the thermal dose is accumulated to beyond a threshold level to result in stasis in the tissue vasculature [38]. The initial increase of blood perfusion in response to heating in normal tissue has been attributed to the host mediated physiological response, in which elevated blood perfusion is used to dissipate heat to bring the temperature back to normal. Tumor tissue vasculature is believed to be generally less well organized and leaky and therefore, lacks such a host mediated physiological response [318-320]. To model the blood perfusion response to heating in normal tissue, Equation 29 can be corrected as follows [38]:

(30)  w = w37C fee

where f(e-Ω) is a correction function. The blood perfusion response in thermal therapy of normal porcine kidney tissue has been reported (see Fig. 15) and can be fit using Equation 30 with the following correction function:

(31)  fxNPKT = 33.6x4+37x351.8x2+32.2x7.5104,x>0.9276.9x4315.3x3+169.2x242.9x+5.1,x0.9

where the subscript NPKT represents normal porcine kidney tissue and x represents exp(-Ω) with Ω given by Equation 6. The two kinetic parameters (ΔE, and A) in Equation 6 for calculating x are 98.73 kJ/mole and 1.81×1013 s-1, respectively [38]. Further studies to quantify the two kinetic parameters and determine the correction function (for normal tissues only) for various types of living systems are necessary for the wide application of this time and temperature dependent model for blood perfusion and therefore, accurate modeling of the thermal response in heating tumor and the surrounding normal living tissue.

4.4. Latent Heat and Its Absorption/Release Pattern

A major component of biological materials is water, which undergoes phase change from liquid to either vapor when boiling occurs as a result of heating or ice when solidification/freezing occurs as a result of cooling. For an extremely slow cooling/heating process, phase change of pure liquid water to either the vapor or ice phase should occur at its equilibrium freezing (~ 0 oC) or boiling (~ 100 oC) point accompanying with a significant amount of release (for freezing) or absorption (for boiling) of latent heat, respectively. Due to the existence of impurities such as salts and cryoprotectants in tissue water, the latent heat release during freezing of biological materials, however, occurs over a temperature range called the mushy zone [321]. For thermal modeling, the upper boundary (θub) of the mushy zone are typically taken as the equilibrium freezing point (θefp) of a solution, which is a colligative quantity and can be estimated based on the osmolality of the solution as follows [322]:

(32)  efp = 1.858Osm

where Osm represents osmolality and the subscript efp represents equilibrium freezing point. Osmolality is the ratio of the total amount (in mole) of particles (ions and non-dissociated molecules) in a solution to the total amount (in kg) of solvent (water) in the solution. For a dilute, aqueous solution at room temperature (when the density of water is ~ 1 kg/L), the value of osmolality is similar to that of osmolarity that is the ratio of the total amount (in mole) of particles in a solution to the total volume (in liter or L) of the solution. For a typical physiological solution with an osmolality of 0.3, the equilibrium freezing point is around -0.56 oC according to the equation. The lower boundary (θlb) of the mushy zone used for thermal modeling varies significantly from studies to studies (typically between -4 and -10 oC and can be as low as -40 oC). Various latent heat release patterns (^) in the mushy zone have been used as summarized in [142]. A simple but quite frequently used linear release pattern of latent heat with regard to temperature is given as follows:

(33)   = ub/lbub,lbub

It has been reported that the magnitude of latent heat (L) has a much more significant impact on thermal modeling than its release pattern [323]. Biological samples also contain lipids, proteins, and other organic and inorganic substances with essentially negligible latent heat release comparing with water during freezing. Moreover, some amount of water may bind with solutes and never become frozen (i.e., the so-called bound water) during cooling. Therefore, the latent heat of a biological sample should be less than that of pure water that is confirmed by recent studies on measuring the latent heat of some biological solutions using DSC [324,325]. In the absence of experimental data, however, the latent heat of a biological sample may be estimated as a first cut by using the total tissue water content and the latent heat of pure water as follows [142]:

(34)  L = Lw

where φ is total water content in a biological sample in weight and the latent heat of solidification for pure water (Lw) is 335 J/g. Some studies approximates the latent heat of biological tissue as the product of the total water content (φ) and the latent heat of a physiological solution (1× PBS, ~ 300 J/g [324,325]), which is much closer to tissue than pure water in terms of chemical composition. The total water content of a number of tissues can be found in the literature [146,326] and/or determined quantitatively by drying (i.e., dry weight vs. wet weight) or magnetic resonance spectrometry.

Very few studies have been reported on the absorption pattern of latent heat during boiling in heating biological materials, probably because boiling is typically avoided in many thermal therapy protocols. The rationale is that boiling off tissue water at a temperature higher than 100 oC leads to tissue carbonization, a type of wound or lesion that is sometimes very difficult to heal in the host. In the event that evaporation does occur, the amount of latent heat release may be calculated similarly to Equation 34 based on the weight fraction of water in the sample. The latent heat of evaporation for pure water is 2,250 J/g. Similarly, the pattern of latent heat absorption may be modeled as being linearly dependent on temperature over a small temperature range (e.g., 99-100 oC) [316].

4.5. Thermophysical Properties

The knowledge of thermophysical properties (mainly thermal conductivity and specific heat) of biological tissues is crucial for accurate modeling of the thermal response in a bioheat transfer process. Many studies have been performed to measure the properties of various biological tissues, and the results from these studies have been reviewed and tabulated in the literature [146,285,326,327]. Nonetheless, data on thermophysical properties of tissue below -30 °C or above 40 °C are still scant. Such data are undoubtedly in demand because these temperatures are typically encountered in thermal destruction and stabilization applications.

The main factors that affect tissue thermophysical properties are temperature and sample composition (mainly water content). Because thermophysical properties of water (> 0 oC) and ice (< 0 oC) over a wide temperature range (from liquid nitrogen to boiling water) are well-known and the total water content of many tissues are available [38,228,321,326,328], the properties of biological samples may be estimated using those of water/ice based on their total water content. By lumping together the effects of all non-aqueous contents, the specific heat (c) of unfrozen biomaterials may be estimated as follows [142]:

(35)  cuf = 0.41cw

where the subscript uf represents unfrozen and the first and second terms in the brackets represent the effect on specific heat of aqueous and non-aqueous contents in tissue, respectively. A similar relationship between water content and thermal conductivity (kT) of unfrozen biomaterials is as follows [142]:

(36)  kT,uf = 0.281kT,w

where the first and second terms in the brackets represent the effect on thermal conductivity of the aqueous and non-aqueous contents in tissue, respectively. One can suggest a similar approach based on water content to estimate the thermophysical properties of frozen tissue. In this case, Equations 35 and 36 can be modified to predict the specific heat and thermal conductivity of frozen tissue as follows [142]:

(37)  cf = cice0.41cw

(38)  kT,f = kT,ice0.281kT,w

where the subscript f represents frozen and the first and second terms in Equations 37 and 38 represent the effect of tissue water and non-aqueous contents on the thermophysical properties, respectively. Because the non-aqueous contents do not undergo phase change during freezing, their contributions to the specific heat and thermal conductivity represented by the second term of Equations 37 and 38 are calculated based on that of water rather than ice. In addition, one might use the thermal conductivity and specific heat of water at 0 °C to approximate that of sub-cooled water below 0 °C as the first cut considering that the latter are still not well-established. Finally, determining the thermophysical properties of subcooled and vitrified water or solutions below their freezing point is crucial for accurate thermal modeling of the bioheat transfer processes in thermal stabilization by vitrification (without ice formation).


Temperatures over a wide range (from cryogenic to hyperthermic regimens) have been utilized to either selectively destroy unwanted or stabilize important biologicals for the treatment of diseases. Thermal destruction can be realized by freezing biologicals to below -20 oC using cryoprobes or heating them to more than 50 oC using RF/microwave, laser, and HIFU as the energy source in a minimally invasive manner. Cryopreservation of biologicals by slow-freezing is associated with injury due to ice formation and freeze concentration induced cell dehydration. Conventional vitrification is associated with cytotoxicity of high concentration of cryoprotectants. Research focus in cryopreservation is shifting to achieve vitrification using a low concentration of cryoprotectants (low-CPA vitrification), which combines the benefits of both slow-freezing and conventional vitrification while avoiding their shortcomings. Eventually, thermal stabilization at ambient temperature (dry or lyopreservation) is the solution to the increasing demand of wide and convenient distribution of biopreserved products for modern cell-based medicine. However, lyopreservation of eukaryotic mammalian cells at ambient temperature is still challenging today. Integrating modern micro and nanotechnology into the field of cell preservation might be the solution to the challenges. Studies have been performed to understand both hyperthermic and cryothermic injury at the molecular, cell, and tissue levels. The mechanism of injury to biologicals of hyperthermic temperatures is better understood than that of cryogenic temperatures. Interestingly, using therapeutic agents such as TNF-α and anticancer drugs to augment thermal injury and the targeted/controlled delivery of the agents to tumor site using nanotechnology are attracting more and more attention. Mechanistic models for hyperthermic injury have been developed with success whereas developing such models for cryothermic injury is more challenging and very much in need. Much more work is still in need for thermal modeling of bioheat transfer. For hyperthermic applications, it is particularly challenging to model the change of blood perfusion with regard to both temperature and time (or the extent of vascular stasis) for various types of tissues. For cryogenic applications, the heat transfer process is affected by additional thermophysical events such as ice formation and glass transition. Therefore, measurement of latent heat (as a result of ice formation and glass transition) and thermophysical properties of various biological materials at cryogenic (either frozen or vitrified) temperatures is of great importance for accurate thermal modeling of the bioheat transfer process. With the development of advanced models strengthened by further experimental investigations, the thermostability of biologicals over a wide temperature range (from that of liquid nitrogen to that of boiling water) can be engineered for either effective destruction of unwanted biomaterials (e.g., tumor) or banking of desired biologicals (e.g., stem cells) for the emerging cell-based medicine.


This work was partially supported by grants from the Wendy Will Case Cancer Fund, National Science Foundation (CBET-1033426), and National Institute of Biomedical Imaging and Bioengineering of the National Institutes of Health (R01EB012108). Special gratitude is due to Dr. John C. Bischof at the University of Minnesota and Dr. Mehmet Toner at Harvard University (Medical School) for their mentorship and guidance during the author’s doctoral and postdoctoral training in their labs where much of the work reviewed in this article was done. The author would like to thank Dr. Alex Fowler at the University of Massachusetts, Dartmouth, for many useful discussions on thermal stabilization of biologicals. Many thanks are due to Dr. Dieter Haemmerich at the Medical University of South Carolina for a critical review of the manuscript.


[1] Lu Gd, Needham J. Celestial Lancets: A History and Rationale of Acupuncture and Moxa 1980.
[2] Sivin N. "Science and Medicine in Imperial China-The State of the Field" J Asian Stud 1988; 47: 41-90.
[3] Unschuld PU. Huang Di nei jing su wen: Nature, Knowledge, Imagery in an Ancient Chinese Medical Text. Berkeley and Los Angeles: University of California Press 2003.
[4] Waddle JA. "Radiofrequency ablation of liver and lung tumors" Radiologic Technology 2006 Sept-Oct;
[5] Duffell E. "Curative power of fever" Lancet 2001; 358: 1276.
[6] Aufderheide AC. The Scientific Study of Mummies. Cambridge, UK: Cambridge University Press 2003.
[7] Biello D. "Cold Comfort: Cool Receptors May Ease Chronic Pain" Sci Am 2006 August 22;
[8] Adams F. "Aphorisms by Hippocrates (Section VII)", 2010 [accessed: April 13, 2010];
[9] Phillips C. "Radiofrequency Ablation Making Inroads as Cancer Treatment" NCI Cancer Bulletin 2005 July 19;
[10] Bommadevara M, Zhu L. "Temperature difference between the body core and arterial blood supplied to the brain during hyperthermia or hypothermia in humans" Biomech Model Mechanobiol 2002; 1: 137-49.
[11] Cabanac M. "Selective brain cooling in humans: "fancy" or fact?" FASEB J 1993; 7: 1143-6. discussion 1146-1147
[12] Storm FK. Hyperthermia in Cancer Therapy. Boston: GK Hall Medical Publisher 1989.
[13] Gautherie M. Biological Basis of Oncologic Thermotherapy. Berlin-Heidelberg: Springer-Verlag 1990.
[14] Hall EJ. "Hyperthermia", in Radiobiology for the Radiobiologist. 4th Edition. Philadelphia. PA: Lippincott Company 1994.
[15] Paliwal BR, Hetzel FW, Dewhirst MW. Biol, Phys and Clini Aspe Hyperthermia. New York: American Institute of Physics 1987.
[16] Field B, Hand JW. Introduction to the Clinical Aspects of Clinical Hyperthermia. Philadelphia, PA: Taylor and Francis 1990.
[17] Desai MM, Gill IS. "Current status of cryoablation and radiofrequency ablation in the management of renal tumors" Curr Opin Urol 2002; 12: 387-93.
[18] Jahangiri M, Weir G, Mandal K, Savelieva I, Camm J. "Current strategies in the management of atrial fibrillation" Ann Thorac Surg 2006; 82: 357-64.
[19] Marberger M. "Energy-based ablative therapy of prostate cancer: high-intensity focused ultrasound and cryoablation" Curr Opin Urol 2007; 17: 194-9.
[20] Savage SJ, Gill IS. "Renal tumor ablation: energy-based technologies" World J Urol 2000; 18: 283-8.
[21] Diller KR, Ryan TP. "Heat transfer in living systems: Current opportunities" J Heat Transfer 1998; 120: 810-29.
[22] Schumann C, Rieder C, Bieberstein J, et al. "State of the Art in Computer-Assisted Planning, Intervention, and Assessment of Liver-Tumor Ablation" Crit Rev Biomed Eng 2010; 38: 31-52.
[23] Bernardo NO, Gill IS. "Laparoscopic partial nephrectomy: current status" Arch Esp Urol 2002; 55: 868-0.
[24] Hall-Craggs MA, Vaidya JS. "Minimally invasive therapy for the treatment of breast tumours" Eur J Radiol 2002; 42: 52-7.
[25] Spaliviero M, Moinzadeh A, Gill IS. "Laparoscopic cryotherapy for renal tumors" Technol Cancer Res Treat 2004; 3: 177-80.
[26] Gage AA, Baust JG. "Cryosurgery for tumors" J Am Coll Surg 2007; 205: 342-56.
[27] Haemmerich D, Laeseke PF. "Thermal tumour ablation: devices, clinical applications and future directions" Int J Hyperthermia 2005; 21: 755-60.
[28] Diederich CJ. "Thermal ablation and high-temperature thermal therapy: overview of technology and clinical implementation" Int J Hyperthermia 2005; 21: 745-53.
[29] Shitzer A. "Cryosurgery: Analysis and Experimentation of Cryoprobes in Phase Changing Media" J Heat Transfer 2011; 133: 011005(011001-011012).
[30] Rewcastle J C, Sandison G A, Saliken J C, Donnelly B J, McKinnon J G. "Considerations during clinical operation of two commercially available cryomachines" J Surg Oncol 1999; 71: 106-1.
[31] Reif F. "Chapter 5 – Simple applications of macroscopic thermodynamics" In: in Fundamentals of Statistical and Thermal Physics. New York: McGraw-Hill 1965.
[32] Stafford RJ, Fuentes D, Elliott AA, Weinberg JS, Ahrar K. "Laser-Induced Thermal Therapy for Tumor Ablation" Crit Rev Biomed Eng 2010; 38: 79-100.
[33] Simpson CR, Kohl M, Essenpreis M, Cope M. "Near-infrared optical properties of ex vivo human skin and subcutaneous tissues measured using the Monte Carlo inversion technique" Phys Med Biol 1998; 43: 2465-78.
[34] Sershen S R, Westcott S L, Halas N J, West J L. "Temperature-sensitive polymer-nanoshell composites for photothermally modulated drug delivery" J Biomed Mater Res 2000; 51: 293-8.
[35] Ter Haar G, Coussios C. "High intensity focused ultrasound: Physical principles and devices" Int J Hyperther 2007; 23: 89-104.
[36] Farny CH, Holt RG, Roy RA. "The Correlation Between Bubble-Enhanced HIFU Heating and Cavitation Power" IEEE Trans Biomed Eng 2010; 57: 175-84.
[37] Coussios CC, Farny CH, Ter Haar G, Roy RA. "Role of acoustic cavitation in the delivery and monitoring of cancer treatment by high-intensity focused ultrasound (HIFU)" Int J Hyperther 2007; 23: 105-20.
[38] He X, Mcgee S, Coad JE, et al. "Investigation of the thermal and tissue injury behaviour in microwave thermal therapy using a porcine kidney model" Int J Hyperther 2004; 20: 567-93.
[39] Zhu L, Xu L X, Chencinski N. "Quantification of the 3-D electromagnetic power absorption rate in tissue during transurethral prostatic microwave thermotherapy using heat transfer model" IEEE Trans Biomed Eng 1998; 45: 1163-72.
[40] Polk C, Postow E. Handbook of Biological Effects of Electromagnetic Fields. 2nd edition. Boca Raton, FL: CRC Press LLC 1995.
[41] Haemmerich D. "Biophysics of Radiofrequency Ablation" Crit Rev Biomed Eng 2010; 38: 53-63.
[42] Brace CL. "Microwave Tissue Ablation: Biophysics, Technology, and Applications" Crit Rev Biomed Eng 2010; 38: 65-78.
[43] Day ES, Morton JG, West JL. "Nanoparticles for thermal cancer therapy" J Biomech Eng 2009; 131: 074001.
[44] Huang X, El-Sayed IH, El-Sayed MA. "Applications of gold nanorods for cancer imaging and photothermal therapy" Methods Mol Biol 2010; 624: 343-57.
[45] Belzer FO, Southard JH. "Organ preservation and transplantation" Prog Clin Biol Res 1986; 224: 291-303.
[46] Maathuis MH, Leuvenink HG, Ploeg RJ. "Perspectives in organ preservation" Transplantation 2007; 83: 1289-98.
[47] Fuller BJ, Lee CY. "Hypothermic perfusion preservation: the future of organ preservation revisited?" Cryobiology 2007; 54: 129-45.
[48] Fahy GM, MacFarlane DR, Angell CA, Meryman HT. "Vitrification as an approach to cryopreservation" Cryobiology 1984; 21: 407-26.
[49] Karlsson J O M, Toner M. "Long-term storage of tissues by cryopreservation: Critical issues" Biomaterials 1996; 17: 243-56.
[50] Mazur P. "Freezing of living cells: mechanisms and implications" Am J Physiol 1984; 247: C125-42.
[51] Rall WF, Fahy GM. "Ice-Free Cryopreservation of Mouse Embryos at -196-Degrees-C by Vitrification" Nature 1985; 313: 573-.
[52] Coger R, Toner M. "Preservation Techniques for Biomaterials" In: Bronzino JD, Ed. in The Biomedical Engineering Handbook. 2nd Edition. Boca Raton: CRC Press LLC 2000.
[53] He X. "Preservation of Embryonic Stem Cells" In: in Embryonic Stem Cells. Vienna, Austria: Intech 2010.
[54] Fahy GM, Wowk B, Wu J, et al. "Cryopreservation of organs by vitrification: perspectives and recent advances" Cryobiology 2004; 48: 157-78.
[55] Fahy GM, Wowk B, Wu J, Paynter S. "Improved vitrification solutions based on the predictability of vitrification solution toxicity" Cryobiology 2004; 48: 22-35.
[56] Fowler A, Toner M. "Cryo-injury and biopreservation" Ann N Y Acad Sci 2005; 1066: 119-35.
[57] Heng BC, Kuleshova LL, Bested SM, Liu H, Cao T. "The cryopreservation of human embryonic stem cells" Biotechnol Appl Biochem 2005; 41: 97-104.
[58] Hunt CJ, Pegg DE, Armitage SE. "Optimising cryopreservation protocols for haematopoietic progenitor cells: A methodological approach for umbilical cord blood" Cryoletters 2006; 27: 73-83.
[59] Chen SU, Lien YR, Cheng YY, Chen HF, Ho HN, Yang YS. "Vitrification of mouse oocytes using closed pulled straws (CPS) achieves a high survival and preserves good patterns of meiotic spindles, compared with conventional straws, open pulled straws (OPS) and grids" Hum Reprod 2001; 16: 2350-6.
[60] Chen SU, Lien YR, Chen HF, Chao KH, Ho HN, Yang YS. "Open pulled straws for vitrification of mature mouse oocytes preserve patterns of meiotic spindles and chromosomes better than conventional straws" Hum Reprod 2000; 15: 2598-603.
[61] Fuller BJ. "Cryoprotectants: The essential antifreezes to protect life in the frozen state" Cryoletters 2004; 25: 375-88.
[62] Beattie GM, Crowe JH, Lopez AD, Cirulli V, Ricordi C, Hayek A. "Trehalose: A cryoprotectant that enhances recovery and preserves function of human pancreatic islets after long-term storage" Diabetes 1997; 46: 519-23.
[63] Wowk B, Leitl E, Rasch CM, Mesbah-Karimi N, Harris SB, Fahy GM. "Vitrification enhancement by synthetic ice blocking agents" Cryobiology 2000; 40: 228-36.
[64] He X, Park EY, Fowler A, Yarmush ML, Toner M. "Vitrification by ultra-fast cooling at a low concentration of cryoprotectants in a quartz micro-capillary: a study using murine embryonic stem cells" Cryobiology 2008; 56: 223-32.
[65] Boutron P. "Comparison with the theory of the kinetics and extent of ice crystallization and of the glass-forming tendency in aqueous cryoprotective solutions" Cryobiology 1986; 23: 88-102.
[66] Berejnov V, Husseini NS, Alsaied OA, Thorne RE. "Effects of cryoprotectant concentration and cooling rate on vitrification of aqueous solutions" J Appl Crystallogr 2006; 39: 244-51.
[67] Karlsson JOM, Cravalho EG, Toner M. "A Model of Diffusion-Limited Ice Growth inside Biological Cells during Freezing" J Appl Phys 1994; 75: 4442-5.
[68] Toner M, Cravalho EG, Karel M. "Thermodynamics and Kinetics of Intracellular Ice Formation during Freezing of Biological Cells" J Appl Phys 1990; 67: 1582-93.
[69] Yang G, Zhang A, Xu LX, He X. "Modeling the cell-type dependence of diffusion-limited intracellular ice nucleation and growth during both vitrification and slow freezing" J Appl Phys 2009; 105: 114701-.
[70] Bruggeller P, Mayer E. "Complete Vitrification in Pure Liquid Water and Dilute Aqueous-Solutions" Nature 1980; 288: 569-71.
[71] Bhat SN, Sharma A, Bhat SV. "Vitrification and glass transition of water: Insights from spin probe ESR" Phys Rev Lett 2005; 95: 2357021-4.
[72] Yavin S, Arav A. "Measurement of essential physical properties of vitrification solutions" Theriogenology 2007; 67: 81-9.
[73] Gardner DK, Sheehan CB, Rienzi L, Katz-Jaffe M, Larman MG. "Analysis of oocyte physiology to improve cryopreservation procedures" Theriogenology 2007; 67: 64-72.
[74] Vajta G, Nagy ZP. "Are programmable freezers still needed in the embryo laboratory? Review on vitrification" Reprod Biomed Online 2006; 12: 779-96.
[75] Han X, Ma H, Jiao A, Critser JK. "Investigations on the heat transport capability of a cryogenic oscillating heat pipe and its application in achieving ultra-fast cooling rates for cell vitrification cryopreservation" Cryobiology 2008; 56: 195-203.
[76] Jiao A, Han X, Critser JK, Ma H. "Numerical investigations of transient heat transfer characteristics and vitrification tendencies in ultra-fast cell cooling processes" Cryobiology 2006; 52: 386-92.
[77] Lee HJ, Elmoazzen H, Wright D, et al. "Ultra-rapid vitrification of mouse oocytes in low cryoprotectant concentrations" Reprod Biomed Online 2010; 20: 201-8.
[78] Zhang W, Yang G, Zhang A, Xu LX, He X. "Preferential vitrification of water in small alginate microcapsules significantly augments cell cryopreservation by vitrification" Biomed Microdevices 2010; 12: 89-96.
[79] Crowe JH, Hoekstra FA, Crowe LM. "Anhydrobiosis" Ann Rev Phys 1992; 54: 579-99.
[80] Browne J, Tunnacliffe A, Burnell A. "Anhydrobiosis: plant desiccation gene found in a nematode" Nature 2002; 416: 38.
[81] Perry R N. "Desiccation survival of parasitic nematodes" Parasitology 1999; 119(Suppl): S19-30.
[82] Clegg J S. "Cryptobiosis--a peculiar state of biological organization" Comp Biochem Physiol B Biochem Mol Biol 2001; 128: 613-24.
[83] Crowe J H, Crowe L M, Tablin F. "Stabilization of dry mammalian cells: Lessons from nature" Integr Comp Biol 2004; 44: 542-2.
[84] Crowe J H, Crowe L M. "Preservation of mammalian cells-learning nature's tricks" Nat Biotechnol 2000; 18: 145-6.
[85] Crowe JH, Cooper AF. "Cryptobiosis" Sci Am 1971; 225: 30-6.
[86] Armstrong WP. "The Power Of Plants-July, 1996" [April 21, 2010];
[87] Crowe JH, Crowe LM, Oliver AE, Tsvetkova N, Wolkers W, Tablin F. "The trehalose myth revisited: introduction to a symposium on stabilization of cells in the dry state" Cryobiology 2001; 43: 89-105.
[88] Crowe JH, Crowe LM, Wolkers WF, et al. "Stabilization of dry mammalian cells: Lessons from nature" Integr Comp Biol 2005; 45: 810-20.
[89] Crowe J H, Crowe L M. "Preservation of mammalian cells-learning nature's tricks" Nat Biotechnol 2000; 18: 145-6.
[90] Wolkers WF, Tablin F, Crowe JH. "From anhydrobiosis to freeze-drying of eukaryotic cells" Comp Biochem Physiol A Mol Integr Physiol 2002; 131: 535-43.
[91] Eroglu A, Russo MJ, Bieganski R, et al. "Intracellular trehalose improves the survival of cryopreserved mammalian cells" Nat Biotechnol 2000; 18: 163-67.
[92] Rey L, May JC. Freeze-drying/lyophilization of pharmaceutical and biological products. New York: Marcel Dekker, Inc 1999.
[93] Felix F. Freeze-drying of pharmaceuticals and biopharmaceuticals. Cambridge, UK: RSC Publishing 2007.
[94] Abdul-Fattah AM, Kalonia DS, Pikal MJ. "The challenge of drying method selection for protein pharmaceuticals: product quality implications" J Pharm Sci 2007; 96: 1886-916.
[95] Nail SL, Jiang S, Chongprasert S, Knopp SA. "Fundamentals of freeze-drying" Pharm Biotechnol 2002; 14: 281-360.
[96] Bellows R J, King C J. "Freeze-drying of aqueous solutions: Maximum allowable operating temperature" Cryobiology 1972; 9: 559-61.
[97] Pikal MJ, Shah S. "The collapse temperature in freeze drying: Dependence on measurement methodology and rate of water removal from the glassy phase" Int J Pharm 1990; 62: 165-86.
[98] Gieseler H, Lee H, Mulherkar B, Pikal MJ. "Applicability of Manomentric Temperature Measurement (MTM) and SMART Freeze-Dryer Technology to Development of an Optimized Freeze-Drying Cycle: Prelim-inary Investigation of Two Amorphous Systems". Proc. 1st European Congress on Life Science Process Technology; Nuremberg, Germany. 2005.
[99] Meister E, Marek E, Gieseler H. "Collapse Temperature Measurement by Freeze-Dry Microscopy and Transferability to Freeze Drying Processes: Influence of Solute Concentration on Collapse Behavior and Effect on Cycle Design". Proc. AAPS Annual Meeting and Exposition; San Antonio (TX), USA. 2006.
[100] Pikal MJ. "Use of laboratory data in freeze drying process design: heat and mass transfer coefficients and the computer simulation of freeze drying" J Parenter Sci Technol 1985; 39: 115-38.
[101] MacKenzie AP. "Basic principles of freeze-drying for pharmaceutics" Bull Parenteral Drug Assoc 1966; 20: 101-29.
[102] Pikal MJ, Shah S, Senior D, Lang JE. "Physical chemistry of freeze-drying: measurement of sublimation rates for frozen aqueous solutions by a microbalance technique" J Pharm Sci 1983; 72: 635-50.
[103] Meister E, Gieseler H. "Evaluation of Collapse Temperatures by Freeze-Dry Microscopy: Impact of Excipient Concentration on Measured Transition and the Overall Dependance on Measurement Methodology" 2006. Proc. 5th World Meeting on Pharmaceutics and Pharmaceutical Technology; Geneva Switzerland. 2006.
[104] Kramer T, Kremer DM, Pikal MJ, Petre WJ, Shalaev EY, Gatlin LA. "A procedure to optimize scale-up for the primary drying phase of lyophilization" J Pharm Sci 2009; 98: 307-18.
[105] Fonseca F, Passot S, Lieben P, Marin M. "Collapse temperature of bacterial suspensions: the effect of cell type and concentration" Cryo Lett 2004; 25: 425-34.
[106] Fonseca F, Passot S, Cunin O, Marin M. "Collapse temperature of freeze-dried Lactobacillus bulgaricus suspensions and protective media" Biotechnol Prog 2004; 20: 229-38.
[107] He XM, Fowler A, Toner M. "Water activity and mobility in solutions of glycerol and small molecular weight sugars: Implication for cryo- and lyopreservation" J Appl Phys 2006; 100: 074702 (074711pp).
[108] He X, Fowler A, Menze M, Hand S, Toner M. "Desiccation kinetics and biothermodynamics of glass forming trehalose solutions in thin films" Ann Biomed Eng 2008; 36: 1428-39.
[109] Hancock BC, Shamblin SL, Zografi G. "Molecular mobility of amorphous pharmaceutical solids below their glass transition temperatures" Pharm Res 1995; 12: 799-806.
[110] Yang G, Gilstrap K, Zhang A, Xu LX, He X. "Collapse temperature of solutions important for lyopreservation of living cells at ambient temperature" Biotechnol Bioeng 2010; 106: 247-59.
[111] Acker JP, Chen T, Fowler A, Toner M. "Engineering desiccation tolerance in mammalian cells: tools and techniques" In: Fuller BJ, Lane L, Benson EE, Eds. in Life in the frozen state. Boca Raton, FL: CRC Press LLC 2004; pp. 563-81.
[112] Chen T, Acker JP, Eroglu A, et al. "Beneficial effect of intracellular trehalose on the membrane integrity of dried mammalian cells" Cryobiology 2001; 43: 168-81.
[113] Eroglu A, Toner M, Toth TL. "Beneficial effect of microinjected trehalose on the cryosurvival of human oocytes" Fertil Steril 2002; 77: 152-8.
[114] Bhowmick P, Eroglu A, Wright DL, Toner M, Toth TL. "Osmometric behavior of mouse oocytes in the presence of different intracellular sugars" Cryobiology 2002; 45: 183-7.
[115] Eroglu A, Lawitts JA, Toner M, Toth TL. "Quantitative microinjection of trehalose into mouse oocytes and zygotes, and its effect on development" Cryobiology 2003; 46: 121-34.
[116] Eroglu A, Elliott G, Wright DL, Toner M, Toth TL. "Progressive elimination of microinjected trehalose during mouse embryonic development" Reprod Biomed Online 2005; 10: 503-10.
[117] Gordon SL, Oppenheimer SR, Mackay AM, Brunnabend J, Puhlev I, Levine F. "Recovery of human mesenchymal stem cells following dehydration and rehydration" Cryobiology 2001; 43: 182-7.
[118] Guo N, Puhlev I, Brown DR, Mansbridge J, Levine F. "Trehalose expression confers desiccation tolerance on human cells" Nat Biotechnol 2000; 18: 168-71.
[119] Puhlev I, Guo N, Brown DR, Levine F. "Desiccation tolerance in human cells" Cryobiology 2001; 42: 207-17.
[120] Acker JP, Lu XM, Young V, et al. "Measurement of trehalose loading of mammalian cells porated with a metal-actuated switchable pore" Biotechnol Bioeng 2003; 82: 525-32.
[121] Elliott GD, Liu XH, Cusick JL, et al. "Trehalose uptake through P2X7 purinergic channels provides dehydration protection" Cryobiology 2006; 52: 114-27.
[122] Liu XH, Aksan A, Menze MA, Hand SC, Toner M. "Trehalose loading through the mitochondrial permeability transition pore enhances desiccation tolerance in rat liver mitochondria" Biochim Biophys Acta 2005; 1717: 21-6.
[123] Shirakashi R, Kostner CM, Muller KJ, Kurschner M, Zimmermann U, Sukhorukov VL. "Intracellular delivery of trehalose into mammalian cells by electropermeabilization" J Membr Biol 2002; 189: 45-54.
[124] Reuss R, Ludwig J, Shirakashi R, et al. "Intracellular delivery of carbohydrates into mammalian cells through swelling-activated pathways" J Membr Biol 2004; 67-81.
[125] Oliver AE, Jamil K, Crowe JH, Tablin F. "Loading MSCs with trehalose by endocytosis" Cell Preserv Technol 2004; 2: 35-49.
[126] Wolkers WF, Looper SA, Fontanilla RA, Tsvetkova NM, Tablin F, Crowe JH. "Temperature dependence of fluid phase endocytosis coincides with membrane properties of pig platelets" Biochim Biophys Acta 2003; 1612: 154-63.
[127] He XM, Amin AA, Fowler A, Toner M. "Thermally induced introduction of trehalose into primary rat hepatocytes" Cell Preserv Technol 2006; 4: 178-87.
[128] Beattie GM, Crowe JH, Lopez AD, Cirulli V, Ricordi C, Hayek A. "Trehalose: a cryoprotectant that enhances recovery and preserves function of human pancreatic islets after long-term storage" Diabetes 1997; 46: 519-23.
[129] Crowe JH. "Trehalose as a "chemical chaperone": fact and fantasy" Adv Exp Med Biol 2007; 594: 143-58.
[130] Kanias T, Acker JP. "Mammalian cell desiccation: facing the challenges" Cell Preserv Technol 2006; 4: 253-76.
[131] Zhang W, Rong J, Wang Q, He X. "The encapsulation and intracellular delivery of trehalose using a thermally responsive nanocapsule" Nanotechnol 2009; 20: 275101.
[132] Holovati JL, Gyongyossy-Issa MI, Acker JP. "Effects of trehalose-loaded liposomes on red blood cell response to freezing and post-thaw membrane quality" Cryobiology 2009; 58: 75-83.
[133] Holovati JL, Acker JP. "Spectrophotometric measurement of intraliposomal trehalose" Cryobiology 2007; 55: 98-107.
[134] Scott KL, Gyongyossy-Issa MIC, Acker JP. "Response of hemotopoietic progenitor cells to trehalose-loaded liposomes (Abstract)" Cryibiology 2006; 53: 380.
[135] Zhang W, Gilstrap K, Wu N, et al. "Synthesis and characterization of thermally responsive pluronic F127-Chitosan nanocapsules for controlled release and intracellular delivery of small molecules" ACS Nano 2010; 4: 6747-59.
[136] de Castro AG, Lapinski J, Tunnacliffe A. "Anhydrobiotic engineering" Nat Biotechnol 2000; 18: 473.
[137] Iturriaga G. "The LEA proteins and trehalose loving couple: a step forward in anhydrobiotic engineering" Biochem J 2008; 410: e1-2.
[138] Hand SC, Jones D, Menze MA, Witt TL. "Life without water: Expression of plant LEA genes by an anhydrobiotic arthropod" J Exp Zool A Ecol Genet Physiol 2007; 307A: 62-.
[139] Huang Z, Tunnacliffe A. "Desiccation response of mammalian cells: anhydrosignaling" Methods Enzymol 2007; 428: 269-77.
[140] Tunnacliffe A, Garcia de Castro A, Manzanera M. "Anhydrobiotic engineering of bacterial and mammalian cells: is intracellular trehalose sufficient?" Cryobiology 2001; 43: 124-32.
[141] Li D, He X. "Desiccation induced structural alterations in a 66-amino acid fragment of an anhydrobiotic nematode late embryogenesis abundant (LEA) protein" Biomacromolecules 2009; 10: 1469-77.
[142] He X, Bischof JC. "Quantification of temperature and injury response in thermal therapy and cryosurgery" Crit Rev Biomed Eng 2003; 31: 355-422.
[143] Bischof JC, He X. "Thermal stability of proteins" Ann N Y Acad Sci 2005; 1066: 12-33.
[144] Bischof JC. "Micro and nanoscale phenomenon in bioheat transfer" Heat Mass Transf 2006; 42: 955-66.
[145] Wright NT. "On a relationship between the Arrhenius parameters from thermal damage studies" J Biomech Eng 2003; 125: 300-4.
[146] Diller KR, Valvano JW, Pearce JA. "Bioheat transfer" In: Kreith F, Ed. in CRC Thermal Engineering Handbook. Boca Raton, FL: CRC Press 2000.
[147] Hoffmann N E, Bischof J C. "The cryobiology of cryosurgical injury" Urology 2002; 60: 40-9.
[148] Pearce J, Thomsen S. "Rate process analysis of thermal damage" In: Welch AJ, van Germert MJC, Eds. in Optical Thermal Response of Laser-Irradiated Tissue. New York: Plenum Press 1995.
[149] Dewhirst MW, Viglianti BL, Lora-Michiels M, Hanson M, Hoopes PJ. "Basic principles of thermal dosimetry and thermal thresholds for tissue damage from hyperthermia" Int J Hyperther 2003; 19: 267-94.
[150] Dewey WC. "Failla memorial lecture. The search for critical cellular targets damaged by heat" Radiat Res 1989; 120: 191-204.
[151] Alberts B, Bray D, Johnson A, et al. Essential Cell Biology: An introduction to the molecular biology of the cell. New York: Garland Publishing, Inc 1998.
[152] Lepock JR. "Involvement of membranes in cellular responses to hyperthermia" Radiat Res 1982; 92: 433-8.
[153] Lepock J R. "Cellular effects of hyperthermia: relevance to the minimum dose for thermal damage" Int J Hyperther 2003; 19: 252-66.
[154] Wimberly BT, Brodersen DE, Clemons WM, et al. "Structure of the 30S ribosomal subunit" Nature 2000; 407: 327-9.
[155] Ritchie KP, Keller BM, Syed KM, Lepock JR. "Hyperthermia (heat shock)-induced protein denaturation in liver, muscle and lens tissue as determined by differential scanning calorimetry" Int J Hyperther 1994; 10: 605-18.
[156] Lepock JR, Frey HE, Ritchie KP. "Protein denaturation in intact hepatocytes and isolated cellular organelles during heat shock" J Cell Biol 1993; 122: 1267-76.
[157] Freeman ML, Borrelli MJ, Syed K, Senisterra G, Stafford DM, Lepock JR. "Characterization of a signal generated by oxidation of protein thiols that activates the heat shock transcription factor" J Cell Physiol 1995; 164: 356-66.
[158] He X, Wolkers WF, Crowe JH, Swanlund DJ, Bischof JC. "In situ thermal denaturation of proteins in dunning AT-1 prostate cancer cells: implication for hyperthermic cell injury" Ann Biomed Eng 2004; 32: 1384-98.
[159] Lepock JR, Frey HE, Bayne H, Markus J. "Relationship of hyperthermia-induced hemolysis of human erythrocytes to the thermal denaturation of membrane proteins" Biochim Biophys Acta 1989; 980: 191-201.
[160] Cress AE, Gerner EW. "Cholesterol levels inversely reflect the thermal sensitivity of mammalian cells in culture" Nature 1980; 283: 677-9.
[161] Palecz D, Grzelinska E, Bartosz G, Leyko W, Moller JV. "Ligand and lipid domain stabilization of a membraneous Ca2+-ATPase during hyperthermia" Biochim Biophys Acta 1988; 937: 23-30.
[162] Anghileri LJ, Robert J. "Effects of hyperthermia and lanthanum on tumor cell leakage" Int J Clin Pharmacol Ther Toxicol 1987; 25: 374-8.
[163] Quinn PJ. "Principles of membrane stability and phase behavior under extreme conditions" J Bioenerg Biomembr 1989; 21: 3-19.
[164] Yatvin MB, Vorpahl JW, Gould MN, Lyte M. "The effects of membrane modification and hyperthermia on the survival of P-388 and V-79 cells" Eur J Cancer Clin Oncol 1983; 19: 1247-53.
[165] Yatvin MB, Cramp WA. "Role of cellular membranes in hyperthermia: some observations and theories reviewed" Int J Hyperther 1993; 9: 165-85.
[166] Rao W, Deng Z-S, Liu J. "A Review of Hyperthermia Combined With Radiotherapy/Chemotherapy on Malignant Tumors" Crit Rev Biomed Eng 2010; 38: 101-6.
[167] Visaria RK, Griffin RJ, Williams BW, et al. "Enhancement of tumor thermal therapy using gold nanoparticle-assisted tumor necrosis factor-alpha delivery" Mol Cancer Ther 2006; 5: 1014-20.
[168] Unezaki S, Maruyama K, Takahashi N, et al. "Enhanced delivery and antitumor activity of doxorubicin using long-circulating thermosensitive liposomes containing amphipathic polyethylene glycol in combination with local hyperthermia" Pharm Res 1994; 11: 1180-5.
[169] Han HD, Choi MS, Hwang T, et al. "Hyperthermia-induced antitumor activity of thermosensitive polymer modified temperature-sensitive liposomes" J Pharm Sci 2006; 95: 1909-7.
[170] Kong G, Dewhirst MW. "Hyperthermia and liposomes" Int J Hyperther 1999; 15: 345-70.
[171] Visaria R, Bischof JC, Loren M, et al. "Nanotherapeutics for enhancing thermal therapy of cancer" Int J Hyperther 2007; 23: 501-11.
[172] Ponce AM, Vujaskovic Z, Yuan F, Needham D, Dewhirst MW. "Hyperthermia mediated liposomal drug delivery" Int J Hyperther 2006; 22: 205-13.
[173] Dewhirst MW. "Hyperthermia and nanotechnology - a note from the Editor-in-chief" Int J Hyperther 2008; 24: 449-50.
[174] Dennis CL, Jackson AJ, Borchers JA, et al. "Nearly complete regression of tumors via collective behavior of magnetic nanoparticles in hyperthermia" Nanotechnol 2009; 20: 395103.
[175] Burke A, Ding X, Singh R, et al. "Long-term survival following a single treatment of kidney tumors with multiwalled carbon nanotubes and near-infrared radiation" Proc Natl Acad Sci USA 2009; 106: 12897-902.
[176] Park H, Yang J, Lee J, Haam S, Choi IH, Yoo KH. "Multifunctional nanoparticles for combined doxorubicin and photothermal treatments" ACS Nano 2009; 3: 2919-6.
[177] Urano M. "Tumor response to hyperthermia" Hyperther in Oncol 1988; 1: 161-200.
[178] Diller KR. "Stress protein expression kinetics" Annu Rev Biomed Eng 2006; 8: 403-24.
[179] Rylander MN, Feng Y, Bass J, Diller KR. "Thermally induced injury and heat-shock protein expression in cells and tissues" Ann N Y Acad Sci 2005; 1066: 222-42.
[180] Rylander MN, Diller KR, Wang S, Aggarwal SJ. "Correlation of HSP70 expression and cell viability following thermal stimulation of bovine aortic endothelial cells" J Biomech Eng 2005; 127: 751-7.
[181] Rylander MN, Feng Y, Bass J, Diller KR. "Heat shock protein expression and injury optimization for laser therapy design" Lasers Surg Med 2007; 39: 731-46.
[182] Song AS, Diller KR. "Modeling heat shock protein expression produced by a heat wrap" J Biomech Eng 2009; 131: 074510.
[183] Wang S, Diller KR, Aggarwal SJ. "Kinetics study of endogenous heat shock protein 70 expression" J Biomech Eng 2003; 125: 794-7.
[184] Wang S, Xie W, Rylander MN, Tucker PW, Aggarwal S, Diller KR. "HSP70 kinetics study by continuous observation of HSP-GFP fusion protein expression on a perfusion heating stage" Biotechnol Bioeng 2008; 99: 146-54.
[185] He X M, Bischof J C. "The kinetics of thermal injury in human renal carcinoma cells" Ann Biomed Eng 2005; 33: 502-10.
[186] Bhowmick S, Coad JE, Swanlund DJ, Bischof JC. "In vitro thermal therapy of AT-1 Dunning prostate tumours" Int J Hyperther 2004; 20: 73-92.
[187] He X, Bhowmick S, Bischof JC. "Thermal therapy in urologic systems: a comparison of arrhenius and thermal isoeffective dose models in predicting hyperthermic injury" J Biomech Eng 2009; 131: 074507.
[188] Shah B, Bhowmick S. "Evaluation of important treatment parameters in supraphysiological thermal therapy of human liver cancer HepG2 cells" Ann Biomed Eng 2006; 34: 1745-57.
[189] Bhowmick S, Hoffmann NE, Bischof JC. "Thermal therapy of prostate tumor tissue in the dorsal skin flap chamber" Microvasc Res 2002; 64: 170-3.
[190] Shen Y, Liu P, Zhang A, Xu LX. "Study on tumor microvasculature damage induced by alternate cooling and heating" Ann Biomed Eng 2008; 36: 1409-9.
[191] Coad JE, Bischof JC. "Histologic differences between cryothermic and hyperthermic therapies" Proc SPIE 2003; 4954: 27-36.
[192] Cotran RS, Kumar V, Robbins SL. Robbins pathologic basis of disease. 5th edition. Philadelphia, PA: W.B. Sanders Company 1994.
[193] Hennings L, Kaufmann Y, Griffin R, et al. "Dead or alive? Autofluorescence distinguishes heat-fixed from viable cells" Int J Hyperther 2009; 25: 355-63.
[194] Toner M. "Nucleation of ice crystals inside biological cells" In: Steponkus P, Ed. Advances in Low-Temperature Biology. London: JAI Press 1993; 2: pp. 1-51.
[195] Bischof JC, Wolkers WF, Tsvetkova NM, Oliver AE, Crowe JH. "Lipid and protein changes due to freezing in Dunning AT-1 cells" Cryobiology 2002; 45: 22-32.
[196] Wolkers WF, Balasubramanian SK, Ongstad EL, Zec HC, Bischof JC. "Effects of freezing on membranes and proteins in LNCaP prostate tumor cells" Biochim Biophys Acta 2007; 1768: 728-36.
[197] Balasubramanian SK, Wolkers WF, Bischof JC. "Membrane hydration correlates to cellular biophysics during freezing in mammalian cells" Biochim Biophys Acta 2009; 1788: 945-53.
[198] Steponkus PL, Lynch DV. "Freeze/thaw-induced destabilization of the plasma membrane and the effects of cold acclimation" J Bioenerg Biomembr 1989; 21: 21-41.
[199] Yamada T, Kuroda K, Jitsuyama Y, Takezawa D, Arakawa K, Fujikawa S. "Roles of the plasma membrane and the cell wall in the responses of plant cells to freezing" Planta 2002; 215: 770-8.
[200] Chao BH, He X, Bischof JC. "Pre-treatment inflammation induced by TNF-alpha augments cryosurgical injury on human prostate cancer" Cryobiology 2004; 49: 10-27.
[201] Goel R, erson K, Slaton J, et al. "Adjuvant approaches to enhance cryosurgery" J Biomech Eng 2009; 131: 074003.
[202] Goel R, Swanlund D, Coad J, Paciotti GF, Bischof JC. "TNF-alpha-based accentuation in cryoinjury--dose, delivery, and response" Mol Cancer Ther 2007; 6: 2039-47.
[203] Han B, Iftekhar A, Bischof JC. "Improved cryosurgery by use of thermophysical and inflammatory adjuvants" Technol Cancer Res Treat 2004; 3: 103-1.
[204] Han B, Swanlund DJ, Bischof JC. "Cryoinjury of MCF-7 human breast cancer cells and inhibition of post-thaw recovery using TNF-alpha" Technol Cancer Res Treat 2007; 6: 625-34.
[205] Jiang J, Goel R, Iftekhar MA, et al. "Tumor necrosis factor-alpha-induced accentuation in cryoinjury: mechanisms in vitro and in vivo" Mol Cancer Ther 2008; 7: 2547-55.
[206] Koushafar H, Rubinsky B. "Effect of antifreeze proteins on frozen primary prostatic adenocarcinoma cells" Urology 1997; 49: 421-5.
[207] Pham L, Dahiya R, Rubinsky B. "An in vivo study of anti-freeze protein adjuvant cryosurgery" Cryobiology 1999; 38: 169-75.
[208] Han B, Bischof JC. "Direct cell injury associated with eutectic crystallization during freezing" Cryobiology 2004; 48: 8-21.
[209] Wang CL, Teo KY, Han B. "An amino acidic adjuvant to augment cryoinjury of MCF-7 breast cancer cells" Cryobiology 2008; 57: 52-9.
[210] Sugiyama Y, Saji S, Miya K, et al. "[Therapeutic effect of multimodal therapy, such as cryosurgery, locoregional immunotherapy and systemic chemotherapy against far advanced breast cancer]" Gan To Kagaku Ryoho 2001; 28: 1616-9.
[211] Kawaguchi Y, Sugiyama Y, Saji S. "Cryoimmunological therapy with local injection of OK-432 against advance or recurrent breast cancer" Gan To Kagaku Ryoho 2003; 30: 1583-86.
[212] Forest V, Peoc'h M, Campos L, Guyotat D, Vergnon JM. "Effects of cryotherapy or chemotherapy on apoptosis in a non-small-cell lung cancer xenografted into SCID mice" Cryobiology 2005; 50: 29-37.
[213] Forest V, Peoc'h M, Campos L, Guyotat D, Vergnon JM. "Benefit of a combined treatment of cryotherapy and chemotherapy on tumour growth and late cryo-induced angiogenesis in a non-small-cell lung cancer model" Lung Cancer 2006; 54: 79-86.
[214] Ikekawa S, Ishihara K, Tanaka S, Ikeda S. "Basic studies of cryochemotherapy in a murine tumor system" Cryobiology 1985; 22: 477-83.
[215] Clarke DM, Baust JM, Van Buskirk RG, Baust JG. "Chemo-cryo combination therapy: an adjunctive model for the treatment of prostate cancer" Cryobiology 2001; 42: 274-85.
[216] Mir LM, Rubinsky B. "Treatment of cancer with cryochemotherapy" Br J Cancer 2002; 86: 1658-60.
[217] Yuan F, Zhou W, Zhang J, et al. "Anticancer drugs are synergistic with freezing in induction of apoptosis in HCC cells" Cryobiology 2008; 57: 60-5.
[218] Clarke DM, Baust JM, Van Buskirk RG, Baust JG. "Addition of anticancer agents enhances freezing-induced prostate cancer cell death: implications of mitochondrial involvement" Cryobiology 2004; 49: 45-61.
[219] Yan J F, Liu J. "Nanocryosurgery and its mechanisms for enhancing freezing efficiency of tumor tissues" Nanomed 2008; 4: 79-87.
[220] Liu J, Deng ZS. "Nano-cryosurgery: advances and challenges" J Nanosci Nanotechnol 2009; 9: 4521-2.
[221] Trad FS, Toner M, Biggers JD. "Effects of cryoprotectants and ice-seeding temperature on intracellular freezing and survival of human oocytes" Hum Reprod 1999; 14: 1569-77.
[222] Yarmush ML, Toner M, Dunn JC, Rotem A, Hubel A, Tompkins RG. "Hepatic tissue engineering. Development of critical technologies" Ann N Y Acad Sci 1992; 665: 238-52.
[223] Balasubramanian S K, Bischof JC, Hubel A. "Water transport and IIF parameters for a connective tissue equivalent" Cryobiology 2006; 52: 62-73.
[224] Balasubramanian SK, Venkatasubramanian RT, Menon A, Bischof JC. "Thermal injury prediction during cryoplasty through in vitro characterization of smooth muscle cell biophysics and viability" Ann Biomed Eng 2008; 36: 86-101.
[225] Irimia D, Karlsson JO. "Kinetics and mechanism of intercellular ice propagation in a micropatterned tissue construct" Biophys J 2002; 82: 1858-68.
[226] Stott SL, Irimia D, Karlsson JO. "Parametric analysis of intercellular ice propagation during cryosurgery, simulated using monte carlo techniques" Technol Cancer Res Treat 2004; 3: 113-23.
[227] Armitage WJ, Juss BK. "Freezing monolayers of cells without gap junctions" Cryobiology 2003; 46: 194-6.
[228] He X, Bischof JC. "Analysis of thermal stress in cryosurgery of kidneys" J Biomech Eng 2005; 127: 656-1.
[229] Rabin Y, Olson P, Taylor MJ, Steif PS, Julian TB, Wolmark N. "Gross damage accumulation on frozen rabbit liver due to mechanical stress at cryogenic temperatures" Cryobiology 1997; 34: 394-405.
[230] Rabin Y, Steif PS. "Analysis of thermal stresses around a cryosurgical probe" Cryobiology 1996; 33: 276-90.
[231] Rabin Y, Steif PS, Taylor MJ, Julian TB, Wolmark N. "An experimental study of the mechanical response of frozen biological tissues at cryogenic temperatures" Cryobiology 1996; 33: 472-82.
[232] Steif PS, Palastro MC, Rabin Y. "Analysis of the effect of partial vitrification on stress development in cryopreserved blood vessels" Med Eng Phys 2007; 29: 661-70.
[233] Steif PS, Palastro MC, Rabin Y. "The effect of temperature gradients on stress development during cryopreservation via vitrification" Cell Preserv Technol 2007; 5: 104-5.
[234] Steif PS, Noday DA, Rabin Y. "Can thermal expansion differences between cryopreserved tissue and cryoprotective agents alone cause cracking? Cryo Lett 2009; 30: 414-21.
[235] Najimi S, Rubinsky B. "Non-invasive detection of thermal stress fractures in frozen biological materials" Cryo Lett 1997; 18: 209-16.
[236] Venkatasubramanian RT, Grassl ED, Barocas VH, Lafontaine D, Bischof JC. "Effects of freezing and cryopreservation on the mechanical properties of arteries" Ann Biomed Eng 2006; 34: 823-32.
[237] Zhang A, Cheng S, Gao D, Xu LX. "Thermal stress study of two different artery cryopreservation methods" Cryo Lett 2005; 26: 113-20.
[238] Schmidlin FR, Rupp CC, Hoffmann NE, et al. "Measurement and prediction of thermal behavior and acute assessment of injury in a pig model of renal cryosurgery" J Endourol 2001; 15: 193-7.
[239] Rupp CC, Hoffmann NE, Schmidlin FR, Swanlund DJ, Bischof JC, Coad JE. "Cryosurgical changes in the porcine kidney: histologic analysis with thermal history correlation" Cryobiology 2002; 45: 167-82.
[240] Hruby G, Edelstein A, Karpf J, et al. "Risk factors associated with renal parenchymal fracture during laparoscopic cryoablation" BJU Int 2008; 102: 723-26.
[241] Janzen NK, Perry KT, Han KR, et al. "The effects of intentional cryoablation and radio frequency ablation of renal tissue involving the collecting system in a porcine model" J Urol 2005; 173: 1368-74.
[242] Rukstalis DB, Khorsandi M, Garcia FU, Hoenig DM, Cohen JK. "Clinical experience with open renal cryoablation" Urology 2001; 57: 34-9.
[243] Nakada SY, Lee FT Jr, Warner T, Chosy SG, Moon TD. "Laparoscopic cryosurgery of the kidney in the porcine model: an acute histological study" Urology 1998; 51: 161-6.
[244] Han B, Miller JD, Jung JK. "Freezing-induced fluid-matrix interaction in poroelastic material" J Biomech Eng 2009; 131: 021002.
[245] Baicu S, Taylor MJ, Chen Z, Rabin Y. "Vitrification of carotid artery segments: An integrated study of thermophysical events and functional recovery toward scale-up for clinical applications" Cell Pre Tech 2006; 4: 236-44.
[246] Wusteman M, Robinson M, Pegg D. "Vitrification of large tissues with dielectric warming: biological problems and some approaches to their solution" Cryobiology 2004; 48: 179-89.
[247] Hoffmann NE, Bischof JC. "Cryosurgery of normal and tumor tissue in the dorsal skin flap chamber: Part II--injury response" J Biomech Eng 2001; 123: 310-6.
[248] Pegg DE. "Viability assays for preserved cells, tissues, and organs" Cryobiology 1989; 26: 212-31.
[249] Sahu SK, Song CW. "Thermal sensitivity and kinetics of thermotolerance in bovine aortic endothelial cells in culture" Int J Hyperther 1991; 7: 103-1.
[250] Gearhart J. "New potential for human embryonic stem cells" Science 1998; 282: 1061-2.
[251] Bhowmick P, Coad JE, Bhowmick S, et al. "In vitro assessment of the efficacy of thermal therapy in human benign prostatic hyperplasia" Int J Hyperther 2004; 20: 421-39.
[252] Cravalho EG, Toner M, Gaylor DC, Lee RC. "Response of cells to supraphysiological temperatures: experimental measurements and kinetic models" In: Lee RC, Cravalho EG, Burke JF, Eds. Electrical Trauma: The Pathophysiology, Manifestations and Clinical Management. Cambridge, UK: Cambridge Univ. Press 1992.
[253] Lee JM, Choi SH, Park HS, et al. "Radiofrequency thermal ablation in canine femur: evaluation of coagulation necrosis reproducibility and MRI-histopathologic correlation" AJR Am J Roentgenol 2005; 185: 661-7.
[254] Donnelly SE, Donnelly BJ, Saliken JC, Raber EL, Vellet AD. "Prostate cancer: gadolinium-enhanced MR imaging at 3 weeks compared with needle biopsy at 6 months after cryoablation" Radiology 2004; 232: 830-33.
[255] Porter CA t, Woodrum DA, Callstrom MR, et al. "MRI after technically successful renal cryoablation: early contrast enhancement as a common finding" AJR Am J Roentgenol 2010; 194: 790-3.
[256] Mala T, Frich L, Aurdal L, et al. "Intraoperative contrast-enhanced MR-imaging as predictor of tissue damage during cryoablation of porcine liver" Magn Reson Imaging 2003; 21: 733-40.
[257] Cheng HL, Purcell CM, Bilbao JM, Plewes B. "Prediction of subtle thermal histopathological change using a novel analysis of Gd-DTPA kinetics" J Magn Reson Imaging 2003; 18: 585-98.
[258] Feng Y, Tinsley Oden J, Rylander MN. "A two-state cell damage model under hyperthermic conditions: theory and in vitro experiments" J Biomech Eng 2008; 130: 041016.
[259] Dienes GJ. "A kinetic model of biological radiation response" Radiat Res 1966; 28: 183-202.
[260] Landry J, Marceau N. "Rate-limiting events in hyperthermic cell killing" Radiat Res 1978; 75: 573-85.
[261] Chang IA. "Considerations for thermal injury analysis for RF ablation devices" Open Biomed Eng J 4: 3-12.
[262] Henriques FC, Moritz AR. "Studies of Thermal Injury: I. The Conduction of Heat to and through Skin and the Temperatures Attained Therein. A Theoretical and an Experimental Investigation" Am J Pathol 1947; 23: 530-49.
[263] Moritz AR, Henriques FC. "Studies of thermal injury: II. The relative importance of time and surface temperature in the causation of cutaneous burns" Am J Pathol 1947; 23: 695-720.
[264] Moritz AR, Henriques FC Jr. "The reciprocal relationship of surface temperature and time in the production of hyperthermic cutaneous injury" Am J Pathol 1947; 23: 897.
[265] Eyring H, Stearn AE. "The application of the theory of absolute reaction rates to proteins" Chemical Reviews 1939; 24: 253-70.
[266] Johnson FH, Eyring H, Stoner BJ. "Temperature" In: in The Theory of Rate Process in Biology and Medicine. New York: John Wiley & Sons, Inc 1974; pp. 549-98.
[267] Sapareto SA, Dewey WC. "Thermal dose determination in cancer therapy" Int J Radiat Oncol Biol Phys 1984; 10: 787-800.
[268] Pearce JA. "Models for Thermal Damage in Tissues: Processes and Applications" Crit Rev Biomed Eng 2010; 38: 1-20.
[269] Payne SJ, Pang T, O'Neill DP. "Mathematical modeling of thermal ablation" Crit Rev Biomed Eng 2010; 38: 21-30.
[270] Pitt RE. "Cryobiological implications of diff erent methods of calculating the chemical potential of water in partially frozen suspending media" Cryo Lett 1990; 11: 227-40.
[271] Hobbs F. Ice physics. London: Oxford University Press 1974.
[272] Toner M, Tompkins RG, Cravalho EG, Yarmush ML. "Transport phenomena during freezing of isolated hepatocytes" Aiche J 1992; 38: 1512-22.
[273] Yang G, Veres M, Szalai G, Zhang A, Xu LX, He X. "Biotransport phenomena in freezing mammalian oocytes" Ann Biomed Eng 2011; 39: 580-91.
[274] Toner M, Cravalho EG, Karel M, Armant DR. "Cryomicroscopic analysis of intracellular ice formation during freezing of mouse oocytes without cryoadditives" Cryobiology 1991; 28: 55-71.
[275] Diller KR. "Bioheat and mass transfer as viewed through a microscope" J Biomech Eng 2005; 127: 67-84.
[276] Chen C, Li WZ. "Diffusion controlled ice growth with soft impingement inside biological cells during freezing" Cryo letters 2008; 29: 371-81.
[277] Zhao G, Luo D, Gao D. "Universal model for intracellular ice formation and its growth" AIChE J 2006; 52: 2596-606.
[278] Karlsson JO, Cravalho EG, Borel Rinkes IH, Tompkins RG, Yarmush ML, Toner M. "Nucleation and growth of ice crystals inside cultured hepatocytes during freezing in the presence of dimethyl sulfoxide" Biophys J 1993; 65: 2524-36.
[279] Mazur P. "Kinetics of water loss from cells at subzero temperatures and the likelihood of intracellular freezing" J Gen Physiol 1963; 47: 347-69.
[280] Levin RL, Cravalho EG, Huggins CE. "A membrane model describing the effect of temperature on the water conductivity of erythrocyte membranes at subzero temperatures" Cryobiology 1976; 13: 415-29.
[281] Devireddy RV, Raha D, Bischof JC. "Measurement of water transport during freezing in cell suspensions using a differential scanning calorimeter" Cryobiology 1998; 36: 124-55.
[282] Devireddy RV, Coad JE, Bischof JC. "Microscopic and calorimetric assessment of freezing processes in uterine fibroid tumor tissue" Cryobiology 2001; 42: 225-43.
[283] Luo D, Han X, He L, et al. "A modified differential scanning calorimetry for determination of cell volumetric change during the freezing process" Cryo Lett 2002; 23: 229-36.
[284] Bischof JC. "Quantitative measurement and prediction of bio-physical response during freezing in tissues" Ann Rev Biomed Eng 2000; 2: 257-88.
[285] Han B, Bischof JC. "Engineering challenges in tissue preservation" Cell Preserv Technol 2004; 2: 91-112.
[286] Acker JP, McGann LE. "Innocuous intracellular ice improves survival of frozen cells" Cell Transplant 2002; 11: 563-71.
[287] Acker JP, McGann LE. "Protective effect of intracellular ice during freezing? Cryobiology 2003; 46: 197-202.
[288] Konstas AA, Neimark MA, Laine AF, Pile-Spellman J. "A theoretical model of selective cooling using intracarotid cold saline infusion in the human brain" J Appl Physiol 2007; 102: 1329-40.
[289] Xu X, Tikuisis P, Giesbrecht G. "A mathematical model for human brain cooling during cold-water near-drowning" J Appl Physiol 1999; 86: 265-72.
[290] King RWP, Trembly BS, Strohbehn JW. "The electromagnetic-field of an insulated antenna in a conducting or dielectric medium" IEEE Trans Microw Theory Tech 1983; 31: 574-83.
[291] Gentili GB, Gori F, Leoncini M. "Electromagnetic and thermal models of a water-cooled dipole radiating in a biological tissue" IEEE Trans Biomed Eng 1991; 38: 98-103.
[292] Enderle J, Blanchard SM, Bronzino J. Intro Biomedl Eng. Burlington, MA: Elsvier Academic Press 2005.
[293] Pennes HH. "Analysis of tissue and arterial blood temperatures in the resting human forearm" J Appl Physiol 1948; 1: 93-122.
[294] Chato J C. "Heat-Transfer to blood-vessels" J Biomech Eng 1980; 102: 110-8.
[295] Chen MM, Holmes KR. "Microvascular contributions in tissue heat transfer" Ann N Y Acad Sci 1980; 335: 137-50.
[296] Weinbaum S, Jiji LM, Lemons DE. "Theory and experiment for the effect of vascular microstructure on surface tissue heat transfer--Part I: Anatomical foundation and model conceptualization" J Biomech Eng 1984; 106: 321-0.
[297] Keller KH, Seiler L. "Analysis of peripheral heat transfer in man" J Appl Physiol 1971; 30: 779.
[298] Mitchell JW, Myers GE. "An analytical model of the counter-current heat exchange phenomena" Biophys J 1968; 8: 897-911.
[299] Mooibroek J, Lagendijk JJ. "A fast and simple algorithm for the calculation of convective heat transfer by large vessels in three-dimensional inhomogeneous tissues" IEEE Trans Biomed Eng 1991; 38: 490-501.
[300] Huang HW, Chen ZP, Roemer RB. "A counter current vascular network model of heat transfer in tissues" J Biomech Eng 1996; 118: 120-9.
[301] Weinbaum S, Jiji LM. "A new simplified bioheat equation for the effect of blood flow on local average tissue temperature" J Biomech Eng 1985; 107: 131-9.
[302] Wissler EH. "Comments on the new bioheat equation proposed by Weinbaum and Jiji" J Biomech Eng 1987; 109: 226-33.
[303] Shrivastava D, McKay B, Roemer RB. "An analytical study of heat transfer in finite tissue with two blood vessels and uniform dirichlet boundary conditions" J Heat Transfer 2005; 127: 179-88.
[304] Shrivastava D, Roemer R. "An analytical study of heat transfer in a finite tissue region with two blood vessels and general Dirichlet boundary conditions" Int J Heat Mass Transf 2005; 48: 4090-102.
[305] Shrivastava D, Roemer RB. "Readdressing the issue of thermally significant blood vessels using a countercurrent vessel network" J Biomech Eng 2006; 128: 210-6.
[306] Arkin H, Xu LX, Holmes KR. "Recent developments in modeling heat transfer in blood perfused tissues" IEEE Trans Biomed Eng 1994; 41: 97-107.
[307] Charney C K. "Mathematical models of bioheat transfer" Advances in Heat Tansfer 1992; 22: 19-155.
[308] Wissler EH. "Pennes' 1948 paper revisited" J Appl Physiol 1998; 85: 35-41.
[309] Nelson DA. "Invited editorial on "Pennes' 1948 paper revisited" J Appl Physiol 1998; 85: 2-3.
[310] Roemer RB. "Engineering aspects of hyperthermia therapy" Ann Rev Biomed Eng 1999; 1: 347-76.
[311] Kolios MC, Worthington AE, Sherar MD, Hunt JW. "Experimental evaluation of two simple thermal models using transient temperature analysis" Phys Med Biol 1998; 43: 3325-40.
[312] Song CW. "Effect of local hyperthermia on blood flow and microenvironment: a review" Cancer Res 1984; 44: 4721s-30s.
[313] Vaupel P. "Pathophysiological mechanisms of hyperthermia in cancer therapy" In: Gautherie M, Ed. in Biological basis of oncologic thermotherapy. Berlin: Springer-Verlag 1990.
[314] Dewhirst M, Sim D, Gross J, MA K. "Effect of heating rate on tumour and normal tissue microcirculatory function" In: Overgaard J, Ed. in Hyperthermic Oncology. London, UK: Taylor & Francis 1984.
[315] Schutt DJ, Haemmerich D. "Effects of variation in perfusion rates and of perfusion models in computational models of radio frequency tumor ablation" Med Phys 2008; 35: 3462-70.
[316] Abraham JP, Sparrow EM. "A thermal-ablation bioheat model including liquid-to-vapor phase change, pressure- and necrosis-dependent perfusion, and moisture-dependent properties" Int J Heat and Mass Transf 2007; 50: 2537-44.
[317] Skinner MG, Iizuka MN, Kolios MC, Sherar MD. "A theoretical comparison of energy sources--microwave, ultrasound and laser--for interstitial thermal therapy" Phys Med Biol 1998; 43: 3535-47.
[318] Brown SL, Hunt JW, Hill RP. "Differential thermal sensitivity of tumour and normal tissue microvascular response during hyperthermia" Int J Hyperther 1992; 8: 501-14.
[319] Dudar TE, Jain RK. "Differential response of normal and tumor microcirculation to hyperthermia" Cancer Res 1984; 44: 605-12.
[320] Acker JC, Dewhirst MW, Honore GM, Samulski TV, Tucker JA, Oleson JR. "Blood perfusion measurements in human tumours: evaluation of laser Doppler methods" Int J Hyperther 1990; 6: 287-304.
[321] Alexiades V, Solomon AD. Mathematical modeling of melting and freezing processes. Washington: Hemisphere Publishing Corp 1993.
[322] Silbey RJ, Alberty RA, Bawendi MG. Physical Chemistry. 4th edition. Hoboken, NJ: John Wiley & Sons, Inc. 2004, 00: 0 2004.
[323] Devireddy RV, Smith DJ, Bischof JC. "Effect of microscale mass transport and phase change on numerical prediction of freezing in biological tissues" J Heat Transfer 2002; 124: 365-74.
[324] Han B, Bischof J C. "Thermodynamic nonequilibrium phase change behavior and thermal properties of biological solutions for cryobiology applications" J Biomech Eng 2004; 126: 196-203.
[325] Han B, Choi JH, Dantzig JA, Bischof JC. "A quantitative analysis on latent heat of an aqueous binary mixture" Cryobiology 2006; 52: 146-51.
[326] Duck FA. Physical properties of tissue: a comprehensive reference book. San Diego, CA: Academic Press 1990.
[327] Choi J, Bischof JC. "Review of biomaterial thermal property measurements in the cryogenic regime and their use for prediction of equilibrium and non-equilibrium freezing applications in cryobiology" Cryobiology 2010; 60: 52-70.
[328] Hoffmann NE, Bischof JC. "Cryosurgery of normal and tumor tissue in the dorsal skin flap chamber: Part I--thermal response" J Biomech Eng 2001; 123: 301-09.