Possible mechanisms and simulation modeling of FLASH radiotherapy

FLASH radiotherapy (FLASH-RT) has great potential to improve patient outcomes. It delivers radiation doses at an ultra-high dose rate (UHDR: ≥ 40 Gy/s) in a single instant or a few pulses. Much higher irradiation doses can be administered to tumors with FLASH-RT than with conventional dose rate (0.01-0.40 Gy/s) radiotherapy. UHDR irradiation can suppress toxicity in normal tissues while sustaining antitumor efficiency, which is referred to as the FLASH effect. However, the mechanisms underlying the effects of the FLASH remain unclear. To clarify these mechanisms, the development of simulation models that can contribute to treatment planning for FLASH-RT is still underway. Previous studies indicated that transient oxygen depletion or augmented reactions between secondary reactive species produced by irradiation may be involved in this process. To discuss the possible mechanisms of the FLASH effect and its clinical potential, we summarized the physicochemical, chemical, and biological perspectives as well as the development of simulation modeling for FLASH-RT.

A step-by-step simulation code for estimating yields of water radiolysis species based on electron track-structure mode in the PHITS code

Objective. Time-dependent yields of chemical products resulting from water radiolysis play a great role in evaluating DNA damage response after exposure to ionizing radiation. Particle and Heavy Ion Transport code System (PHITS) is a general-purpose Monte Carlo simulation code for radiation transport, which simulates atomic interactions originating from discrete energy levels of ionizations and electronic excitations as well as molecular excitations as physical stages. However, no chemical code for simulating water radiolysis products exists in the PHITS package.Approach.Here, we developed a chemical simulation code dedicated to the PHITS code, hereafter calledPHITS-Chemcode, which enables the calculation of theGvalues of water radiolysis species (•OH, eaq-, H2, H2O2etc) by electron beams.Main results.The estimatedGvalues during 1 μs are in agreement with the experimental ones and other simulations. ThisPHITS-Chemcode also simulates the radiolysis in the presence of OH radical scavengers, such as tris(hydroxymethyl)aminomethane and dimethyl sulfoxide. Thank to this feature, the contributions of direct and indirect effects on DNA damage induction under various scavenging capacities can be analyzed.Significance.This chemical code coupled with PHITS could contribute to elucidating the mechanism of radiation effects by connecting physical, physicochemical, and chemical processes.

Modeling for predicting survival fraction of cells after ultra-high dose rate irradiation

Objective. FLASH radiotherapy (FLASH-RT) with ultra-high dose rate (UHDR) irradiation (i.e. > 40 Gy s-1) spares the function of normal tissues while preserving antitumor efficacy, known as the FLASH effect. The biological effects after conventional dose rate-radiotherapy (CONV-RT) with ≤0.1 Gy s-1have been well modeled by considering microdosimetry and DNA repair processes, meanwhile modeling of radiosensitivities under UHDR irradiation is insufficient. Here, we developed anintegrated microdosimetric-kinetic(IMK)model for UHDR-irradiationenabling the prediction of surviving fraction after UHDR irradiation.Approach.TheIMK model for UHDR-irradiationconsiders the initial DNA damage yields by the modification of indirect effects under UHDR compared to CONV dose rate. The developed model is based on the linear-quadratic (LQ) nature with the dose and dose square coefficients, considering the reduction of DNA damage yields as a function of dose rate.Main results.The estimate by the developed model could successfully reproduce thein vitroexperimental dose-response curve for various cell line types and dose rates.Significance.The developed model would be useful for predicting the biological effects under the UHDR irradiation.

First-principles simulation of an ejected electron produced by monochromatic deposition energy to water at the femtosecond order

This study uses a time-dependent first-principles simulation code to investigate the transient dynamics of an ejected electron produced in the monochromatic deposition energy from 11 to 19 eV in water. The energy deposition forms a three-body single spur comprising a hydroxyl radical (OH˙), hydronium ion (H3O+), and hydrated electron (eaq-). The earliest formation involves electron thermalization and delocalization dominated by the molecular excitation of water. Our simulation results show that the transient electron dynamics primarily depends on the amount of deposition energy to water; the thermalization time varies from 200 to 500 fs, and the delocalization varies from 3 to 10 nm in this energy range. These features are crucial for determining the earliest single-spur formation and facilitating a sequential simulation from an energy deposition to a chemical reaction in water photolysis or radiolysis. The spur radius obtained from the simulation correlates reasonably with the experimental-based estimations. Our results should provide universalistic insights for analysing ultrafast phenomena dominated by the molecular excitation of water in the femtosecond order.

Improvement of the hybrid approach between Monte Carlo simulation and analytical function for calculating microdosimetric probability densities in macroscopic matter

OBJECTIVE

Estimate of the probability density of the microdosimetric quantities in macroscopic matter is indispensable for applying the concept of microdosimetry to medical physics and radiological protection. Particle and Heavy Ion Transport code System (PHITS) enables estimating the microdosimetric probability densities due to its unique hybrid modality between the Monte Carlo and analytical approaches called the microdosimetric function. It can convert the deposition energies calculated by the macroscopic Monte Carlo radiation transport simulation to microdosimetric probability densities in water using an analytical function based on the track-structure simulations.

METHODOLOGY

In this study, we improved this function using the latest track-structure simulation codes implemented in PHITS. The improved function is capable of calculating the probability densities of not only the conventional microdosimetric quantities such as lineal energy but also the numbers of ionization events occurred in a target site, the so-called ionization cluster size distribution, for arbitrary site diameters from 3 nm to 1 μm.

MAIN RESULTS

The accuracy of the improved function was well verified by comparing the microdosimetric probability densities measured by tissue-equivalent proportional counters with the corresponding data calculated in this study. Test calculations for clonogenic cell survival using the improved function coupled with the modified microdosimetric kinetic model suggested a slight increase of its relative biological effectiveness compared with our previous estimations. As a new application of the improved function, we calculated the relative biological effectiveness of the single-strand break and double-strand break yields for proton irradiations using the updated PHITS coupled with the simplified DNA damage estimation model, and confirmed its equivalence in accuracy and its superiority in computational time compared to our previously proposed method based on the track-structure simulation.

SIGNIFICANCE

From these features, we concluded that the improved function could expand the application fields of PHITS by bridging the gap between microdosimetry and macrodosimetry.

Modelling oxygen effects on the in- and out-of-field radiosensitivity of cells exposed to intensity-modulated radiation fields

OBJECTIVE

The delivery of intensity-modulated radiation fields has improved the conformity of dose to tumour targets during radiotherapy (RT). Previously, it has been shown that intercellular communication between cells positioned in- and outside of the radiation field impacts cellular radiosensitivity under hypoxic and normoxic conditions. However, the mechanism of intercellular communication in hypoxia remains to be fully understood. In this study, the cell-killing effects of intercellular communication in hypoxia were modelled in an effort to better understand the underlying mechanisms of response. 
Approach: By irradiating a 50% area of the culture dish (half-field exposure), experimental dose-response curves for cell survival and residual DNA double-strand breaks (DSBs) were generated in prostate (DU145) and non-small cell lung cancer (H1299) cells. The oxygen enhancement ratio (OER) was determined from early DSB yields (corresponding to relative direct damage) and used to model the in- and out-of-field radiosensitivity. 
Main results: The developed integrated microdosimetric-kinetic (IMK) model successfully predicted the experimental dose responses for survival and lethal lesions, and provides a mechanistic interpretation that the probability of hits for releasing cell-killing signals is dependent on oxygen. This experimental and modelling study also suggests that residual DSBs correspond to logarithmic survival fraction (meaning lethal lesions) for in- and out-of-field cells. Our data suggest that the OER value determined using uniform-field exposure can be applied to predict the in- and out-of-field radiosensitivity of cells following exposure to intensity modulated beams. 
Significance: The developed IMK model facilitates a more precise understanding of intercellular signalling following exposure to intensity-modulated radiation fields.&#xD.

Initial yield of hydrated electron production from water radiolysis based on first-principles calculation

Many scientific insights into water radiolysis have been applied for developing life science, including radiation-induced phenomena, such as DNA damage and mutation induction or carcinogenesis. However, the generation mechanism of free radicals due to radiolysis remains to be fully understood. Consequently, we have encountered a crucial problem in that the initial yields connecting radiation physics to chemistry must be parameterized. We have been challenged in the development of a simulation tool that can unravel the initial free radical yields, from physical interaction by radiation. The presented code enables the first-principles calculation of low energy secondary electrons resulting from the ionization, in which the secondary electron dynamics are simulated while considering dominant collision and polarization effects in water. In this study, using this code, we predicted the yield ratio between ionization and electronic excitation from a delocalization distribution of secondary electrons. The simulation result presented a theoretical initial yield of hydrated electrons. In radiation physics, the initial yield predicted from parameter analysis of radiolysis experiments in radiation chemistry was successfully reproduced. Our simulation code helps realize a reasonable spatiotemporal connection from radiation physics to chemistry, which would contribute to providing new scientific insights for precise understanding of underlying mechanisms of DNA damage induction.

Translational study for stereotactic body radiotherapy against non-small cell lung cancer, including oligometastases, considering cancer stem-like cells enable predicting clinical outcome from in vitro data

BACKGROUND

Curative effects of stereotactic body radiotherapy (SBRT) for non-small cell lung cancer (NSCLC) have been evaluated using various biophysical models. Because such model parameters are empirically determined based on clinical experience, there is a large gap between in vitro and clinical studies. In this study, considering the heterogeneous cell population, we performed a translational study to realize the possible linkage based on a modeling approach.

METHODS

We modeled cell-killing and tumor control probability (TCP) considering two populations: progeny and cancer stem-like cells. The model parameters were determined from in vitro survival data of A549 and EBC-1 cells. Based on the cellular parameters, we predicted TCP and compared it with the corresponding clinical data from 553 patients collected at Hirosaki University Hospital.

RESULTS

Using an all-in-one developed model, the so-called integrated microdosimetric-kinetic (IMK) model, we successfully reproduced both in vitro survival after acute irradiation and the 3-year TCP with various fractionation schemes (6-10 Gy per fraction). From the conventional prediction without considering cancer stem cells (CSCs), this study revealed that radioresistant CSCs play a key role in the linkage between in vitro and clinical outcomes.

CONCLUSIONS

This modeling study provides a possible generalized biophysical model that enables precise estimation of SBRT worldwide.

An Analytical Method for Quantifying the Yields of DNA Double-Strand Breaks Coupled with Strand Breaks by γ-H2AX Focus Formation Assay Based on Track-Structure Simulation

Complex DNA double-strand break (DSB), which is defined as a DSB coupled with additional strand breaks within 10 bp in this study, induced after ionizing radiation or X-rays, is recognized as fatal damage which can induce cell death with a certain probability. In general, a DSB site inside the nucleus of live cells can be experimentally detected using the γ-H2AX focus formation assay. DSB complexity is believed to be detected by analyzing the focus size using such an assay. However, the relationship between focus size and DSB complexity remains uncertain. In this study, using Monte Carlo (MC) track-structure simulation codes, i.e., an in-house WLTrack code and a Particle and Heavy Ion Transport code System (PHITS), we developed an analytical method for qualifying the DSB complexity induced by photon irradiation from the microscopic image of γ-H2AX foci. First, assuming that events (i.e., ionization and excitation) potentially induce DNA strand breaks, we scored the number of events in a water cube (5.03 × 5.03 × 5.03 nm³) along electron tracks. Second, we obtained the relationship between the number of events and the foci size experimentally measured by the γ-H2AX focus formation assay. Third, using this relationship, we evaluated the degree of DSB complexity induced after photon irradiation for various X-ray spectra using the foci size, and the experimental DSB complexity was compared to the results estimated by the well-verified DNA damage estimation model in the PHITS code. The number of events in a water cube was found to be proportional to foci size, suggesting that the number of events intrinsically related to DSB complexity at the DNA scale. The developed method was applicable to focus data measured for various X-ray spectral situations (i.e., diagnostic kV X-rays and therapeutic MV X-rays). This method would contribute to a precise understanding of the early biological impacts of photon irradiation by means of the γ-H2AX focus formation assay.

Application of a simple DNA damage model developed for electrons to proton irradiation

Proton beam therapy allows irradiating tumor volumes with reduced side effects on normal tissues with respect to conventional x-ray radiotherapy. Biological effects such as cell killing after proton beam irradiations depend on the proton kinetic energy, which is intrinsically related to early DNA damage induction. As such, DNA damage estimation based on Monte Carlo simulations is a research topic of worldwide interest. Such simulation is a mean of investigating the mechanisms of DNA strand break formations. However, past modellings considering chemical processes and DNA structures require long calculation times. Particle and heavy ion transport system (PHITS) is one of the general-purpose Monte Carlo codes that can simulate track structure of protons, meanwhile cannot handle radical dynamics simulation in liquid water. It also includes a simple model enabling the efficient estimation of DNA damage yields only from the spatial distribution of ionizations and excitations without DNA geometry, which was originally developed for electron track-structure simulations. In this study, we investigated the potential application of the model to protons without any modification. The yields of single-strand breaks, double-strand breaks (DSBs) and the complex DSBs were assessed as functions of the proton kinetic energy. The PHITS-based estimation showed that the DSB yields increased as the linear energy transfer (LET) increased, and reproduced the experimental and simulated yields of various DNA damage types induced by protons with LET up to about 30 keV μm−1. These results suggest that the current DNA damage model implemented in PHITS is sufficient for estimating DNA lesion yields induced after protons irradiation except at very low energies (below 1 MeV). This model contributes to evaluating early biological impacts in radiation therapy.

Impact of the Lorentz force on electron track structure and early DNA damage yields in magnetic resonance-guided radiotherapy

Magnetic resonance-guided radiotherapy (MRgRT) has been developed and installed in recent decades for external radiotherapy in several clinical facilities. Lorentz forces modulate dose distribution by charged particles in MRgRT; however, the impact of Lorentz forces on low-energy electron track structure and early DNA damage induction remain unclear. In this study, we estimated features of electron track structure and biological effects in a static magnetic field (SMF) using a general-purpose Monte Carlo code, particle and heavy ion transport code system (PHITS) that enables us to simulate low-energy electrons down to 1 meV by track-structure mode. The macroscopic dose distributions by electrons above approximately 300 keV initial energy in liquid water are changed by both perpendicular and parallel SMFs against the incident direction, indicating that the Lorentz force plays an important role in calculating dose within tumours. Meanwhile, DNA damage estimation based on the spatial patterns of atomic interactions indicates that the initial yield of DNA double-strand breaks (DSBs) is independent of the SMF intensity. The DSB induction is predominantly attributed to the secondary electrons below a few tens of eV, of which energy deposition patterns are not considerably affected by the Lorentz force. Our simulation study suggests that treatment planning for MRgRT can be made with consideration of only changed dose distribution.

Inflammatory Signaling and DNA Damage Responses after Local Exposure to an Insoluble Radioactive Microparticle

Simple Summary

A cesium-bearing microparticle (Cs-BMP) is an insoluble radioactive microparticle possessing high specific radioactivity, which was discovered after the incident at the Fukushima nuclear power plant. Due to their insoluble nature, such Cs-BMPs are assumed to adhere in the long term to normal tissue, leading to chronic local exposure. However, radiation risk due to the intake of internal exposure to radioactive cesium is conventionally estimated from the organ dose given by uniform exposure to soluble cesium. As such, it is critical to clarify the normal tissue effects posed by heterogeneous exposure to Cs-BMPs. This in vitro study reports on the relationship between the inflammatory responses and DNA damage induction during local exposure to a Cs-BMP.

Abstract

Cesium-bearing microparticles (Cs-BMPs) can reach the human respiratory system after inhalation, resulting in chronic local internal exposure. We previously investigated the spatial distribution of DNA damage induced in areas around a Cs-BMP; however, the biological impacts have not been fully clarified due to the limited amount of data. Here, we investigated the inflammatory signaling and DNA damage responses after local exposure to a Cs-BMP in vitro. We used two normal human lung cell lines, i.e., lung fibroblast cells (WI-38) and bronchial epithelial cells (HBEC3-KT). After 24 h exposure to a Cs-BMP, inflammation was evaluated by immunofluorescent staining for nuclear factor κB (NF-κB) p65 and cyclooxygenase 2 (COX-2). The number of DNA double-strand breaks (DSBs) was also detected by means of phospholylated histone H2AX (γ-H2AX) focus formation assay. Cs-BMP exposure significantly increased NF-κB p65 and COX-2 expressions, which were related to the number of γ-H2AX foci in the cell nuclei. Compared to the uniform (external) exposure to ¹³⁷Cs γ-rays, NF-κB tended to be more activated in the cells proximal to the Cs-BMP, while both NF-κB p65 and COX-2 were significantly activated in the distal cells. Experiments with chemical inhibitors for NF-κB p65 and COX-2 suggested the involvement of such inflammatory responses both in the reduced radiosensitivity of the cells proximal to Cs-BMP and the enhanced radiosensitivity of the cells distal from Cs-BMP. The data show that local exposure to Cs-BMP leads to biological effects modified by the NF-κB pathway, suggesting that the radiation risk for Cs-BMP exposure can differ from that estimated based on conventional uniform exposure to normal tissues.

Track-structure modes in particle and heavy ion transport code system (PHITS): application to radiobiological research

PURPOSE

In radiation physics, Monte Carlo radiation transport simulations are powerful tools to evaluate the cellular responses after irradiation. When investigating such radiation-induced biological effects, it is essential to perform track structure simulations by explicitly considering each atomic interaction in liquid water at the sub-cellular and DNA scales. The Particle and Heavy-Ion Transport code System (PHITS) is a Monte Carlo code which enables to calculate track structure at DNA scale by employing the track-structure modes for electrons, protons and carbon ions. In this paper, we review the recent development status and future prospects of the track-structure modes in the PHITS code.

CONCLUSIONS

To date, the physical features of these modes have been verified using the available experimental data and Monte Carlo simulation results reported in literature. These track-structure modes can be used for calculating microdosimetric distributions to estimate cell survival and for estimating initial DNA damage yields. The use of PHITS track-structure mode is expected not only to clarify the underlying mechanisms of radiation effects but also to predict curative effects in radiation therapy. The results of PHITS simulations coupled with biophysical models will contribute to the radiobiological studies by precisely predicting radiation-induced biological effects based on the Monte Carlo approach.

4-Methylumbelliferone administration enhances radiosensitivity of human fibrosarcoma by intercellular communication

Hyaluronan synthesis inhibitor 4-methylumbelliferone (4-MU) is a candidate of radiosensitizers which enables both anti-tumour and anti-metastasis effects in X-ray therapy. The curative effects under such 4-MU administration have been investigated in vitro; however, the radiosensitizing mechanisms remain unclear. Here, we investigated the radiosensitizing effects under 4-MU treatment from cell experiments and model estimations. We generated experimental surviving fractions of human fibrosarcoma cells (HT1080) after 4-MU treatment combined with X-ray irradiation. Meanwhilst, we also modelled the pharmacological effects of 4-MU treatment and theoretically analyzed the synergetic effects between 4-MU treatment and X-ray irradiation. The results show that the enhancement of cell killing by 4-MU treatment is the greatest in the intermediate dose range of around 4 Gy, which can be reproduced by considering intercellular communication (so called non-targeted effects) through the model analysis. As supposed to be the involvement of intercellular communication in radiosensitization, the oxidative stress level associated with reactive oxygen species (ROS), which leads to DNA damage induction, is significantly higher by the combination of 4-MU treatment and irradiation than only by X-ray irradiation, and the radiosensitization by 4-MU can be suppressed by the ROS inhibitors. These findings suggest that the synergetic effects between 4-MU treatment and irradiation are predominantly attributed to intercellular communication and provide more efficient tumour control than conventional X-ray therapy.

Oxygen enhancement ratios of cancer cells after exposure to intensity modulated x-ray fields: DNA damage and cell survival

Hypoxic cancer cells within solid tumours show radio-resistance, leading to malignant progression in fractionated radiotherapy. When prescribing dose to tumours under heterogeneous oxygen pressure with intensity-modulated radiation fields, intercellular signalling could have an impact on radiosensitivity between in-field and out-of-field (OF) cells. However, the impact of hypoxia on radio-sensitivity under modulated radiation intensity remains to be fully clarified. Here, we investigate the impact of hypoxia on in-field and OF radio-sensitivities using two types of cancer cells, DU145 and H1299. Using a nBIONIX hypoxic culture kit and a shielding technique to irradiate 50% of a cell culture flask, oxygen enhancement ratios for double-strand breaks (DSB) and cell death endpoints were determined. Thesein vitromeasurements indicate that hypoxia impacts OF cells, although the hypoxic impacts on OF cells for cell survival were dose-dependent and smaller compared to those for in-field and uniformly irradiated cells. These decreased radio-sensitivities of OF cells were shown as a consistent tendency for both DSB and cell death endpoints, suggesting that radiation-induced intercellular communication is of importance in advanced radiotherapy dose-distributions such as with intensity-modulated radiotherapy.

Verification of KURBUC-based ion track structure mode for proton and carbon ions in the PHITS code

The particle and heavy ion transport code system (PHITS) is a general-purpose Monte Carlo radiation transport simulation code. It has the ability to handle diverse particle types over a wide range of energy. The latest PHITS development enables the generation of track structure for proton and carbon ions (¹H⁺, ¹²C⁶⁺) based on the algorithms in the KURBUC code, which is considered as one of the most verified track-structure codes worldwide. This ion track-structure mode is referred to as the PHITS-KURBUC mode. In this study, the range, radial dose distributions, and microdosimetric distributions were calculated using the PHITS-KURBUC mode. Subsequently, they were compared with the corresponding data obtained from the original KURBUC and from other studies. These comparative studies confirm the successful inclusion of the KURBUC code in the PHITS code. As results of the synergistic effect between the macroscopic and microscopic radiation transport codes, this implementation enabled the detailed calculation of the microdosimetric and nanodosimetric quantities under complex radiation fields, such as proton beam therapy with the spread-out Bragg peak.

Development of a new microdosimetric biological weighting function for the RBE10 assessment in case of the V79 cell line exposed to ions from 1H to 238U

An improved biological weighting function (IBWF) is proposed to phenomenologically relate microdosimetric lineal energy probability density distributions with the relative biological effectiveness (RBE) for the in vitro clonogenic cell survival (surviving fraction = 10%) of the most commonly used mammalian cell line, i.e. the Chinese hamster lung fibroblasts (V79). The IBWF, intended as a simple and robust tool for a fast RBE assessment to compare different exposure conditions in particle therapy beams, was determined through an iterative global-fitting process aimed to minimize the average relative deviation between RBE calculations and literature in vitro data in case of exposure to various types of ions from ¹H to ²³⁸U. By using a single particle- and energy- independent function, it was possible to establish an univocal correlation between lineal energy and clonogenic cell survival for particles spanning over an unrestricted linear energy transfer range of almost five orders of magnitude (0.2 keV µm^-1 to 15 000 keV µm^-1 in liquid water). The average deviation between IBWF-derived RBE values and the published in vitro data was ∼14%. The IBWF results were also compared with corresponding calculations (in vitro RBE₁₀ for the V79 cell line) performed using the modified microdosimetric kinetic model (modified MKM). Furthermore, RBE values computed with the reference biological weighting function (BWF) for the in vivo early intestine tolerance in mice were included for comparison and to further explore potential correlations between the BWF results and the in vitro RBE as reported in previous studies. The results suggest that the modified MKM possess limitations in reproducing the experimental in vitro RBE₁₀ for the V79 cell line in case of ions heavier than ²⁰Ne. Furthermore, due to the different modelled endpoint, marked deviations were found between the RBE values assessed using the reference BWF and the IBWF for ions heavier than ²H. Finally, the IBWF was unchangingly applied to calculate RBE values by processing lineal energy density distributions experimentally measured with eight different microdosimeters in 19 ¹H and ¹²C beams at ten different facilities (eight clinical and two research ones). Despite the differences between the detectors, irradiation facilities, beam profiles (pristine or spread out Bragg peak), maximum beam energy, beam delivery (passive or active scanning), energy degradation system (water, PMMA, polyamide or low-density polyethylene), the obtained IBWF-based RBE trends were found to be in good agreement with the corresponding ones in case of computer-simulated microdosimetric spectra (average relative deviation equal to 0.8% and 5.7% for ¹H and ¹²C ions respectively).

Implications of radiation microdosimetry for accelerator-based boron neutron capture therapy: a radiobiological perspective

Boron neutron capture therapy (BNCT) has great potential to selectively destroy cancer cells while sparing surrounding normal cells. The basic concept of BNCT was developed in the 1930s, but it has not yet been commonly used in clinical practice, even though there is now a large number of experimental and translational studies demonstrating its marked therapeutic potential. With the development of neutron accelerators that can be installed in medical institutions, accelerator-based BNCT is expected to become available at several medical institutes around the world in the near future. In this commentary, from the point of view of radiation microdosimetry, we discuss the biological effects of BNCT, especially the underlying mechanisms of compound biological effectiveness. Radiobiological perspectives provide insight into the effectiveness of BNCT in creating a synergy effect in the field of clinical oncology.

A Model for Estimating Dose-Rate Effects on Cell-Killing of Human Melanoma after Boron Neutron Capture Therapy

Boron neutron capture therapy (BNCT) is a type of radiation therapy for eradicating tumor cells through a ¹⁰B(n,α)⁷Li reaction in the presence of ¹⁰B in cancer cells. When delivering a high absorbed dose to cancer cells using BNCT, both the timeline of ¹⁰B concentrations and the relative long dose-delivery time compared to photon therapy must be considered. Changes in radiosensitivity during such a long dose-delivery time can reduce the probability of tumor control; however, such changes have not yet been evaluated. Here, we propose an improved integrated microdosimetric-kinetic model that accounts for changes in microdosimetric quantities and dose rates depending on the ¹⁰B concentration and investigate the cell recovery (dose-rate effects) of melanoma during BNCT irradiation. The integrated microdosimetric–kinetic model used in this study considers both sub-lethal damage repair and changes in microdosimetric quantities during irradiation. The model, coupled with the Monte Carlo track structure simulation code of the Particle and Heavy Ion Transport code System, shows good agreement with in vitro experimental data for acute exposure to ⁶⁰Co γ-rays, thermal neutrons, and BNCT with ¹⁰B concentrations of 10 ppm. This indicates that microdosimetric quantities are important parameters for predicting dose-response curves for cell survival under BNCT irradiations. Furthermore, the model estimation at the endpoint of the mean activation dose exhibits a reduced impact of cell recovery during BNCT irradiations with high linear energy transfer (LET) compared to ⁶⁰Co γ-rays irradiation with low LET. Throughout this study, we discuss the advantages of BNCT for enhancing the killing of cancer cells with a reduced dose-rate dependency. If the neutron spectrum and the timelines for drug and dose delivery are provided, the present model will make it possible to predict radiosensitivity for more realistic dose-delivery schemes in BNCT irradiations.

A theoretical cell-killing model to evaluate oxygen enhancement ratios at DNA damage and cell survival endpoints in radiation therapy

Radio-resistance induced under low oxygen pressure plays an important role in malignant progression in fractionated radiotherapy. For the general approach to predict cell killing under hypoxia, cell-killing models (e.g. the Linear-Quadratic model) have to be fitted to in vitro experimental survival data for both normoxia and hypoxia to obtain the oxygen enhancement ratio (OER). In such a case, model parameters for every oxygen condition needs to be considered by model-fitting approaches. This is inefficient for fractionated irradiation planning. Here, we present an efficient model for fractionated radiotherapy the integrated microdosimetric-kinetic model including cell-cycle distribution and the OER at DNA double-strand break endpoint (OER_DSB). The cell survival curves described by this model can reproduce the in vitro experimental survival data for both acute and chronic low oxygen concentrations. The OER_DSB used for calculating cell survival agrees well with experimental DSB ratio of normoxia to hypoxia. The important parameters of the model are oxygen pressure and cell-cycle distribution, which enables us to predict cell survival probabilities under chronic hypoxia and chronic anoxia. This work provides biological effective dose (BED) under various oxygen conditions including its uncertainty, which can contribute to creating fractionated regimens for multi-fractionated radiotherapy. If the oxygen concentration in a tumor can be quantified by medical imaging, the present model will make it possible to estimate the cell-killing and BED under hypoxia in more realistic intravital situations.

A Simplified Cluster Analysis of Electron Track Structure for Estimating Complex DNA Damage Yields

Complex DNA damage, defined as at least two vicinal lesions within 10-20 base pairs (bp), induced after exposure to ionizing radiation, is recognized as fatal damage to human tissue. Due to the difficulty of directly measuring the aggregation of DNA damage at the nano-meter scale, many cluster analyses of inelastic interactions based on Monte Carlo simulation for radiation track structure in liquid water have been conducted to evaluate DNA damage. Meanwhile, the experimental technique to detect complex DNA damage has evolved in recent decades, so both approaches with simulation and experiment get used for investigating complex DNA damage. During this study, we propose a simplified cluster analysis of ionization and electronic excitation events within 10 bp based on track structure for estimating complex DNA damage yields for electron and X-ray irradiations. We then compare the computational results with the experimental complex DNA damage coupled with base damage (BD) measured by enzymatic cleavage and atomic force microscopy (AFM). The computational results agree well with experimental fractions of complex damage yields, i.e., single and double strand breaks (SSBs, DSBs) and complex BD, when the yield ratio of BD/SSB is assumed to be 1.3. Considering the comparison of complex DSB yields, i.e., DSB + BD and DSB + 2BD, between simulation and experimental data, we find that the aggregation degree of the events along electron tracks reflects the complexity of induced DNA damage, showing 43.5% of DSB induced after 70 kVp X-ray irradiation can be classified as a complex form coupled with BD. The present simulation enables us to quantify the type of complex damage which cannot be measured through in vitro experiments and helps us to interpret the experimental detection efficiency for complex BD measured by AFM. This simple model for estimating complex DNA damage yields contributes to the precise understanding of the DNA damage complexity induced after X-ray and electron irradiations.

Intensity Modulated Radiation Fields Induce Protective Effects and Reduce Importance of Dose-Rate Effects

In advanced radiotherapy, intensity modulated radiation fields and complex dose-delivery are utilized to prescribe higher doses to tumours. Here, we investigated the impact of modulated radiation fields on radio-sensitivity and cell recovery during dose delivery. We generated experimental survival data after single-dose, split-dose and fractionated irradiation in normal human skin fibroblast cells (AGO1522) and human prostate cancer cells (DU145). The dose was delivered to either 50% of the area of a T25 flask containing the cells (half-field) or 100% of the flask (uniform-field). We also modelled the impact of dose-rate effects and intercellular signalling on cell-killing. Applying the model to the survival data, it is found that (i) in-field cell survival under half-field exposure is higher than uniform-field exposure for the same delivered dose; (ii) the importance of sub-lethal damage repair (SLDR) in AGO1522 cells is reduced under half-field exposure; (iii) the yield of initial DNA lesions measured with half-field exposure is smaller than that with uniform-field exposure. These results suggest that increased cell survival under half-field exposure is predominantly attributed not to rescue effects (increased SLDR) but protective effects (reduced induction of initial DNA lesions). In support of these protective effects, the reduced DNA damage leads to modulation of cell-cycle dynamics, i.e., less G₁ arrest 6 h after irradiation. These findings provide a new understanding of the impact of dose-rate effects and protective effects measured after modulated field irradiation.

Analysis of the high-dose-range radioresistance of prostate cancer cells, including cancer stem cells, based on a stochastic model

In radiotherapy, cancer stem cells (CSCs) are well recognized as one of the radioresistant cell types. Even in a small subpopulation, CSCs may have an influence on tumor control probability, represented by cell killing after irradiation. However, the relationship between the percentage content of CSCs and the cell survival dose-response curve has not yet been quantitatively clarified. In this study, we developed a cell-killing model for two cell populations (CSCs and progeny cells) to predict the surviving fractions, and compared it with the conventional linear-quadratic (LQ) model. Three prostate cancer cell lines (DU145, PC3 and LNCaP) were exposed to X-rays at doses ranging from 0 to 10 Gy. After the irradiation, we performed clonogenic survival assays to generate the cell survival curves, and carried out flow-cytometric analyses to estimate the percentage content of CSCs for each cell line. The cell survival curves for DU145 cells and PC3 cells seemed not to follow the conventional LQ model in the high dose range (>8 Gy). However, the outputs of the developed model agreed better with the experimental cell survival curves than those of the LQ model. The percentage content of CSCs predicted by the developed model was almost coincident with the measured percentage content for both DU145 cells and PC3 cells. The experiments and model analyses indicate that a small subpopulation of radioresistant CSCs has lower radiosensitivity in the high-dose range, which may lessen the clinical outcome for patients with prostate cancer after high-dose radiation therapy.

Investigation of dose-rate effects and cell-cycle distribution under protracted exposure to ionizing radiation for various dose-rates

During exposure to ionizing radiation, sub-lethal damage repair (SLDR) competes with DNA damage induction in cultured cells. By virtue of SLDR, cell survival increases with decrease of dose-rate, so-called dose-rate effects (DREs). Here, we focused on a wide dose-rate range and investigated the change of cell-cycle distribution during X-ray protracted exposure and dose-response curves via hybrid analysis with a combination of in vitro experiments and mathematical modelling. In the course of flow-cytometric cell-cycle analysis and clonogenic assays, we found the following responses in CHO-K1 cells: (1) The fraction of cells in S phase gradually increases during 6 h exposure at 3.0 Gy/h, which leads to radio-resistance. (2) Slight cell accumulation in S and G2/M phases is observed after exposure at 6.0 Gy/h for more than 10 hours. This suggests that an increase of SLDR rate for cells in S phase during irradiation may be a reproducible factor to describe changes in the dose-response curve at dose-rates of 3.0 and 6.0 Gy/h. By re-evaluating cell survival for various dose-rates of 0.186-60.0 Gy/h considering experimental-based DNA content and SLDR, it is suggested that the change of S phase fraction during irradiation modulates the dose-response curve and is possibly responsible for some inverse DREs.

Investigation of dose-rate effects and cell-cycle distribution under protracted exposure to ionizing radiation for various dose-rates

During exposure to ionizing radiation, sub-lethal damage repair (SLDR) competes with DNA damage induction in cultured cells. By virtue of SLDR, cell survival increases with decrease of dose-rate, so-called dose-rate effects (DREs). Here, we focused on a wide dose-rate range and investigated the change of cell-cycle distribution during X-ray protracted exposure and dose-response curves via hybrid analysis with a combination of in vitro experiments and mathematical modelling. In the course of flow-cytometric cell-cycle analysis and clonogenic assays, we found the following responses in CHO-K1 cells: (1) The fraction of cells in S phase gradually increases during 6 h exposure at 3.0 Gy/h, which leads to radio-resistance. (2) Slight cell accumulation in S and G₂/M phases is observed after exposure at 6.0 Gy/h for more than 10 hours. This suggests that an increase of SLDR rate for cells in S phase during irradiation may be a reproducible factor to describe changes in the dose-response curve at dose-rates of 3.0 and 6.0 Gy/h. By re-evaluating cell survival for various dose-rates of 0.186-60.0 Gy/h considering experimental-based DNA content and SLDR, it is suggested that the change of S phase fraction during irradiation modulates the dose-response curve and is possibly responsible for some inverse DREs.

Estimation of the radiation-induced DNA double-strand breaks number by considering cell cycle and absorbed dose per cell nucleus

DNA double-strand breaks (DSBs) are thought to be the main cause of cell death after irradiation. In this study, we estimated the probability distribution of the number of DSBs per cell nucleus by considering the DNA amount in a cell nucleus (which depends on the cell cycle) and the statistical variation in the energy imparted to the cell nucleus by X-ray irradiation. The probability estimation of DSB induction was made following these procedures: (i) making use of the Chinese Hamster Ovary (CHO)-K1 cell line as the target example, the amounts of DNA per nucleus in the logarithmic and the plateau phases of the growth curve were measured by flow cytometry with propidium iodide (PI) dyeing; (ii) the probability distribution of the DSB number per cell nucleus for each phase after irradiation with 1.0 Gy of 200 kVp X-rays was measured by means of γ-H2AX immunofluorescent staining; (iii) the distribution of the cell-specific energy deposition via secondary electrons produced by the incident X-rays was calculated by WLTrack (in-house Monte Carlo code); (iv) according to a mathematical model for estimating the DSB number per nucleus, we deduced the induction probability density of DSBs based on the measured DNA amount (depending on the cell cycle) and the calculated dose per nucleus. The model exhibited DSB induction probabilities in good agreement with the experimental results for the two phases, suggesting that the DNA amount (depending on the cell cycle) and the statistical variation in the local energy deposition are essential for estimating the DSB induction probability after X-ray exposure.

Integrated Modelling of Cell Responses after Irradiation for DNA-Targeted Effects and Non-Targeted Effects

Intercellular communication after ionizing radiation exposure, so-called non-targeted effects (NTEs), reduces cell survival. Here we describe an integrated cell-killing model considering NTEs and DNA damage along radiation particle tracks, known as DNA-targeted effects (TEs) based on repair kinetics of DNA damage. The proposed model was applied to a series of experimental data, i.e., signal concentration, DNA damage kinetics, cell survival curve and medium transfer bystander effects (MTBEs). To reproduce the experimental data, the model considers the following assumptions: (i) the linear-quadratic (LQ) function as absorbed dose to express the hit probability to emit cell-killing signals, (ii) the potentially repair of DNA lesions induced by NTEs, and (iii) lower efficiency of repair for the damage in NTEs than that in TEs. By comparing the model results with experimental data, we found that signal-induced DNA damage and lower repair efficiency in non-hit cells are responsible for NTE-related repair kinetics of DNA damage, cell survival curve with low-dose hyper-radiosensitivity (HRS) and MTBEs. From the standpoint of modelling, the integrated cell-killing model with the LQ relation and a different repair function for NTEs provide a reasonable signal-emission probability and a new estimation of low-dose HRS linked to DNA repair efficiency.

Markov chain Monte Carlo analysis for the selection of a cell-killing model under high-dose-rate irradiation

PURPOSE

High-dose-rate irradiation with 6 MV linac x rays is a wide-spread means to treat cancer tissue in radiotherapy. The treatment planning relies on a mathematical description of surviving fraction (SF), such as the linear-quadratic model (LQM) formula. However, even in the case of high-dose-rate treatment, the repair kinetics of DNA damage during dose-delivery time plays a function in predicting the dose-SF relation. This may call the SF model selection into question when considering the dose-delivery time or dose-rate effects (DREs) in radiotherapy and in vitro cell experiments. In this study, we demonstrate the importance of dose-delivery time at high-dose-rate irradiations used in radiotherapy by means of Bayesian estimation.

METHODS

To evaluate the model selection for SF, three types of models, the LQM and two microdosimetric-kinetic models with and without DREs (MKM_DR and MKM) were applied to describe in vitroSF data (our work and references). The parameters in each model were evaluated by a Markov chain Monte Carlo (MCMC) simulation.

RESULTS

The MCMC analysis shows that the cell survival curve by the MKM_DR fits the experimental data the best in terms of the deviance information criterion (DIC). In the fractionated regimen with 30 fractions to a total dose of 60 Gy, the final cell survival estimated by the MKM_DR was higher than that by the LQM. This suggests that additional fractions are required for attaining the total dose equivalent to yield the same effect as the conventional regimen using the LQM in fractionated radiotherapy.

CONCLUSIONS

Damage repair during dose-delivery time plays a key role in precisely estimating cell survival even at a high dose rate in radiotherapy. Consequently, it was suggested that the cell-killing model without repair factor during a short dose-delivery time may overestimate actual cell killing in fractionated radiotherapy.

Inorganic polyphosphate enhances radio-sensitivity in a human non-small cell lung cancer cell line, H1299

Inorganic polyphosphate is a linear polymer containing tens to hundreds of orthophosphate residues linked by high-energy phosphoanhydride bonds. Polyphosphate has been recognized as a potent anti-metastasis reagent. However, the molecular mechanism underlying polyphosphate action on cancer cells is poorly understood. In this study, we investigated the involvement of polyphosphate in radio-sensitivity using a human non-small cell lung cancer cell line, H1299. We found that polyphosphate treatment decreases cellular adenosine triphosphate levels, suggesting a disruption of energy metabolism. We also found that the induction of DNA double-strand breaks was enhanced in polyphosphate-treated cells after X-ray irradiation and colony formation assay revealed that cell survival decreased compared with that of the control groups. These findings suggest that polyphosphate is a promising radio-sensitizer for cancer cells. Therefore, we hypothesized that polyphosphate treatment disrupts adenosine triphosphate-mediated energy transfer for cellular survival and DNA repair, thereby reducing the cellular capability to resist X-ray irradiation.

Modeling cell survival and change in amount of DNA during protracted irradiation

Hyper-radiosensitivity (HRS) is a well-known bioresponse under low-dose or low-dose-rate exposures. Although disorder of the DNA repair function, non-targeted effects and accumulation of cells in G2 have been experimentally observed, the mechanism for inducing HRS by long-term irradiation is still unclear. On the basis of biological experiments and a theoretical study, we have shown that change in the amount of DNA associated with accumulation of cells in G2 enhances radiosensitivity. To demonstrate continuous irradiation with 250 kVp X-rays, we adopted a fractionated regimen of 0.186 or 1.00 Gy per fraction at intervals of 1 h (i.e. 0.186 Gy/h, 1.00 Gy/h on average) to Chinese Hamster Ovary (CHO)-K1 cells. The change in the amount of DNA during irradiation was quantified by flow cytometric analysis with propidium iodide (PI). Concurrently, we attempted a theoretical evaluation of the DNA damage by using a microdosimetric-kinetic (MK) model that was modified to incorporate the change in the amount of DNA. Our experimental results showed that the fraction of the cells in G2/M phase increased by 6.7% with 0.186 Gy/h and by 22.1% with 1.00 Gy/h after the 12th irradiation. The MK model considering the change in amount of DNA during the irradiation exhibited a higher radiosensitivity at a high dose range, which could account for the experimental clonogenic survival. The theoretical results suggest that HRS in the high dose range is associated with an increase in the total amount of DNA during irradiation.

Educational Activity for the Radiation Emergency System in the Northern Part of Japan: Meeting Report on "The 3rd Educational Symposium on Radiation and Health (ESRAH) by Young Scientists in 2016"

In the northern part of Japan, close cooperation is essential in preparing for any possible emergency response to radiation accidents because several facilities, such as the Low-Level Radioactive Waste Disposal Facility, the MOX Fuel Fabrication Plant and the Vitrified Waste Storage Center, exist in Rokkasho Village (Aomori Prefecture). After the accident at Fukushima Daiichi Nuclear Power Plant in 2011, special attention should be given to the relationship between radiation and human health, as well as establishing a system for managing with a radiation emergency. In the area of Hokkaido and Aomori prefectures in Japan, since 2008 an exchange meeting between Hokkaido University and Hirosaki University has been held every year to have opportunities to discuss radiation effects on human health and to collect the latest news on monitoring environmental radiation. This meeting was elevated to an international meeting in 2014 titled "Educational Symposium on Radiation and Health (ESRAH) by Young Scientists". The 3rd ESRAH meeting took place in 2016, with the theme "Investigating Radiation Impact on the Environmental and Health". Here we report the meeting findings on the continuing educational efforts after the Fukushima incident, what was accomplished in terms of building a community educational approaches, and future goals.

Evaluation of the cell survival curve under radiation exposure based on the kinetics of lesions in relation to dose-delivery time

We have investigated the dose rate effects on cell damage caused by photon-beam irradiation. During a relatively long dose-delivery time with a low dose rate, lesions created in cells may undergo some reactions, such as DNA repair. In order to investigate these reactions quantitatively, we adopted the microdosimetric-kinetic (MK) model and deduced a cell surviving fraction (SF) formula for continuous irradiation. This model enabled us to estimate the SF from dose and dose rate. The parameters in the MK model were determined so as to generate the SF, and we attempted to evaluate the dose rate effects on the SF. To deduce the cell-specific parameters in the SF formula, including the dose rate, we performed a split-dose experiment and a single-dose experiment with a constant dose-delivery time (10 min) (to retain the condition for equivalent behavior of cell lesions) by means of a clonogenic assay. Then, using the MK model parameters, the SFs were reproduced for a variety of dose rates (1.0, 0.31, 0.18, 0.025 and 0.0031 Gy/min) and were compared with reported experimental data. The SF curves predicted by the MK model agreed well with the experimental data, suggesting that the dose rate effects appear in the kinetics of cell lesions during the dose-delivery time. From fitting the analysis of the model formula to the experimental data, it was shown that the MK model could illustrate the characteristics of log-SF in a rectilinear form at a high dose range with a relatively low dose rate.

Quantitative estimation of DNA damage by photon irradiation based on the microdosimetric-kinetic model

The microdosimetric-kinetic (MK) model is one of the models that can describe the fraction of cells surviving after exposure to ionizing radiation. In the MK model, there are specific parameters, k and yD, where k is an inherent parameter to represent the number of potentially lethal lesions (PLLs) and yD indicates the dose-mean lineal energy in keV/μm. Assuming the PLLs to be DNA double-strand breaks (DSBs), the rate equations are derived for evaluating the DSB number in the cell nucleus. In this study, we estimated the ratio of DSBs for two types of photon irradiation (6 MV and 200 kVp X-rays) in Chinese hamster ovary (CHO-K1) cells and human non-small cell lung cancer (H1299) cells by observing the surviving fraction. The estimated ratio was then compared with the ratio of γ-H2AX foci using immunofluorescent staining. For making a comparison of the number of DSBs among a variety of radiation energy cases, we next utilized the survival data in the literature for both cells exposed to other photon types, such as (60)Co γ-rays, (137)Cs γ-rays and 100 kVp X-rays. The ratio of DSBs based on the MK model with conventional data was consistent with the ratio of γ-H2AX foci numbers, confirming that the γ-H2AX focus is indicative of DSBs. It was also shown that the larger yD is, the larger the DSB number is. These results suggest that k and yD represent the characteristics of the surviving fraction and the biological effects for photon irradiation.

Enhancement of macrosphelide-induced apoptosis by mild hyperthermia

Hyperthermia is a useful adjunct in cancer therapy as it can increase the effectiveness and decrease the toxicity of currently available cancer treatments such as chemotherapy and radiation. In the present study, we investigated whether 41 degrees C hyperthermia (mild HT) for 20 min can enhance macrosphelide (MS5)-induced apoptosis in human lymphoma U937 cells. Our results revealed that, compared with MS5 (5 microM) and mild HT alone, the combined treatment exhibited significant enhancement in apoptosis at 6 h, which was evaluated by observing morphological changes and DNA fragmentation. Marked increase in the reactive oxygen species (ROS) generation was observed immediately after the combined treatment. Significant increase in Fas externalization, caspase-8 and caspase-3 activation, and loss of mitochondrial membrane potential (MMP) was found after the combined treatment compared with MS5 and mild HT alone. Moreover, this combination can also alter the expression of apoptosis-related proteins as evident by the cleavage of Bid and down-regulation of Bcl-2 while no change in the expression of Bax was observed. Furthermore, an immediate rise in the intracellular calcium ion ([Ca(2+)]i) concentration was observed after the combined treatment, which continuously increased in a time-dependent manner. In addition, mild HT treatment alone also increases [Ca(2+)]i concentration without inducing apoptosis. Our data indicate that early increase in ROS generation is mainly responsible for the enhancement of apoptosis after the combined treatment.

MAS-NMR spectroscopy studies in the setting reaction of glass ionomer cements

OBJECTIVES

The main objective is the characterisation of the setting reaction in glass ionomer cements based on experimental ionomer glasses with different fluorine content and a commercial glass ionomer cement liquid by using 13C CP/MAS-NMR, 29Si, 27Al and 31P MAS-NMR spectroscopy in order to receive information specifically about the cross-linking process.

METHODS

Different fluorine containing glass compositions based on 4.5SiO2-3Al2O3-1.5P2O5-(5-z)CaO-zCaF(2) where z=0-3, were mixed with a commercially available polymer liquid to form glass ionomer cements. The cements were subjected to 27Al, 13C CP/MAS, 29Si, and 31P MAS-NMR analysis.

RESULTS

The 27Al spectra showed clearly the formation of six-fold coordinate Al(VI), that may crosslink the carboxyl groups in the poly-acid molecules. A shift towards to more positive values of the carboxyl peak in the 13C CP/MAS-NMR spectra showed clearly the proton dissociation of the carboxyl groups. A shift towards more negative values was observed in the 29Si MAS-NMR spectra, suggesting formation of hydrated silica gel and consequently formation of additional Si-O-Si bonds. 31P MAS-NMR spectra also reflected changes in the coordination state around a PO4(3-) tetrahedron. Increasing the fluorine content of the glasses resulted generally in increased reactivity during setting, due to promoting cross-linking and repolymerisation of the silicate phase, followed by clear changes in the MAS-NMR spectra.

CONCLUSIONS

The cross-linking process during the setting reaction of glass ionomer cements can be followed by MAS-NMR spectroscopy observing the conversion of Al(IV) to Al(VI). The acid base setting reaction is completed in 1 day and no further significant changes in the MAS-NMR spectra can be observed. Further study is required in order to understand the role of phosphorus.

Release profile of antimicrobial agents from alpha-tricalcium phosphate cement

A new method for treating carious dentine with alpha-tricalcium phosphate (alpha-TCP) dental cement containing antimicrobial agents has been recently introduced. However, the release behavior of antimicrobial agents from this cement has not yet been clarified. The aim of this study is therefore to examine the release profile of the antimicrobial agents from the alpha-TCP cement. Three kinds of antimicrobial agents (metronidazole, cefaclor and ciprofloxacin) were added to two commercially available alpha-TCP cements (new apatite liner type I and type II). The set cements were then immersed in water at 37 degrees C and the released antimicrobial agents and Ca ion were determined at regular intervals for three months. In addition, scanning electron microscopic observations were conducted before and after immersion for three months. The release profile of the cements containing antimicrobial agents varied depending on the types of antimicrobial agents. The incorporation of antimicrobial agents affected the setting reaction of the cements. The release behavior of the drugs also varied depending on the types of the cements. The differences in the release profile between type I and type II cements reflected the structures and compositions of their matrices.

Composite resin restoration and postoperative sensitivity: clinical follow-up in an undergraduate program

OBJECTIVES

To analyze the relationship between the cavity depth and liners with postoperative sensitivity of resin composite restorations.

METHODS

A clinical follow-up was conducted on 319 resin composite restorations made in the final year of an undergraduate program over a 3-year period. Along with the analyses of cavity type, cavity depth, type of pulpal protection and the materials used, the postoperative sensitivity was also examined on each restoration.

RESULTS

Thirty-nine percent of the restorations had no protective layer (Group 1). As the depth of the prepared cavities increased, the restorations received one of the three pulpal protection methods; a calcium hydroxide base (Group 2), glass ionomer cement (Group 3), or protection with a calcium hydroxide base in combination with glass ionomer cement (Group 4). The incidence of postoperative sensitivity showed no significant difference among Groups 1, 2 and 3, but was significantly lower in Group 1 than in Group 4. The restorations made in shallow and medium depth cavities demonstrated significantly less-postoperative sensitivity than those made in deep cavities. The newer generation dentine-bonding agents showed a significantly lower incidence of postoperative sensitivity than the early generation group.

CONCLUSIONS

Postoperative sensitivity in resin composite restorations was not related to the absence of protective layers but increased with the depth of cavities restored with the resin composite. The type of dentine-bonding agents could also be responsible for postoperative sensitivity.

Radical scavenging by N-aminoazaaromatics

N-Aminoazaaromatics were found to react with nitric oxide in the presence of oxygen to afford deaminated products in high yields. The reaction proceeded almost instantaneously in various solvents including water, and one to two equivalent of NO was consumed depending upon the amount of oxygen coexisted, and 1 equivalent of N2O was released in the reaction. In addition, N-aminoazoles were deaminated by potassium superoxide to give parent azoles in good yields. Two equivalents of superoxide was consumed, and about half equivalents of both nitrite and nitrate ion were released. The results demonstrated that N-aminoazoles have ability to protect the biological system against the oxidation promoted by radicals such as nitrogen oxides and superoxide.

Total syntheses of novel cytocidal beta-carboline alkaloids, oxopropalines D and G

A new type of beta-carboline nucleus, N-methoxymethyl-4-methyl-beta-carboline (4) was synthesized by thermal electrocyclic reaction of a 1-azahexatriene system, involving the indole 2,3-bond. The key compound N-methoxymethyl-1-methoxycarbonyl-4-methyl-beta-carboline (2) was then prepared in a four-step sequence. The total synthesis of oxopropaline G (1e) was achieved from this key compound in four steps. Furthermore, the enantioselective total syntheses of (+)-oxopropaline D (1c) and its enantiomer were also achieved by application of the Sharpless oxidation-procedure in nine steps from 2.

Formulation of glutaraldehyde disinfectant for alginate impressions

The effect of buffer agents incorporated in glutaraldehyde disinfectants on the surface quality of dental models was examined by the measurement of surface roughness, X-ray diffraction analysis and SEM observation. Seven experimental glutaraldehyde disinfectants were prepared using two buffer agents, potassium acetate alone or potassium acetate and sodium hydrogen carbonate in combination. Four kinds of sulfate--zinc, calcium, potassium and magnesium sulfate--were added to these disinfectants in order to accelerate the hydration of calcium sulfate hemihydrate. The impressions treated with the experimental disinfectants for 1 h produced stone surfaces which had significantly lower surface roughness values than those treated with the commercial disinfectants (p < 0.05). The X-ray diffraction analysis and SEM observation showed that these superior surfaces were produced as a result of significant reductions in the amount of residual calcium sulfate hemihydrate. Replacement of buffer agents in commercial glutaraldehyde disinfectants with chemicals such as those studied in the present study will improve the surface quality of dental stone.

The relationship between water absorption characteristics and the mechanical strength of resin-modified glass-ionomer cements in long-term water storage

The purpose of this study is to elucidate the water absorption characteristics of resin-modified glass-ionomer cements and to also investigate the relationship between the characteristics and mechanical strength after long-term water storage. The mechanism of water diffusion in these cements is also discussed. Water absorption was measured using a gravimetric analysis for 12 m, while the diffusion coefficient was calculated using Fick's law of diffusion. Water solubility was determined based on the weight of the residue in the immersed water. The compressive and diametral tensile strength were measured at 1, 2, 6, and 12 m. A correlation was observed between the diffusion coefficient and equilibrium water uptake, which thus suggests the water in the cements to diffuse through micro-voids in accordance with the 'Free volumetric theory'. A correlation was seen between the solubility and diffusion coefficient of the cements. The deterioration ratio, defined as the ratio of the strength at 12 m versus that at 1 m, was also calculated. Finally, a negative correlation was observed between the deterioration ratio of the compressive strength and the diffusion coefficients of the cements.

Structure of bioactive glass and its application to glass ionomer cement

We prepared a new glass ionomer cement using bioactive CaO-P2O5-SiO2(-MgO) glass and investigated its setting process using FT-IR and MAS NMR analyses. The compressive strengths of the cements depended on the glass composition and a maximum strength of 33.3 +/- 4.7 MPa was obtained using cement with the glass composition of MgO:4.6, CaO:44.9, SiO2:34.2 and P2O5:16.3% in weight. FT-IR analysis showed that the COOH group in the polyacrylic acid decreased and carboxylate ion (COO-Ca2+) increased after the setting reaction. A broad signal appeared around -82 ppm in 29Si MAS-NMR spectra of the glass and a new signal corresponding to hydrated silica gel formation appeared around -102 and -111 ppm after setting. This suggests that Ca2+ was released from the glass powder to form carboxylate salt and that a degree of polymerization in the silicate network increased. The setting mechanism of the cement was found to be essentially the same as in conventional glass ionomer cement.

Effect of powder grinding on hydroxyapatite formation in a polymeric calcium phosphate cement prepared from tetracalcium phosphate and poly(methyl vinyl ether-maleic acid)

The primary aim of this study was to determine if cements based on poly(methyl vinyl ether-maleic acid) (PMVE-Ma) and tetracalcium phosphate resulted in hydroxyapatite formation. In addition, the mechanical strength of this type of polymeric calcium phosphate cement was evaluated. Cements were prepared by mixing, in a powder/liquid mass ratio of 3.0, an aqueous solution of PMVE-Ma (mass fraction = 25%) and tetracalcium phosphate powders ground for various periods of time. The tetracalcium phosphate powders and set cements were characterized by means of X-ray powder diffraction and scanning electron microscopy. Mechanical strengths of the cements were tested 24 h after mixing. Prolonged grinding of tetracalcium phosphate powder decreased particle size and/or crystallite size and increased lattice distortion. This enhanced the reactivity of the tetracalcium phosphate powder and elevated the extent of crosslinking between PMVE-Ma molecules, resulting in improved mechanical strength. Hydroxyapatite formation was detected in the cement prepared with the most finely ground tetracalcium phosphate powder. The conversion of residual tetracalcium phosphate particles to more thermodynamically stable hydroxyapatite crystals will reduce the solubility of the polymeric cement and increase its biocompatibility.

Effect of fluoride on apatite formation from Ca4(PO4)2O in 0.1 mol L(-1) KH2PO4

The effect of fluoride on the hydrolysis of tetracalcium phosphate (TTCP; Ca4(PO4)2O) was investigated in 0.1 mol l(-1)KH2PO4 containing 0-83 mol l(-1) KF. Characterization of the final apatite phase formed by the hydrolysis was made with X-ray diffraction and SEM. The initial pH was between 4.5 and 5.4, depending on the solutions, and the pH rapidly increased and was kept constant between 7.3 and 6.5. An increase in KF concentration tended to lower the pH in the final stage of hydrolysis. The calcium concentration was considerably lower than the phosphorus concentration throughout the reaction. The fluoride concentration decreased shortly after the start of hydrolysis. The hydrolysis of TTCP in 0.1 mol l(-1) KH2PO4 proceeded to form hydroxyapatite via DCPD when the KF concentration was low. The hydrolysis product was a calcium-deficient non-stoichiometric hydroxyapatite with a Ca/P ratio of about 1.5. With an increase in the KF concentration in the 0.1 mol l(-1) KH2PO4 solution, TTCP directly transformed into hydroxyapatite containing F- ions or fluorapatite and with improved crystallinity. The addition of fluoride in the solution initially accelerated the formation of apatite. However, the layer of newly formed apatite adhering to the TTCP particles retarded TTCP dissolution; as a result, hydrolysis was delayed. IR analysis showed that the apatite phase contained HPO2-4 ions in the structure. The formula for the hydrolysis product of TTCP in the presence of fluoride can be expressed as follows: Ca10-x(HPO4)x(PO4)10-x(OH)2-x-yFy.

Polymeric calcium phosphate cements derived from poly(methyl vinyl ether-maleic acid)

OBJECTIVES

The purpose of this study was to assess the feasibility of forming polymeric calcium phosphate cements from a mixed powder of dicalcium phosphate/tetracalcium phosphate or only tetracalcium phosphate and poly(methyl vinyl ether-maleic acid) (PMVE-Ma), and to study their setting reaction.

METHODS

The setting reaction process of the polymeric cements was evaluated by mechanical strength tests, infrared spectroscopy and x-ray diffraction analysis and compared with that of a water-setting calcium phosphate cement. The mechanical strength data were analyzed using ANOVA and Scheffé's multiple comparisons test.

RESULTS

Cements prepared from the mixed powder and 25-30 wt% aqueous solutions of PMVE-Ma had high mechanical strength after 24 h storage in distilled water at 37 degrees C. The hardening mechanism depended on an acid-base reaction between the carboxyl groups of PMVE-Ma and the mixed powder, especially its tetracalcium phosphate component. The formation of hydroxyapatite in the polymeric calcium phosphate cement was not detected and is apparently inhibited as a result of the competing reaction of PMVE-Ma with the mixed powder.

SIGNIFICANCE

The cement-forming reaction was significantly faster than that of a water-setting calcium phosphate cement and slower than that observed with the mixed powder and polyacids such as poly(acrylic acid). The characteristics of the polymeric cements suggest that the materials may be useful in cavity lining or endodontic sealing.

Effect of abietic acid and poly(methyl methacrylate) on the dissolution process of zinc oxide-eugenol cement

The effects of abietic acid and poly(methyl methacrylate) (PMMA) on the dissolution process of zinc oxide-eugenol cement were studied in 0.01 M acetic acid buffer solution (pH 4). Abietic acid markedly reduced the elution of zinc from the cement, because of the matrix stabilization by the reaction between the acid and zinc oxide. The dissolution process of the cement without additives was controlled by surface reaction. With the cement containing abietic acid or PMMA, the diffusion of eluted species through the cement matrix participated in the dissolution process.

[Studies on as-triazine derivatives. XIX. Synthesis of 2,3-diarylpyrazine 2,3-diarylpyridine derivatives as blood platelet aggregation inhibitors]

4,5-Diphenyl-2-ethoxypyrimidine (1), 3,4-diphenyl-6-ethoxypyridazine (2) and 2,3-diphenyl-5-ethoxypyrazine (3) were evaluated for inhibitory activity towards arachidonic acid-induced aggregation of rabbit blood platelet in vitro. 2,3-Diphenyl-5-ethoxypyrazine (3) exhibited significant inhibitory activity. Thus, various 5-substituted 2,3-bis(4-methoxyphenyl)pyrazines were synthesized by the nucleophilic substitution reaction of 5-chloro-2,3-bis(4-methoxyphenyl)pyrazine (9). In a similar manner, substituted 2,3-bis(4-methoxyphenyl)pyridines were prepared from 2,3-bis(4-methoxyphenyl)-6-methylsulfonylpyridine (17), which was synthesized by the cycloaddition retro Diels-Alder reaction of 5,6-bis(4-methoxyphenyl)-3-methylsulfonyl-1,2,4-triazine (16) with norbornadiene. Among the compounds prepared, 6-isopropoxy-2,3-bis(4-methoxyphenyl)-pyrazine (10f) showed the most potent inhibitory activity, which was more than the activity of anitrazafen[5,6-bis(4-methoxyphenyl)-3-methyl-1,2,4-triazine.