Comparative Bioactivities of Chemically Modified Fucoidan and λ-Carrageenan toward Cells Encapsulated in Covalently Cross-Linked Hydrogels

The sulfated marine polysaccharides, fucoidan and λ-carrageenan, are known to possess anti-inflammatory, immunomodulatory, and cellular protective properties. Although they hold considerable promise for tissue engineering constructs, their covalent cross-linking in hydrogels and comparative bioactivities to cells are absent from the literature. Thus, fucoidan and λ-carrageenan were modified with methacrylate groups and were covalently cross-linked with the synthetic polymer poly(vinyl alcohol)-methacrylate (PVA-MA) to form 20 wt % biosynthetic hydrogels. Identical degrees of methacrylation were confirmed by 1H NMR, and covalent conjugation was determined by using a colorimetric 1,9-dimethyl-methylene blue (DMMB) assay. Pancreatic beta cells were encapsulated in the hydrogels, followed by culturing in the 3D environment for a prolonged period of 32 days and evaluation of the cellular functionality by live/dead, adenosine 5'-triphosphate (ATP) level, and insulin secretion. The results confirmed that fucoidan and λ-carrageenan exhibited ∼12% methacrylate substitution, which generated hydrogels with stable conjugation of the polysaccharides with PVA-MA. The cells encapsulated in the PVA-fucoidan hydrogels demonstrated consistently high ATP levels over the culture period. Furthermore, only cells in the PVA-fucoidan hydrogels retained glucose responsiveness, demonstrating comparatively higher insulin secretion in response to glucose. In contrast, cells in the PVA-λ-carrageenan and the PVA control hydrogels lost all glucose responsiveness. The present work confirms the superior effects of chemically modified fucoidan over λ-carrageenan on pancreatic beta cell survival and function in covalently cross-linked hydrogels, thereby illustrating the importance of differential polysaccharide structural features on their biological effects.

Simulating impaired left ventricular-arterial coupling in aging and disease: a systematic review

Aortic stenosis, hypertension, and left ventricular hypertrophy often coexist in the elderly, causing a detrimental mismatch in coupling between the heart and vasculature known as ventricular-vascular (VA) coupling. Impaired left VA coupling, a critical aspect of cardiovascular dysfunction in aging and disease, poses significant challenges for optimal cardiovascular performance. This systematic review aims to assess the impact of simulating and studying this coupling through computational models. By conducting a comprehensive analysis of 34 relevant articles obtained from esteemed databases such as Web of Science, Scopus, and PubMed until July 14, 2022, we explore various modeling techniques and simulation approaches employed to unravel the complex mechanisms underlying this impairment. Our review highlights the essential role of computational models in providing detailed insights beyond clinical observations, enabling a deeper understanding of the cardiovascular system. By elucidating the existing models of the heart (3D, 2D, and 0D), cardiac valves, and blood vessels (3D, 1D, and 0D), as well as discussing mechanical boundary conditions, model parameterization and validation, coupling approaches, computer resources and diverse applications, we establish a comprehensive overview of the field. The descriptions as well as the pros and cons on the choices of different dimensionality in heart, valve, and circulation are provided. Crucially, we emphasize the significance of evaluating heart-vessel interaction in pathological conditions and propose future research directions, such as the development of fully coupled personalized multidimensional models, integration of deep learning techniques, and comprehensive assessment of confounding effects on biomarkers.

Neural activity of retinal ganglion cells under continuous, dynamically-modulated high frequency electrical stimulation

Objective.Current retinal prosthetics are limited in their ability to precisely control firing patterns of functionally distinct retinal ganglion cell (RGC) types. The aim of this study was to characterise RGC responses to continuous, kilohertz-frequency-varying stimulation to assess its utility in controlling RGC activity.Approach.We usedin vitropatch-clamp experiments to assess electrically-evoked ON and OFF RGC responses to frequency-varying pulse train sequences. In each sequence, the stimulation amplitude was kept constant while the stimulation frequency (0.5-10 kHz) was changed every 40 ms, in either a linearly increasing, linearly decreasing or randomised manner. The stimulation amplitude across sequences was increased from 10 to 300µA.Main results.We found that continuous stimulation without rest periods caused complex and irreproducible stimulus-response relationships, primarily due to strong stimulus-induced response adaptation and influence of the preceding stimulus frequency on the response to a subsequent stimulus. In addition, ON and OFF populations showed different sensitivities to continuous, frequency-varying pulse trains, with OFF cells generally exhibiting more dependency on frequency changes within a sequence. Finally, the ability to maintain spiking behaviour to continuous stimulation in RGCs significantly reduced over longer stimulation durations irrespective of the frequency order.Significance.This study represents an important step in advancing and understanding the utility of continuous frequency modulation in controlling functionally distinct RGCs. Our results indicate that continuous, kHz-frequency-varying stimulation sequences provide very limited control of RGC firing patterns due to inter-dependency between adjacent frequencies and generally, different RGC types do not display different frequency preferences under such stimulation conditions. For future stimulation strategies using kHz frequencies, careful consideration must be given to design appropriate pauses in stimulation, stimulation frequency order and the length of continuous stimulation duration.

The direct influence of retinal degeneration on electrical stimulation efficacy: Significant implications for retinal prostheses

Photoreceptor loss and inner retinal network remodeling severely impacts the ability of retinal prosthetic devices to create artificial vision. We developed a computational model of a degenerating retina based on rodent data and tested its response to retinal electrical stimulation. This model includes detailed network connectivity and diverse neural intrinsic properties, capable of exploring how the degenerated retina influences the performance of electrical stimulation during the degeneration process. Our model suggests the possibility of quantitatively modulating retinal ON and OFF pathways between phase II and III of retinal degeneration without requiring any differences between ON and OFF RGC intrinsic cellular properties. The model also provided insights about how remodeling events influence stage-dependent differential electrical responses of ON and OFF pathways.Clinical Relevance-This data-driven model can guide future development of retinal prostheses and stimulation strategies that may benefit patients at different stages of retinal disease progression, particularly in the early and mid-stages, thus increasing their global acceptance.

Modulating functionally-distinct vagus nerve fibers using microelectrodes and kilohertz frequency electrical stimulation

Modulation of functionally distinct nerve fibers with bioelectronic devices provides a therapeutic opportunity for various diseases. In this study, we began by developing a computational model including four major subtypes of myelinated fibers and one unmyelinated fiber. Second, we used an intrafascicular electrode to perform kHz-frequency electric stimulation to preferentially modulate a population of fibers. Our model suggests that fiber physical properties and electrode-to-fascicle distance severely impacts stimulus-response relationships. Large diameter fibers (Aα- and Aβ-) were only minimally influenced by the fascicle size and electrode location, while smaller diameter fibers (Aδ-, B- and C-) indicated a stronger dependency.Clinical Relevance- Our findings support the possibility of selectively modulating functionally-distinct nerve fibers using electrical stimulation in a small, localized region. Our model provides an effective tool to design next-generation implantable devices and therapeutic stimulation strategies toward minimizing off-target effects.

Cells-in-Touch: 3D Printing in Reconstruction and Modelling of Microscopic Biological Geometries for Education and Future Research Applications

Additive manufacturing (3D printing) and computer-aided design (CAD) still have limited uptake in biomedical and bioengineering research and education, despite the significant potential of these technologies. The utility of organ-scale 3D-printed models of living structures is widely appreciated, while the workflows for microscopy data translation into tactile accessible replicas are not well developed yet. Here, we demonstrate an accessible and reproducible CAD-based methodology for generating 3D-printed scalable models of human cells cultured in vitro and imaged using conventional scanning confocal microscopy with fused deposition modeling (FDM) 3D printing. We termed this technology CiTo-3DP (Cells-in-Touch for 3D Printing). As a proof-of-concept, we created dismountable CiTo-3DP models of human epithelial, mesenchymal, and neural cells by using selectively stained nuclei and cytoskeletal components. We also provide educational and research context for the presented cellular models. In the future, the CiTo-3DP approach can be adapted to different imaging and 3D printing modalities and comprehensively present various cell types, subcellular structures, and extracellular matrices. The resulting CAD and 3D printed models could be used for a broad spectrum of education and research applications.

Repeatability of brain phase-based magnetic resonance electric properties tomography methods and effect of compressed SENSE and RF shimming

Magnetic resonance electrical properties tomography (MREPT) is an emerging imaging modality to noninvasively measure tissue conductivity and permittivity. Implementation of MREPT in the clinic requires repeatable measurements at a short scan time and an appropriate protocol. The aim of this study was to investigate the repeatability of conductivity measurements using phase-based MREPT and the effects of compressed SENSE (CS), and RF shimming on the precision of conductivity measurements. Conductivity measurements using turbo spin echo (TSE) and three-dimensional balanced fast field echo (bFFE) with CS factors were repeatable. Conductivity measurement using bFFE phase showed smaller mean and variance that those measured by TSE. The conductivity measurements using bFFE showed minimal deviation with CS factors up to 8, with deviation increasing at CS factors > 8. Subcortical structures produced less consistent measurements than cortical parcellations at higher CS factors. RF shimming using full slice coverage 2D dual refocusing echo acquisition mode (DREAM) and full coverage 3D dual TR approaches further improved measurement precision. BFFE is a more optimal sequence than TSE for phase-based MREPT in brain. Depending on the area of the brain being measured, the scan can be safely accelerated with compressed SENSE without sacrifice of precision, offering the potential to employ MREPT in clinical research and applications. RF shimming with better field mapping further improves precision of the conductivity measures.

Simulating the impact of photoreceptor loss and inner retinal network changes on electrical activity of the retina

Objective.A major reason for poor visual outcomes provided by existing retinal prostheses is the limited knowledge of the impact of photoreceptor loss on retinal remodelling and its subsequent impact on neural responses to electrical stimulation. Computational network models of the neural retina assist in the understanding of normal retinal function but can be also useful for investigating diseased retinal responses to electrical stimulation.Approach.We developed and validated a biophysically detailed discrete neuronal network model of the retina in the software package NEURON. The model includes rod and cone photoreceptors, ON and OFF bipolar cell pathways, amacrine and horizontal cells and finally, ON and OFF retinal ganglion cells with detailed network connectivity and neural intrinsic properties. By accurately controlling the network parameters, we simulated the impact of varying levels of degeneration on retinal electrical function.Main results.Our model was able to reproduce characteristic monophasic and biphasic oscillatory patterns seen in ON and OFF neurons during retinal degeneration (RD). Oscillatory activity occurred at 3 Hz with partial photoreceptor loss and at 6 Hz when all photoreceptor input to the retina was removed. Oscillations were found to gradually weaken, then disappear when synapses and gap junctions were destroyed in the inner retina. Without requiring any changes to intrinsic cellular properties of individual inner retinal neurons, our results suggest that changes in connectivity alone were sufficient to give rise to neural oscillations during photoreceptor degeneration, and significant network connectivity destruction in the inner retina terminated the oscillations.Significance.Our results provide a platform for further understanding physiological retinal changes with progressive photoreceptor and inner RD. Furthermore, our model can be used to guide future stimulation strategies for retinal prostheses to benefit patients at different stages of disease progression, particularly in the early and mid-stages of RD.

kHz-frequency electrical stimulation selectively activates small, unmyelinated vagus afferents

BACKGROUND

Vagal reflexes regulate homeostasis in visceral organs and systems through afferent and efferent neurons and nerve fibers. Small, unmyelinated, C-type afferents comprise over 80% of fibers in the vagus and form the sensory arc of autonomic reflexes of the gut, lungs, heart and vessels and the immune system. Selective bioelectronic activation of C-afferents could be used to mechanistically study and treat diseases of peripheral organs in which vagal reflexes are involved, but it has not been achieved.

METHODS

We stimulated the vagus in rats and mice using trains of kHz-frequency stimuli. Stimulation effects were assessed using neuronal c-Fos expression, physiological and nerve fiber responses, optogenetic and computational methods.

RESULTS

Intermittent kHz stimulation for 30 min activates specific motor and, preferentially, sensory vagus neurons in the brainstem. At sufficiently high frequencies (>5 kHz) and at intensities within a specific range (7-10 times activation threshold, T, in rats; 15-25 × T in mice), C-afferents are activated, whereas larger, A- and B-fibers, are blocked. This was determined by measuring fiber-specific acute physiological responses to kHz stimulus trains, and by assessing fiber excitability around kHz stimulus trains through compound action potentials evoked by probing pulses. Aspects of selective activation of C-afferents are explained in computational models of nerve fibers by how fiber size and myelin shape the response of sodium channels to kHz-frequency stimuli.

CONCLUSION

kHz stimulation is a neuromodulation strategy to robustly and selectively activate vagal C-afferents implicated in physiological homeostasis and disease, over larger vagal fibers.

Slotted Monopole Patch Antenna for Microwave-Based Head Imaging Applications

A modified monopole patch antenna for microwave-based hemorrhagic or ischemic stroke recognition is presented in this article. The designed antenna is fabricated on a cost-effective FR-4 lossy material with a 0.02 loss tangent and 4.4 dielectric constant. Its overall dimensions are 0.32 λ × 0.28 λ × 0.007 λ, where λ is the lower bandwidth 1.3 GHz frequency wavelength. An inset feeding approach is utilized to feed the antenna to reduce the input impedance (z = voltage/current). A total bandwidth (below -10 dB) of 2.4 GHz (1.3-3.7 GHz) is achieved with an effective peak gain of over 6 dBi and an efficiency of over 90%. A time-domain analysis confirms that the antenna produces minimal signal distortion. Simulated and experimental findings share a lot of similarities. Brain tissue is penetrated by the antenna to a satisfactory degree, while still exhibiting a safe specific absorption rate (SAR). The maximum SAR value measured for the head model is constrained to be equal to or below 0.1409 W/kg over the entire usable frequency band. Evaluation of theoretical and experimental evidence indicates the intended antenna is appropriate for Microwave Imaging (MWI) applications.

A varying-radius cable equation for the modelling of impulse propagation in excitable fibres

In this study, we present a varying-radius cable equation for nerve fibres taking into account the varying diameter along the neuronal segments. Finite element neuronal models utilising the classical (fixed-radius) and varying-radius cable formulations were compared using simple and realistic morphologies under intra- and extracellular electrical stimulation protocols. We found that the use of the classical cable equation to model intracellular neural electrical stimulation exhibited an error of 17% in a passive resistive cable model with abrupt change in radius from 1 to 2 μm, when compared to the known analytical solution and varying-radius cable formulation. This error was observed to increase substantially using more realistic neuron morphologies and branching structures. In the case of extracellular stimulation however, the difference between the classical and varying-radius formulations was less pronounced, but we expect this difference will increase under more complex stimulation paradigms such as high-frequency stimulation. We conclude that for computational neuroscience applications, it is essential to use the varying-radius cable equation for accurate prediction of neuronal responses under electrical stimulation.

Fluid structure computational model of simulating mitral valve motion in a contracting left ventricle

BACKGROUND

Fluid structure interaction simulations h hold promise in studying normal and abnormal cardiac function, including the effect of fluid dynamics on mitral valve (MV) leaflet motion. The goal of this study was to develop a 3D fluid structure interaction computational model to simulate bileaflet MV when interacting with blood motion in left ventricle (LV).

METHODS

The model consists of ideal geometric-shaped MV leaflets and the LV, with MV dimensions based on human anatomical measurements. An experimentally-based hyperelastic isotropic material was used to model the mechanical behaviour of the MV leaflets, with chordae tendineae and papillary muscle tips also incorporated. LV myocardial tissue was prescribed using a transverse isotropic hyperelastic formulation. Incompressible Navier-Stokes fluid formulations were used to govern the blood motion, and the Arbitrary Lagrangian Eulerian (ALE) method was employed to determine the mesh deformation of the fluid and solid domains due to trans-valvular pressure on MV boundaries and the resulting leaflet movement.

RESULTS

The LV-MV generic model was able to reproduce physiological MV leaflet opening and closing profiles resulting from the time-varying atrial and ventricular pressures, as well as simulating normal and prolapsed MV states. Additionally, the model was able to simulate blood flow patterns after insertion of a prosthetic MV with and without left ventricular outflow tract flow obstruction. In the MV-LV normal model, the regurgitant blood flow fraction was 10.1 %, with no abnormality in cardiac function according to the mitral regurgitation severity grades reported by the American Society of Echocardiography.

CONCLUSION

Our simulation approach provides insights into intraventricular fluid dynamics in a contracting LV with normal and prolapsed MV function, as well as aiding in the understanding of possible complications after transcatheter MV implantation prior to clinical trials.

An artificial intelligence powered platform for auto-analyses of spine alignment irrespective of image quality with prospective validation

BACKGROUND

Assessment of spine alignment is crucial in the management of scoliosis, but current auto-analysis of spine alignment suffers from low accuracy. We aim to develop and validate a hybrid model named SpineHRNet+, which integrates artificial intelligence (AI) and rule-based methods to improve auto-alignment reliability and interpretability.

METHODS

From December 2019 to November 2020, 1,542 consecutive patients with scoliosis attending two local scoliosis clinics (The Duchess of Kent Children's Hospital at Sandy Bay in Hong Kong; Queen Mary Hospital in Pok Fu Lam on Hong Kong Island) were recruited. The biplanar radiographs of each patient were collected with our medical machine EOS™. The collected radiographs were recaptured using smartphones or screenshots, with deidentified images securely stored. Manually labelled landmarks and alignment parameters by a spine surgeon were considered as ground truth (GT). The data were split 8:2 to train and internally test SpineHRNet+, respectively. This was followed by a prospective validation on another 337 patients. Quantitative analyses of landmark predictions were conducted, and reliabilities of auto-alignment were assessed using linear regression and Bland-Altman plots. Deformity severity and sagittal abnormality classifications were evaluated by confusion matrices.

FINDINGS

SpineHRNet+ achieved accurate landmark detection with mean Euclidean distance errors of 2·78 and 5·52 pixels on posteroanterior and lateral radiographs, respectively. The mean angle errors between predictions and GT were 3·18° and 6·32° coronally and sagittally. All predicted alignments were strongly correlated with GT (p < 0·001, R2 > 0·97), with minimal overall difference visualised via Bland-Altman plots. For curve detections, 95·7% sensitivity and 88·1% specificity was achieved, and for severity classification, 88·6-90·8% sensitivity was obtained. For sagittal abnormalities, greater than 85·2-88·9% specificity and sensitivity were achieved.

INTERPRETATION

The auto-analysis provided by SpineHRNet+ was reliable and continuous and it might offer the potential to assist clinical work and facilitate large-scale clinical studies.

FUNDING

RGC Research Impact Fund (R5017-18F), Innovation and Technology Fund (ITS/404/18), and the AOSpine East Asia Fund (AOSEA(R)2019-06).

Secondary flow in bifurcations - Important effects of curvature, bifurcation angle and stents

Coronary bifurcations have complex flow patterns including secondary flow zones and helical flow, which directly affect pathophysiological mechanisms such as the development of atherosclerosis. The objective of this study was to generate insights into the effects of curvature, bifurcation angle and the presence of stents on flow patterns and resulting haemodynamics in coronary left main bifurcations. The blood flow and associated metrics were modelled in both idealised and patient-specific bifurcations with varying curvature and bifurcation angles with and without stents, resulting in a total of 128 geometries considered. The results showed that larger curvature of bifurcating vessels has a significant influence on secondary flow, especially with distance to the bifurcation region, causing a skew, spin and asymmetry of Dean vortices, an increase in helical flow intensity with symmetry loss, and a decrease in adversely low time-average wall shear stress (TAWSS). Generally, asymmetric flow patterns coincided with adversely low TAWSS regions. In identical stented geometries, the presence of the stents induced local recirculation immediately adjacent to the stent struts, thus generating adversely low TAWSS in these areas, with some effect on the overall secondary flow. Overall, the effect of stents outweighed the effect of curvature and BA. This new knowledge contributes to a better understanding of the joint effects of curvature, bifurcation angle, and stents on flow patterns and haemodynamics in coronary bifurcations.

Selective myelinated nerve fiber stimulation via temporal interfering electric fields

We have investigated selective electrical stimulation of myelinated nerve fibers using a computational model of temporal interfering (TI) fields. The model consists of two groups of electrodes placed on the outer bundle surface, each group stimulated at a different frequency. We manipulated the stimulus waveform, magnitude and frequency of short-duration stimuli (70ms), and investigated fiber-specific stimulus-elicited compound action potentials. Results show that under 100Hz & 200Hz TI stimulation with 0.6mA total current shared by the electrodes, continuous action potentials were generated in deeper nerve fibers, and that the firing region was steerable by changing individual electrode currents. This study provides a promising platform for non-invasive nerve bundle stimulation by TI fields.

Effects of modifying the electrode placement and pulse width on cognitive side effects with unilateral ECT: A pilot randomised controlled study with computational modelling

BACKGROUND

The electrode placement and pulse width for electroconvulsive therapy (ECT) are important treatment parameters associated with ECT related retrograde memory side-effects. Modification of these parameters with right unilateral (RUL) ECT may have utility for further reducing these side-effects.

OBJECTIVE

This study explored use of the frontoparietal (FP) placement for reducing retrograde memory side effects with ECT. We hypothesised that superior retrograde memory outcomes would occur with FP compared to temporoparietal (TP) placement and with ultrabrief (UB: 0.3 ms) compared to brief pulse (BP: 1.0 ms) width ECT.

METHODS

In this randomised cross-over, double-blinded study, participants received a single treatment of BP TP, BP FP, UB TP and UB FP ECT. Neuropsychological testing was conducted prior to and immediately following each treatment. Computational modelling was conducted to explore associations between E-fields in regions-of-interest associated with memory.

RESULTS

Nine participants completed the study. The FP placement was not superior to TP for retrograde memory outcomes. For both electrode placements UB pulse width was associated with significantly better visual retrograde memory compared to BP (p < .05). With TP ECT, higher E-fields in regions-of-interest were significantly associated with greater visual retrograde memory side-effects (hippocampi: r = -0.77, p = .04; inferior frontal gyri: r = -0.92, p < .01; middle frontal gyri: r = -0.84, p = .02).

CONCLUSIONS

Modification of pulse-width had greater effects than electrode placement for reducing retrograde memory side-effects with RUL ECT. Preliminary findings suggested that higher E-fields may be associated with greater cognitive side-effects with ECT.

The role of regional myocardial topography post-myocardial infarction on infarct extension

Infarct extension involves necrosis of healthy myocardium in the border zone (BZ), progressively enlarging the infarct zone (IZ) and recruiting the remote zone (RZ) into the BZ, eventually leading to heart failure. The mechanisms underlying infarct extension remain unclear, but myocyte stretching has been suggested as the most likely cause. Using human patient-specific left-ventricular (LV) numerical simulations established from cardiac magnetic resonance imaging (MRI) of myocardial infarction (MI) patients, the correlation between infarct extension and regional mechanics abnormality was investigated by analysing the fibre stress-strain loops (FSSLs). FSSL abnormality was characterised using the directional regional external work (DREW) index, which measures FSSL area and loop direction. Sensitivity studies were also performed to investigate the effect of infarct stiffness on regional myocardial mechanics and potential for infarct extension. We found that infarct extension was correlated to severely abnormal FSSL in the form of counter-clockwise loop at the RZ close to the infarct, as indicated by negative DREW values. In regions demonstrating negative DREW values, we observed substantial fibre stretching in the isovolumic relaxation (IVR) phase accompanied by a reduced rate of systolic shortening. Such stretching in IVR phase in part of the RZ was due to its inability to withstand the high LV pressure that was still present and possibly caused by regional myocardial stiffness inhomogeneity. Further analysis revealed that the occurrence of severely abnormal FSSL due to IVR fibre stretching near the RZ-BZ boundary was due to a large amount of surrounding infarcted tissue, or an excessively stiff IZ.

Multiphysics computational modelling of the cardiac ventricles

Development of cardiac multiphysics models has progressed significantly over the decades and simulations combining multiple physics interactions have become increasingly common. In this review, we summarise the progress in this field focusing on various approaches of integrating ventricular structures. electrophysiological properties, myocardial mechanics, as well as incorporating blood hemodynamics and the circulatory system. Common coupling approaches are discussed and compared, including the advantages and shortcomings of each. Currently used strategies for patient-specific implementations are highlighted and potential future improvements considered.

Fucoidan- and carrageenan-based biosynthetic poly(vinyl alcohol) hydrogels for controlled permeation

Since the permeation of the inflammatory cytokines into hydrogel scaffolds has been shown to cause dysfunction of encapsulated cells, appropriate design strategies to circumvent this are essential. In the present work, it was hypothesized that highly crosslinked PVA-fucoidan and PVA-carrageenan hydrogels can control permeation of the trefoil-shaped inflammatory cytokine IL-1β while allowing the permeation of the globular protein albumin. PVA, fucoidan, and carrageenans were functionalized with methacrylate groups and the functionalized polymers were co-crosslinked by UV photopolymerization. The resultant hydrogels were characterized physicochemically and the release of fucoidan and carrageenans was quantified by developing a colorimetric assay, which was validated by XPS analysis. The permeability characteristics of the hydrogels were evaluated using bovine serum albumin (BSA), IgG, and IL-1β. The results demonstrated an increase in hydrogel swelling through the incorporation of the polysaccharides with minimal overall mass loss. The release studies showed hydrogel stability, where the formulations exhibited ~43% retention of fucoidan and ~60-80% retention of carrageenans consistently up to 7 days. The permeation data revealed very low permeation of IgG and IL-1β through the hydrogels, with <1% permeation after 24 h, while allowing >6% permeation of BSA. These data indicate that such hydrogels can be used as the basis for cytokine-protective implantable devices for clinical applications.

The left anterior right temporal (LART) placement for electroconvulsive therapy: A computational modelling study

Electrode placement in electroconvulsive therapy (ECT) has a major impact on treatment efficacy and cognitive side effects. Left Anterior Right Temporal (LART) is a lesser utilised bilateral montage which may produce more optimal clinical outcomes relative to standard bitemporal ECT. In this study we used computational modelling to explore how stimulation effects from LART and two novel variants (LART - F3 and LART - Frontal) compared to the more common bilateral placements of bitemporal and bifrontal ECT. High resolution finite element human head models were generated from MRI scans of three subjects with Major Depressive Disorder. Differences in regional stimulation were examined through parametric tests for regions of interest and subtraction maps. Compared to bitemporal ECT, LART - Original resulted in significantly greater stimulation of the left cingulate gyrus (hypothesised to be associated with treatment efficacy), and relatively reduced stimulation of the bilateral hippocampi (potentially associated with cognitive side effects). No additional clinical benefit was suggested with the novel LART placements compared to the original LART. The original LART placement is a promising montage for further clinical investigation.

A Computational Model of Functionally-distinct Cervical Vagus Nerve Fibers

Cervical vagus nerve stimulation (VNS) is a neuromodulation therapy used in the treatment of several chronic disorders. In order to maximize the therapeutic effectiveness of VNS, it has become increasingly important to deliver fiber-specific neurostimulation, so that undesired effects can be minimized. Assessing the activation of different vagal fiber types through electrical stimulation is therefore essential for developing fiber-selective VNS therapies. Towards this goal, we conducted in silico investigations using a generic model of functionally distinct nerve fibers and clinically relevant cuff electrodes using COMSOL. Our model is constrained by histological observations from rat cervical vagus nerves and its outputs are validated against averaged compound nerve action potentials (CNAPs) obtained from rat vagus nerve recordings. We propose this model as an effective tool to design fiber-specific stimulation protocols before testing them in experimental animals.

Creation of virtual channels in the retina using synchronous and asynchronous stimulation - a modelling study

OBJECTIVE

Implantable retinal prostheses aim to provide artificial vision to those suffering from retinal degenerative diseases by electrically stimulating the remaining retinal neurons using a multi-electrode array. The spatial resolution of these devices can be improved by creation of so-called virtual channels (VCs) that are commonly achieved through synchronized stimulation of multiple electrodes. It is largely unclear though if VCs can be created using asynchronous stimulation, which was the primary aim of this study.

APPROACH

A computational model of multi-layered retina and epi-retinal dual-electrode stimulation was developed to simulate the neural activity of populations of retinal ganglion cells (RGCs) using the VC strategy under both synchronous and asynchronous stimulation conditions.

MAIN RESULTS

Our simulation suggests that VCs can be created using asynchronous stimulation. VC performance under both synchronous and asynchronous stimulation conditions can be improved by optimizing stimulation parameters such as current intensity, current ratio (α) between two electrodes, electrode spacing and the stimulation waveform. In particular, two VC performance measures; (1) linear displacement of the centroid of RGC activation, and (2) the RGC activation size consistency as a function of different current ratios α, have comparable performance under asynchronous and synchronous stimulation with appropriately selected stimulation parameters.

SIGNIFICANCE

Our findings support the possibility of creating VCs in the retina under both synchronous and asynchronous stimulation conditions. The results provide theoretical evidence for future retinal prosthesis designs with higher spatial resolution and power efficiency whilst reducing the number of current sources required to achieve these outcomes.

Towards Controlling Functionally-Distinct Retinal Ganglion Cells In Degenerate Retina

Present retinal neuroprostheses have limited performance capabilities due to indiscriminate activation of different neural pathways. Based on our success in differentially activating ON and OFF cells using high frequency stimuli in a healthy retina, in this study we explored whether we could achieve similar differential activation between these two cell types but in degenerate retina. We found that after blocking the synaptic network, ON retinal ganglion cells (RGCs) could be differentially activated at higher frequencies (4 - 6 kHz) and amplitudes (200 - 240 µA), and OFF RGCs at relatively lower amplitudes (80 - 160 µA) across all tested frequencies. We further found that both cell types could be controlled by quickly modulating the frequency using short stimulation bursts. This work takes us one step closer to reducing the likelihood of indiscriminate activation of RGCs by accurately controlling the activation of functionally-distinct neural pathways.

Abstract 472: What electric field strength is necessary for maximum tumor-treating fields efficacy

Why Is 4 V/cm Electric Field Strength Necessary for Maximum Tumor-Treating Fields Efficacy? INTRODUCTION: Tumor Treating Fields (TTFields) extend overall survival in glioblastoma patients. Threshold field strength for TTFields is typically cited as 2V/cm. Yet plots of cell death vs. field strength in vitro indicate a power law relationship asymptotic to 4 V/cm for 100% efficacy. The mechanisms of action (MoA) of TTFields underlying this relationship are currently under investigation. TTFields are suspected of disrupting critical mitotic processes performed by large polar cellular molecules such as microtubules (MT). METHODS: Using the power law relationship for the energy (intensity, I) carried by an electric field and its amplitude (Eq. 1), we calculated energy vs. TTFields amplitude. (1) RESULTS: We suggest Eq. 1 is a likely source of the relationship between TTFields' amplitude and cell death proportion that has been demonstrated empirically. Further, Table 1 shows energy transmitted per unit MT surface area per TTFields cycle, specifically, the energy per MT dimer band around an MT helix per TTFields period of 5 µs.

TTFields Field Strength (V/cm)Radiant Flux Density (W/nm2)Energy/(Cycle-Area) (J/nm2)Energy/MT Dimer/Cycle (J)Energy/MT Dimer Band/Cycle (J)113.3 × 10−186.63 × 10−234.01 × 10−224.17 × 10−20253.1 × 10−1826.6 × 10−231.60 × 10−211.67 × 10−193119 × 10−1859.7 × 10−233.61 × 10−213.75 × 10−194212 × 10−18106 × 10−236.42 × 10−216.67 × 10−195332 × 10−18166 × 10−231.00 × 10−201.04 × 10−186478 × 10−18239 × 10−231.44 × 10−201.50 × 10−18

A metric for disruption of cellular processes is 1-2 orders of magnitude greater than the thermal background energy in the cell, kT = 4.3 × 10−21 J. Table 1 shows that the energy absorbed by a MT dimer is insufficient to disrupt cellular processes at 2 V/cm, while due to energy increasing with square of field strength, at 5 V/cm it may be disruptive. For an entire dimer band around the MT helix, 2 V/cm appears to be sufficient, while 4 V/cm seems amply strong to disrupt MT function, in accordance with in vitro experiments. CONCLUSION: Assuming the power law of energy vs. electric field strength and the resulting energy absorbed by microtubule components, we found a correlation between typical TTFields' amplitudes and estimated disruption thresholds of MT function. TTFields MoA. 4 V/cm, not 2V/cm, should be the target field strength to realize the full efficacy of TTFields.

Citation Format: Kristen W. Carlson, Jack A. Tuszynski, Socrates Dokos, Nirmal Paudel, Zeev Bomzon. What electric field strength is necessary for maximum tumor-treating fields efficacy [abstract]. In: Proceedings of the Annual Meeting of the American Association for Cancer Research 2020; 2020 Apr 27-28 and Jun 22-24. Philadelphia (PA): AACR; Cancer Res 2020;80(16 Suppl):Abstract nr 472.

MRI-SegFlow: a novel unsupervised deep learning pipeline enabling accurate vertebral segmentation of MRI images

Most deep learning based vertebral segmentation methods require laborious manual labelling tasks. We aim to establish an unsupervised deep learning pipeline for vertebral segmentation of MR images. We integrate the sub-optimal segmentation results produced by a rule-based method with a unique voting mechanism to provide supervision in the training process for the deep learning model. Preliminary validation shows a high segmentation accuracy achieved by our method without relying on any manual labelling.The clinical relevance of this study is that it provides an efficient vertebral segmentation method with high accuracy. Potential applications are in automated pathology detection and vertebral 3D reconstructions for biomechanical simulations and 3D printing, facilitating clinical decision making, surgical planning and tissue engineering.

Optimizing Stimulation Strategies for Retinal Electrical Stimulation: a Modelling Study

In this research, a continuum multi-compartmental model of retinal electrical stimulation was utilized to find the best strategy for activating retinal ganglion cells (RGCs). Two types of return electrodes configuration placed suprachoroidally were used: monopolar and hexapolar. The current was delivered either simultaneously or sequentially with two kinds of waveforms: biphasic symmetric charge-balanced cathodic and anodic first pulses. Our results revealed there is no significant difference in current threshold between single monopolar and hexapolar stimulation regardless of the applied current stimulus waveform. Moreover, sequential stimulation for both monopolar or hexapolar was more effective in reducing current threshold than simultaneous stimulation when biphasic cathodic first pulses were used. Concurrent monopolar stimulation was significant in reducing the current threshold compared to single monopolar whereas concurrent hexapolar did not alter the current threshold. Overall, concurrent monopolar stimulation was efficacious in reducing current threshold regardless of the stimulus waveforms and sequential stimulation was more useful only with biphasic cathodic first pulses.

XLIF interbody cage reduces stress and strain of fixation in spinal reconstructive surgery in comparison with TLIF cage with bilateral or unilateral fixation: a computational analysis

In treating recalcitrant low back pain, extreme lateral lumbar interbody fusion (XLIF) with a large cage is reported to have better stability compared to approach of transforaminal lumbar interbody fusion (TLIF) using a small cage. In addition, bilateral pedicle screw fixation (PSF) in comparison with unilateral fixation achieved no inferior fusion rate, but with a significant reduction in operation time and blood loss. The aim of the study was to understand the mechanism underpinning the stability of lumbar interbody fusion using different cage sizes with unilateral or bilateral PSF. A computer model of human lumbar vertebrae L4 and L5 with implants was reconstructed based on CT scans and simulated in Ansys Workbench. Simulation results demonstrated that for either XLIF or TLIF cages, the maximum values of rod stress were comparable with bilateral and unilateral PSF. However, the stability was considerably reduced with unilateral PSF for TLIF due to significantly increased facet joint strain for TLIF; whereas for XLIF with left unilateral PSF, the max facet joint strain was comparable to bilateral PSF, possibly due to facet tropism of this specific subject.

The role of end-diastolic myocardial fibre stretch on infarct extension

Myocardial infarct extension, a process involving the enlargement of infarct and border zone, leads to progressive degeneration of left ventricular (LV) function and eventually gives rise to heart failure. Despite carrying a high risk, the causation of infarct extension is still a subject of much speculation. In this study, patient-specific LV models were developed to investigate the correlation between infarct extension and impaired regional mechanics. Subsequently, sensitivity analysis was performed to examine the causal factors responsible for the impaired regional mechanics observed in regions surrounding the infarct and border zone. From our simulations, fibre strain, fibre stress and fibre stress-strain loop (FSSL) were the key biomechanical variables affected in these regions. Among these variables, only FSSL was correlated with infarct extension, as reflected in its work density dissipation (WDD) index value, with high WDD indices recorded at regions with infarct extension. Impaired FSSL is caused by inadequate contraction force generation during the isovolumic contraction and ejection phases. Our further analysis revealed that the inadequacy in contraction force generation is not necessarily due to impaired myocardial intrinsic contractility, but at least in part, due to inadequate muscle fibre stretch at end-diastole, which depresses the ability of myocardium to generate adequate contraction force in the subsequent systole (according to the Frank-Starling law). Moreover, an excessively stiff infarct may cause its neighbouring myocardium to be understretched at end-diastole, subsequently depressing the systolic contractile force of the neighbouring myocardium, which was found to be correlated with infarct extension.

Computational evaluation of an extra-aortic elastic-wrap applied to simulated aging anisotropic human aorta models

Structural changes occurring to the aortic wall can result in vascular stiffening. This is represented by a loss of vascular compliance during pulsatile flow, resulting in increased systolic and pulse blood pressure, particularly in populations aged 50 and over. Aortic stiffness is thought to be permanent and an active de-stiffening strategy is yet to be developed. Extra aortic elastic wrapping has been proposed as a surgical technique to boost aortic distensibility and treat hypertension in the elderly. Previously, in-vivo and in-vitro testing have suggested a pulse-pressure reduction potential of elastic wrapping in the stiffened aortas. Herein, we explore the feasibility of elastic aortic wrapping to improve simulated aortic compliance across the age span. Detailed computational studies of the anisotropic aortic wall mechanics, using data from human subjects, were performed, evaluating key performance properties for the interaction between the aortic wall and elastic aortic wrap procedure. Main determinants of the procedure’s efficiency are identified using a pre-defined aortic stiffness and wrap elasticity. Finite element analysis predicts that segmental aortic distensibility can be increased if elastic wrapping is applied to a simulated stiff aorta. Elastic aortic wrapping is calculated to have little impact on the compliance of an initially distensible aorta.

SCIDOT-52. IN SITU DELIVERY OF TUMOR TREATING FIELDS THERAPY IN GLIOBLASTOMA PATIENTS

The efficacy of Tumor-Treating Fields (TTFields) rests on the result of a large-scale clinical trial that demonstrated an increase in the survival of newly diagnosed glioblastoma patients when combined with temozolomide chemotherapy. Overall survival now extends to over 60 months in some of our patients when dexamethasone, which we suspected of interference with TTFields effects, is replaced with celecoxib to control tumor-associated inflammation. The transcranial method of delivering TTFields has not changed in light of ongoing advances in deep brain stimulation (DBS) and transcranial electric stimulation (TES), notably that the resistivity of the skull is the principle obstacle to placing therapeutic electric field strength of 2 V/cm into target tumor sites and variation in skull thickness is the main difference in TES efficiency across individuals. Realistic human head finite element modeling (FEM) predicted that surgical craniectomy beneath TTFields’ electrodes would enhance field strength at target tumor sites. Here we show that 2 V/cm can be reliably delivered to tumor sites using minimally-invasive DBS cylindrical leads or ribbon electrode arrays, pre- or post-resection. Two objections arise to the in situ method: 1) Will TTFields stimulate axons in situ? 2) Will field strength exceed safety limits for cell damage? Neural stimulation modeling and experiments show that TTFields’ frequency of 200 kHz, 1–3 orders of magnitude higher than ion channel time constants, is too high to stimulate them. Furthermore, 2 V/cm is well below cell damage limits of 700 V/mm. Thus we propose a new delivery method to improve tumor control in glioblastoma patients and to provide valuable information on TTFields’ effects via cell studies using in situ electric fields at 200 kHz.

A systematic review and computational modelling analysis of unilateral montages in electroconvulsive therapy

OBJECTIVE

To examine the clinical outcomes of ECT unilateral placements compared in prior studies and apply insights from computational modelling to understand differences between placements.

METHODS

PubMed, Embase, Scopus and PsycINFO and reference lists were systematically searched for studies of depressed patients where two unilateral placements were compared and clinical outcomes were reported. Computational modelling was done to generate electric field maps for each unilateral placement identified in the systematic review.

RESULTS

A total of 29 studies met criteria for inclusion. Eight studies reported efficacy outcomes and 23 studies reported cognitive outcomes. Most studies found no significant difference in efficacy between right unilateral (RUL) and left unilateral (LUL) ECT, and no difference was found between temporo-parietal and fronto-temporal ECT. For the majority of studies, RUL placements had better verbal anterograde memory outcomes compared with the LUL placements. There was some evidence suggestive of cognitive advantages for fronto-frontal and fronto-parietal placements relative to temporo-parietal ECT.

CONCLUSIONS

For efficacy, studies mainly focused on the comparison of right vs. left hemispheric stimulation, with the available evidence suggesting no substantive difference. RUL placements tended to have better verbal anterograde memory outcomes relative to LUL placements, though limited differences were found between the RUL placements.

Computational Models And Tools For Developing Sophisticated Stimulation Strategies For Retinal Neuroprostheses

Improvements to the efficacy of retinal neuroprostheses can be achieved by developing more sophisticated neural stimulation strategies to enable selective or preferential activation of specific retinal ganglion cells (RGCs). Computational models are particularly well suited for these investigations. The electric field can be accurately described by mathematical formalisms, and the population-based neural responses to the electrical stimulation can be investigated at resolutions well beyond those achievable by current state-of-the-art biological techniques. In this study, we used a biophysically-and morphologically-detailed RGC model to explore the ability of high frequency electrical stimulation (HFS) to preferentially activate ON and OFF RGC subtypes, the two major information pathways of the retina. The performance of a wide range of electrical stimulation amplitudes (0 - $100~\mu \mathbf {A}$) and frequencies (1 - 10 kHz) on functionally-distinct RGC responses were evaluated. We found that ON RGCs could be preferentially activated at relatively higher stimulation amplitudes $( > 50 {\mu } \mathrm {A})$ and frequencies $( >2$ kHz) while OFF RGCs were activated by lower stimulation amplitudes (10 to $50 {\mu } \mathrm {A})$ across all tested frequencies. These stimuli also show great promise in eliciting RGC responses that parallel RGC encoding: one RGC type exhibited an increase in spiking activity during electrical stimulation whilst another exhibited decreased spiking activity, given the same stimulation parameters.

Computational Modelling of the Effect of Infarct Stiffness on Regional Myocardial Mechanics

Ventricular remodeling after myocardial infarction increases the rate of mortality and is highly associated with the extent of infarct transmurality. It is hypothesized that infarct stiffness alters regional mechanics and affects the likelihood of human ventricular remodeling. However, this is yet to be studied in detail. In this paper, we present simulations from an actively-contracting left ventricular model to investigate the effects of transmural infarct stiffness on myofiber regional mechanics. Results show that higher infarct stiffness reduces systolic stress at the infarct and border zones, minimizing infarct bulging but increasing the diastolic stress at the endocardial border zone. Determining a proper amount of infarct stiffness is required to achieve a balanced regional mechanics across the cardiac cycle that may be useful in therapy, such as myocardial hydrogel injection to adjust its stiffness and reduce stress to prevent ventricular remodeling.

Impact of myocardial infarction on intraventricular vortex and flow energetics assessed using computational simulations

Flow energetics have been proposed as early indicators of progressive left ventricular (LV) functional impairment in patients with myocardial infarction (MI), but its correlation with individual MI parameters has not been fully explored. Using electro-fluid-structure interaction LV models, this study investigated the correlation between four MI parameters: infarct size, infarct multiplicity, regional enhancement of contractility at the viable myocardium area (RECVM), and LV mechanical dyssynchrony (LVMD) with intraventricular vortex and flow energetics. In LV with small infarcts, our results showed that infarct appearance amplified the energy dissipation index (DI), where substantial viscous energy loss was observed in areas with high flow velocity and near the infarct-vortex interface. The LV with small multiple infarcts and RECVM showed remarkable DI increment during systole and diastole. In correlation analysis, the systolic kinetic energy fluctuation index (E') was positively related to ejection fraction (EF) (R²  = 0.982) but negatively correlated with diastolic E' (R²  = 0.970). Diastolic E' was inversely correlated with vortex kinetic energy (R²  = 0.960) and vortex depth (R²  = 0.876). We showed an excessive systolic DI could differentiate infarcted LV with normal EF from healthy LV. Strong flow acceleration, LVMD, and vortex-infarct interactions were predominant factors that induced excessive DI in infarcted LVs. Instead of causing undesired flow turbulence, high systolic E' suggested the existence of energetic flow acceleration, while high diastolic E' implied an inefficient diastolic filling. Thus, systolic E' is not a suitable early indicator for progressive LV dysfunction in MI patients, while diastolic E' may be a useful index to indicate diastolic impairment in these patients.

Mediating Retinal Ganglion Cell Spike Rates Using High-Frequency Electrical Stimulation

Recent retinal studies have directed more attention to sophisticated stimulation strategies based on high-frequency (>1.0 kHz) electrical stimulation (HFS). In these studies, each retinal ganglion cell (RGC) type demonstrated a characteristic stimulus-strength-dependent response to HFS, offering the intriguing possibility of focally targeting retinal neurons to provide useful visual information by retinal prosthetics. Ionic mechanisms are known to affect the responses of electrogenic cells during electrical stimulation. However, how these mechanisms affect RGC responses is not well understood at present, particularly when applying HFS. Here, we investigate this issue via an in silico model of the RGC. We calibrate and validate the model using an in vitro retinal preparation. An RGC model based on accurate biophysics and realistic representation of cell morphology, was used to investigate how RGCs respond to HFS. The model was able to closely replicate the stimulus-strength-dependent suppression of RGC action potentials observed experimentally. Our results suggest that spike inhibition during HFS is due to local membrane hyperpolarization caused by outward membrane currents near the stimulus electrode. In addition, the extent of HFS-induced inhibition can be largely altered by the intrinsic properties of the inward sodium current. Finally, stimulus-strength-dependent suppression can be modulated by a wide range of stimulation frequencies, under generalized electrode placement conditions. In vitro experiments verified the computational modeling data. This modeling and experimental approach can be extended to further our understanding on the effects of novel stimulus strategies by simulating RGC stimulus-response profiles over a wider range of stimulation frequencies and electrode locations than have previously been explored.

A Multi-Domain Continuum Model of Electrical Stimulation of Healthy and Degenerate Retina

A continuum multi-domain model of electrical stimulation of the retina is presented and validated against retinal ganglion cell (RGC) excitation thresholds reported in a recently published in vitro experimental study. We applied our model to investigate the response of the RGC layer to electrical stimulation during mid-to-late stage retinal degeneration for both epiretinal and suprachoroidal configurations. Interestingly, our model predicted that suprachoroidal stimulation of the degenerate retina required increased current thresholds, mainly because of the presence of the glial scar layer. In contrast, epiretinal stimulation thresholds were almost similar for both healthy and degenerate models. The latter finding implies that there is no influence of the glial scar layer on epiretinal stimulation current thresholds.

A Multiphysics Biventricular Cardiac Model: Simulations With a Left-Ventricular Assist Device

Computational models have become essential in predicting medical device efficacy prior to clinical studies. To investigate the performance of a left-ventricular assist device (LVAD), a fully-coupled cardiac fluid-electromechanics finite element model was developed, incorporating electrical activation, passive and active myocardial mechanics, as well as blood hemodynamics solved simultaneously in an idealized biventricular geometry. Electrical activation was initiated using a simplified Purkinje network with one-way coupling to the surrounding myocardium. Phenomenological action potential and excitation-contraction equations were adapted to trigger myocardial contraction. Action potential propagation was formulated within a material frame to emulate gap junction-controlled propagation, such that the activation sequence was independent of myocardial deformation. Passive cardiac mechanics were governed by a transverse isotropic hyperelastic constitutive formulation. Blood velocity and pressure were determined by the incompressible Navier-Stokes formulations with a closed-loop Windkessel circuit governing the circulatory load. To investigate heart-LVAD interaction, we reduced the left ventricular (LV) contraction stress to mimic a failing heart, and inserted a LVAD cannula at the LV apex with continuous flow governing the outflow rate. A proportional controller was implemented to determine the pump motor voltage whilst maintaining pump motor speed. Following LVAD insertion, the model revealed a change in the LV pressure-volume loop shape from rectangular to triangular. At higher pump speeds, aortic ejection ceased and the LV decompressed to smaller end diastolic volumes. After multiple cycles, the LV cavity gradually collapsed along with a drop in pump motor current. The model was therefore able to predict ventricular collapse, indicating its utility for future development of control algorithms and pre-clinical testing of LVADs to avoid LV collapse in recipients.

A continuum model of electrical stimulation of multi-compartmental retinal ganglion cells

A continuum multi-domain model of electrical stimulation of the retina is presented. Each point in the retinal ganglion cell layer could be thought of as representing a single cell, whose biophysics is described using a four-compartment formulation incorporating varying ion channel expressions in the soma, axon initial segment, dendrites and axon. Our continuum model was validated against a discrete morphologically-realistic OFF RGC model, using intra- and extra-cellular electrical stimulation scenarios. Simulations from the continuum model reproduced the same results as that of the discrete model. Our continuum model is the first multi-domain model to represent all main RGC compartments, not just the soma. Moreover, we demonstrated that this model allows the investigation of axonal activation which has been observed to influence the perception of phosphenes.

Abstract 3209: Numerical modeling of intracellular mechanisms in tumor-treating fields

Tumor Treating Fields (TTFields) are 100-500 kHz electric fields with intensities of about 1-4 V/cm. They are known to exert an antimitotic effect on cancer cells, most likely by exerting forces on highly polar tubulin dimers, thereby disrupting spindle formation. Calculations show that TTFields-tubulin interaction energy is negligible compared to the thermal energy in the cell (1). Therefore, this interaction is unlikely to disrupt cellular function. Conductivity of polymerized tubulin, microtubules (MTs), was measured to be 20 S/m, which is 400 times higher than that of the ambient cytosol (0.05 S/m) (2). Thus when TTFields penetrate the cytosol, they may induce electric currents along MTs that are strong enough to disrupt key cellular functions. In particular, if the power (energy per unit time) deposited by these currents is on par with that the power consumed by the molecular motor kinesin, then TTFields may disrupt the forces needed for cell division, thereby disrupting mitosis. To test this hypothesis, we performed numerical simulations evaluating the magnitude of the electric current along MTs exposed to TTFields at 200 kHz. Based on studies of MTs and their microenvironment, we model the MT as a layered cylindrical structure (1): Innermost is the lumen (15 nm in thickness), surrounded by 13 strands of alpha-beta tubulin dimers linked in a helix (4.5 nm). C-termini extend out from the helix with a thickness of 3.5 nm. MTs carry net negative charge; thus they are surrounded by a counter-ion layer (2 nm), and an outer nonconductive Bjerrum layer (3 nm). We built a finite element model in COMSOL Multiphysics (tm) incorporating these layers and examined the current density induced in each layer by TTFields for MTs varying in length from 1-10 µm within an ambient 200 kHz AC electric field of 1-4 V/cm. Our model shows that MTs act as electrical "shunts" that conduct electric current within them. The highest current flows through the counter-ion layer surrounding the C-termini. The current density in this layer exceeds the level likely to disrupt the motor protein kinesin "walk" along the C-termini. The current density is highest when both the field and the MTs are aligned with the cell axis, in accord with in vitro experiments (3). Our model is consistent with the hypothesis that when cells are exposed to TTFields, MTs act as cables carrying high-density electric currents strong enough to disrupt the function of molecular motors, ultimately disrupting mitosis.

References:

1. Tuszynski JA et al. An overview of sub-cellular mechanisms involved in the action of TTFields. Int J Environ Res Public Health 2016.

2. Santelices IB et al. Response to alternating electric fields of tubulin dimers and microtubule ensembles in electrolytic solutions. Sci Rep 2017.

3. Kirson ED et al. Alternating electric fields arrest cell proliferation in animal tumor models and human brain tumors. Proc Natl Acad Sci U S A 2007.

Citation Format: Kristen W. Carlson, Jack A. Tuszynski, Socrates Dokos. Numerical modeling of intracellular mechanisms in tumor-treating fields [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2018; 2018 Apr 14-18; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2018;78(13 Suppl):Abstract nr 3209.

Differential electrical responses in retinal ganglion cell subtypes: effects of synaptic blockade and stimulating electrode location

OBJECTIVE

Visual prostheses have shown promising results in restoring visual perception to blind patients. The ability to differentially activate retinal ganglion cell (RGC) subtypes could further improve the efficacy of these medical devices.

APPROACH

Using whole-cell patch clamp, we investigated membrane potential differences between ON and OFF RGCs in the mouse retina when their synaptic inputs were blocked by synaptic blockers, and examined the differences in stimulation thresholds under such conditions. By injecting intracellular current, we further confirmed the relationship between RGC stimulation thresholds and resting membrane potentials (RMPs). In addition, we investigated the effects of stimulating electrode location on the differences in stimulation thresholds between ON and OFF RGCs.

MAIN RESULTS

With synaptic blockade, ON RGCs became significantly more hyperpolarized (from  -61.8  ±  1.4 mV to  -70.8  ±  1.6 mV), while OFF RGCs depolarized slightly (from  -60.5  ±  0.7 mV to  -58.6  ±  0.9 mV). RGC stimulation thresholds were negatively correlated with their RMPs (Pearson r value:  -0.5154; p-value: 0.0042). Thus, depriving ON RGCs of synaptic inputs significantly increased their thresholds (from 14.7  ±  1.3 µA to 22.3  ±  2.1 µA) over those of OFF RGCs (from 13.2  ±  0.7 µA to 13.1  ±  1.1 µA). However, with control solution, ON and OFF RGC stimulation thresholds were not significantly different. Finally, placement of the stimulating electrode away from the axon enhanced differences in stimulation thresholds between ON and OFF RGCs, facilitating preferential activation of OFF RGCs.

SIGNIFICANCE

Since ON and OFF RGCs have antagonistic responses to natural light, achieving differential RGC activation could convey more natural visual information, leading to better visual prosthesis outcomes.

Closed-Loop Efficient Searching of Optimal Electrical Stimulation Parameters for Preferential Excitation of Retinal Ganglion Cells

The ability for visual prostheses to preferentially activate functionally-distinct retinal ganglion cells (RGCs) is important for improving visual perception. This study investigates the use of high frequency stimulation (HFS) to elicit RGC activation, using a closed-loop algorithm to search for optimal stimulation parameters for preferential ON and OFF RGC activation, resembling natural physiological neural encoding in response to visual stimuli. We evaluated the performance of a wide range of electrical stimulation amplitudes and frequencies on RGC responses in vitro using murine retinal preparations. It was possible to preferentially excite either ON or OFF RGCs by adjusting amplitudes and frequencies in HFS. ON RGCs can be preferentially activated at relatively higher stimulation amplitudes (>150 μA) and frequencies (2-6.25 kHz) while OFF RGCs are activated by lower stimulation amplitudes (40-90 μA) across all tested frequencies (1-6.25 kHz). These stimuli also showed great promise in eliciting RGC responses that parallel natural RGC encoding: ON RGCs exhibited an increase in spiking activity during electrical stimulation while OFF RGCs exhibited decreased spiking activity, given the same stimulation amplitude. In conjunction with the in vitro studies, in silico simulations indicated that optimal HFS parameters could be rapidly identified in practice, whilst sampling spiking activity of relevant neuronal subtypes. This closed-loop approach represents a step forward in modulating stimulation parameters to achieve appropriate neural encoding in retinal prostheses, advancing control over RGC subtypes activated by electrical stimulation.