Search

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. [ABSTRACT FROM AUTHOR]

The neural interface is a critical factor in governing efficient and safe charge transfer between a stimulating electrode and biological tissue. The interface plays a crucial role in the efficacy of electric stimulation in chronic implants and both electromechanical properties and biological properties shape this. In the case of cochlear implants, it has long been recognized that neurotrophins can stimulate growth of the target auditory nerve fibers into a favorable apposition with the electrode array, and recently such arrays have been re-purposed to enable electrotransfer (electroporation)-based neurotrophin gene augmentation to improve the "bionic ear." For both this acute bionic array-directed electroporation and for chronic conventional cochlear implant arrays, the electric fields generated in target tissue during pulse delivery are central to efficacy, but are challenging to map. We present a computational model for predicting electric fields generated by array-driven DNA electrotransfer in the cochlea. The anatomically realistic model geometry was reconstructed from magnetic resonance images of the guinea pig cochlea and an eight-channel electrode array embedded within this geometry. The model incorporates a description of both Faradaic and non-Faradaic mechanisms occurring at the electrode-electrolyte interface with frequency dependency optimized to match experimental impedance spectrometry measurements. Our simulations predict that a tandem electrode configuration with four ganged cathodes and four ganged anodes produces three to fourfold larger area in target tissue where the electric field is within the range for successful gene transfer compared to an alternate paired anode-cathode electrode configuration. These findings matched in vivo transfection efficacy of a green fluorescent protein (GFP) reporter following array-driven electrotransfer of the reporter-encoding plasmid DNA. This confirms utility of the developed model as a tool to optimize the safety and efficacy of electrotransfer protocols for delivery of neurotrophin growth factors, with the ultimate aim of using gene augmentation approaches to improve the characteristics of the electrode-neural interfaces in chronically implanted neurostimulation devices. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]