Development and evaluation of a computerized algorithm for the interpretation of pulmonary function tests.

Pulmonary function tests (PFTs) are usually interpreted by clinicians using rule-based strategies and pattern recognition. The interpretation, however, has variabilities due to patient and interpreter errors. Most PFTs have recognizable patterns that can be categorized into specific physiological defects. In this study, we developed a computerized algorithm using the python package (pdfplumber) and validated against clinicians' interpretation. We downloaded PFT reports in the electronic medical record system that were in PDF format. We digitized the flow volume loop (FVL) and extracted numeric values from the reports. The algorithm used FEV1/FVC<0.7 for obstruction, TLC<80%pred for restriction and <80% or >120%pred for abnormal DLCO. The algorithm also used a small airway disease index (SADI) to quantify late expiratory flattening of the FVL to assess small airway dysfunction. We devised keywords for the python Natural Language Processing (NLP) package (spaCy) to identify obstruction, restriction, abnormal DLCO and small airway dysfunction in the reports. The algorithm was compared to clinicians' interpretation in 6,889 PFTs done between March 1st, 2018, and September 30th, 2020. The agreement rates (Cohen's kappa) for obstruction, restriction and abnormal DLCO were 94.4% (0.868), 99.0% (0.979) and 87.9% (0.750) respectively. In 4,711 PFTs with FEV1/FVC≥0.7, the algorithm identified 190 tests with SADI < lower limit of normal (LLN), suggesting small airway dysfunction. Of these, the clinicians (67.9%) also flagged 129 tests. When SADI was ≥ LLN, no clinician's reports indicated small airway dysfunction. Our results showed the computerized algorithm agreed with clinicians' interpretation in approximately 90% of the tests and provided a sensitive objective measure for assessing small airway dysfunction. The algorithm can improve efficiency and consistency and decrease human errors in PFT interpretation. The computerized algorithm works directly on PFT reports in PDF format and can be adapted to incorporate a different interpretation strategy and platform.

Prokaryotic voltage-gated sodium channels are more effective than endogenous Nav1.5 channels in rescuing cardiac action potential conduction: an in silico study.

Methods to augment Na+ current in cardiomyocytes hold potential for the treatment of various cardiac arrhythmias involving conduction slowing. Because the gene coding cardiac Na+ channel (Nav1.5) is too large to fit in a single adeno-associated virus (AAV) vector, new gene therapies are being developed to enhance endogenous Nav1.5 current (by overexpression of chaperon molecules or use of multiple AAV vectors) or to exogenously introduce prokaryotic voltage-gated Na+ channels (BacNav) whose gene size is significantly smaller than that of the Nav1.5. In this study, based on experimental measurements in heterologous expression systems, we developed an improved computational model of the BacNav channel, NavSheP D60A. We then compared in silico how NavSheP D60A expression vs. Nav1.5 augmentation affects the electrophysiology of cardiac tissue. We found that the incorporation of BacNav channels in both adult guinea pig and human cardiomyocyte models increased their excitability and reduced action potential duration. When compared with equivalent augmentation of Nav1.5 current in simulated settings of reduced tissue excitability, the addition of the BacNav current was superior in improving the safety of conduction under conditions of current source-load mismatch, reducing the vulnerability to unidirectional conduction block during premature pacing, preventing the instability and breakup of spiral waves, and normalizing the conduction and ECG in Brugada syndrome tissues with mutated Nav1.5. Overall, our studies show that compared with a potential enhancement of the endogenous Nav1.5 current, expression of the BacNav channels with their slower inactivation kinetics can provide greater anti-arrhythmic benefits in hearts with compromised action potential conduction.NEW & NOTEWORTHY Slow action potential conduction is a common cause of various cardiac arrhythmias; yet, current pharmacotherapies cannot augment cardiac conduction. This in silico study compared the efficacy of recently proposed antiarrhythmic gene therapy approaches that increase peak sodium current in cardiomyocytes. When compared with the augmentation of endogenous sodium current, expression of slower-inactivating bacterial sodium channels was superior in preventing conduction block and arrhythmia induction. These results further the promise of antiarrhythmic gene therapies targeting sodium channels.

Rule-based definition of muscle bundles in patient-specific models of the left atrium.

Atrial fibrillation (AF) is the most common arrhythmia encountered clinically, and as the population ages, its prevalence is increasing. Although the CHA2DS2- VASc score is the most used risk-stratification system for stroke risk in AF, it lacks personalization. Patient-specific computer models of the atria can facilitate personalized risk assessment and treatment planning. However, a challenge faced in creating such models is the complexity of the atrial muscle arrangement and its influence on the atrial fiber architecture. This work proposes a semi-automated rule-based algorithm to generate the local fiber orientation in the left atrium (LA). We use the solutions of several harmonic equations to decompose the LA anatomy into subregions. Solution gradients define a two-layer fiber field in each subregion. The robustness of our approach is demonstrated by recreating the fiber orientation on nine models of the LA obtained from AF patients who underwent WATCHMAN device implantation. This cohort of patients encompasses a variety of morphology variants of the left atrium, both in terms of the left atrial appendages (LAAs) and the number of pulmonary veins (PVs). We test the fiber construction algorithm by performing electrophysiology (EP) simulations. Furthermore, this study is the first to compare its results with other rule-based algorithms for the LA fiber architecture definition available in the literature. This analysis suggests that a multi-layer fiber architecture is important to capture complex electrical activation patterns. A notable advantage of our approach is the ability to reconstruct the main LA fiber bundles in a variety of morphologies while solving for a small number of harmonic fields, leading to a comparatively straightforward and reproducible approach.

Development and Evaluation of a Small Airway Disease Index Derived From Modeling the Late-Expiratory Flattening of the Flow-Volume Loop.

Excessive decrease in the flow of the late expiratory portion of a flow volume loop (FVL) or "flattening", reflects small airway dysfunction. The assessment of the flattening is currently determined by visual inspection by the pulmonary function test (PFT) interpreters and is highly variable. In this study, we developed an objective measure to quantify the flattening. We downloaded 172 PFT reports in PDF format from the electronic medical records and digitized and extracted the expiratory portion of the FVL. We located point A (the point of the peak expiratory flow), point B (the point corresponding to 75% of the expiratory vital capacity), and point C (the end of the expiratory portion of the FVL intersecting with the x-axis). We did a linear fitting to the A-B segment and the B-C segment. We calculated: 1) the AB-BC angle (∠ABC), 2) BC-x-axis angle (∠BCX), and 3) the log ratio of the BC slope over the vertical distance between point A and x-axis [log (BC/A-x)]. We asked an expert pulmonologist to assess the FVLs and separated the 172 PFTs into the flattening and the non-flattening groups. We defined the cutoff value as the mean minus one standard deviation using data from the non-flattening group. ∠ABC had the best concordance rate of 80.2% with a cutoff value of 149.7°. We then asked eight pulmonologists to evaluate the flattening with and without ∠ABC in another 168 PFTs. The Fleiss' kappa was 0.320 (lower and upper confidence intervals [CIs]: 0.293 and 0.348 respectively) without ∠ABC and increased to 0.522 (lower and upper CIs: 0.494 and 0.550) with ∠ABC. There were 147 CT scans performed within 6 months of the 172 PFTs. Twenty-six of 55 PFTs (47.3%) with ∠ABC <149.7° had CT scans showing small airway disease patterns while 44 of 92 PFTs (47.8%) with ∠ABC ≥149.7° had no CT evidence of small airway disease. We concluded that ∠ABC improved the inter-rater agreement on the presence of the late expiratory flattening in FVL. It could be a useful addition to the assessment of small airway disease in the PFT interpretation algorithm and reporting.

Engineered bacterial voltage-gated sodium channel platform for cardiac gene therapy.

Therapies for cardiac arrhythmias could greatly benefit from approaches to enhance electrical excitability and action potential conduction in the heart by stably overexpressing mammalian voltage-gated sodium channels. However, the large size of these channels precludes their incorporation into therapeutic viral vectors. Here, we report a platform utilizing small-size, codon-optimized engineered prokaryotic sodium channels (BacNav) driven by muscle-specific promoters that significantly enhance excitability and conduction in rat and human cardiomyocytes in vitro and adult cardiac tissues from multiple species in silico. We also show that the expression of BacNav significantly reduces occurrence of conduction block and reentrant arrhythmias in fibrotic cardiac cultures. Moreover, functional BacNav channels are stably expressed in healthy mouse hearts six weeks following intravenous injection of self-complementary adeno-associated virus (scAAV) without causing any adverse effects on cardiac electrophysiology. The large diversity of prokaryotic sodium channels and experimental-computational platform reported in this study should facilitate the development and evaluation of BacNav-based gene therapies for cardiac conduction disorders.

Application of Machine Learning in Pulmonary Function Assessment Where Are We Now and Where Are We Going?

Analysis of pulmonary function tests (PFTs) is an area where machine learning (ML) may benefit clinicians, researchers, and the patients. PFT measures spirometry, lung volumes, and carbon monoxide diffusion capacity of the lung (DLCO). The results are usually interpreted by the clinicians using discrete numeric data according to published guidelines. PFT interpretations by clinicians, however, are known to have inter-rater variability and the inaccuracy can impact patient care. This variability may be caused by unfamiliarity of the guidelines, lack of training, inadequate understanding of lung physiology, or simply mental lapses. A rules-based automated interpretation system can recapitulate expert's pattern recognition capability and decrease errors. ML can also be used to analyze continuous data or the graphics, including the flow-volume loop, the DLCO and the nitrogen washout curves. These analyses can discover novel physiological biomarkers. In the era of wearables and telehealth, particularly with the COVID-19 pandemic restricting PFTs to be done in the clinical laboratories, ML can also be used to combine mobile spirometry results with an individual's clinical profile to deliver precision medicine. There are, however, hurdles in the development and commercialization of the ML-assisted PFT interpretation programs, including the need for high quality representative data, the existence of different formats for data acquisition and sharing in PFT software by different vendors, and the need for collaboration amongst clinicians, biomedical engineers, and information technologists. Hurdles notwithstanding, the new developments would represent significant advances that could be the future of PFT, the oldest test still in use in clinical medicine.

High-Resolution Measurement of Local Activation Time Differences From Bipolar Electrogram Amplitude.

Localized changes in myocardial conduction velocity (CV) are pro-arrhythmic, but high-resolution mapping of local CV is not yet possible during clinical electrophysiology procedures. This is in part because measurement of local CV at small spatial scales (1 mm) requires accurate annotation of local activation time (LAT) differences with very high temporal resolution (≤1 ms), beyond that of standard clinical methods. We sought to develop a method for high-resolution measurement of LAT differences and validate against existing techniques. First, we use a simplified theoretical model to identify a quantitative relationship between the LAT difference of a pair of electrodes and the peak amplitude of the bipolar EGM measured between them. This allows LAT differences to be calculated from bipolar EGM peak amplitude, by a novel "Determination of EGM Latencies by Transformation of Amplitude" (DELTA) method. Next, we use simulated EGMs from a computational model to validate this method. With 1 kHz sampling, LAT differences less than 4 ms were more accurately measured with DELTA than by standard LAT annotation (mean error 3.8% vs. 22.9%). In a 1-dimensional and a 2-dimension model, CV calculations were more accurate using LAT differences found by the DELTA method than by standard LAT annotation (by unipolar dV/dt timing). DELTA-derived LAT differences were more accurate than standard LAT annotation in simulated complex fractionated EGMs from a model incorporating fibrosis. Finally, we validated the DELTA method in vivo using 18,740 bipolar EGMs recorded from the left atrium of 10 atrial fibrillation patients undergoing catheter ablation. Using clinical EGMs, there was agreement in LAT differences found by DELTA, standard LAT annotation, and unipolar waveform cross-correlation. These results demonstrate an underlying relationship between a bipolar EGM's peak amplitude and the activation time difference between its two electrodes. Our computational modeling and clinical results suggest this relationship can be leveraged clinically to improve measurement accuracy for small LAT differences, which may improve CV measurement at small spatial scales.

Excitation properties of computational models of unmyelinated peripheral axons.

Biophysically based computational models of nerve fibers are important tools for designing electrical stimulation therapies, investigating drugs that affect ion channels, and studying diseases that affect neurons. Although peripheral nerves are primarily composed of unmyelinated axons (i.e., C-fibers), most modeling efforts focused on myelinated axons. We implemented the single-compartment model of vagal afferents from Schild et al. (1994) (Schild JH, Clark JW, Hay M, Mendelowitz D, Andresen MC, Kunze DL. J Neurophysiol 71: 2338-2358, 1994) and extended the model into a multicompartment axon, presenting the first cable model of a C-fiber vagal afferent. We also implemented the updated parameters from the Schild and Kunze (1997) model (Schild JH, Kunze DL. J Neurophysiol 78: 3198-3209, 1997). We compared the responses of these novel models with those of three published models of unmyelinated axons (Rattay F, Aberham M. IEEE Trans Biomed Eng 40: 1201-1209, 1993; Sundt D, Gamper N, Jaffe DB. J Neurophysiol 114: 3140-3153, 2015; Tigerholm J, Petersson ME, Obreja O, Lampert A, Carr R, Schmelz M, Fransén E. J Neurophysiol 111: 1721-1735, 2014) and with experimental data from single-fiber recordings. Comparing the two models by Schild et al. (1994, 1997) revealed that differences in rest potential and action potential shape were driven by changes in maximum conductances rather than changes in sodium channel dynamics. Comparing the five model axons, the conduction speeds and strength-duration responses were largely within expected ranges, but none of the models captured the experimental threshold recovery cycle-including a complete absence of late subnormality in the models-and their action potential shapes varied dramatically. The Tigerholm et al. (2014) model best reproduced the experimental data, but these modeling efforts make clear that additional data are needed to parameterize and validate future models of autonomic C-fibers.NEW & NOTEWORTHY Peripheral nerves are primarily composed of unmyelinated axons, and there is growing interest in electrical stimulation of the autonomic nervous system to treat various diseases. We present the first cable model of an unmyelinated vagal nerve fiber and compare its ion channel isoforms and conduction responses with other published models of unmyelinated axons, establishing important tools for advancing modeling of autonomic nerves.

Composite Backward Differentiation Formula for the Bidomain Equations.

The bidomain equations have been widely used to model the electrical activity of cardiac tissue. While it is well-known that implicit methods have much better stability than explicit methods, implicit methods usually require the solution of a very large nonlinear system of equations at each timestep which is computationally prohibitive. In this work, we present two fully implicit time integration methods for the bidomain equations: the backward Euler method and a second-order one-step two-stage composite backward differentiation formula (CBDF2) which is an L-stable time integration method. Using the backward Euler method as fundamental building blocks, the CBDF2 scheme is easily implementable. After solving the nonlinear system resulting from application of the above two fully implicit schemes by a nonlinear elimination method, the obtained nonlinear global system has a much smaller size, whose Jacobian is symmetric and possibly positive definite. Thus, the residual equation of the approximate Newton approach for the global system can be efficiently solved by standard optimal solvers. As an alternative, we point out that the above two implicit methods combined with operator splittings can also efficiently solve the bidomain equations. Numerical results show that the CBDF2 scheme is an efficient time integration method while achieving high stability and accuracy.

Modeling the effects of sinusoidal stimulation and synaptic plasticity on linked neural oscillators.

The brain exhibits intrinsic oscillatory behavior, which plays a vital role in communication and information processing. Abnormalities in brain rhythms have been linked to numerous disorders, including depression and schizophrenia. Rhythmic electrical stimulation (e.g., transcranial magnetic stimulation and transcranial alternating current stimulation) has been used to modulate these oscillations and produce lasting changes in neural activity. In this computational study, we investigate the combined effects of sinusoidal stimulation and synaptic plasticity on model networks comprised of simple, tunable four-neuron oscillators. While not intended to model a specific brain circuit, this idealization was created to provide some intuition on how electrical modulation can induce plastic changes in the oscillatory state. Linked pairs of oscillators were stimulated with sinusoidal current, and their behavior was measured as a function of their intrinsic frequencies, inter-oscillator synaptic strengths, and stimulus strength and frequency. Under certain stimulus conditions, sinusoidal current can disrupt the network's natural firing patterns. Synaptic plasticity can induce weight imbalances that permanently change the characteristic firing behavior of the network. Grids of 100 oscillators with random frequencies were also subjected to a wide array of stimulus conditions. The characteristics of the post-stimulus network activity depend heavily on the stimulus frequency and amplitude as well as the initial strength of inter-oscillator connections. Synchronization arises at the network level from complex patterns of activity propagation, which are enhanced or disrupted by different stimuli. The findings may prove important to the design of novel neuromodulation treatments and techniques seeking to affect oscillatory activity in the brain.

Mechanism and magnitude of bipolar electrogram directional sensitivity: Characterizing underlying determinants of bipolar amplitude.

BACKGROUND: The amplitude of bipolar electrograms (EGMs) is directionally sensitive, decreasing when measured from electrode pairs oriented oblique to a propagating wavefront. OBJECTIVE: The purpose of this study was to use a computational model and clinical data to establish the mechanism and magnitude of directional sensitivity. METHODS: Simulated EGMs were created using a computational model with electrode pairs rotated relative to a passing wavefront. A clinical database of 18,740 EGMs with varying electrode separation and orientations was recorded from the left atrium of 10 patients with atrial fibrillation during pacing. For each EGM, the angle of incidence between the electrodes and the wavefront was measured using local conduction velocity (CV) mapping. RESULTS: A theoretical model was derived describing the effect of the changing angle of incidence, electrode spacing, and CV on the local activation time difference between a pair of electrodes. Model predictions were validated using simulated and clinical EGMs. Bipolar amplitude measured by an electrode pair is decreased (directionally sensitive) at angles of incidence resulting in local activation time differences shorter than unipolar downstroke duration. Directional sensitivity increases with closer electrode spacing, faster CV, and longer unipolar EGM duration. For narrowly spaced electrode pairs (<5 mm), it is predicted at all orientations. CONCLUSION: Directional sensitivity occurs because bipolar amplitude is reduced when the component unipolar EGMs overlap, such that neither electrode is "indifferent." At the electrode spacing of clinical catheters, this is predicted to occur regardless of catheter orientation. This suggests that bipolar directional sensitivity can be lessened but not overcome by recently introduced catheters with additional rotated electrode pairs.

Diabetes Mobile Care: Aggregating and Visualizing Data from Multiple Mobile Health Technologies.

As the appeal and use of mobile health (mHealth) technologies continues to grow, where does mHealth fit into clinical practice? This article explores the approach and obstacles encountered when integrating mHealth data into existing clinical frameworks and explores data visualization design tradeoffs. Specifically, this paper discusses the successes and challenges that arose when using commercial mHealth technologies, synthesizing multiple mHealth device data, and tailoring visualizations based on iterative feedback from type II diabetes mellitus patients. This research aims to influence the development of patient portals within electronic health records by understanding and addressing the challenges involved in acquiring, interpreting, and displaying this data set. In particular, we need to ensure that the presentation of these data is accessible and understandable by diverse populations.

Rhythmic modulation of thalamic oscillations depends on intrinsic cellular dynamics.

OBJECTIVE: Rhythmic brain stimulation has emerged as a powerful tool to modulate cognition and to target pathological oscillations related to neurological and psychiatric disorders. However, we lack a systematic understanding of how periodic stimulation interacts with endogenous neural activity as a function of the brain state and target. APPROACH: To address this critical issue, we applied periodic stimulation to a unified biophysical thalamic network model that generates multiple distinct oscillations, and examined thoroughly the impact of rhythmic stimulation on different oscillatory states. MAIN RESULTS: We found that rhythmic perturbation induces four basic response mechanisms: entrainment, acceleration, resonance and suppression. Importantly, the appearance and expression of these mechanisms depend highly on the intrinsic cellular dynamics in each state. Specifically, the low-threshold bursting of thalamocortical cells (TCs) in delta (δ) oscillation renders the network relatively insensitive to entrainment; the high-threshold bursting of TCs in alpha (α) oscillation leads to widespread oscillation suppression while the tonic spiking of TC cells in gamma (γ) oscillation results in prominent entrainment and resonance. In addition, we observed entrainment discontinuity during α oscillation that is mediated by firing pattern switching of high-threshold bursting TC cells. Furthermore, we demonstrate that direct excitatory stimulation of the lateral geniculate nucleus (LGN) entrains thalamic oscillations via an asymmetric Arnold tongue that favors higher frequency entrainment and resonance, while stimulation of the inhibitory circuit, the reticular nucleus, induces much weaker and more symmetric entrainment and resonance. These results support the notion that rhythmic stimulation engages brain oscillations in a state- and target-dependent manner. SIGNIFICANCE: Overall, our study provides, for the first time, insights into how the biophysics of thalamic oscillations guide the emergence of complex, state-dependent mechanisms of target engagement, which can be leveraged for the future rational design of novel therapeutic stimulation modalities.

Muscle Thickness and Curvature Influence Atrial Conduction Velocities.

Electroanatomical mapping is currently used to provide clinicians with information about the electrophysiological state of the heart and to guide interventions like ablation. These maps can be used to identify ectopic triggers of an arrhythmia such as atrial fibrillation (AF) or changes in the conduction velocity (CV) that have been associated with poor cell to cell coupling or fibrosis. Unfortunately, many factors are known to affect CV, including membrane excitability, pacing rate, wavefront curvature, and bath loading, making interpretation challenging. In this work, we show how endocardial conduction velocities are also affected by the geometrical factors of muscle thickness and wall curvature. Using an idealized three-dimensional strand, we show that transverse conductivities and boundary conditions can slow down or speed up signal propagation, depending on the curvature of the muscle tissue. In fact, a planar wavefront that is parallel to a straight line normal to the mid-surface does not remain normal to the mid-surface in a curved domain. We further demonstrate that the conclusions drawn from the idealized test case can be used to explain spatial changes in conduction velocities in a patient-specific reconstruction of the left atrial posterior wall. The simulations suggest that the widespread assumption of treating atrial muscle as a two-dimensional manifold for electrophysiological simulations will not accurately represent the endocardial conduction velocities in regions of the heart thicker than 0.5 mm with significant wall curvature.

Microheterogeneity-induced conduction slowing and wavefront collisions govern macroscopic conduction behavior: A computational and experimental study.

The incidence of cardiac arrhythmias is known to be associated with tissue heterogeneities including fibrosis. However, the impact of microscopic structural heterogeneities on conduction in excitable tissues remains poorly understood. In this study, we investigated how acellular microheterogeneities affect macroscopic conduction under conditions of normal and reduced excitability by utilizing a novel platform of paired in vitro and in silico studies to examine the mechanisms of conduction. Regular patterns of nonconductive micro-obstacles were created in confluent monolayers of the previously described engineered-excitable Ex293 cell line. Increasing the relative ratio of obstacle size to intra-obstacle strand width resulted in significant conduction slowing up to 23.6% and a significant increase in wavefront curvature anisotropy, a measure of spatial variation in wavefront shape. Changes in bulk electrical conductivity and in path tortuosity were insufficient to explain these observed macroscopic changes. Rather, microscale behaviors including local conduction slowing due to microscale branching, and conduction acceleration due to wavefront merging were shown to contribute to macroscopic phenomena. Conditions of reduced excitability led to further conduction slowing and a reversal of wavefront curvature anisotropy due to spatially non-uniform effects on microscopic slowing and acceleration. This unique experimental and computation platform provided critical mechanistic insights in the impact of microscopic heterogeneities on macroscopic conduction, pertinent to settings of fibrotic heart disease.

A Scalable Weight-Free Learning Algorithm for Regulatory Control of Cell Activity in Spiking Neuronal Networks.

Recent developments in neural stimulation and recording technologies are providing scientists with the ability of recording and controlling the activity of individual neurons in vitro or in vivo, with very high spatial and temporal resolution. Tools such as optogenetics, for example, are having a significant impact in the neuroscience field by delivering optical firing control with the precision and spatiotemporal resolution required for investigating information processing and plasticity in biological brains. While a number of training algorithms have been developed to date for spiking neural network (SNN) models of biological neuronal circuits, exiting methods rely on learning rules that adjust the synaptic strengths (or weights) directly, in order to obtain the desired network-level (or functional-level) performance. As such, they are not applicable to modifying plasticity in biological neuronal circuits, in which synaptic strengths only change as a result of pre- and post-synaptic neuron firings or biological mechanisms beyond our control. This paper presents a weight-free training algorithm that relies solely on adjusting the spatiotemporal delivery of neuron firings in order to optimize the network performance. The proposed weight-free algorithm does not require any knowledge of the SNN model or its plasticity mechanisms. As a result, this training approach is potentially realizable in vitro or in vivo via neural stimulation and recording technologies, such as optogenetics and multielectrode arrays, and could be utilized to control plasticity at multiple scales of biological neuronal circuits. The approach is demonstrated by training SNNs with hundreds of units to control a virtual insect navigating in an unknown environment.

Unified thalamic model generates multiple distinct oscillations with state-dependent entrainment by stimulation.

The thalamus plays a critical role in the genesis of thalamocortical oscillations, yet the underlying mechanisms remain elusive. To understand whether the isolated thalamus can generate multiple distinct oscillations, we developed a biophysical thalamic model to test the hypothesis that generation of and transition between distinct thalamic oscillations can be explained as a function of neuromodulation by acetylcholine (ACh) and norepinephrine (NE) and afferent synaptic excitation. Indeed, the model exhibited four distinct thalamic rhythms (delta, sleep spindle, alpha and gamma oscillations) that span the physiological states corresponding to different arousal levels from deep sleep to focused attention. Our simulation results indicate that generation of these distinct thalamic oscillations is a result of both intrinsic oscillatory cellular properties and specific network connectivity patterns. We then systematically varied the ACh/NE and input levels to generate a complete map of the different oscillatory states and their transitions. Lastly, we applied periodic stimulation to the thalamic network and found that entrainment of thalamic oscillations is highly state-dependent. Our results support the hypothesis that ACh/NE modulation and afferent excitation define thalamic oscillatory states and their response to brain stimulation. Our model proposes a broader and more central role of the thalamus in the genesis of multiple distinct thalamo-cortical rhythms than previously assumed.

Modeling dynamics in diseased cardiac tissue: Impact of model choice.

Cardiac arrhythmias have been traditionally simulated using continuous models that assume tissue homogeneity and use a relatively large spatial discretization. However, it is believed that the tissue fibrosis and collagen deposition, which occur on a micron-level, are critical factors in arrhythmogenesis in diseased tissues. Consequently, it remains unclear how well continuous models, which use averaged electrical properties, are able to accurately capture complex conduction behaviors such as re-entry in fibrotic tissues. The objective of this study was to compare re-entrant behavior in discrete microstructural models of fibrosis and in two types of equivalent continuous models, a homogenous continuous model and a hybrid continuous model with distinct heterogeneities. In the discrete model, increasing levels of tissue fibrosis lead to a substantial increase in the re-entrant cycle length which is inadequately reflected in the homogenous continuous models. These cycle length increases appear to be primarily due to increases in the tip path length and to altered restitution behavior, and suggest that it is critical to consider the discrete effects of fibrosis on conduction when studying arrhythmogenesis in fibrotic myocardium. Hybrid models are able to accurately capture some aspects of re-entry and, if carefully tuned, may provide a framework for simulating conduction in diseased tissues with both accuracy and efficiency.