The invention generally relates to anatomically accurate brain phantoms and, more particularly, unique brain phantoms, methods of making brain phantoms, and methods of using brain phantoms. In some embodiments brain phantoms are tailored for transcranial magnetic stimulation (TMS) simulation or other non-invasive neuromodulation techniques such as transcranial direct current stimulation (tDCS) or invasive neuromodulation techniques such as deep brain stimulation (DBS) or combination of invasive and non-invasive neuromodulation techniques.
Transcranial Magnetic Stimulation (TMS) is a non-invasive technique used for treatment and diagnosis of many neurological conditions and diseases. However, the experimental measurement of induced electric fields in the brain tissues is not well established or understood due to non-availability of anatomically realistic head/brain phantoms. The lack of anatomically realistic brain phantoms has made the experimental verification of induced electric fields in the brain tissues an impediment to the development of new treatment protocols.
Transcranial magnetic stimulation (TMS) is a non-invasive treatment for neurological and psychiatric disorders, but the strength of the induced electric field decreases with distance from the TMS coil and the precise nature of induced electrical and magnetic fields may not be safely measured in living patients being treated with TMS. Deep brain effects of TMS are still unclear and not well-studied.
An aspect of some embodiments is a 3-D anatomically realistic brain phantom developed using segmentation of MRIs, construction of volumetric brain models, 3-D printing, casting using composite polymer that mimic the brain conductivities, or a combination thereof. Exemplary phantoms may be used for the purpose of evaluation of neuromodulation techniques such as Transcranial Magnetic Stimulation (TMS) or other non-invasive neuromodulation techniques such as transcranial direct current stimulation (tDCS) or invasive neuromodulation techniques such as DBS or combination of invasive and non-invasive neuromodulation techniques. The phantoms may also be used in the evaluation of quality assurance and quality control (QA/QT) of neuroimaging modalities like magnetic resonance imaging (MRI).
An aspect of some embodiments is the enablement of the professional in the field of brain modulation and treatment to test and preform actual brain stimulations that are accurate and matches the clinical setting of the of the treatments of TMS, tDCS, DBS or combination of them. There are currently no brain heterogenous phantoms that can experimentally verify TMS and tDCS parameters. Embodiments herein comprise brain phantoms and experimental verification of TMS, tDCS, DBS or combination treatment parameters with such brain phantoms.
An exemplary method is capable of producing patient specific three dimensional anatomically realistic head and brain models from MRI data. It may be used to test the safety of neuroimaging and invasive and non-invasive neuromodulation procedures prior to them being performed on the actual patient.
Exemplary embodiments may involve a 3-D anatomically realistic brain phantom that mimics the electrical conduction and mechanical stiffness of the brain. An exemplary phantom is producible using MRI images, software for brain tissue segmentation and image reconstruction, a 3-D printer, and polymer with conductive fillers. Varied loading of fillers may be used to differentiate different types of tissue in the brain or in the head.
In some embodiments, brain tissues of phantoms are divided into cerebrospinal fluid (CSF), white matter (WM), grey matter (GM), ventricles (containing CSF), and cerebellum. In an exemplary phantom according to the present invention, shells are printed for each tissue layer/structure of the brain and head. After printing the shells are filled with a conductive material (e.g., polydimethyl-siloxane (PDMS) or silicon with nanoparticles) that is capable of mimicking the conductive properties of different brain tissues.
In contrast to commercially available concentric spherical phantoms, exemplary embodiments disclosed herein involve producing anatomically accurate phantoms quickly and economically through 3D printing. These phantoms can also be 3D printed for specific patients using their MRI data which can be valuable for complex brain stimulation procedures such as deep brain stimulation (DBS).
In some embodiments, diffusion tensor imaging (DTI) is used to construct 3D fiber tract models usable for either finite element analysis for TMS simulation as well as for integration with physical brain phantoms for real world experimentation (as opposed to a computer simulation). The induced electric field throughout a fiber tract can be assessed, and parts of the deep brain are identifiable as receiving stimulation when the outer cortex of the brain is stimulated.
The patent or application file contains at least one drawing executed in color. Copies of this patent or patent application publication with color drawing(s) will be provided by the Office upon request and payment of the necessary fee.
For some exemplary phantoms according to the present invention, shells are 3D printed for each tissue layer of the brain. After printing, these shells are filled with a conductive material (e.g., silicon or PDMS with nanoparticles) that is configured to mimic the conductive properties of brain tissue. Different conductive material compositions may be used for different tissues (e.g., white matter vs. grey matter vs. CNF; brain matter vs. bone vs. skin).
At block 503 the MRI brain images are segmented and reconstructed using computer program tools. At the time of this application MRI data may have .nii file extensions, for example, and require conversion for 3D modeling prior to supplying the data to a 3D printer. Software suitable for block 503 at the time includes FreeSurfer, FSL, and simNIBS. The output of such programs may be imported to another program (e.g., Meshmixer) to actually make the shells.
Blocks 503 to 508 walk through the production, use, and disposal of shells. In other words, these steps detail both mold making, casting, and mold removal. At block 503 shells are actually manufactured in a tangible form. 3D printing is an exemplary means for producing the shells in a cost effective manner. An exemplary material for the 3D printing process is an ABS (Acrylonitrile butadiene styrene) material. 3D printing may require the printing of supporting structures which do not actually have any anatomical analog. In such case these supporting structures may be removed at block 505 by, for example, chemically dissolving the parts (e.g., with acetone for ABS). Next the conductive “tissue” material is poured between shells (e.g., between and upper and lower shell pair, and/or between an inner and outer shell pair) at block 506, and permitted to cure at block 507. The curing process may involve time during which the chemical composition of the “tissue” material reacts and sets. The curing process may involve exposing the “tissue” material to some form of electromagnetic radiation that triggers curing or just keeping the material at environmental conditions for a finite duration of time.
After the conductive “tissue” material cures, the shells and the conductive material are placed in an appropriate chemical bath (e.g., Acetone) to dissolve all remaining shell material (e.g., ABS) at block 508, leaving only the cast “tissue” material for the phantom behind. The mold-making and casting of blocks 504 to 508 are repeated for subsequent parts. As will be discussed in greater detail below, for some tissue structures a prior casting may be used in place of one or more shells. As a result some tissue structures of the phantom are produced using two or more shells, some with only one shell, and some without any shells. Advantages of this approach are many. Fewer shells means less 3D printing which means lower costs of production. Using a prior casting of an existing part as the “mold” for the next part also means the two tissues will intimately share a boundary and reduce or avoid the possibility of gaps between phantom layers which could negatively affect the conductive behavior across the material-to-material boundary. At the conclusion of process 500 all shells have been removed and a multi-layered brain phantom remains and is ready for use.
The end result of the process illustrated by
At stage 601, inner and outer shells 601′ for grey matter (GM) are 3D printed and any supporting structures related to the printing process (and not anatomy) are removed. At stage 602 conductive material 602′ is poured in between the inner and outer shells. The conductive material cures. At stage 603 the shells are removed. For example, the shells and conductive material are bathed in acetone (or other suitable dissolving agent) to dissolve the shell material, leaving only the cured grey matter material 603′ remaining. The cured grey matter material is configured to be anatomically accurate with (e.g., mimic) grey matter found in naturally occurring grey matter of the brain.
At stage 605, a shell 605′ is 3D printed for the ventricles and any supporting structures related to the printing process removed. At stage 606 conductive material 606′ is poured inside the shell. The conductive material cures. At stage 607 the shell and conductive material are bathed in dissolving agent to dissolve the shell material, leaving only the cured ventricle material 607′. The cured ventricle material is configured to be anatomically accurate with (e.g., mimic) CSF ordinarily found in naturally occurring ventricles of the brain.
At stage 609, the cured and shell-less grey matter material from stage 603 and ventricle material from stage 607 are arranged together with anatomically correct three-dimensional spacing and orientation. The ventricle material is placed in the proper position inside the grey matter (GM) material. At stage 610 conductive material 610′ simulating white matter (WM) is poured in the gaps between the ventricle material and GM and cured. Stage 610 differs notably from stages 602 and 606 in that no shells are necessary or indeed used according to some exemplary embodiments. Instead, the existing ventricle and GM castings serve as the mold for defining the boundaries of the WM, just as in real anatomy. The elimination of any shells to make the WM reduces both time and money costs involved in 3D printing and provides excellent interfaces among the grey matter, white matter, and ventricle material.
At stage 612 a shell 612′ is 3D printed for the cerebellum and any supporting structures related to the printing process removed. At stage 613 conductive material 613′ is poured inside the shell. The conductive material cures. At stage 614 the shell and conductive material are bathed in dissolving agent to dissolve the shell material, leaving only the cured cerebellum material 614′. The cured cerebellum material is configured to be anatomically accurate with (e.g., mimic) a cerebellum found in a naturally occurring brain.
At stage 616 a shell 616′ is 3D printed for the cerebrospinal fluid (CSF) layer that envelopes the brain. At stage 617, the shell, the cerebellum, and the already assembled ventricle/WM/GM pieces are arranged together with anatomically correct three-dimensional spacing and orientation. One or more reference frames may be used to position the cerebellum in the assembly with the right location and orientation. For example, a 3D Cartesian table may be used to accurately assembly the brain regions. At stage 618 conductive material 618′ which mimics CSF properties is poured into the spaces within the shell and cured. At stage 619 the shell is dissolved, leaving a fully cured and assembled unitary brain phantom 620 from the CSF layer and deeper. The cured CSF material 619′ is configured to be anatomically accurate with (e.g., mimic) CSF ordinarily found in a naturally occurring brain. In some embodiments the brain phantom 620 produced at stage 619 may serve as an end product of the manufacturing process. The phantom 620 is suitable for use in simulations, tests, and experimentation relating to open surgery on the brain (e.g., where the skin and bone are removed).
At stage 626 a shell 626′ is 3D printed which defines the outer boundary of the skin surrounding the skull. At stage 627 the phantom from stage 624 is arranged within the shell 626′. At stage 628 the space between the shell and bone layer is filled with material 628′ configured to mimic skin tissue and cured. At stage 629 the shell is removed. The cured skin material 629′ is configured to be anatomically accurate with (e.g., mimic) skin ordinarily found on a naturally occurring head. The final result is a unitary brain phantom 600 which resembles the brain of the original MRI data and has distinguishable parts including skin, bone, CSF, cerebellum, grey matter, white matter, and ventricles (containing CSF).
In the preceding descriptions, removal of shells after material cast within the shell is cured has generally been described as performed via dissolving in a chemical bath. Alternative shell removal techniques may also be used in embodiments. For instance, shells may in some cases be broken and the resulting fragments removed (without any dissolving necessary).
An exemplary phantom is usable for the purpose of evaluation of the neuromodulation such as transcranial magnetic stimulation (TMS). It enables the professional in the field of the brain modulation and treatment to test and perform actual brain stimulations on the phantom that are accurate and match the clinical setting of the of TMS treatment. Prior to the instant invention, no brain phantoms existed to the knowledge of the inventor which were capable of experimentally verifying TMS parameters. An exemplary phantom is examinable under different TMS parameters and suitable for comparison with FEM modelling of induced electric field and magnetic fields in different tissues of the brain. Microelectrodes may be placed at different locations/depths on the phantom to measure the current I and resistance Ω. Since the phantom exhibits same electrical properties of the brain, close readings to actual TMS procedures may be achieved.
At block 803, data is collected as the medical procedure is performed. If the procedure is a neuromodulation procedure, electrical or magnetic properties in the brain phantom may be obtained during the neuromodulation procedure (e.g., see
In connection with the selection process in block 804, the data collected at block 803 may be compared to reference data, and the selection may be based on the comparison. For example, particular stimulatory conditions may be desired at a particular location within a patient's brain. The desired stimulatory conditions may be saved as reference data. TMS may then be performed using initial settings and the stimulatory response measured. The TMS device settings can be adjusted until the desired stimulatory conditions within the brain phantom are reached. The final TMS settings may then be used for actually treating the patient.
At block 805 the subsequent neuromodulation or neuroimaging procedure is actually performed on the patient using the selected parameters.
Exemplary methods like method 800 may be useful for personalized medicine for any patient, human or animal. Method 800 is especially advantageous for patients with unique or abnormal conditions which set them apart from the anatomy typical of most patients. For example, some patients may already have a brain implant such as a DBS device. In such a case any effects of the presence of the existing implant on neuroimaging or neuromodulation may be unknown or uncertain. Method 800 permits a safe and reliable means for assessing how such a procedure may go using an anatomically accurate phantom mimicking the patient's brain and containing a copy of the patient's implant to select the parameters of performing the procedure. The procedure may then be performed on the actual patient with a greater degree of certainty and safety.
In the above descriptions for manufacturing brain phantoms, the materials used for the phantom layers are generally described as conductive materials. Addressing the materials directly, an exemplary conductive material is a silicon or silicone based compound (a compound containing silicon, Si) or PDMS with one or more of graphite, multi walled or single walled carbon nanotubes (MWCNT/SWCNT), and silver nanoparticles and nanowires that is capable of mimicking the electrical conductive properties of different brain tissues based on the respective amounts of these constituents. Some embodiments use specific ratios of all three additives in the silicon base. The silicon base may be, for example, silicon polymer polydimethylsiloxane PDM. The conductivity of layers of an exemplary phantom may be in the range between 0.2-3.0 Sm−1. For the skin layer the conductivity range may be lower, e.g., as low as 0.1 Sm−1. In some other embodiments the layers may each be in the range of 0.2-1.8 Sm−1. In a particular example, the electrical conductivity of different brain tissue that was matched in a phantom was as follows: ventricles & CSF=1.77 Sm−1, GM=0.23 Sm−1, WM=0.24 Sm−1, and cerebellum=0.65 Sm−1.
An “anatomically accurate” brain phantom mimics the brain of a living organism, e.g. a mammalian brain (e.g., a human brain). Anatomically accurate may mean the three dimensional geometry (e.g., sizes, relative sizes, dimensions, relative dimensions, locations or positions, relative locations or positions, etc.) of the phantom matches or substantially matches the three dimensional geometry of a real brain (e.g., an actual mammalian brain). Anatomically accurate may mean one or more electrical properties (e.g., electrical conductivity) of the brain phantom match or substantially match one or more electrical properties of a real brain (e.g., an actual mammalian brain). Anatomically accurate may mean one or more material properties (e.g., mass density, viscosity, etc.) of the brain phantom match or substantially match one or more material properties of a real brain (e.g., an actual mammalian brain). A brain phantom may match or substantially match a real brain if at least one layer/structure of the brain phantom matches or substantially matches the corresponding real brain structure. A brain phantom may match or substantially match a real brain only if all the layers/structures of the brain phantom match or substantially match the real brain. Table 1 below presents exemplary but non-limiting material properties which may be used in a computer simulation or physical brain phantom which is anatomically accurate.
The published literature on the conductivity of healthy adult brain's white matter, grey matter, CSF, skull and skin varies significantly. However, the most trusted values of conductivities for grey and white matter fall in the range of 0.1 to 0.5 S/m. Therefore an anatomically accurate brain phantom may be produced with the conductivities of the grey matter and white matter in the range of 0.1 to 0.5 Sm−1. Different conductivities may be used for different structures/layers/regions of the brain phantom. To achieve different conductivities, different composite polymers may be prepared and used. For example, exemplary brain phantoms or structures/layers thereof may comprise a composite polymer of a silicon-based compound (e.g., PDMS) and carbon nanotubes (in particular multi-walled carbon nanotubes, MWCNTs) with the conductivity/resistivity varied among the structures/layers by variable wt % of the MWCNTs. Table 2 presents the relationship between resistivity and composition of MWCNTs in PDMS.
Resistivity of 300-400 ohms/cm corresponds to 0.3-0.5 S/m (an exemplary target value range for WM and GM). An exemplary composition used for Example 1 below is 11.5 wt % of MWCNTs composition for GM and WM.
Layers or structures of exemplary brain phantoms may have electrical conductivities varied with respect to one another by varying one or more of the materials or compositional ratios with respect to the other layers/structures. For example, different layers or structures may be configured to have different electrical conductivities based on nanotubes of different lengths in one layer versus another layer (e.g., shorter in one layer versus longer in another layer). Different layers or structures may be configured to have different electrical conductivities based on different materials for the nanotubes in one layer versus another layer (e.g., carbon versus silver). Different layers or structures may be configured to have different electrical conductivities based on different types of nanotubes in one layer versus another layer (e.g., single walled versus multi walled nanotubes). In some exemplary embodiments which use a combination of PDMS with MWCNTs, the wt % of CNT in PDMS may be between 10.5 and 15.3 with a resistivity between 1000 and 35 ohm/cm. In some exemplary embodiments which use a combination of PDMS with MWCNTs, the wt % of CNT in PDMS may be 10.5 to 12.5 with a resistivity between 1000 and 300 ohm/cm. In some exemplary embodiments which use a combination of PDMS with MWCNTs, the wt % of CNT in PDMS may be about 11.5 (e.g., 11.5±0.5) with a resistivity of or about 500 ohm/cm. Other exemplary specifications may be used in other embodiments.
In some embodiments, brain phantoms may comprise integrated fiber tracts. An important feature in the brain that is normally (and disadvantageously) ignored by researchers while calculating induced electric fields is fiber tracts due to their anatomical complexities and small dimensions. The fiber tracts are part of the white matter with high conductivity and impart anisotropy to the white matter. In some embodiments, fiber tracts may be made an integral part of a brain phantom.
Following is an example approach for integrating one or more fiber tracts into a brain phantom. High resolution diffusion tensor imaging (DTI) data is collected. Then fiber tracts are extracted from the DTI data using, for example, a graphical toolbox ExploreDTI. Then fiber tracts are constructed using whole brain tractography. Tracts are drawn from a seed region of interest (ROI) such as M1 with parameters such as the following: a seed fractional anisotropy threshold of 0.2, minimum fiber length of 50 mm, and angle threshold of 30 degrees.
The 3D fiber (CAD) models can be exported (e.g., as .STL files) for finite element analysis in a simulation program such as Sim4Life.
The 3D fiber (CAD) models may also be integrated directly into a brain phantom.
Integration of fiber tracts to a brain phantom may involve different or modified manufacturing processes than described above. Integration of fiber tracts in brain phantom needs conductive polymers that have less viscosity. Preparation of soft polymers with high conductivity is a complex process with many variables. The volume fraction/loading factor of conductive materials in a polymer composite increases the viscosity of the polymer and makes it difficult to conform to complex structures, such as fiber tracts, for white and grey matter of the brain. The conductivity of the polymer composite is sensitive to the length of carbon nanotubes or the nanowires used. In some embodiments, carbon nanotubes longer than 50 μm are used to make the hardness of the polymer composite a suitable level. Phantoms may be fabricated using the PDMS/CNT composite in a 3D Bioplotter® which is a multi-nozzle, 3D printer of UV curable soft materials. With a 3D printer like the 3D Bioplotter®, the use of shells and casting may be avoided (although the cost of production may be substantially increased). A 3D printer like the Bioplotter® is used or usable to print complex shapes to simultaneously fabricate all the regions of the phantom (e.g., the parts described above), including the fiber tracts. Electrical and magnetic fields may be induced in the resulting brain phantom as previously described in connection with
Where computer software is discussed herein, it should be understood that such software may be embodied in computer readable instructions which may be provided to one or more processors of a general purpose computer, special purpose computer, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions/acts specified in the description above, in one of the flowcharts, and/or in one or more block diagram blocks. These computer readable program instructions may also be stored in a computer readable storage medium that can direct a computer, a programmable data processing apparatus, and/or other devices to function in a particular manner, such that the computer readable storage medium having instructions stored therein comprises an article of manufacture including instructions which implement aspects of the function/act specified in the flowcharts and/or block diagrams. Embodiments herein may comprise one or more computers, one or more processors, one or more computer readable storage media, and/or appropriate input/output devices therefore, as well as additional supporting hardware as necessary.
Unless the context indicates otherwise, block diagrams and flowcharts are exemplary and may involve fewer or greater number of blocks and/or a different order of items or steps. In some embodiments elements or steps may be concurrent, combined, or otherwise organized differently than is depicted or described.
A phantom was produced by 3-D printed shells for each tissue layer of the brain. Brain tissues were divided into cerebrospinal fluid (CSF), white matter (WM), grey matter (GM), ventricles, and cerebellum. These layers were made into shells defining their geometric spatial boundaries and 3D printed. The shells were then filled with a conductive material (silicon polymer polydimethylsiloxane PDM with electrically conductive filler multi walled carbon nanotubes MWCNT) to impart electrical conduction to the brain phantom. Then, the shells were broken or dissolved to finally produce the brain phantom. The electrical conductivity of the brain phantom tissue was is in the range of 0.4-1.0 Sm−1. The phantom was then examined under different TMS parameters and compared with FEM modelling of induced electric and magnetic fields in the brain.
TMS device coils were positioned on the brain phantom and an oscilloscope probe was placed just underneath the surface of the phantom in order to measure the voltage (phantom probe). Another probe was placed at the same distance (from the coils) as the first probe but outside the phantom so as to measure the voltage induced on the probe just from the TMS coils (reference probe). Then, the magnetic field from the coils was applied. The process was repeated for four distances total: 1, 2, 3, and 4 cm. At each distance measurements were repeated at four different power intensities: 25, 50, 75, and 100%. The brain phantom and experimental set up corresponds with
Next comparing both graphs, experimental (
DTI data was collected from a single subject, healthy 30 year old male patient on a GE 3T HDx scanner with 60 directions. This data was visualized in the graphical toolbox ExploreDTI. Fiber tracts were constructed using whole brain tractography with the following parameters: seed fractional anisotropy threshold of 0.2, minimum fiber length of 50 mm, angle threshold of 30 degrees, and a step size of 1. Tracts were calculated and drawn from a seed ROI drawn on the z-plane around the M1 region. The coordinates of each individual tract were imported into SolidWorks as a curve which could then be extruded to a solid, 3D object with a diameter of 0.25 mm. These 3D fiber models were exported as STL files for finite element analysis to simulate TMS alongside 3D brain tissue models.
The ability to 3D model fiber tracts allows for a more comprehensive study of the deep brain effects of cortical stimulation. Tractography constructs the pathways of fiber tracts, allowing for visualization of connecting regions of the brain. Finite element analysis calculates the magnitude of the electric field at any point on the tract. Knowing this can help determine if and where TMS can be used to stimulate the deep parts of the brain by stimulating the cortical regions and improve therapies for disorders affecting the deep brain.
While exemplary embodiments of the present invention have been disclosed herein, one skilled in the art will recognize that various changes and modifications may be made without departing from the scope of the invention as defined by the following claims.
This application claims the benefit of U.S. Provisional Patent Application No. 62/546,810, filed Aug. 17, 2017, the complete contents of which are herein incorporated by reference.
Number | Date | Country | |
---|---|---|---|
62546810 | Aug 2017 | US |