WO2022014761A1 - 흡수에너지 기반 전기장 암치료 계획 시스템 및 방법 - Google Patents
흡수에너지 기반 전기장 암치료 계획 시스템 및 방법 Download PDFInfo
- Publication number
- WO2022014761A1 WO2022014761A1 PCT/KR2020/009527 KR2020009527W WO2022014761A1 WO 2022014761 A1 WO2022014761 A1 WO 2022014761A1 KR 2020009527 W KR2020009527 W KR 2020009527W WO 2022014761 A1 WO2022014761 A1 WO 2022014761A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- electric field
- dose
- prescription
- tumor
- cancer treatment
- Prior art date
Links
- 230000005684 electric field Effects 0.000 title claims abstract description 166
- 238000011282 treatment Methods 0.000 title claims abstract description 129
- 206010028980 Neoplasm Diseases 0.000 title claims abstract description 116
- 238000000034 method Methods 0.000 title claims abstract description 50
- 201000011510 cancer Diseases 0.000 title claims abstract description 45
- 210000000056 organ Anatomy 0.000 claims abstract description 27
- 230000000704 physical effect Effects 0.000 claims abstract description 7
- 238000010521 absorption reaction Methods 0.000 claims description 17
- 230000004071 biological effect Effects 0.000 claims description 4
- 238000002591 computed tomography Methods 0.000 claims description 4
- 239000000463 material Substances 0.000 claims description 4
- 238000002595 magnetic resonance imaging Methods 0.000 claims description 2
- 238000005457 optimization Methods 0.000 abstract description 15
- 230000004663 cell proliferation Effects 0.000 description 38
- 230000002401 inhibitory effect Effects 0.000 description 36
- 230000000694 effects Effects 0.000 description 16
- 238000010586 diagram Methods 0.000 description 13
- 230000008569 process Effects 0.000 description 12
- 208000005017 glioblastoma Diseases 0.000 description 9
- 238000002560 therapeutic procedure Methods 0.000 description 8
- 230000004083 survival effect Effects 0.000 description 5
- 230000001225 therapeutic effect Effects 0.000 description 5
- 230000006870 function Effects 0.000 description 4
- 238000011160 research Methods 0.000 description 4
- 238000002512 chemotherapy Methods 0.000 description 3
- 238000004070 electrodeposition Methods 0.000 description 3
- 208000003174 Brain Neoplasms Diseases 0.000 description 2
- 206010027476 Metastases Diseases 0.000 description 2
- 206010061902 Pancreatic neoplasm Diseases 0.000 description 2
- BPEGJWRSRHCHSN-UHFFFAOYSA-N Temozolomide Chemical compound O=C1N(C)N=NC2=C(C(N)=O)N=CN21 BPEGJWRSRHCHSN-UHFFFAOYSA-N 0.000 description 2
- 230000006907 apoptotic process Effects 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 238000012423 maintenance Methods 0.000 description 2
- 208000015486 malignant pancreatic neoplasm Diseases 0.000 description 2
- 230000009401 metastasis Effects 0.000 description 2
- 210000000920 organ at risk Anatomy 0.000 description 2
- 201000002528 pancreatic cancer Diseases 0.000 description 2
- 208000008443 pancreatic carcinoma Diseases 0.000 description 2
- 238000004393 prognosis Methods 0.000 description 2
- 230000000306 recurrent effect Effects 0.000 description 2
- 229960004964 temozolomide Drugs 0.000 description 2
- 201000010915 Glioblastoma multiforme Diseases 0.000 description 1
- 206010018338 Glioma Diseases 0.000 description 1
- 241001465754 Metazoa Species 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000011394 anticancer treatment Methods 0.000 description 1
- 230000004611 cancer cell death Effects 0.000 description 1
- 230000034994 death Effects 0.000 description 1
- 238000013135 deep learning Methods 0.000 description 1
- 230000003111 delayed effect Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000004980 dosimetry Methods 0.000 description 1
- 230000005672 electromagnetic field Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 208000029824 high grade glioma Diseases 0.000 description 1
- 230000002700 inhibitory effect on cancer Effects 0.000 description 1
- 210000004072 lung Anatomy 0.000 description 1
- 201000011614 malignant glioma Diseases 0.000 description 1
- 230000001394 metastastic effect Effects 0.000 description 1
- 206010061289 metastatic neoplasm Diseases 0.000 description 1
- 231100000957 no side effect Toxicity 0.000 description 1
- 238000011275 oncology therapy Methods 0.000 description 1
- 230000001151 other effect Effects 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 238000001959 radiotherapy Methods 0.000 description 1
- 230000010076 replication Effects 0.000 description 1
- 238000012552 review Methods 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000003325 tomography Methods 0.000 description 1
- 230000035899 viability Effects 0.000 description 1
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61N—ELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
- A61N1/00—Electrotherapy; Circuits therefor
- A61N1/18—Applying electric currents by contact electrodes
- A61N1/32—Applying electric currents by contact electrodes alternating or intermittent currents
- A61N1/36—Applying electric currents by contact electrodes alternating or intermittent currents for stimulation
- A61N1/36002—Cancer treatment, e.g. tumour
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/0033—Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room
- A61B5/0036—Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room including treatment, e.g., using an implantable medical device, ablating, ventilating
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/05—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves
- A61B5/053—Measuring electrical impedance or conductance of a portion of the body
- A61B5/0537—Measuring body composition by impedance, e.g. tissue hydration or fat content
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/05—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves
- A61B5/055—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves involving electronic [EMR] or nuclear [NMR] magnetic resonance, e.g. magnetic resonance imaging
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/48—Other medical applications
- A61B5/4836—Diagnosis combined with treatment in closed-loop systems or methods
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7235—Details of waveform analysis
- A61B5/7264—Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/03—Computed tomography [CT]
- A61B6/032—Transmission computed tomography [CT]
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
- A61B6/5211—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
- A61B6/5217—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data extracting a diagnostic or physiological parameter from medical diagnostic data
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61N—ELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
- A61N1/00—Electrotherapy; Circuits therefor
- A61N1/02—Details
- A61N1/04—Electrodes
- A61N1/0404—Electrodes for external use
- A61N1/0472—Structure-related aspects
- A61N1/0476—Array electrodes (including any electrode arrangement with more than one electrode for at least one of the polarities)
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61N—ELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
- A61N1/00—Electrotherapy; Circuits therefor
- A61N1/02—Details
- A61N1/08—Arrangements or circuits for monitoring, protecting, controlling or indicating
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61N—ELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
- A61N1/00—Electrotherapy; Circuits therefor
- A61N1/18—Applying electric currents by contact electrodes
- A61N1/32—Applying electric currents by contact electrodes alternating or intermittent currents
- A61N1/36—Applying electric currents by contact electrodes alternating or intermittent currents for stimulation
- A61N1/36014—External stimulators, e.g. with patch electrodes
- A61N1/3603—Control systems
- A61N1/36034—Control systems specified by the stimulation parameters
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61N—ELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
- A61N1/00—Electrotherapy; Circuits therefor
- A61N1/40—Applying electric fields by inductive or capacitive coupling ; Applying radio-frequency signals
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/11—Region-based segmentation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10072—Tomographic images
- G06T2207/10081—Computed x-ray tomography [CT]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10072—Tomographic images
- G06T2207/10088—Magnetic resonance imaging [MRI]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20081—Training; Learning
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30096—Tumor; Lesion
Definitions
- the present application relates to an absorption energy-based electric field cancer treatment planning system and method, and more specifically, to each tissue based on the strength of the electric field given to each tissue in the human body, the tissue density, the tissue electrical conductivity, and the electric field application time. Calculate the given physical absorbed energy, determine the patient's prescription dose based on this, and establish a treatment plan so that the determined prescription dose is best applied to the tumor and minimized to the surrounding normal tissues. It relates to an absorbed energy-based electric field cancer treatment planning system and method for performing treatment based on the present invention.
- the first strength of electric field therapy is that it is known that the electric field has a great effect only on dividing cells, so it will intensively affect cancer cells that divide faster than normal cells. The point is that it becomes (refer to Non-Patent Document 2).
- Non-Patent Document 3 the side effects of electric field treatment were significantly lower in 7 items out of a total of 9 comparison items of side effects between chemotherapy and electric field treatment, and in 2 items, anticancer Treatment and electric field treatment showed almost equal side effects.
- the third strength is that when the electric field is applied extensively, including the treatment area, it is expected that it will have a therapeutic effect even on very fine tumors that are not visible on medical images such as CT.
- an electric field is applied around the tumor, an electric field that cannot be ignored is transmitted not only to the tumor but also to the tumor area. Therefore, it is expected that the cancer metastasis probability can also be dramatically reduced (see Non-Patent Document 5).
- Non-Patent Document 6 discloses a method for treating cancer diseases.
- electric field cancer treatment has been approved by the US FDA for recurrent glioblastoma cancer in 2011 and for the first diagnosed glioblastoma in 2015. In Japan, it has been approved for treatment of patients with recurrent glioblastoma.
- the number of treated patients is increasing rapidly every year, showing a 50-fold increase from 152 in 2014 to 8813 in 2018 (see Non-Patent Document 6).
- the effect of electric field treatment on cancer cells differs depending on the strength of the electric field applied to the tissue and the application time of the electric field. Specifically, the greater the intensity and application time of the electric field, the greater the apoptosis and division inhibitory effect on cancer cells, so it can be said that the strength and application time of the electric field are in a proportional relationship with the apoptosis effect of the electric field (see Non-Patent Documents 7 and 8). ).
- an appropriate electric field strength or treatment time is not determined to treat a specific tumor, and treatment is being performed with the maximum electric field strength and maximum treatment time within a range where there are no side effects.
- electric field therapy does not have a concept of a dose generally used in radiation therapy or chemotherapy, and it can be said that a prescription dose that is applied differently depending on the type or condition of the tumor is not applied.
- the power loss density which is the unit of power loss per unit volume in electric field treatment
- the unit of power loss density can be regarded as the standard of almost the only prescription suggested so far for electric field therapy, there is a fatal limit to using this unit as a dose unit for electric field therapy.
- the degree of cancer cell death is proportional to the electric field application time
- the power loss density does not include items related to the total treatment time that the patient needs to be treated (i.e., the electric field application time). and treatment cannot be optimized (see Non-Patent Document 8).
- the power loss density may be suitable for use as a factor predicting the prognosis of electric field treatment, but there is a limit to using it as a dose unit for electric field treatment. There is a need for a unit of dose.
- Non-Patent Document 1 Eilon D. Kirson et al, disruption of cancer cell replication by alternating electric fields, cancer research, 64, 3288-3295 (2004)
- Non-Patent Document 2 Miklos Pless, Uri Weinberg, tumor treating fields: concept, evidence, future, Expert Opinion, 20(8), 1099-1106 (2011)
- Non-Patent Document 3 Angela M. Davies et al, Tumor treating fields: a new frontier in cancer therapy, Annals of the New York academy of sciences, 1291, 86-95 (2013)
- Non-Patent Document 4 Stupp et al, Effect of Tumor-Treating Fields Plus Maintenance Temozolomide vs Maintenance Temozolomide Alone on Survival in Patients With Glioblastoma: A Randomized Clinical Trial, Journal of the American Medical Association, 318(23), 2306-2316 (2017)
- Non-Patent Document 5 Eilon D. Kirson et al. Alternating electric fields (TTFields) inhibit metastatic spread of solid tumors to the lungs, Clin Exp Metastasis 26, 633-640 (2009)
- Non-Patent Document 6 Novocure Corporate Presentation (https://3sj0u94bgxp33grbz1fkt62h-wpengine.netdna-ssl.com/wp-content/uploads/2019/2017201905_NVCR_Corporate_Presentation_vFF.pdf)
- Non-Patent Document 7 Eilon D. Kirson et al, alternating electric fields arrest cell proliferation in animal tumor models and human brain tumors, PNAS, 104(24), 10152-10157 (2007)
- Non-Patent Document 8 Yunhui Jo et al, Effectiveness of a Fractionated Therapy Scheme in Tumor Treating Fields Therapy, Technology in Cancer Research & Treatment, 18, 1-10 (2019)
- Non-Patent Document 9 Denise Fabian et al, Treatment of Glioblastoma (GBM) with the Addition of Tumor-Treating Fields (TTF): A Review, Cancers, 11, 174 (2019)
- Non-Patent Document 10 Matthew T. Ballo et al, Correlation of Tumor Treating Fields Dosimetry to Survival Outcomes in Newly Diagnosed Glioblastoma: A Large-Scale Numerical Simulation-Based Analysis of Data from the Phase 3 EF-14 Randomized Trial, Int J Radiation Oncol Biol Phys, Vol. 104, No. 5, pp. 1106-1113 (2019)
- Non-Patent Document 11 Dimitris J. Panagopoulos et al, Evaluation of Specific Absorption Rate as a Dosimetric Quantity for Electromagnetic Fields Bioeffects, PLoS One, 8(6), e62663 (2013)
- Non-Patent Document 12 Stefano Mandija et al. Opening a new window on MR-based Electrical Properties Tomography with deep learning, Scientific Reports, 9, 8895 (2019)
- a treatment plan is established so that the prescribed dose can be applied to the tumor as best as possible and given to the surrounding normal tissue as a minimum, and the parameters derived from this are applied to the electric field treatment device. It is required to set the standard of dose for electric field treatment so that the treatment effect can be maximized and side effects can be minimized by allowing the treatment to be performed in conjunction with the
- an embodiment of the present invention provides an absorption energy-based electric field cancer treatment planning system.
- the absorption energy-based electric field cancer treatment planning system includes: an image classification unit for classifying organs and tumors in a medical image of a patient including organs and tumors; a material property information setting unit configured to set material property information of each area classified by the image classification unit; a prescription information determining unit for setting a prescription method in consideration of the input tumor type and tumor status information, and determining prescription information including prescription dose, number of treatments, treatment time, and treatment frequency; In consideration of the position of each region classified by the image classification unit, the number of electrodes used for electric field cancer treatment, the position of the electrode, and the electric field strength for each electrode are initially set, and the dose to calculate the dose distribution in the body based on the initial setting output unit; And by changing at least one of the number of electrodes, the position of the electrodes, and the electric field strength for each electrode, it may include a dose optimizing unit optimizing the dose so that the prescribed dose is given to the tumor as much as possible and to the surrounding normal tissue at a minimum.
- another embodiment of the present invention provides a method for planning an electric field cancer treatment based on absorbed energy.
- the absorbed energy-based electric field cancer treatment planning method includes: segmenting a region of interest including a tumor and an organ in a patient's medical image; Setting a prescription method in any one of the absorbed energy density-based or absorbed specific energy-based; determining prescription information including a prescription dose, a number of treatments, a treatment time, and a treatment frequency according to a set prescription method; calculating a dose distribution in the body based on the number of electrodes initially set in consideration of the tumor and organs, the positions of the electrodes, and the electric field strength for each electrode; and evaluating and optimizing the treatment plan based on the calculated body dose distribution.
- the dose of electric field treatment can be set based on absorbed energy that is substantially proportional to the cell proliferation inhibitory effect, and based on this, It becomes possible to determine the applicable clinical prescription dose. Through this, it becomes possible to establish and evaluate a systematic treatment plan, which was not possible due to the absence of a dose unit in electric field therapy until now.
- the prescribed dose is delivered to the tumor as much as possible and at the same time to the surrounding normal tissues to the minimum. treatment can be performed.
- 1 is a view showing the cell proliferation inhibitory effect according to the electric field strength when an electric field is applied.
- FIG. 2 is a view showing the cell proliferation inhibitory effect according to the power loss density and the electric field application time when the electric field is applied.
- 3 is a view showing the cell proliferation inhibitory effect when different power loss densities and different electric field application times are applied when an electric field is applied.
- FIG. 4 is a diagram illustrating a process of calculating the absorbed energy density or absorbed specific energy distribution in the body using a medical image of a patient during electric field treatment.
- FIG. 6 is a diagram illustrating a process of optimizing a dose in the body using an absorbed energy-based dose unit when an electric field treatment plan is established.
- FIG. 7 is a diagram illustrating the results of dose optimization performed using absorbed energy-based dose units during electric field treatment assuming an organ at risk composed of virtual tumor tissue and normal cells in the human body model phantom.
- FIG. 8 is a diagram showing the dependence of the electric field frequency of the cell proliferation inhibitory effect upon application of the electric field.
- FIG. 9 is a flowchart of an electric field cancer treatment planning method according to an embodiment of the present invention.
- FIG. 10 is a block diagram of an electric field cancer treatment planning system according to another embodiment of the present invention.
- FIGS. 1 to 3 after applying an electric field with an intensity of ⁇ 1 V/cm and a frequency of 150 kHz to the malignant glioma cell line (U373), the viability of the cells is measured through cell counting. Compared to the control group, the inhibitory effect of cell proliferation of the cell line (experimental group) to which the electric field was applied was expressed as a percentage compared to the control group. This result is an average value after performing the experiment three times, and the percentage of the cell proliferation inhibitory effect compared to the control group was calculated by the following Equation 1.
- the cell proliferation inhibitory effect is 30% compared to the control group in the experimental group.
- FIG. 1 is a view showing the cell proliferation inhibitory effect according to the electric field strength when an electric field is applied, (a) shows the configuration of the experimental equipment for verifying the cell proliferation inhibitory effect by the 150 kHz alternating current electric field, (b) is the electric field strength ( Electric field intensity) is a graph showing the cell proliferation inhibitory effect.
- FIG. 2 is a view showing the power loss density when an electric field is applied and the cell proliferation inhibitory effect according to the electric field application time
- (a) is a diagram showing the relationship between the power loss density and the cell proliferation inhibitory effect by the 150 kHz alternating current electric field A graph
- (b) is a graph showing the relationship between the electric field application time (Treatment duration) and the cell proliferation inhibitory effect when the power loss density is fixed to 11.5 mW/cm 3 .
- the cell proliferation inhibitory effect upon application of an electric field is proportional to the power loss density (Power density).
- the power loss density is defined as 1/2 ⁇ E 2 ( ⁇ : electrical conductivity, E: electric field strength), and the curve shown in (a) is the result of fitting the data to the first-order curve.
- the cell proliferation inhibitory effect upon application of the electric field is proportional to the time of application of the electric field.
- FIG. 3 is a view showing the cell proliferation inhibitory effect when different power loss densities and different electric field application times are applied when an electric field is applied, (a) is a power loss density of 6.5 mW/cm 3 and 11.5 mW/cm 3 , respectively In the case of , it is a graph showing the cell proliferation inhibitory effect according to the electric field application time, and (b) is a graph showing the cell proliferation inhibitory effect according to the absorbed energy per unit volume or per unit mass.
- Equation 3 (b) shows the relationship between the cell proliferation inhibitory effect and absorbed energy density (ie, absorbed energy per unit volume, absorbed energy density) or absorbed specific energy (ie, absorbed energy per unit mass, specific energy absorbed).
- the absorbed energy density and the absorption specific energy can be calculated according to the following Equations 2 and 3, respectively.
- Absorbed energy density power loss density ⁇ electric field application time
- Absorption specific energy (power loss density ⁇ electric field application time) ⁇ (mass density)
- the cell proliferation inhibitory effect of the cell line to which the electric field is applied is proportional to the absorbed energy density or the absorbed specific energy.
- the absorbed energy delivered to the cell that is, absorbed energy density or absorbed specific energy, is a concept that includes the electric field application time (s) proportional to the cell proliferation inhibitory effect along with the power loss density. Compared to that, it can be said to be a more practical and reasonable standard for quantifying the cell proliferation inhibitory effect.
- the absorbed energy density and the absorbed specific energy are classified depending on whether the mass per unit volume is considered, and both can be used as a criterion for quantifying the cell proliferation inhibitory effect when an electric field is applied.
- FIG. 4 is a diagram illustrating a process of calculating the absorbed energy density or absorbed specific energy distribution in the body using a medical image of a patient during electric field treatment.
- treatment planning conditions eg, electrode position, electric field strength for each electrode, treatment time, etc.
- the three-dimensional absorbed energy density distribution in the human body may be calculated according to Equation 2 above (S44).
- (c) denotes the absorption of energy density (J / cm 3) in the case where the mass density of the tumor 1.2 g / cm 3 distributed
- (d) is absorbed in the case of the mass density of the tumor 1.2 g / cm 3 shows the distribution of specific energy (J/g)
- (e) shows the distribution of absorbed energy density (J/cm 3 ) when the mass density of the tumor is 0.6 g/cm 3
- (f) shows the distribution of the mass density of the tumor The absorption specific energy (J/g) distribution in the case of 0.6 g/cm 3 is shown.
- FIG. 6 is a diagram illustrating a process of optimizing a dose in the body using an absorbed energy-based dose unit when an electric field treatment plan is established.
- the prescribed dose may be any one of absorbed energy density and absorbed specific energy.
- the electrode position, the electric field strength for each electrode, and the treatment time are initially set (S63), based on this, the body electric field dose is calculated (S64), and the calculated body electric field dose distribution is evaluated (S65), and a predetermined standard It can be determined whether or not is satisfied (S66).
- the treatment can be performed so that the prescribed dose (ie, absorbed energy density or absorbed specific energy) is best given to the tumor and minimized to surrounding normal tissues (S67).
- the prescribed dose ie, absorbed energy density or absorbed specific energy
- FIG. 7 is a diagram showing the results of dose optimization performed using an absorbed energy-based dose unit during electric field treatment assuming an important organ (organ at risk) consisting of a virtual tumor tissue and normal cells in the human body model phantom, (a) shows the phantom containing the tumor tissue and important organs, (b) shows the dose distribution before dose optimization, (c) shows the dose distribution after dose optimization, (d) shows the An electric field cancer treatment system associated with an electric field cancer treatment planning system according to an embodiment is shown.
- organ organ at risk
- a virtual U-shaped tumor and a spherical important organ which is a normal tissue, were set in the human body model phantom, and then a 200 kHz AC electric field was applied.
- FIG. 8 is a diagram showing the dependence of the cell proliferation inhibitory effect on the electric field frequency when the electric field is applied.
- CFPAC-I pancreatic cancer cell lines
- HPAF-II pancreatic cancer cell lines
- the cell proliferation inhibitory effect which varies depending on the frequency of application, is expressed as a percentage compared to the 150 kHz frequency. That is, in FIG. 8 , since the strength of the electric field and the application time of the electric field are the same and only the frequency is different, the absorbed energy transferred to the cells is all the same, and the biological effect according to the frequency is shown.
- cell proliferation is most effectively inhibited at 150 kHz when the same amount of absorbed energy is delivered to the cells of the two cell lines used in the experiment.
- the cell proliferation inhibitory effect falls to about 60 to 70% level compared to the case of treatment with a frequency of 150 kHz.
- This frequency-dependent phenomenon is a general phenomenon that occurs not only in pancreatic cancer cells but also in all types of cells (see Non-Patent Document 7).
- the frequency (ie, reference frequency) that is the standard for the prescribed dose it is preferable to set the frequency (ie, reference frequency) that is the standard for the prescribed dose together, and it is preferable to use a frequency other than the reference frequency.
- a frequency weighting factor it is necessary to revise the prescribed dose by applying a frequency weighting factor accordingly.
- the prescribed dose at the reference frequency is 100J/g and the frequency weight of the used frequency is 2.0
- the prescribed dose should be set by modifying it to 200J/g.
- the frequency weight for the used frequency may be defined according to Equation 4 below, and the prescribed dose at the used frequency may be modified according to Equation 5 below.
- Frequency weight (frequency used) cell proliferation inhibitory effect (reference frequency) / cell proliferation inhibitory effect (used frequency)
- Prescribed dose (used frequency) Prescribed dose (reference frequency) ⁇ frequency weight
- FIG. 9 is a flowchart of an electric field cancer treatment planning method according to an embodiment of the present invention.
- a medical image (eg, 3D medical image data) of a patient is imported into the treatment planning system ( S91 ), and regions of interest such as a target tumor and important organs can be divided ( S92 ).
- the prescription method is set in either the absorbed energy density-based or the absorbed specific energy-based method (S93, S94), and prescription information including the prescribed dose, number of treatments, treatment time, treatment frequency, etc. according to the set prescription method can be determined (S95).
- the number of electrodes, the position of the electrodes, and the electric field strength for each electrode are initially set in consideration of the location of the tumor and important organs around the tumor, and then the dose distribution in the body can be calculated based on this (S96).
- the treatment plan may be evaluated based on the calculated body dose distribution, and the treatment plan optimization process may be performed as described above with reference to FIG. 6 ( S97 ).
- the parameters derived from the optimized treatment plan may be calculated (S98), and the calculated parameters may be applied to the treatment system to perform electric field treatment (S99).
- FIG. 10 is a block diagram of an electric field cancer treatment planning system according to another embodiment of the present invention.
- the electric field cancer treatment planning system 1000 includes an image classification unit 1100 , a physical property information setting unit 1200 , a prescription information determining unit 1300 , and a dose calculating unit. (1400) and may be configured to include a dose optimization unit (1500).
- the image classification unit 1100 may receive a medical image of a patient including organs and tumors and classify organs and tumors included in the medical images.
- the image classification unit 1100 classifies a tumor and each organ from a three-dimensional medical image of a patient including a plurality of organs and tumors, and reconstructs a three-dimensional image for each classified region to form a tumor and a plurality of tumors. It is possible to distinguish the positional relationship between organs.
- the physical property information setting unit 1200 may set physical property information of each area classified by the image classifying unit 1100 (eg, a tumor, an organ, and other human tissues).
- the physical property information may include electrical conductivity and mass density of each tissue in the human body.
- the electrical conductivity may be set based on magnetic resonance imaging (MRI) or a previously established electrical conductivity DB for each human tissue.
- the mass density may be set based on computed tomography (CT) or a pre-established mass density DB for each human tissue.
- the prescription information determiner 1300 may set a prescription method in consideration of the input tumor type and tumor status information, and may determine prescription information including a prescription dose, the number of treatments, a treatment time, a treatment frequency, and the like.
- the prescription method may be any one of the absorbed energy density-based method and the absorbed specific energy-based method, and accordingly, the unit of the prescribed dose may be determined as any one of the absorbed energy density and the absorbed specific energy.
- the prescribed dose can be determined by the absorbed energy delivered to the tissue with respect to a preset reference frequency, and when the frequency used (ie, treatment frequency) is different from the reference frequency, the frequency weight to reflect the biological effect according to the frequency It can be applied to modify the prescribed dose.
- the frequency used ie, treatment frequency
- the dose calculating unit 1400 considers the positions of each region classified by the image classification unit 1100 (eg, tumors, organs, and other human tissues), the number of electrodes, the positions of the electrodes, and the electric field strength for each electrode (that is, the voltage applied to the electrode) is initially set, and the dose distribution in the body can be calculated based on this.
- the dose distribution in the body may be an absorbed energy density distribution or an absorbed specific energy distribution according to a set prescription method, and these may be calculated according to Equations 2 and 3, respectively.
- the dose optimization unit 1500 changes at least one of the number of electrodes of the electrode pad provided in the electric field cancer treatment device, the position of the electrodes, and the electric field strength for each electrode so that the prescribed dose is given to the tumor as much as possible and to the surrounding normal tissue at the minimum. Thus, the dose can be optimized.
- the dose optimization unit 1500 determines the size of the electrode pad by calculating the number of electrodes and the positions of the electrodes in consideration of the dose distribution in the body calculated by the dose calculating unit 1400, and determines the arrangement of the electrodes.
- the dose may be optimized by repeating the process of changing the electric field strength for each electrode (ie, the voltage applied to the electrode).
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- General Health & Medical Sciences (AREA)
- Animal Behavior & Ethology (AREA)
- Veterinary Medicine (AREA)
- Biomedical Technology (AREA)
- Public Health (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Physics & Mathematics (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Biophysics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Pathology (AREA)
- High Energy & Nuclear Physics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Artificial Intelligence (AREA)
- Physiology (AREA)
- Oncology (AREA)
- Hospice & Palliative Care (AREA)
- Optics & Photonics (AREA)
- Theoretical Computer Science (AREA)
- Neurology (AREA)
- Evolutionary Computation (AREA)
- Fuzzy Systems (AREA)
- Mathematical Physics (AREA)
- Psychiatry (AREA)
- Signal Processing (AREA)
- Pulmonology (AREA)
- General Physics & Mathematics (AREA)
- Radiation-Therapy Devices (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
- Medical Treatment And Welfare Office Work (AREA)
- Electrotherapy Devices (AREA)
- Medicines That Contain Protein Lipid Enzymes And Other Medicines (AREA)
Abstract
본 발명의 일 실시예에 따른 흡수에너지 기반 전기장 암치료 계획 시스템은, 장기와 종양이 포함된 환자의 의료영상에서 장기와 종양을 분류하는 영상 분류부; 상기 영상 분류부에 의해 분류된 각 영역의 물성 정보를 설정하는 물성 정보 설정부; 입력된 종양 종류 및 종양 상태 정보를 고려하여 처방 방법을 설정하고, 처방선량, 치료횟수, 치료시간 및 치료 주파수를 포함하는 처방 정보를 결정하는 처방 정보 결정부; 상기 영상 분류부에 의해 분류된 각 영역의 위치를 고려하여 전기장 암치료를 위해 사용되는 전극의 개수, 전극의 위치 및 전극별 전기장 세기를 초기 설정하고 초기 설정을 기초로 체내 선량 분포를 계산하는 선량 산출부; 및 상기 전극의 개수, 전극의 위치 및 전극별 전기장 세기 중 적어도 하나를 변경하여 상기 처방선량이 종양에 최대한 부여되고 주변의 정상조직에는 최소한으로 부여되도록 선량을 최적화하는 선량 최적화부를 포함할 수 있다.
Description
본 출원은 흡수에너지 기반 전기장 암치료 계획 시스템 및 방법에 관한 것으로, 보다 상세하게는 인체 내 각 조직에 부여된 전기장의 세기, 조직의 밀도, 조직의 전기전도도 및 전기장 인가 시간을 기반으로 각 조직에 부여된 물리적인 흡수에너지를 계산하고 이를 기반으로 환자의 처방선량(prescription dose)을 결정하며 결정된 처방선량이 종양에 최대한 잘 부여되고 주변의 정상조직에는 최소한으로 부여될 수 있도록 치료계획을 설립하고 이를 기초로 치료를 수행하는 흡수에너지 기반 전기장 암치료 계획 시스템 및 방법에 관한 것이다.
2000년대 초, 이스라엘 생물리학 교수인 Yoram Palti는 분열하는 암세포에 교류전기장을 걸어주면 분열이 지연되거나 사멸하는 현상을 처음으로 발견했으며, 2004년에는 Cancer Research 저널에 세계 최초로 전기장 암치료 효과에 대한 연구결과를 발표하였다(비특허문헌 1 참조). 그 후 전기장 암치료에 대한 여러 연구논문들이 발표되었으며 암치료 학계에서는 전기장 치료의 3가지 강점으로 인해 큰 주목을 하고 있다.
전기장 치료의 첫 번째 강점은 전기장은 분열하는 세포에만 큰 영향을 미치는 것으로 알려져 있기 때문에 정상세포에 비해 분열속도가 빠른 암세포에 집중적으로 영향을 미치게 되고, 부작용도 기존 치료법과 비교했을 때 매우 미미할 것으로 예상된다는 점이다(비특허문헌 2 참조). 실제로 2013년에 발표된 논문(비특허문헌 3 참조)에 의하면 항암치료와 전기장 치료의 부작용 비교 항목 총 9개 중 7개의 항목에서 전기장 치료의 부작용이 월등히 낮게 나타났으며, 2개의 항목에 대해서는 항암치료와 전기장 치료가 거의 동등한 수준의 부작용을 보였다.
두 번째로, 전기장을 이용한 암치료는 아직 초기 단계임에도 불구하고 치료효능 측면에서 기존의 치료법에 비해 더 좋은 결과를 보여주고 있다. 그 예로, 난치암 중 하나인 악성 교모세포종(glioblastoma multiforme, GBM) 환자가 항암치료만 받은 경우 Progression Free Survival(PFS)과 Overall Survival(OS) 및 2년 이상 생존확률이 각각 4.0개월, 16.0개월, 31%로 나타난 반면, 전기장 치료가 추가된 경우에는 6.7개월, 20.9개월, 43%로 나타나, 기존 치료법에 비해 약 1.7, 1.3, 1.4배 더 좋은 결과를 보여주었다(비특허문헌 4).
마지막으로 세 번째 강점은 전기장이 치료부위를 포함하여 광범위하게 적용될 경우 CT 등의 의료영상에서 보이지 않는 매우 미세한 종양에 대해서도 치료 효과가 있을 것으로 기대된다는 점이다. 종양을 중심으로 전기장을 걸어주게 되면 종양뿐만 아니라 종양주변에도 무시할 수 없는 전기장이 전달되기 때문에, 종양주변에 존재하나 육안으로 확인 불가능할 정도의 미세한 크기의 종양에도 전기장이 영향을 미쳐 암세포의 분열을 억제시킬 수 있으며 이에 따라 암 전이 확률도 획기적으로 낮출 수 있을 것으로 기대된다(비특허문헌 5 참조).
현재 전기장 암치료법은 2011년 재발한 교모세포종 암에 대해, 2015년에는 최초 진단된 교모세포종에 대해 미국 FDA 승인을 받았으며, 유럽에서는 CE 마크를 획득하여 미국, 독일, 스위스 등지에서 약 1000 여개의 병원을 통해 시행되고 있으며, 일본에서도 재발한 교모세포종 환자에 대하여 치료 허가를 받은 상태이다. 또한 치료 받은 환자 수는 매년 급격히 증가하고 있어 2014년에 152명에서 2018년에는 8813명으로 50배 이상의 증가율을 보여주고 있다(비특허문헌 6 참조).
하지만 현재의 전기장 암치료는 치료효과를 극대화하기 위해 몇 가지 해결해야 할 과제가 있다.
현재까지의 연구 결과에 따르면, 암 세포에 대한 전기장 치료의 효과는 조직에 인가된 전기장의 세기와 전기장의 인가 시간에 따라 차이를 나타낸다. 구체적으로, 전기장의 세기 및 인가 시간이 크면 클수록 암세포에 대한 사멸 및 분열 억제 효과가 커서 전기장의 세기와 인가 시간이 전기장의 세포 사멸 효과와 비례 관계에 있다고 할 수 있다(비특허문헌 7 및 8 참조).
현재 상용화된 전기장 암치료 시스템의 경우, 사용하는 전기장의 세기는 피부에 부작용이 없는 범위에서 최대한의 전기장을 인가하는 것으로 제한되고, 치료 시간은 환자가 수용할 수 있는 한 거의 하루 종일(18-24시간/일) 전기장을 인가하여 치료를 진행하고 있는 실정이다(비특허문헌 7 및 9 참조).
즉, 현재의 전기장 치료는 특정 종양을 치료하기 위해 적절한 전기장의 세기나 치료 시간이 정해져 있지 않으며 부작용이 없는 범위에서 최대한의 전기장의 세기와 최대한의 치료 시간으로 치료가 진행되고 있다. 다시 말해, 전기장 치료는 방사선 치료나 항암치료에서 일반적으로 사용하는 선량(dose)이라는 개념이 없으며 종양의 종류나 상태에 따라 다르게 적용되는 처방선량(Prescription dose)도 적용되지 않는 치료법이라고 할 수 있다.
그러나, 이처럼 처방선량 없이 치료를 진행할 경우, 어느 정도 세기의 전기장을 얼마 동안 종양에 인가해야 치료 효과를 나타낼 수 있는지, 전기장의 세기, 전기장 인가 시간, 전기장에 의해 전달된 에너지, 전기장과 치료 시간의 함수 등 에서 치료 효과와 가장 잘 부합하는 항목이 무엇인지 알 수 없으며, 그 기준도 제시할 수 없다. 다시 말해, 단순히 전기장의 세기와 전기장 인가 시간이 치료 효과와 비례 관계임은 짐작할 수 있으나, 그 관계가 선형적(linear)인지 아니면 좀 더 복잡한 관계식을 갖는지는 정확하게 알 수 없기 때문에 전기장 치료에서 치료 계획을 최적화할 수 없으며 치료도 효과적으로 진행할 수 없다.
340명의 뇌종양 환자에 대한 전기장 치료를 수행한 최근 임상 연구 결과에 의하면, 전기장 치료에서 단위 부피당 전력손실의 단위인 전력손실밀도(power loss density)가 종양에서 0.77 mW/cm3 이상일 경우가 그렇지 않은 경우에 비해서 현저히 좋은 예후를 보이고 있다(비특허문헌 10 참조). 비록 전력손실밀도라는 단위가 전기장 치료에 대해 현재까지 제시된 거의 유일한 처방의 기준으로 간주될 수는 있다고 할지라도 이 단위를 전기장 치료의 선량 단위로 사용하기에는 치명적인 한계가 있다. 즉, 암세포의 사멸 정도가 전기장 인가 시간에 비례함에도 불구하고 전력손실밀도에는 환자가 치료 받아야 할 총 치료 시간(즉, 전기장 인가 시간)과 연관되는 항목이 포함되어 있지 않아서 이를 기반으로는 처방을 할 수 없으며 치료도 최적화할 수 없다(비특허문헌 8 참조).
또한, 인체 내에는 다양한 암 종이 발생할 수 있으며 다양한 암 종을 구성하는 각기 다른 세포들은 전기장의 주파수에 따라 각각 상이하게 반응한다(비특허문헌 7 참조). 따라서 인가된 전기장이 세포에 미치는 영향은 전기장의 세기와 인가된 시간에 의해 결정되는 물리량뿐만 아니라 인가된 주파수의 값에 따라 상이하게 반응하는 생물학적인 효과에 따라 결정된다고 할 수 있다.
이러한 이유들로 인해 전력손실밀도는 전기장 치료의 예후를 예측하는 인자로 사용하기에는 적합할 수 있으나 전기장 치료의 선량 단위로 사용하기에는 한계가 있으며 전기장의 인가 시간과 주파수가 추가적으로 고려된 새로운 개념의 전기장 치료 선량의 단위가 필요한 실정이다.
(비특허문헌 1) Eilon D. Kirson et al, disruption of cancer cell replication by alternating electric fields, cancer research, 64, 3288-3295 (2004)
(비특허문헌 2) Miklos Pless, Uri Weinberg, tumor treating fields: concept, evidence, future, Expert Opinion, 20(8), 1099-1106 (2011)
(비특허문헌 3) Angela M. Davies et al, Tumor treating fields: a new frontier in cancer therapy, Annals of the New York academy of sciences, 1291, 86-95 (2013)
(비특허문헌 4) Stupp et al, Effect of Tumor-Treating Fields Plus Maintenance Temozolomide vs Maintenance Temozolomide Alone on Survival in Patients With Glioblastoma: A Randomized Clinical Trial, Journal of the American Medical Association, 318(23), 2306-2316 (2017)
(비특허문헌 5) Eilon D. Kirson et al. Alternating electric fields (TTFields) inhibit metastatic spread of solid tumors to the lungs, Clin Exp Metastasis 26, 633-640 (2009)
(비특허문헌 6) Novocure Corporate Presentation (https://3sj0u94bgxp33grbz1fkt62h-wpengine.netdna-ssl.com/wp-content/uploads/2019/05/201905_NVCR_Corporate_Presentation_vFF.pdf)
(비특허문헌 7) Eilon D. Kirson et al, alternating electric fields arrest cell proliferation in animal tumor models and human brain tumors, PNAS, 104(24), 10152-10157 (2007)
(비특허문헌 8) Yunhui Jo et al, Effectiveness of a Fractionated Therapy Scheme in Tumor Treating Fields Therapy, Technology in Cancer Research & Treatment,18, 1-10 (2019)
(비특허문헌 9) Denise Fabian et al, Treatment of Glioblastoma (GBM) with the Addition of Tumor-Treating Fields (TTF): A Review, Cancers, 11, 174 (2019)
(비특허문헌 10) Matthew T. Ballo et al, Correlation of Tumor Treating Fields Dosimetry to Survival Outcomes in Newly Diagnosed Glioblastoma: A Large-Scale Numerical Simulation-Based Analysis of Data from the Phase 3 EF-14 Randomized Trial, Int J Radiation Oncol Biol Phys, Vol. 104, No. 5, pp. 1106-1113 (2019)
(비특허문헌 11) Dimitris J. Panagopoulos et al, Evaluation of Specific Absorption Rate as a Dosimetric Quantity for Electromagnetic Fields Bioeffects, PLoS One, 8(6), e62663 (2013)
(비특허문헌 12) Stefano Mandija et al. Opening a new window on MR-based Electrical Properties Tomography with deep learning, Scientific Reports, 9, 8895 (2019)
따라서, 당해 기술분야에서는 전기장 치료 시에 사용 가능한 선량 단위를 도입하여 처방선량이 종양에 최대한 잘 인가되고 주변의 정상조직에는 최소한으로 부여될 수 있도록 치료계획을 설립하고 이를 통해 도출된 파라미터들을 전기장 치료기와 연계하여 치료를 수행하도록 함으로써 치료 효과는 극대화하고 부작용은 최소화할 수 있도록 전기장 치료에 대한 선량의 기준을 설정하는 것이 요구되고 있다.
본 발명의 목적은 이상에서 언급한 목적으로 제한되지 않으며, 언급되지 않은 또 다른 목적들은 아래의 기재로부터 본 발명이 속하는 통상의 지식을 가진 자에게 명확히 이해될 수 있을 것이다.
상기 과제를 해결하기 위해서, 본 발명의 일 실시예는 흡수에너지 기반 전기장 암치료 계획 시스템을 제공한다.
상기 흡수에너지 기반 전기장 암치료 계획 시스템은, 장기와 종양이 포함된 환자의 의료영상에서 장기와 종양을 분류하는 영상 분류부; 상기 영상 분류부에 의해 분류된 각 영역의 물성 정보를 설정하는 물성 정보 설정부; 입력된 종양 종류 및 종양 상태 정보를 고려하여 처방 방법을 설정하고, 처방선량, 치료횟수, 치료시간 및 치료 주파수를 포함하는 처방 정보를 결정하는 처방 정보 결정부; 상기 영상 분류부에 의해 분류된 각 영역의 위치를 고려하여 전기장 암치료를 위해 사용되는 전극의 개수, 전극의 위치 및 전극별 전기장 세기를 초기 설정하고 초기 설정을 기초로 체내 선량 분포를 계산하는 선량 산출부; 및 상기 전극의 개수, 전극의 위치 및 전극별 전기장 세기 중 적어도 하나를 변경하여 상기 처방선량이 종양에 최대한 부여되고 주변의 정상조직에는 최소한으로 부여되도록 선량을 최적화하는 선량 최적화부를 포함할 수 있다.
또한, 본 발명의 다른 실시예는 흡수에너지 기반 전기장 암치료 계획 방법을 제공한다.
상기 흡수에너지 기반 전기장 암치료 계획 방법은, 환자의 의료영상에서 종양 및 장기를 포함하는 관심 영역을 분할하는 단계; 흡수에너지밀도 기반 또는 흡수비에너지 기반 중 어느 하나의 방식으로 처방 방법을 설정하는 단계; 설정된 처방 방법에 따라 처방선량, 치료횟수, 치료시간 및 치료 주파수를 포함하는 처방 정보를 결정하는 단계; 종양 및 장기를 고려하여 초기 설정한 전극의 개수, 전극의 위치 및 전극별 전기장 세기를 기초로 체내 선량 분포를 계산하는 단계; 및 계산된 체내 선량 분포를 기반으로 치료계획을 평가 및 최적화하는 단계를 포함할 수 있다.
덧붙여 상기한 과제의 해결수단은, 본 발명의 특징을 모두 열거한 것이 아니다. 본 발명의 다양한 특징과 그에 따른 장점과 효과는 아래의 구체적인 실시형태를 참조하여 보다 상세하게 이해될 수 있을 것이다.
본 발명의 일 실시예에 따르면, 분열하는 세포에 전기장이 인가될 때 세포증식 억제효과와 실질적으로 비례 관계에 있는 흡수에너지를 기반으로 전기장 치료의 선량을 설정할 수 있고, 이를 기반으로 각 종양 치료에 적용 가능한 임상적인 처방선량을 결정할 수 있게 된다. 이를 통해, 현재까지 전기장 치료에서 선량 단위의 부존재로 인하여 불가능했던 체계적인 치료계획의 수립 및 평가가 가능해진다.
더 나아가, 처방선량이 인체 내에서 최적화(Optimization)될 수 있도록 치료계획을 수립하고 이를 통해 도출된 파라미터들을 전기장 치료기에 적용함으로써 처방선량을 종양에 최대한 전달함과 동시에 주변 정상조직에는 최소한으로 전달하도록 치료를 수행할 수 있게 된다.
또한, 임상결과에 대한 통계적인 분석을 진행할 수 있으며, 이를 통해 환자 치료를 위한 합리적인 처방선량 수립이 가능해져서 전기장 암치료의 치료 효과를 극대화할 수 있게 된다.
본 발명의 효과는 이상에서 언급한 효과로 제한되지 않으며, 언급되지 않은 또 다른 효과들은 아래의 기재로부터 본 발명이 속하는 통상의 지식을 가진 자에게 명확히 이해될 수 있을 것이다.
도 1은 전기장 인가 시 전기장 세기에 따른 세포증식 억제효과를 도시하는 도면이다.
도 2는 전기장 인가 시 전력손실밀도와 전기장 인가 시간에 따른 세포증식 억제효과를 도시하는 도면이다.
도 3은 전기장 인가 시 서로 다른 전력손실밀도와 서로 다른 전기장 인가 시간을 적용한 경우 세포증식 억제효과를 도시하는 도면이다.
도 4는 전기장 치료 시 환자의 의료영상을 이용하여 체내 흡수에너지밀도 또는 흡수비에너지 분포를 계산하는 과정을 도시하는 도면이다.
도 5는 인체 모형 팬텀에서 정상 조직의 질량밀도(mass density)를 1 g/cm3로 가정하고 종양(tumor)의 질량밀도가 1.2 g/cm3인 경우와 0.6 g/cm3인 경우에 대해 200kHz의 전기장을 인가한 경우의 흡수에너지밀도와 흡수비에너지의 분포를 도시하는 도면이다.
도 6은 전기장 치료계획 수립 시 흡수에너지 기반 선량 단위를 이용한 체내 선량 최적화 과정을 도시하는 도면이다.
도 7은 인체 모형 팬텀에 가상의 종양조직과 정상세포로 이루어진 중요장기(organ at risk)를 가정하고 전기장 치료 시 흡수에너지 기반 선량 단위를 이용하여 진행된 선량 최적화 결과를 도시하는 도면이다.
도 8은 전기장 인가 시 세포증식 억제효과의 전기장 주파수 의존도를 도시하는 도면이다.
도 9는 본 발명의 일 실시예에 따른 전기장 암치료 계획 방법의 흐름도이다.
도 10은 본 발명의 다른 실시예에 따른 전기장 암치료 계획 시스템의 구성도이다.
이하, 첨부된 도면을 참조하여 본 발명이 속하는 기술분야에서 통상의 지식을 가진 자가 본 발명을 용이하게 실시할 수 있도록 바람직한 실시예를 상세히 설명한다. 다만, 본 발명의 바람직한 실시예를 상세하게 설명함에 있어, 관련된 공지 기능 또는 구성에 대한 구체적인 설명이 본 발명의 요지를 불필요하게 흐릴 수 있다고 판단되는 경우에는 그 상세한 설명을 생략한다. 또한, 유사한 기능 및 작용을 하는 부분에 대해서는 도면 전체에 걸쳐 동일한 부호를 사용한다.
덧붙여, 명세서 전체에서, 어떤 부분이 다른 부분과 '연결'되어 있다고 할 때, 이는 '직접적으로 연결'되어 있는 경우뿐만 아니라, 그 중간에 다른 소자를 사이에 두고 '간접적으로 연결'되어 있는 경우도 포함한다. 또한, 어떤 구성요소를 '포함'한다는 것은, 특별히 반대되는 기재가 없는 한 다른 구성요소를 제외하는 것이 아니라 다른 구성요소를 더 포함할 수 있다는 것을 의미한다.
우선, 본 발명의 실시예에 따른 흡수에너지 기반 전기장 암치료 계획 시스템 및 방법에 대해 설명하기에 앞서, 세포증식 억제효과를 정량화하기 위한 기준에 대해 구체적으로 설명한다.
이를 위해, 도 1 내지 도 3에서는, 악성신경교종 세포주(U373)에 세기 ~1 V/cm, 주파수 150 kHz의 전기장을 인가한 후 세포의 생존능력(viability)을 세포 카운팅(cell counting)을 통해 대조군과 비교하여 전기장을 인가한 세포주(실험군)의 세포증식 억제효과(inhibitory effect of cell proliferation)를 대조군 대비 백분율로 나타내었다. 이 결과는 3회 실험을 진행한 후 평균값을 표시한 것으로, 세포증식 억제효과의 대조군 대비 백분율은 하기의 수학식 1에 의해 산출하였다.
[수학식 1]
세포증식 억제효과 (%) = (1 -실험군의 세포 수/대조군 세포 수) x 100
예를 들어, 전기장이 인가된 실험군의 생존한 세포 수가 대조군 대비 70% 라면 실험군의 대조군 대비 세포증식 억제효과는 30%가 되는 것이다.
도 1은 전기장 인가 시 전기장 세기에 따른 세포증식 억제효과를 도시하는 도면으로, (a)는 150kHz 교류 전기장에 의한 세포증식 억제효과 검증을 위한 실험 장비의 구성을 나타내고, (b)는 전기장 세기(Electric field intensity)에 따른 세포증식 억제효과를 나타내는 그래프이다.
도 1의 (b)를 참조하면, 전기장 인가 시 세포증식 억제효과는 전기장 세기와 비례 관계임을 알 수 있다. 그러나, 세포증식 억제효과와 전기장 세기는 단순 선형(linear) 관계라기 보다는 비선형적인 관계에 가까움을 알 수 있다. (b)에 도시된 곡선은 데이터를 2차 곡선에 피팅한 결과로서, 결정계수 R2 = 0.9916이다.
도 2는 전기장 인가 시 전력손실밀도와 전기장 인가 시간에 따른 세포증식 억제효과를 도시하는 도면으로, (a)는 150kHz 교류 전기장에 의한 전력손실밀도(power density)와 세포증식 억제효과의 관계를 나타내는 그래프이고, (b)는 전력손실밀도가 11.5mW/cm3로 고정된 경우의 전기장 인가 시간(Treatment duration)과 세포증식 억제효과의 관계를 나타내는 그래프이다.
도 2의 (a)를 참조하면, 전기장 인가 시 세포증식 억제효과는 전력손실밀도(Power density)와 비례 관계임을 알 수 있다. 여기서, 전력손실밀도는 1/2 σE2 (σ: 전기전도도, E: 전기장 세기)로 정의되며, (a)에 도시된 곡선은 데이터를 1차 곡선에 피팅한 결과이다.
도 1의 (b)와 도 2의 (a)를 종합적으로 고려하면, 전기장 인가 시 세포증식 억제효과는 전기장 세기보다는 전기장 세기의 제곱에 직접적으로 비례하고, 이에 따라 전력손실밀도에 좀 더 선형적으로 비례함을 알 수 있다(R2 = 0.9759).
또한, 도 2의 (b)를 참조하면, 전기장 인가 시 세포증식 억제효과는 전기장 인가 시간과 비례함을 알 수 있다.
도 3은 전기장 인가 시 서로 다른 전력손실밀도와 서로 다른 전기장 인가 시간을 적용한 경우 세포증식 억제효과를 도시하는 도면으로, (a)는 전력손실밀도가 각각 6.5mW/cm3 및 11.5mW/cm3인 경우 전기장 인가 시간에 따른 세포증식 억제효과를 나타내는 그래프이고, (b)는 단위부피당 또는 단위질량당 흡수에너지에 따른 세포증식 억제효과를 나타내는 그래프이다.
도 3의 (a)를 참조하면, 전기장 인가 시간이 동일한 경우에 전력손실밀도가 큰 경우 보다 높은 세포증식 억제효과를 얻을 수 있음을 알 수 있다.
도 3의 (b)는 세포증식 억제효과와 흡수에너지밀도(즉, 단위부피당 흡수에너지, absorbed energy density) 또는 흡수비에너지(즉, 단위질량당 흡수에너지, specific energy absorbed)와의 관계를 나타내는 것으로, 여기서 흡수에너지밀도 및 흡수비에너지는 각각 하기의 수학식 2 및 수학식 3에 따라 산출할 수 있다.
[수학식 2]
흡수에너지밀도 = 전력손실밀도 × 전기장 인가 시간
[수학식 3]
흡수비에너지 = (전력손실밀도 × 전기장 인가 시간) ÷ (질량밀도)
도 3의 (b)를 참조하면, 전기장 인가 시 전기장을 인가한 세포주의 세포증식 억제효과는 흡수에너지밀도 또는 흡수비에너지와 비례함을 알 수 있다. 세포에 전달된 흡수에너지, 즉 흡수에너지밀도 또는 흡수비에너지는 전력손실밀도와 더불어 세포증식 억제효과와 비례하는 전기장 인가 시간(s)이 포함된 개념으로 전기장 인가 시간을 고려하지 않은 전력손실밀도에 비해 세포증식 억제효과를 정량화하는데 보다 실질적이고 합리적인 기준이라 할 수 있다. 여기서, 흡수에너지밀도와 흡수비에너지는 단위부피당 질량이 고려되었는지 여부에 따라 구분되는 것으로, 양자 모두 전기장 인가 시 세포증식 억제효과를 정량화하는 기준으로 사용될 수 있다.
도 4는 전기장 치료 시 환자의 의료영상을 이용하여 체내 흡수에너지밀도 또는 흡수비에너지 분포를 계산하는 과정을 도시하는 도면이다.
도 4를 참조하면, 우선 환자의 의료영상을 가져와서 관심 영역(Region of interest)을 분할한 후(S41), 치료계획 조건(예를 들어, 전극 위치, 전극별 전기장 세기, 치료 시간 등)을 설정할 수 있다(S42).
이후, 흡수에너지밀도 분포를 계산하고자 하는 경우, 인체 각 조직의 전기전도도 정보를 획득한 후(S43), 상기의 수학식 2에 따라 인체 내 3차원 흡수에너지밀도 분포를 계산할 수 있다(S44).
한편, 흡수비에너지 분포를 계산하고자 하는 경우, 인체 각 조직의 전기전도도 정보를 획득하고(S45), 상기의 수학식 2에 따라 인체 내 3차원 흡수에너지밀도 분포를 계산한 후(S46), 체내 각 조직의 질량 밀도 정보를 획득하고(S47), 흡수에너지밀도를 체내 각 조직의 질량 밀도로 나누어 줌으로써 인체 내 3차원 흡수비에너지 분포를 계산할 수 있다(S48).
도 5는 인체 모형 팬텀에서 정상 조직의 질량밀도(mass density)를 1 g/cm3로 가정하고 종양(tumor)의 질량밀도가 1.2 g/cm3인 경우와 0.6 g/cm3인 경우에 대해 200kHz의 전기장을 인가한 경우의 흡수에너지밀도와 흡수비에너지의 분포를 도시하는 도면으로, (a) 및 (b)는 각각 인체 모형 팬텀의 관상면(coronal plane) 및 축면(axial plane)를 도시하고, (c)는 종양의 질량밀도가 1.2 g/cm3인 경우의 흡수에너지밀도(J/cm3) 분포를 나타내고, (d)는 종양의 질량밀도가 1.2 g/cm3인 경우의 흡수비에너지(J/g) 분포를 나타내며, (e)는 종양의 질량밀도가 0.6 g/cm3인 경우의 흡수에너지밀도(J/cm3) 분포를 나타내고, (f)는 종양의 질량밀도가 0.6 g/cm3인 경우의 흡수비에너지(J/g) 분포를 나타낸다.
도 5의 (c) 및 (d)를 참조하면, 종양의 질량밀도가 1.2 g/cm3인 경우에 주변 정상 조직의 질량밀도인 1 g/cm3과 차이가 작아서 종양(즉, 중앙의 원)에서의 흡수에너지밀도의 분포와 흡수비에너지의 분포가 거의 유사함을 알 수 있다.
반면, 도 5의 (e) 및 (f)를 참조하면, 종양의 질량밀도가 0.6 g/cm3인 경우에는 주변 정상 조직의 질량밀도인 1 g/cm3과 차이가 상대적으로 커서 종양(즉, 중앙의 원)에서의 흡수에너지밀도의 분포와 흡수비에너지의 분포가 확연히 달라짐을 알 수 있다.
이를 토대로, 종양과 정상 조직의 질량밀도의 차이가 클수록 흡수에너지밀도의 분포와 흡수비에너지의 분포가 크게 상이함을 알 수 있으며, 이는 흡수에너지밀도와 흡수비에너지 중 어느 하나를 선택하여 선량의 단위로 일관성 있게 사용할 필요가 있음을 시사한다.
도 6은 전기장 치료계획 수립 시 흡수에너지 기반 선량 단위를 이용한 체내 선량 최적화 과정을 도시하는 도면이다.
도 6을 참조하면, 우선 환자의 의료영상을 가져와서 관심 영역을 분할한 후(S61), 처방선량, 분할 치료 횟수(예를 들어, 1일 치료 시간 및 총 치료 일수 등), 치료 주파수 등을 포함하는 처방 정보를 결정할 수 있다(S62). 여기서, 처방선량은 흡수에너지밀도와 흡수비에너지 중 어느 하나일 수 있다.
이후, 전극 위치, 전극별 전기장 세기 및 치료 시간을 초기 설정한 후(S63), 이를 기반으로 체내 전기장 선량을 계산하고(S64), 계산된 체내 전기장 선량 분포를 평가하여(S65), 기 정해진 기준을 만족하는지 여부를 판단할 수 있다(S66).
계산된 체내 전기장 선량 분포가 기 정해진 기준을 만족하지 못하는 경우, 전극 위치 및 전극별 전기장 세기 중 적어도 하나를 변경한 후(S68), 체내 전기장 선량 분포가 기 정해진 기준을 만족할 수 있도록 하는 최적화 과정을 반복하여 수행할 수 있다(S64 내지 S66).
이러한 최적화 과정을 통해 처방선량(즉, 흡수에너지밀도 또는 흡수비에너지)이 종양에 최대한 잘 부여되고 주변의 정상조직에는 최소한으로 부여되도록 치료를 수행할 수 있다(S67).
도 7은 인체 모형 팬텀에 가상의 종양조직과 정상세포로 이루어진 중요장기(organ at risk)를 가정하고 전기장 치료 시 흡수에너지 기반 선량 단위를 이용하여 진행된 선량 최적화 결과를 도시하는 도면으로, (a)는 종양조직과 중요장기가 포함된 팬텀을 도시하고, (b)는 선량 최적화 이전의 선량분포를 도시하고, (c)는 선량 최적화 이후의 선량분포를 도시하며, (d)는 본 발명의 일 실시예에 따른 전기장 암치료 계획 시스템과 연계된 전기장 암치료 시스템을 도시한다.
구체적으로, 도 7의 (a)에 도시된 바와 같이 인체 모형 팬텀 내에 가상의 ㄷ자형의 종양과 정상 조직인 구 모양의 중요장기를 설정한 후, 200kHz 교류 전기장을 인가하였다.
도 7의 (b) 및 (c)는 각각 도 6에 도시된 체내 선량 최적화 과정을 진행하기 전 및 진행한 후에 종양과 중요장기의 선량분포를 도시하는 것으로, 최적화 과정을 진행하기 전에는 종양과 중요장기에 인가되는 흡수비에너지는 동일하게 ~ 20 J/g인 반면, 최적화 과정을 진행한 후에는 종양과 중요장기에 인가되는 흡수비에너지는 각각 20 J/g과 10~15 J/g으로, 최적화 과정을 통해 중요장기에 인가되는 흡수비에너지를 최소화할 수 있음을 확인할 수 있다.
도 8은 전기장 인가 시 세포증식 억제효과의 전기장 주파수 의존도를 도시하는 도면으로, 췌장암 세포주(CFPAC-I, HPAF-II)에 전기장의 세기가 약 1.2V/cm인 다양한 주파수의 전기장을 48시간 동안 인가하고 주파수에 따라 달라지는 세포증식 억제효과를 150kHz 주파수 대비 백분율로 나타낸 것이다. 즉, 도 8은 전기장의 세기와 전기장 인가 시간은 동일하고 주파수만 상이하므로, 세포에 전달되는 흡수에너지는 모두 동일하고, 주파수에 따른 생물학적 효과를 나타내는 것이다.
도 8을 참조하면, 실험에 사용된 두 개의 세포주는 동일한 양의 흡수에너지가 세포에 전달될 경우 150kHz에서 가장 효과적으로 세포증식이 억제됨을 알 수 있다. 예를 들어, 100kHz의 주파수로 치료할 경우, 150kHz의 주파수로 치료한 경우에 비해 세포증식 억제효과가 약 60 ~ 70% 수준으로 떨어진다. 이러한 주파수 의존 현상은 췌장암 세포뿐만 아니라 모든 종류의 세포에서 일어나는 일반적인 현상이다 (비특허문헌 7 참조).
따라서, 특정 암 종에 대해 흡수에너지를 기반으로 처방선량을 결정하는 경우, 그 처방선량의 기준이 되는 주파수(즉, 기준 주파수)를 함께 설정하는 것이 바람직하며, 기준 주파수가 아닌 다른 주파수를 사용하는 경우 그에 따른 주파수 가중치(Frequency weighting factor)를 적용하여 처방선량을 수정할 필요가 있다. 예를 들어, 기준 주파수에서 처방선량이 100J/g이고 사용 주파수의 주파수 가중치가 2.0이라 가정하면, 처방선량은 200J/g으로 수정하여 설정되어야 한다. 여기서, 사용 주파수에 대한 주파수 가중치는 하기의 수학식 4에 따라 정의될 수 있으며, 사용 주파수에서의 처방선량은 하기의 수학식 5에 따라 수정될 수 있다.
[수학식 4]
주파수 가중치(사용 주파수)= 세포증식 억제효과(기준 주파수)/세포증식 억제효과(사용 주파수)
[수학식 5]
처방선량(사용 주파수) = 처방선량(기준 주파수) × 주파수 가중치
도 9는 본 발명의 일 실시예에 따른 전기장 암치료 계획 방법의 흐름도이다.
도 9를 참조하면, 우선 치료 계획 시스템에 환자의 의료영상(예를 들어, 3차원 의료영상 데이터)을 가져와서(S91) 대상 종양과 중요장기 등의 관심 영역을 분할할 수 있다(S92).
이후, 흡수에너지밀도 기반 또는 흡수비에너지 기반 중 어느 하나의 방식으로 처방 방법을 설정하고(S93, S94), 설정된 처방 방법에 따라 처방선량, 치료횟수, 치료시간, 치료 주파수 등을 포함하는 처방 정보를 결정할 수 있다(S95).
이후, 종양의 위치 및 종양 주변의 중요장기들을 고려하여 전극의 개수, 전극의 위치 및 전극별 전기장 세기를 초기 설정한 후 이를 기초로 체내 선량 분포를 계산할 수 있다(S96).
이후, 계산된 체내 선량 분포를 기반으로 치료계획을 평가하고 도 6을 참조하여 상술한 바와 같이 치료계획 최적화 과정을 수행할 수 있다(S97).
이후, 최적화된 치료계획을 도출한 파라미터들을 산출하고(S98), 산출된 파라미터들을 치료 시스템에 적용하여 전기장 치료를 수행하도록 할 수 있다(S99).
상술한 바와 같은 본 발명의 실시예에 따르면, 특별한 치료 계획 시스템 없이 행해지고 있는 기존의 전기장 치료 방법에 비해 보다 합리적이고 효율적인 치료를 환자에게 제공할 수 있으며, 더 나아가 실제 환자마다 개별 맞춤형 치료를 가능케 하고 이를 이용하여 세포증식 억제효과를 극대화할 수 있다.
도 10은 본 발명의 다른 실시예에 따른 전기장 암치료 계획 시스템의 구성도이다.
도 10을 참조하면, 본 발명의 다른 실시예에 따른 전기장 암치료 계획 시스템(1000)은, 영상 분류부(1100), 물성 정보 설정부(1200), 처방 정보 결정부(1300), 선량 산출부(1400) 및 선량 최적화부(1500)를 포함하여 구성될 수 있다.
영상 분류부(1100)는 장기와 종양이 포함된 환자의 의료영상을 입력받고 의료영상에 포함된 장기와 종양을 분류할 수 있다.
일 예에 따르면, 영상 분류부(1100)는 복수의 장기와 종양이 포함된 환자의 3차원 의료영상으로부터 종양과 각 장기들을 분류하고 분류된 각 영역에 대해 3차원 영상으로 재구성하여 종양과 복수의 장기들 간의 위치 관계를 구분하도록 할 수 있다.
물성 정보 설정부(1200)는 영상 분류부(1100)에 의해 분류된 각 영역(예를 들어, 종양, 장기 및 이 밖의 인체 조직)의 물성 정보를 설정할 수 있다.
여기서, 물성 정보는 인체 내 각 조직의 전기전도도와 질량밀도를 포함할 수 있다. 전기전도도는 자기공명영상(MRI) 또는 기 구축된 인체 조직별 전기전도도 DB를 기반으로 설정될 수 있다. 또한, 질량밀도는 컴퓨터단층영상(CT) 또는 기 구축된 인체 조직별 질량밀도 DB를 기반으로 설정될 수 있다.
처방 정보 결정부(1300)는 입력된 종양 종류 및 종양 상태 정보를 고려하여 처방 방법을 설정하고, 처방선량, 치료횟수, 치료시간, 치료 주파수 등을 포함하는 처방 정보를 결정할 수 있다. 여기서, 처방 방법은 흡수에너지밀도 기반 방식 및 흡수비에너지 기반 방식 중 어느 하나일 수 있으며, 이에 따라 처방선량의 단위가 흡수에너지밀도와 흡수비에너지 중 어느 하나로 결정될 수 있다.
또한, 처방선량은 기 설정된 기준 주파수에 대해 조직에 전달된 흡수에너지에 의해 결정될 수 있으며, 사용 주파수(즉, 치료 주파수)가 기준 주파수와 상이할 경우 주파수에 따른 생물학적 효과를 반영하기 위한 주파수 가중치를 적용하여 처방선량을 수정할 수 있다.
선량 산출부(1400)는 영상 분류부(1100)에 의해 분류된 각 영역(예를 들어, 종양, 장기 및 이 밖의 인체 조직)의 위치를 고려하여 전극의 개수, 전극의 위치 및 전극별 전기장 세기(즉 전극에 인가되는 전압)를 초기 설정하고 이를 기초로 체내 선량 분포를 계산할 수 있다. 여기서, 체내 선량 분포는 설정된 처방 방법에 따라 흡수에너지밀도 분포 또는 흡수비에너지 분포일 수 있으며, 이들은 각각 상술한 수학식 2 및 수학식 3에 따라 계산될 수 있다.
선량 최적화부(1500)는 처방선량이 종양에 최대한 부여되고 주변의 정상조직에는 최소한으로 부여되도록 전기장 암치료기에 구비되는 전극패드의 전극의 개수, 전극의 위치, 전극별 전기장 세기 중 적어도 하나를 변경하여 선량을 최적화할 수 있다.
일 예에 따르면, 선량 최적화부(1500)는 선량 산출부(1400)에 의해 계산된 체내 선량 분포를 고려하여 전극의 개수 및 전극의 위치를 산정하여 전극패드의 크기를 결정하고, 전극의 배열을 변경하거나, 전극별 전기장 세기(즉 전극에 인가되는 전압)를 변경하는 과정을 반복하여 선량을 최적화할 수 있다.
본 발명은 전술한 실시예 및 첨부된 도면에 의해 한정되는 것이 아니다. 본 발명이 속하는 기술분야에서 통상의 지식을 가진 자에게 있어, 본 발명의 기술적 사상을 벗어나지 않는 범위 내에서 본 발명에 따른 구성요소를 치환, 변형 및 변경할 수 있다는 것이 명백할 것이다.
Claims (10)
- 장기와 종양이 포함된 환자의 의료영상에서 장기와 종양을 분류하는 영상 분류부;상기 영상 분류부에 의해 분류된 각 영역의 물성 정보를 설정하는 물성 정보 설정부;입력된 종양 종류 및 종양 상태 정보를 고려하여 처방 방법을 설정하고, 처방선량, 치료횟수, 치료시간 및 치료 주파수를 포함하는 처방 정보를 결정하는 처방 정보 결정부;상기 영상 분류부에 의해 분류된 각 영역의 위치를 고려하여 전기장 암치료를 위해 사용되는 전극의 개수, 전극의 위치 및 전극별 전기장 세기를 초기 설정하고 초기 설정을 기초로 체내 선량 분포를 계산하는 선량 산출부; 및상기 전극의 개수, 전극의 위치 및 전극별 전기장 세기 중 적어도 하나를 변경하여 상기 처방선량이 종양에 최대한 부여되고 주변의 정상조직에는 최소한으로 부여되도록 선량을 최적화하는 선량 최적화부를 포함하는 흡수에너지 기반 전기장 암치료 계획 시스템.
- 제 1 항에 있어서,상기 물성 정보는 인체 내 각 조직의 전기전도도와 질량밀도를 포함하는 것을 특징으로 하는 흡수에너지 기반 전기장 암치료 계획 시스템.
- 제 2 항에 있어서,상기 전기전도도는 자기공명영상(MRI) 또는 기 구축된 인체 조직별 전기전도도 DB를 기반으로 설정되는 것을 특징으로 하는 흡수에너지 기반 전기장 암치료 계획 시스템.
- 제 2 항에 있어서,상기 질량밀도는 컴퓨터단층영상(CT) 또는 기 구축된 인체 조직별 질량밀도 DB를 기반으로 설정되는 것을 특징으로 하는 흡수에너지 기반 전기장 암치료 계획 시스템.
- 제 1 항에 있어서,상기 처방 방법은 흡수에너지밀도 기반 방식 및 흡수비에너지 기반 방식 중 어느 하나이며, 설정된 처방 방법에 따라 상기 처방선량의 단위가 흡수에너지밀도 및 흡수비에너지 중 어느 하나로 결정되는 것을 특징으로 하는 흡수에너지 기반 전기장 암치료 계획 시스템.
- 제 5 항에 있어서,상기 흡수에너지밀도는 하기의 수학식흡수에너지밀도 = 전력손실밀도 × 전기장 인가 시간에 따라 계산하는 것을 특징으로 하는 흡수에너지 기반 전기장 암치료 계획 시스템.
- 제 5 항에 있어서,상기 흡수비에너지는 하기의 수학식흡수비에너지 = (전력손실밀도 × 전기장 인가 시간) ÷ (질량밀도)에 따라 계산하는 것을 특징으로 하는 흡수에너지 기반 전기장 암치료 계획 시스템.
- 제 1 항에 있어서,상기 처방 정보 결정부는 기 설정된 기준 주파수에 대해 조직에 전달된 흡수에너지에 의해 상기 처방선량을 결정하는 것을 특징으로 하는 흡수에너지 기반 전기장 암치료 계획 시스템.
- 제 8 항에 있어서,상기 처방 정보 결정부는 상기 치료 주파수가 상기 기준 주파수와 상이할 경우 주파수에 따른 생물학적 효과를 반영하기 위한 주파수 가중치를 적용하여 상기 처방선량을 수정하는 것을 특징으로 하는 흡수에너지 기반 전기장 암치료 계획 시스템.
- 환자의 의료영상에서 종양 및 장기를 포함하는 관심 영역을 분할하는 단계;흡수에너지밀도 기반 또는 흡수비에너지 기반 중 어느 하나의 방식으로 처방 방법을 설정하는 단계;설정된 처방 방법에 따라 처방선량, 치료횟수, 치료시간 및 치료 주파수를 포함하는 처방 정보를 결정하는 단계;종양 및 장기를 고려하여 초기 설정한 전극의 개수, 전극의 위치 및 전극별 전기장 세기를 기초로 체내 선량 분포를 계산하는 단계; 및계산된 체내 선량 분포를 기반으로 치료계획을 평가 및 최적화하는 단계를 포함하는 흡수에너지 기반 전기장 암치료 계획 방법.
Priority Applications (4)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202080004054.9A CN114269274B (zh) | 2020-07-16 | 2020-07-20 | 基于吸收能量的电场癌治疗计划系统及方法 |
EP20829173.2A EP3967365A4 (en) | 2020-07-16 | 2020-07-20 | ABSORPTION ENERGY ELECTRIC FIELD CANCER TREATMENT PLANNING SYSTEM AND METHOD |
JP2021502968A JP7459047B2 (ja) | 2020-07-16 | 2020-07-20 | システム及び方法 |
US17/259,492 US12011589B2 (en) | 2020-07-16 | 2020-07-20 | System and method for planning electric field cancer treatment based on absorbed energy |
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
KR1020200088365A KR102490645B1 (ko) | 2020-07-16 | 2020-07-16 | 흡수에너지 기반 전기장 암치료 계획 시스템 및 방법 |
KR10-2020-0088365 | 2020-07-16 |
Publications (1)
Publication Number | Publication Date |
---|---|
WO2022014761A1 true WO2022014761A1 (ko) | 2022-01-20 |
Family
ID=79554379
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/KR2020/009527 WO2022014761A1 (ko) | 2020-07-16 | 2020-07-20 | 흡수에너지 기반 전기장 암치료 계획 시스템 및 방법 |
Country Status (6)
Country | Link |
---|---|
US (1) | US12011589B2 (ko) |
EP (1) | EP3967365A4 (ko) |
JP (1) | JP7459047B2 (ko) |
KR (2) | KR102490645B1 (ko) |
CN (1) | CN114269274B (ko) |
WO (1) | WO2022014761A1 (ko) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20220241603A1 (en) * | 2021-02-03 | 2022-08-04 | Novocure Gmbh | Varying Parameters of Tumor Treating Fields (TTFields) Treatment to Overcome Treatment Resistance |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20240038558A (ko) | 2022-09-16 | 2024-03-25 | 주식회사 필드큐어 | 충격량 기반 처방을 이용한 암치료 계획 결정 시스템 및 암치료 계획 평가 시스템 |
EP4338789A1 (en) * | 2022-09-16 | 2024-03-20 | Fieldcure Co., Ltd. | System for determining and evaluating tumor treatment plan using dose based on impulse |
WO2024090939A1 (ko) * | 2022-10-25 | 2024-05-02 | 주식회사 필드큐어 | 전기장 전달 시스템 및 방법 |
CN116258024B (zh) * | 2023-05-16 | 2023-08-29 | 赛福凯尔(绍兴)医疗科技有限公司 | 目标区域的电场可视化方法、装置和计算机设备 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040176804A1 (en) * | 2000-02-17 | 2004-09-09 | Yoram Palti | Apparatus and method for optimizing tumor treatment efficiency by electric fields |
KR20160076703A (ko) * | 2014-12-23 | 2016-07-01 | 재단법인 아산사회복지재단 | 자기장을 이용한 체내 점막조직 선량 제어 광자빔 방사선 치료장치 |
KR101721798B1 (ko) * | 2016-04-26 | 2017-03-30 | 서울대학교산학협력단 | 방사선 치료 계획의 변경 시점 평가 방법 |
WO2017072706A1 (en) * | 2015-10-28 | 2017-05-04 | Zeev Bomzon | Ttfield treatment with optimization of electrode positions on the head based on mri-based conductivity measurements |
KR20200004228A (ko) * | 2018-07-03 | 2020-01-13 | 고려대학교 산학협력단 | 최적화 알고리즘을 이용한 전기장 암치료장치 및 방법 |
Family Cites Families (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6477426B1 (en) * | 2000-06-20 | 2002-11-05 | Celsion Corporation | System and method for heating the prostate gland to treat and prevent the growth and spread of prostate tumors |
US20080214933A1 (en) * | 2005-08-25 | 2008-09-04 | Koninklijke Philips Electronics, N.V. | Image-Based Planning Methods and Apparatus for Targeted Therapy |
JP5485153B2 (ja) * | 2007-08-14 | 2014-05-07 | ノボキュア リミテッド | 電界による寄生生物治療 |
US11272979B2 (en) * | 2008-04-29 | 2022-03-15 | Virginia Tech Intellectual Properties, Inc. | System and method for estimating tissue heating of a target ablation zone for electrical-energy based therapies |
US8992517B2 (en) | 2008-04-29 | 2015-03-31 | Virginia Tech Intellectual Properties Inc. | Irreversible electroporation to treat aberrant cell masses |
US9938498B2 (en) * | 2012-05-07 | 2018-04-10 | Nkmax Co., Ltd. | Method for the induction and expansion of natural killer cells derived from peripheral blood mononuclear cells |
US9655669B2 (en) * | 2013-05-06 | 2017-05-23 | Novocure Limited | Optimizing treatment using TTFields by changing the frequency during the course of long term tumor treatment |
KR101619131B1 (ko) * | 2013-12-20 | 2016-05-10 | 고려대학교 산학협력단 | 방사선 치료 계획 추천 방법 및 시스템 |
CN112807074A (zh) * | 2014-05-12 | 2021-05-18 | 弗吉尼亚暨州立大学知识产权公司 | 电穿孔系统 |
JP6879996B2 (ja) * | 2015-03-26 | 2021-06-02 | オンコセック メディカル インコーポレイテッド | 悪性腫瘍の治療方法 |
US11446487B2 (en) | 2016-09-23 | 2022-09-20 | Beth Israel Deaconess Medical Center, Inc. | System and methods for cancer treatment using alternating electric fields |
TWI632893B (zh) * | 2017-05-16 | 2018-08-21 | 國立陽明大學 | 偵測與分析消化道黏膜組織之方法及系統 |
CN107545137A (zh) * | 2017-08-16 | 2018-01-05 | 强深智能医疗科技(昆山)有限公司 | 肿瘤放射治疗计划智能优化方法 |
US20200368525A1 (en) * | 2017-11-17 | 2020-11-26 | Abbvie Inc. | Methods of Treating Glioblastoma |
JP2021513390A (ja) | 2018-02-02 | 2021-05-27 | ユニバーシティー ヘルス ネットワーク | 腫瘍の視覚化および除去のためのデバイス、システム、および方法 |
KR20210018228A (ko) * | 2018-05-02 | 2021-02-17 | 온코섹 메디컬 인코포레이티드 | 전기천공 시스템, 방법, 및 장치 |
WO2020009306A1 (ko) * | 2018-07-03 | 2020-01-09 | 고려대학교 산학협력단 | 최적화 알고리즘을 이용한 전기장 암치료장치 및 방법 |
EP3909015B1 (en) * | 2019-01-08 | 2024-08-14 | Novocure GmbH | Evaluating quality of segmentation of an image into different types of tissue for planning treatment using tumor treating fields (ttfields) |
EP4014201A1 (en) * | 2019-08-14 | 2022-06-22 | F. Hoffmann-La Roche AG | Three dimensional object segmentation of medical images localized with object detection |
-
2020
- 2020-07-16 KR KR1020200088365A patent/KR102490645B1/ko active IP Right Grant
- 2020-07-20 WO PCT/KR2020/009527 patent/WO2022014761A1/ko unknown
- 2020-07-20 JP JP2021502968A patent/JP7459047B2/ja active Active
- 2020-07-20 EP EP20829173.2A patent/EP3967365A4/en active Pending
- 2020-07-20 CN CN202080004054.9A patent/CN114269274B/zh active Active
- 2020-07-20 US US17/259,492 patent/US12011589B2/en active Active
-
2022
- 2022-11-07 KR KR1020220146753A patent/KR20220152989A/ko active Application Filing
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040176804A1 (en) * | 2000-02-17 | 2004-09-09 | Yoram Palti | Apparatus and method for optimizing tumor treatment efficiency by electric fields |
KR20160076703A (ko) * | 2014-12-23 | 2016-07-01 | 재단법인 아산사회복지재단 | 자기장을 이용한 체내 점막조직 선량 제어 광자빔 방사선 치료장치 |
WO2017072706A1 (en) * | 2015-10-28 | 2017-05-04 | Zeev Bomzon | Ttfield treatment with optimization of electrode positions on the head based on mri-based conductivity measurements |
KR101721798B1 (ko) * | 2016-04-26 | 2017-03-30 | 서울대학교산학협력단 | 방사선 치료 계획의 변경 시점 평가 방법 |
KR20200004228A (ko) * | 2018-07-03 | 2020-01-13 | 고려대학교 산학협력단 | 최적화 알고리즘을 이용한 전기장 암치료장치 및 방법 |
Non-Patent Citations (13)
Title |
---|
ANGELA M. DAVIES ET AL.: "Tumor treating fields: a new frontier in cancer therapy", ANNALS OF THE NEW YORK ACADEMY OF SCIENCES, vol. 1291, 2013, pages 86 - 95 |
DENISE FABIAN ET AL.: "Treatment of Glioblastoma (GBM) with the Addition of Tumor-Treating Fields (TTF): A Review", CANCERS, vol. 11, 2019, pages 174 |
DIMITRIS J. PANAGOPOULOS ET AL.: "Evaluation of Specific Absorption Rate as a Dosimetric Quantity for Electromagnetic Fields Bioeffects", PLOS ONE, vol. 8, no. 6, 2013, pages e62663 |
EILON D. KIRSON ET AL.: "Alternating electric fields (TTFields) inhibit metastatic spread of solid tumors to the lungs", CLIN EXP METASTASIS, vol. 26, 2009, pages 633 - 640, XP019747316, DOI: 10.1007/s10585-009-9262-y |
EILON D. KIRSON ET AL.: "alternating electric fields arrest cell proliferation in animal tumor models and human brain tumors", PNAS, vol. 104, no. 24, 2007, pages 10152 - 10157, XP002730220, DOI: 10.1073/pnas.0702916104 |
EILON D. KIRSON ET AL.: "disruption of cancer cell replication by alternating electric fields", CANCER RESEARCH, vol. 64, 2004, pages 3288 - 3295, XP007907575, DOI: 10.1158/0008-5472.CAN-04-0083 |
HERSHKOVICH H.S.; URMAN N.; NAVEH A.; LEVI S.; BOMZON Z.: "Power Density Loss and Related Measures can be used to Quantify the Dose of Tumor Treating Fields (TTFields)", INTERNATIONAL JOURNAL OF RADIATION: ONCOLOGY BIOLOGY PHYSICS., vol. 102, no. 3, 20 October 2018 (2018-10-20), USA , XP085516268, ISSN: 0360-3016, DOI: 10.1016/j.ijrobp.2018.07.1496 * |
MATTHEW T. BALLO ET AL.: "Correlation of Tumor Treating Fields Dosimetry to Survival Outcomes in Newly Diagnosed Glioblastoma: A Large-Scale Numerical Simulation-Based Analysis of Data from the Phase 3 EF-14 Randomized Trial", INT J RADIATION ONCOL BIOL PHYS, vol. 104, no. 5, 2019, pages 1106 - 1113, XP085741418, DOI: 10.1016/j.ijrobp.2019.04.008 |
MIKLOS PLESSURI WEINBERG: "tumor treating fields: concept, evidence, future", EXPERT OPINION, vol. 20, no. 8, 2011, pages 1099 - 1106, XP055374864, DOI: 10.1517/13543784.2011.583236 |
See also references of EP3967365A4 |
STEFANO MANDIJA ET AL.: "Opening a new window on MR-based Electrical Properties Tomography with deep learning", SCIENTIFIC REPORTS, vol. 9, 2019, pages 8895, XP055656496, DOI: 10.1038/s41598-019-45382-x |
STUPP ET AL.: "Effect of Tumor-Treating Fields Plus Maintenance Temozolomide vs Maintenance Temozolomide Alone on Survival in Patients With Glioblastoma: A Randomized Clinical Trial", JOURNAL OF THE AMERICAN MEDICAL ASSOCIATION, vol. 318, no. 23, 2017, pages 2306 - 2316 |
YUNHUI JO ET AL.: "Effectiveness of a Fractionated Therapy Scheme in Tumor Treating Fields Therapy", TECHNOLOGY IN CANCER RESEARCH & TREATMENT, vol. 18, 2019, pages 1 - 10 |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20220241603A1 (en) * | 2021-02-03 | 2022-08-04 | Novocure Gmbh | Varying Parameters of Tumor Treating Fields (TTFields) Treatment to Overcome Treatment Resistance |
Also Published As
Publication number | Publication date |
---|---|
KR20220152989A (ko) | 2022-11-17 |
CN114269274A (zh) | 2022-04-01 |
EP3967365A4 (en) | 2022-12-28 |
KR20220009764A (ko) | 2022-01-25 |
JP7459047B2 (ja) | 2024-04-01 |
KR102490645B1 (ko) | 2023-01-25 |
JP2022544720A (ja) | 2022-10-21 |
US12011589B2 (en) | 2024-06-18 |
EP3967365A1 (en) | 2022-03-16 |
US20220184391A1 (en) | 2022-06-16 |
CN114269274B (zh) | 2024-06-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
WO2022014761A1 (ko) | 흡수에너지 기반 전기장 암치료 계획 시스템 및 방법 | |
KR102533703B1 (ko) | 흡수에너지 기반 전기장 암치료 계획 시스템 및 방법 | |
Yepes et al. | Fixed-versus variable-RBE computations for intensity modulated proton therapy | |
Descovich et al. | Improving plan quality and consistency by standardization of dose constraints in prostate cancer patients treated with CyberKnife | |
Bijdekerke et al. | TomoTherapy: implications on daily workload and scheduling patients | |
WO2012036351A1 (ko) | 방사선치료계획 평가장치 및 평가방법 | |
Haciislamoglu et al. | Effect of dose constraint on the thyroid gland during locoregional intensity‐modulated radiotherapy in breast cancer patients | |
Lin et al. | Right arcuate fasciculus as outcome predictor after low‐frequency repetitive transcranial magnetic stimulation in nonfluent aphasic stroke | |
Kendall et al. | Comparison of volumetric modulated arc therapy and intensity modulated radiation therapy for whole brain hippocampal sparing treatment plans based on radiobiological modeling | |
Liu et al. | Effects of error on dose of target region and organs at risk in treating nasopharynx cancer with intensity modulated radiation therapy | |
Kuck et al. | Modeling trans-spinal direct current stimulation in the presence of spinal implants | |
Lu et al. | Active and adequate exposure of the facial nerve and chorda tympani nerve to improve the safety of cochlear implantation | |
WO2022265232A1 (ko) | 전기요법에서의 전극 어레이 구조 최적화 방법 및 시스템 | |
Gentilal et al. | The Impact of Scalp's Temperature on the Choice of the Best Layout for TTFields Treatment | |
WO2023054879A1 (ko) | 전기 치료를 위한 전극 시스템, 이를 포함하는 전기치료기 및 전기 치료를 위한 전극 시스템의 최적화 방법 | |
Suzuki et al. | Impact of lifetime attributable risk of radiation-induced secondary cancer in proton craniospinal irradiation with vertebral-body-sparing for young pediatric patients with medulloblastoma | |
Li et al. | Quantifying the dosimetric effects of neck contour changes and setup errors on the spinal cord in patients with nasopharyngeal carcinoma: Establishing a rapid estimation method | |
Suvarnamani et al. | Analyze of the Model for Cancer Transmission | |
Wang et al. | A novel approach to embed eye shields in customized bolus on nasal dorsum treatment for electron radiotherapy | |
Nanda et al. | A Dosimetric Comparison of Volumetric-modulated Arc Therapy and IMRT for Cochlea-sparing Radiation Therapy in Locally Advanced Nasopharyngeal Cancer | |
US20240096474A1 (en) | System for determining and evaluating tumor treatment plan using dose based on impulse | |
Liorni et al. | Computational assessment of pregnant woman models exposed to uniform elf-magnetic fields: Compliance with the European current exposure regulations for the general public and occupational exposures at 50 Hz | |
Mikic et al. | Standardizing Skullremodeling Surgery and Electrode Array Layout to Improve Tumor Treating Fields Using Computational Head Modeling and Finite Element Methods | |
Mikic et al. | Tumor-treating fields dosimetry in glioblastoma: Insights into treatment planning, optimization, and dose–response relationships | |
Ryan et al. | The impact of fluorine-18-fluoroethyltyrosine positron emission tomography scan timing on radiotherapy planning in newly diagnosed patients with glioblastoma |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
ENP | Entry into the national phase |
Ref document number: 2021502968 Country of ref document: JP Kind code of ref document: A |
|
ENP | Entry into the national phase |
Ref document number: 2020829173 Country of ref document: EP Effective date: 20210108 |
|
NENP | Non-entry into the national phase |
Ref country code: DE |