Influence of occlusal thickness on the fracture resistance of chairside milled lithium disilicate posterior full-coverage single-unit prostheses containing virgilite: A comparative in vitro study.

PURPOSE: To evaluate the fracture resistance of chairside computer-aided design and computer-aided manufacturing (CAD-CAM) lithium disilicate mandibular posterior crowns with virgilite of different occlusal thicknesses and compare them to traditional lithium disilicate crowns. MATERIALS AND METHODS: Seventy-five chairside CAD-CAM crowns were fabricated for mandibular right first molars, 60 from novel lithium disilicate with virgilite (CEREC Tessera, Dentsply Sirona), and 15 from traditional lithium disilicate (e.max CAD, Ivoclar Vivadent). These crowns were distributed across five groups based on occlusal thickness and material: Group 1 featured CEREC Tessera crowns with 0.8 mm thickness, Group 2 had 1.0 mm thickness, Group 3 had 1.2 mm thickness, Group 4 with 1.5 mm thickness, and Group 5 included e.max CAD crowns with 1.0 mm thickness. These crowns were luted onto 3D-printed resin dies using Multilink Automix resin cement (Ivoclar Vivadent). Subsequently, they underwent cyclic loading (2,000,000 cycles at 1 Hz with a 275 N force) and loading until fracture. Scanning electron microscopy (SEM) assessed the fractured specimens. Statistical analysis involved one-way ANOVA and the Kruskal-Wallis Test (α = 0.05). RESULTS: Fracture resistance varied significantly (<0.001) across mandibular molar crowns fabricated from chairside CAD-CAM lithium disilicate containing virgilite, particularly between crowns with 0.8 mm and those with 1.2 and 1.5 mm occlusal thickness. However, no significant differences were found when comparing crowns with 1, 1.2, and 1.5 mm thicknesses. CEREC Tessera crowns with 1.5 mm thickness exhibited the highest resistance (2119 N/mm2), followed by those with 1.2 mm (1982 N/mm2), 1.0 mm (1763 N/mm2), and 0.8 mm (1144 N/mm2) thickness, whereas e.max CAD crowns with 1.0 mm occlusal thickness displayed the lowest resistance (814 N/mm2). CONCLUSIONS: The relationship between thickness and fracture resistance in the virgilite lithium disilicate full-coverage crowns was directly proportional, indicating that increased thickness corresponded to higher fracture resistance. No significant differences were noted among crowns with thicknesses ranging from 1 to 1.5 mm. This novel ceramic exhibited superior fracture resistance compared to traditional lithium disilicate.

Sensitivity Analysis of a Left Ventricle Model in the Context of Intraventricular Dyssynchrony.

The objective of the current study was to propose a sensitivity analysis of a 3D left ventricle model in order to assess the influence of parameters on myocardial mechanical dispersion. A finite element model of LV electro-mechanical activity was proposed and a screening method was used to evaluate the sensitivity of model parameters on the standard deviation of time to peak strain. Results highlight the importance of propagation parameters associated with septal and lateral segments activation. Simulated curves were compared to myocardial strains, obtained from echocardiography of one healthy subject and one patient diagnosed with intraventricular dyssynchrony and coronary artery disease. Results show a close match between simulation and clinical strains and illustrate the model ability to reproduce myocardial strains in the context of intraventricular dyssynchrony.

Recursive Model Identification for the Evaluation of Baroreflex Sensitivity.

A method for the recursive identification of physiological models of the cardiovascular baroreflex is proposed and applied to the time-varying analysis of vagal and sympathetic activities. The proposed method was evaluated with data from five newborn lambs, which were acquired during injection of vasodilator and vasoconstrictors and the results show a close match between experimental and simulated signals. The model-based estimation of vagal and sympathetic contributions were consistent with physiological knowledge and the obtained estimators of vagal and sympathetic activities were compared to traditional markers associated with baroreflex sensitivity. High correlations were observed between traditional markers and model-based indices.

Parametric analysis of an integrated cardio-respiratory model in preterm newborns during apnea.

The analysis of the complex interactions involved in the acute physiological response to apnea-bradycardia events in preterm newborns remains a challenging task. This paper presents a novel integrated model of cardio-respiratory interactions, adapted to preterm newborns. A sensitivity analysis, based Morris' screening method, was applied to study the effects of physiological parameters on heart rate and desaturation, during the simulation of a 15-seconds apnea-bradycardia episode. The most sensitive parameters are associated with fundamental, integrative physiological mechanisms involving: (i) respiratory mechanics (intermediate airways and lung compliance), (ii) fraction of inspired oxygen, (iii) metabolic rates (oxygen consumption rate), (iv) heart rate regulation and (v) chemoreflex (gain). Results highlight the relevant influence of physiological variables, involved in preterm apnea-bradycardia events.

Mathematical modeling of respiratory system mechanics in the newborn lamb.

In this paper, a mathematical model of the respiratory mechanics is used to reproduce experimental signal waveforms acquired from three newborn lambs. As the main challenge is to determine specific lamb parameters, a sensitivity analysis has been realized to find the most influent parameters, which are identified using an evolutionary algorithm. Results show a close match between experimental and simulated pressure and flow waveforms obtained during spontaneous ventilation and pleural pressure variations acquired during the application of positive pressure, since root mean square errors equal to 0.0119, 0.0052 and 0.0094. The identified parameters were discussed in light of previous knowledge of respiratory mechanics in the newborn.

Fracture Resistance of Chairside CAD/CAM Molar Crowns Fabricated with Different Lithium Disilicate Ceramic Materials

PURPOSE: To compare the fracture resistance of five different groups of chairside CAD/CAM molar crowns fabricated from various lithium disilicate ceramic materials (LDC): one conventional precrystallized CAD/CAM LDC, two novel precrystallized LDCs, and one fully crystallized LDC tested both with and without optional sintering. MATERIALS AND METHODS: A total of 60 chairside CAD/CAM lithium disilicate molar crowns (n = 12 per group) with 1.5-mm occlusal thickness and a 1.0-mm chamfer finish were designed and fabricated with a chairside CAD/CAM system (CEREC, Dentsply Sirona). The restorations were divided into five groups: (1) IPS e.max CAD; (2) Amber Mill; (3) Straumann n!ce; (4) Straumann n!ce with optional sintering; and (5) Supreme CAD. Restorations were cemented using conventional resin luting cement and primer system to 3D-printed resin dies. Bonded restorations were loaded for 100,000 cycles with 275-N force, and the load at break (LB) and peak load (PL) until fracture were measured. SEM images of fracture surfaces on the printed dies were obtained. RESULTS: Fracture resistance was significantly different depending on the material. Supreme CAD showed the highest fracture resistance (LB: 1,557.2 N; PL: 1,785.8 N), followed by Amber Mill (LB: 1,393.0 N; PL: 1,604.2 N) and IPS e.max CAD (LB: 1,315.7 N; PL: 1,461.9 N). Straumann n!ce without (LB: 862.4 N; PL: 942.9 N) and with the optional sintering (LB: 490.4 N; PL: 541.0 N) showed significantly lower fracture resistance than the others. CONCLUSION: The fracture resistance of chairside CAD/CAM lithium disilicate molar crowns varied depending on the material, and the novel materials did not perform as well as the conventional equivalents. Fully crystallized lithium disilicate ceramic block materials showed lower fracture resistance than precrystallized counterparts and should be used with caution in the clinic, especially with optional sintering.

Fracture Resistance of Chairside CAD/CAM Molar Crowns Fabricated with Different Lithium Disilicate Ceramic Materials.

PURPOSE: To compare the fracture resistance of five different groups of chairside CAD/CAM molar crowns fabricated from various lithium disilicate ceramic materials (LDC): one conventional precrystallized CAD/CAM LDC, two novel precrystallized LDCs, and one fully crystallized LDC tested both with and without optional sintering. MATERIALS AND METHODS: A total of 60 chairside CAD/CAM lithium disilicate molar crowns (n = 12 per group) with 1.5-mm occlusal thickness and a 1.0-mm chamfer finish were designed and fabricated with a chairside CAD/CAM system (CEREC, Dentsply Sirona). The restorations were divided into five groups: (1) IPS e.max CAD; (2) Amber Mill; (3) Straumann n!ce; (4) Straumann n!ce with optional sintering; and (5) Supreme CAD. Restorations were cemented using conventional resin luting cement and primer system to 3D-printed resin dies. Bonded restorations were loaded for 100,000 cycles with 275-N force, and the load at break (LB) and peak load (PL) until fracture were measured. SEM images of fracture surfaces on the printed dies were obtained. RESULTS: Fracture resistance was significantly different depending on the material. Supreme CAD showed the highest fracture resistance (LB: 1,557.2 N; PL: 1,785.8 N), followed by Amber Mill (LB: 1,393.0 N; PL: 1,604.2 N) and IPS e.max CAD (LB: 1,315.7 N; PL: 1,461.9 N). Straumann n!ce without (LB: 862.4 N; PL: 942.9 N) and with the optional sintering (LB: 490.4 N; PL: 541.0 N) showed significantly lower fracture resistance than the others. CONCLUSION: The fracture resistance of chairside CAD/CAM lithium disilicate molar crowns varied depending on the material, and the novel materials did not perform as well as the conventional equivalents. Fully crystallized lithium disilicate ceramic block materials showed lower fracture resistance than precrystallized counterparts and should be used with caution in the clinic, especially with optional sintering.

Apnea bradycardia detection based on new coupled hidden semi Markov model.

In this paper, a method for apnea bradycardia detection in preterm infants is presented based on coupled hidden semi Markov model (CHSMM). CHSMM is a generalization of hidden Markov models (HMM) used for modeling mutual interactions among different observations of a stochastic process through using finite number of hidden states with corresponding resting time. We introduce a new set of equations for CHSMM to be integrated in a detection algorithm. The detection algorithm was evaluated on a simulated data to detect a specific dynamic and on a clinical dataset of electrocardiogram signals collected from preterm infants for early detection of apnea bradycardia episodes. For simulated data, the proposed algorithm was able to detect the desired dynamic with sensitivity of 96.67% and specificity of 98.98%. Furthermore, the method detected the apnea bradycardia episodes with 94.87% sensitivity and 96.52% specificity with mean time delay of 0.73 s. The results show that the algorithm based on CHSMM is a robust tool for monitoring of preterm infants in detecting apnea bradycardia episodes. Graphical Abstract Apnea Bradycardia detection using Coupled hidden semi Markov Model from electrocardiography. In this model, a sequence of hidden states is assigned to each observation based on the effects of previous states of all observations.

Detection of Apnea Bradycardia from ECG Signals of Preterm Infants Using Layered Hidden Markov Model.

Apnea-bradycardia (AB) is a common complication in prematurely born infants, which is associated with reduced survival and neurodevelopmental outcomes. Thus, early detection or predication of AB episodes is critical for initiating preventive interventions. To develop automatic real-time operating systems for early detection of AB, recent advances in signal processing can be employed. Hidden Markov Models (HMM) are probabilistic models with the ability of learning different dynamics of the real time-series such as clinical recordings. In this study, a hierarchy of HMMs named as layered HMM was presented to detect AB episodes from pre-processed single-channel Electrocardiography (ECG). For training the hierarchical structure, RR interval, and width of QRS complex were extracted from ECG as observations. The recordings of 32 premature infants with median 31.2 (29.7, 31.9) weeks of gestation were used for this study. The performance of the proposed layered HMM was evaluated in detecting AB. The best average accuracy of 97.14 ± 0.31% with detection delay of - 5.05 ± 0.41 s was achieved. The results show that layered structure can improve the performance of the detection system in early detecting of AB episodes. Such system can be incorporated for more robust long-term monitoring of preterm infants.

Early bradycardia detection and therapeutic interventions in preterm infant monitoring.

In very preterm infants, cardio-respiratory events and associated hypoxemia occurring during early postnatal life have been associated with risks of retinopathy, growth alteration and neurodevelopment impairment. These events are commonly detected by continuous cardio-respiratory monitoring in neonatal intensive care units (NICU), through the associated bradycardia. NICU nurse interventions are mainly triggered by these alarms. In this work, we acquired data from 52 preterm infants during NICU monitoring, in order to propose an early bradycardia detector which is based on a decentralized fusion of three detectors. The main objective is to improve automatic detection under real-life conditions without altering performance with respect to that of a monitor commonly used in NICU. We used heart rate lower than 80 bpm during at least 10 sec to define bradycardia. With this definition we observed a high rate of false alarms (64%) in real-life and that 29% of the relevant alarms were not followed by manual interventions. Concerning the proposed detection method, when compared to current monitors, it provided a significant decrease of the detection delay of 2.9 seconds, without alteration of the sensitivity (97.6% vs 95.2%) and false alarm rate (63.7% vs 64.1%). We expect that such an early detection will improve the response of the newborn to the intervention and allow for the development of new automatic therapeutic strategies which could complement manual intervention and decrease the sepsis risk.

Echocardiographic view and feature selection for the estimation of the response to CRT.

Cardiac resynchronization therapy (CRT) is an implant-based therapy applied to patients with a specific heart failure (HF) profile. The identification of patients that may benefit from CRT is a challenging task and the application of current guidelines still induce a non-responder rate of about 30%. Several studies have shown that the assessment of left ventricular (LV) mechanics by speckle tracking echocardiography can provide useful information for CRT patient selection. A comprehensive evaluation of LV mechanics is normally performed using three different echocardioraphic views: 4, 3 or 2-chamber views. The aim of this study is to estimate the relative importance of strain-based features extracted from these three views, for the estimation of CRT response. Several features were extracted from the longitudinal strain curves of 130 patients and different methods of feature selection (out-of-bag random forest, wrapping and filtering) have been applied. Results show that more than 50% of the 20 most important features are calculated from the 4-chamber view. Although features from the 2- and 3-chamber views are less represented in the most important features, some of the former have been identified to provide complementary information. A thorough analysis and interpretation of the most informative features is also provided, as a first step towards the construction of a machine-learning chain for an improved selection of CRT candidates.

Prediction of response to cardiac resynchronization therapy using a multi-feature learning method.

We hypothesized that a multiparametric evaluation, based on the combination of electrocardiographic and echocardiographic parameters, could enhance the appraisal of the likelihood of reverse remodeling and prognosis of favorable clinical evolution to improve the response of cardiac resynchronization therapy (CRT). Three hundred and twenty-three heart failure patients were retrospectively included in this multicenter study. 221 patients (68%) were responders, defined by a decrease in left ventricle end-systolic volume ≥15% at the 6-month follow-up. In addition, strain data coming from echocardiography were analyzed with custom-made signal processing methods. Integrals of regional longitudinal strain signals from the beginning of the cardiac cycle to strain peak and to the instant of aortic valve closure were analyzed. QRS duration, septal flash and different other features manually extracted were also included in the analysis. The random forest (RF) method was applied to analyze the relative feature importance, to select the most significant features and to build an ensemble classifier with the objective of predicting response to CRT. The set of most significant features was composed of Septal Flash, E, E/A, E/EA, QRS, left ventricular end-diastolic volume and eight features extracted from strain curves. A Monte Carlo cross-validation method with 100 runs was applied, using, in each run, different random sets of 80% of patients for training and 20% for testing. Results show a mean area under the curve (AUC) of 0.809 with a standard deviation of 0.05. A multiparametric approach using a combination of echo-based parameters of left ventricular dyssynchrony and QRS duration helped to improve the prediction of the response to cardiac resynchronization therapy.

Symmetry of horizontal and sagittal condylar path angles: an in vivo study.

The specific aims of the study were to determine (1) the mean condylar path inclination (CP) and Bennett angles (BA) in occlusally healthy dentate subjects and (2) whether there was right-left asymmetry of CP and BA. Eccentric movements of three mm and five mm from reference positions were recorded for 45 subjects using a Cadiax system. The mean right CP angles: 48.8 +/- 8 degrees (three mm) and 48.2 +/- 7 degrees (five mm); the mean left CP angles: 48.8 +/- 8 degrees (three mm) and 49.5 +/- 70 (5mm) were steeper than suggested means for setting articulators. The mean right BA: 7.7 +/- 3 degrees (3mm) and 7.9 +/- 4 degrees (five mm); and the mean left BA: 8.1 +/- 4 degrees (three mm) and 8.5 +/- 4 degrees (five mm) were lower than suggested means for setting articulators. Paired-t-tests showed no significant differences between right-left CP and right-left BA. These results suggest that current recommended average settings for semi-adjustable articulators for dentate individuals need to be reassessed.

Model-based estimation of left ventricular pressure and myocardial work in aortic stenosis.

This paper proposes a model-based estimation of left ventricular (LV) pressure for the evaluation of constructive and wasted myocardial work of patients with aortic stenosis (AS). A model of the cardiovascular system is proposed, including descriptions of i) cardiac electrical activity, ii) elastance-based cardiac cavities, iii) systemic and pulmonary circulations and iv) heart valves. After a sensitivity analysis of model parameters, an identification strategy was implemented using a Monte-Carlo cross-validation approach. Parameter identification procedure consists in two steps for the estimation of LV pressures: step 1) from invasive, intraventricular measurements and step 2) from non-invasive data. The proposed approach was validated on data obtained from 12 patients with AS. The total relative errors between estimated and measured pressures were on average 11.9% and 12.27% and mean R2 were equal to 0.96 and 0.91, respectively for steps 1 and 2 of parameter identification strategy. Using LV pressures obtained from non-invasive measurements (step 2) and patient-specific simulations, Global Constructive (GCW), Wasted (GWW) myocardial Work and Global Work Efficiency (GWE) parameters were calculated. Correlations between measures and model-based estimations were 0.88, 0.80, 0.91 respectively for GCW, GWW and GWE. The main contributions concern the proposal of the parameter identification procedure, applied on an integrated cardiovascular model, able to reproduce LV pressure specifically to each AS patient, by non-invasive procedures, as well as a new method for the non-invasive estimation of constructive, wasted myocardial work and work efficiency in AS.

Dynamic changes in ventricular depolarization during exercise in patients with Brugada syndrome.

Brugada syndrome (BS) is a genetic pathological condition associated with a high risk for sudden cardiac death (SCD). Ventricular depolarization disorders have been suggested as a potential electrophysiological mechanism associated with high SCD risk on patients with BS. This paper aims to characterize the dynamic changes of ventricular depolarization observed during physical exercise in symptomatic and asymptomatic BS patients. To this end, cardiac ventricular depolarization features were automatically extracted from 12-lead ECG recordings acquired during standardized exercise stress test in 110 BS patients, of whom 25 were symptomatic. Conventional parameters were evaluated, including QRS duration, R and S wave amplitudes ([Formula: see text], [Formula: see text]), as well as QRS morphological features, such as up-stroke and down-stroke slopes of the R and S waves ([Formula: see text], [Formula: see text] and [Formula: see text]). The effects of physical exercise and recovery on the dynamics of these markers were assessed in both BS populations. Features showing significantly different dynamics between the studied groups were used alone and in combination with the clinical characteristics of the patients in a logistic regression analysis. Results show larger changes in the second half of the QRS complex through [Formula: see text] and [Formula: see text] measured in the right precordial leads for asymptomatic patients, especially during recovery, when the vagal tone is more pronounced. Multivariate analysis involving both types of features resulted in a reduced model of three relevant features ([Formula: see text] in lead V2, Sex and heart rate recovery, HRR), which achieved a suitable discrimination performance between groups; sensitivity = 80% and specificity = 75% (AUC = 83%). However, after controlling the model for possible confounding factors, only one feature ([Formula: see text]) remained meaningful. This adjusted model significantly improved the overall discrimination performance by up to: sensitivity = 84% and specificity = 100% (AUC = 94%). The study highlights the importance of physical exercise test to unmask differentiated behaviors between symptomatic and asymptomatic BS patients through depolarization dynamic analysis. This analysis together with the obtained model may help to identify asymptomatic patients at low or high risk of future cardiac events, but it should be confirmed by further prospective studies.

A tissue-level electromechanical model of the left ventricle: application to the analysis of intraventricular pressure.

The ventricular pressure profile is characteristic of the cardiac contraction progress and is useful to evaluate the cardiac performance. In this contribution, a tissue-level electromechanical model of the left ventricle is proposed, to assist the interpretation of left ventricular pressure waveforms. The left ventricle has been modeled as an ellipsoid composed of twelve mechano-hydraulic sub-systems. The asynchronous contraction of these twelve myocardial segments has been represented in order to reproduce a realistic pressure profiles. To take into account the different energy domains involved, the tissue-level scale and to facilitate the building of a modular model, multiple formalisms have been used: Bond Graph formalism for the mechano-hydraulic aspects and cellular automata for the electrical activation. An experimental protocol has been defined to acquire ventricular pressure signals from three pigs, with different afterload conditions. Evolutionary Algorithms have been used to identify the model parameters in order to minimize the error between experimental and simulated ventricular pressure signals. Simulation results show that the model is able to reproduce experimental ventricular pressure. In addition, electro-mechanical activation times have been determined in the identification process. For example, the maximum electrical activation time is reached, respectively, 96.5, 139.3 and 131.5 ms for the first, second, and third pigs. These preliminary results are encouraging for the application of the model on non-invasive data like ECG, arterial pressure or myocardial strain.

Differences in Brugada Syndrome Patients Through Ventricular Repolarization Analysis During Sleep.

Brugada syndrome (BS) is a genetic pathology that might cause sudden cardiac death (SCD) in patients with a structurally normal heart. Repolarization disorders have been postulated as a potential substrate for triggering cardiac arrhythmia in BS, that usually occur at rest or during sleep. In this paper, we have characterized ventricular repolarization markers during sleep on patients suffering from BS. To this end, standard 12-lead ECG recordings were analyzed in a population of 110 BS patients (25 symptomatic). The QT and the T-wave peak to T-wave end intervals (respectively QT and Tpe) were assessed from lead V5. The linear relationship between these markers and the instantaneous heart rate period (RR interval) are determined during each hour and for the whole sleep period. From the models obtained, corrected QT and Tpe measures were then estimated for each patient at 60 beats/min (QT60 and Tpe60) and at the mean heart rate observed during the involved time interval (QTHR and TpeHR). Results show larger values for symptomatic patients in all markers, with significant differences with respect to the asymptomatic group in the case of Tpe (Tpe60: p = 0.0012; TpeHR: p = 0.0014). Moreover, the temporal profiles of these markers reveal major differences among BS subgroups during the last 3 hours of sleep, where symptomatic patients presented increased QT60/HR (p = 0.01) and Tpe60/HR (p <; 0.001), as compared to the initial sleep hours. We conclude that BS patients present different repolarization properties according to their symptomatology, especially during the final stage of sleep.

Improving methodology in heart rate variability analysis for the premature infants: Impact of the time length.

BACKGROUND: Heart rate variability (HRV) has been emerging in neonatal medicine. It may help for the early diagnosis of pathology and estimation of autonomous maturation. There is a lack of standardization and automation in the selection of the sequences to analyze and some features have not been explored in this specific population. The main objective of this study was to analyze the impact of the time length of the sequences on the estimation of linear and non-linear HRV features, including horizontal visibility graphs (HVG). METHODS: HRV features were repeatedly measured with linear and non-linear methods on 2-, 5-, 10-minute sequences selected from the longest 15-min sequence and recorded on a weekly basis in 39 infants less than 31 weeks at birth. The associations between HRV measurements were analyzed through principal component analysis and k-means clustering. The effects of the time lengths on HRV measurements and post-menstrual age (PMA) were analyzed by linear mixed effect model for repeated measures. RESULTS: The domains of analysis were concordant for their descriptive parameters of short (rMSSD, SD1 and HF) and long-term (SD, SD2 and LF) variability. α1 was correlated with the LF/HF and SD2/SD1. DC and AC were correlated with short-term variability estimates and significantly increased with GA and PMA. Shortening the windows of analysis increased the random measurement error for all the features and increased the bias for all but short term features and HVGs. CONCLUSION: The linear and non-linear measurements of HRV are correlated each other. Shortening the windows of analysis increased the random error for all the features and increased the bias for all but short term features and HVGs. Short-term HRV can be an index for evaluating the maturation in whatever sequence length.

Quantification of Neural Conduction Block on the Rat Sciatic Nerve based on EMG Response.

Neural conduction block performed by balanced-charge kilohertz frequency alternating currents (KHFAC) has been identified as a potential technique for therapy delivery in different clinical setups. The underlying mechanisms that contribute to this phenomenon have been studied through computational models and animal experiments. However, the optimal stimulation parameters to achieve axonal conduction block are difficult to define, since they depend on the species, the nerve being targeted, as well as the technical and experimental setup. This study proposes an experimental setup along with an original data processing approach for the quantification of the effectiveness of neural conduction block. Experiments were performed on the sciatic nerve of two Sprague-Dawley rats, by evaluating different groups of stimulation parameters with varying amplitudes and frequencies, ranging from 1 to 10 mA and from 2 to 10 kHz, respectively. Results suggest that the effectiveness of axonal conduction block strongly depends on the selection of the stimulation parameters. In this work, more effective blockages were achieved for frequencies around 4 kHz and within an approximate amplitude range of 2 to 8 mA.

Recursive model identification for the analysis of the autonomic response to exercise testing in Brugada syndrome.

This paper proposes the integration and analysis of a closed-loop model of the baroreflex and cardiovascular systems, focused on a time-varying estimation of the autonomic modulation of heart rate in Brugada syndrome (BS), during exercise and subsequent recovery. Patient-specific models of 44 BS patients at different levels of risk (symptomatic and asymptomatic) were identified through a recursive evolutionary algorithm. After parameter identification, a close match between experimental and simulated signals (mean error = 0.81%) was observed. The model-based estimation of vagal and sympathetic contributions were consistent with physiological knowledge, enabling to observe the expected autonomic changes induced by exercise testing. In particular, symptomatic patients presented a significantly higher parasympathetic activity during exercise, and an autonomic imbalance was observed in these patients at peak effort and during post-exercise recovery. A higher vagal modulation during exercise, as well as an increasing parasympathetic activity at peak effort and a decreasing vagal contribution during post-exercise recovery could be related with symptoms and, thus, with a worse prognosis in BS. This work proposes the first evaluation of the sympathetic and parasympathetic responses to exercise testing in patients suffering from BS, through the recursive identification of computational models; highlighting important trends of clinical relevance that provide new insights into the underlying autonomic mechanisms regulating the cardiovascular system in BS. The joint analysis of the extracted autonomic parameters and classic electrophysiological markers could improve BS risk stratification.