Publications

Hybrid free-standing biomimetic materials are developed by integrating the VDAC36 β-barrel protein into robust and flexible three-layered polymer nanomembranes. The first and third layers are prepared by spin-coating a mixture of poly(lactic acid) (PLA) and poly(vinyl alcohol) (PVA). PVA nanofeatures are transformed into controlled nanoperforations by solvent-etching. The two nanoperforated PLA layers are separated by an electroactive layer, which is successfully electropolymerized by introducing a conducting sacrificial substrate under the first PLA nanosheet. Finally, the nanomaterial is consolidated by immobilizing the VDAC36 protein, active as an ion channel, into the nanoperforations of the upper layer. The integration of the protein causes a significant reduction of the material resistance, which decreases from 21.9 to 3.9 kΩ cm2. Electrochemical impedance spectroscopy studies using inorganic ions and molecular metabolites (i.e.l-lysine and ATP) not only reveal that the hybrid films behave as electrochemical supercapacitors but also indicate the most appropriate conditions to obtain selective responses against molecular ions as a function of their charge. The combination of polymers and proteins is promising for the development of new devices for engineering, biotechnological and biomedical applications.

JTD Keywords: channels, evolution, filter, Outer-membrane proteins

This paper describes the design of a linear phase low-pass differentiator filter with a finite impulse response (FIR) for extracting transient features of gas sensor signals (the so-called “bouts”). The detection of these bouts is relevant for estimating the distance of a gas source in a turbulent plume. Our current proposal addresses the shortcomings of previous ‘bout’ estimation methods, namely: (i) they were based in non-causal digital filters precluding real time operation, (ii) they used non-linear phase filters leading to waveform distortions and (iii) the smoothing action was achieved by two filters in cascade, precluding an easy tuning of filter performance. The presented method is based on a low-pass FIR differentiator, plus proper post-processing, allowing easy algorithmic implementation for real-time robotic exploration. Linear phase filters preserve signal waveform in the bandpass region for maximum reliability concerning both ‘bout’ detection and amplitude estimation. As a case study, we apply the proposed filter to predict the source distance from recordings obtained with metal oxide (MOX) gas sensors in a wind tunnel. We first perform a joint optimization of the cut-off frequency of the filter and the bout amplitude threshold, for different wind speeds, uncovering interesting relationships between these two parameters. We demonstrate that certain combinations of parameters can reduce the prediction error to 8 cm (in a distance range of 1.45 m) improving previously reported performances in the same dataset by a factor of 2.5. These results are benchmarked against traditional source distance estimators such as the mean, variance and maximum of the response. We also study how the length of the measurement window affects the performance of different signal features, and how to select the filter parameters to make the predictive models more robust to changes in wind speed. Finally, we provide a MATLAB implementation of the bout detection algorithm and all analysis code used in this study.

JTD Keywords: Gas sensors, Differentiator, Low pass filter, Metal oxide semiconductor, MOX sensors, Signal processing, Feature extraction, Gas source localization, Robotics

Neural respiratory drive as measured by the electromyography allows the study of the imbalance between the load on respiratory muscles and its capacity. Surface respiratory electromyography (sEMG) is a non-invasive tool used for indirectly assessment of NRD. It also provides a way to evaluate the level and pattern of respiratory muscle activation. The prevalence of electrocardiographic activity (ECG) in respiratory sEMG signals hinders its proper evaluation. Moreover, the occurrence of abnormal heartbeats or cardiac arrhythmias in respiratory sEMG measures can make even more challenging the NRD estimation. Respiratory sEMG can be evaluated using the fixed sample entropy (fSampEn), a technique which is less affected by cardiac artefacts. The aim of this work was to investigate the performance of the fSampEn, the root mean square (RMS) and the average rectified value (ARV) on respiratory sEMG signals with supraventricular arrhythmias (SVA) for NRD estimation. fSampEn, ARV and RMS parameters increased as the inspiratory load increased during the test. fSampEn was less influenced by ECG with SVAs for the NRD estimation showing a greater response to respiratory sEMG, reflected with a higher percentage increase with increasing load (228 % total increase, compared to 142 % and 135 % for ARV and RMS, respectively).

JTD Keywords: Electrocardiography, Muscles, Electrodes, Estimation, Band-pass filters, Electromyography, Heart beat

To optimize long-term nocturnal non-invasive ventilation in patients with chronic obstructive pulmonary disease, surface diaphragm electromyography (EMGdi) might be helpful to detect patient-ventilator asynchrony. However, visual analysis is labor-intensive and EMGdi is heavily corrupted by electrocardiographic (ECG) activity. Therefore, we developed an automatic method to detect inspiratory onset from EMGdi envelope using fixed sample entropy (fSE) and a dynamic threshold based on kernel density estimation (KDE). Moreover, we combined fSE with adaptive filtering techniques to reduce ECG interference and improve onset detection. The performance of EMGdi envelopes extracted by applying fSE and fSE with adaptive filtering was compared to the root mean square (RMS)-based envelope provided by the EMG acquisition device. Automatic onset detection accuracy, using these three envelopes, was evaluated through the root mean square error (RMSE) between the automatic and mean visual onsets (made by two observers). The fSE-based method provided lower RMSE, which was reduced from 298 ms to 264 ms when combined with adaptive filtering, compared to 301 ms provided by the RMS-based method. The RMSE was negatively correlated with the proposed EMGdi quality indices. Following further validation, fSE with KDE, combined with adaptive filtering when dealing with low quality EMGdi, indicates promise for detecting the neural onset of respiratory drive.

JTD Keywords: Fixed sample entropy, Adaptive filtering, Root mean square, Diaphragm electromyography, Non-invasive mechanical ventilation, Chronic obstructive pulmonary disease

Obstructive sleep apnea (OSA) is a common disorder with a low diagnosis ratio, leaving many patients undiagnosed and untreated. In the last decades, accelerometry has been found to be a feasible solution to obtain respiratory activity and a potential tool to monitor OSA. On the other hand, many smartphone-based systems have already been developed to propose solutions for OSA monitoring and treatment. The objective of this work was to develop an automatic event detector based on smartphone accelerometry and pulse oximetry, and to assess its ability to detect thoracic movements. It was validated with a commercial OSA monitoring system at home. Results of this preliminary pilot study showed that the proposed event detector for accelerometry signals is a feasible tool to detect abnormal respiratory events, such as apneas and hypopneas, and has potential to be included in smartphone-based systems for OSA assessment.

JTD Keywords: Sleep apnea, Detectors, Pulse oximetry, Monitoring, Manuals, Band-pass filters, Pulse oximeter

The electrical activity of the diaphragm measured by surface electromyography (sEMGdi) provides indirect information on neural respiratory drive. Moreover, it allows evaluating the ventilatory pattern from the onset and offset (ntoff) estimation of the neural inspiratory time. sEMGdi amplitude variation was quantified using the fixed sample entropy (fSampEn), a less sensitive method to the interference from cardiac activity. The detection of the ntoff is controversial, since it is located in an intermediate point between the maximum value and the cessation of sEMGdi inspiratory activity, evaluated by the fSampEn. In this work ntoff detection has been analyzed using thresholds between 40% and 100 % of the fSampEn peak. Furthermore, fSampEn was evaluated analyzing the r parameter from 0.05 to 0.6, using a m equal to 1 and a sliding window size equal to 250 ms. The ntoff has been compared to the offset time (toff) obtained from the airflow during a controlled respiratory protocol varying the fractional inspiratory time from 0.54 to 0.18 whilst the respiratory rate was constant at 16 bpm. Results show that the optimal threshold values were between 66.0 % to 77.0 % of the fSampEn peak value. r values between 0.25 to 0.50 were found suitable to be used with the fSampEn.

JTD Keywords: Protocols, Low pass filters, Electrodes, Standards, Band-pass filters, Muscles, Cutoff frequency

EMG signals reflect the neuromuscular activation patterns related to the execution of a certain movement or task. In this work, we focus on reaching and grasping (R&G) movements in rats. Our objective is to develop an automatic algorithm to detect the onsets and offsets of muscle activity and use it to study muscle latencies in R&G maneuvers. We had a dataset of intramuscular EMG signals containing 51 R&G attempts from 2 different animals. Simultaneous video recordings were used for segmentation and comparison. We developed an automatic onset/offset detector based on the ratio of local maxima of Teager-Kaiser Energy (TKE). Then, we applied it to compute muscle latencies and other features related to the muscle activation pattern during R&G cycles. The automatic onsets that we found were consistent with visual inspection and video labels. Despite the variability between attempts and animals, the two rats shared a sequential pattern of muscle activations. Statistical tests confirmed the differences between the latencies of the studied muscles during R&G tasks. This work provides an automatic tool to detect EMG onsets and offsets and conducts a preliminary characterization of muscle activation during R&G movements in rats. This kind of approaches and data processing algorithms can facilitate the studies on upper limb motor control and motor impairment after spinal cord injury or stroke.

JTD Keywords: Muscles, Electromyography, Rats, Low pass filters, Microsoft Windows, Band-pass filters

Measuring diaphragmatic electromyography (EMGdi) provides an indirect quantification of neural respiratory drive and allows the delimitation of diaphragm neural activation and deactivation during inspiration. EMGdi recordings have been incorporated in novel modes of assisted mechanical ventilation, such as neurally adjusted ventilatory assist (NAVA), to trigger and cycle-off the ventilator. The EMGdi signal improves the assistance delivered by more conventional ventilatory modes, in which the ventilator is synchronized with the patient employing a pneumatic triggering. In this work, we evaluate the time delay between the onset and offset of inspiratory activity estimated from EMGdi and three respiratory mechanical signals: the respiratory flow (FL), the transdiaphragmatic pressure (Pdi) and the diaphragm length (Ldi) signals. To this purpose, these signals were acquired in three mongrel dogs surgically instrumented under general anesthesia. Onsets and offsets were estimated manually and by automatic algorithms on these signals. The highest delays were obtained between EMGdi and FL (100 ms) while the lowest delays were obtained between EMGdi and Pdi (8 ms). Moreover, differences between manual and automatic estimations showed a mean absolute error lower than 45 ms. In conclusion, our study points out that both EMGdi and Pdi signals detect the onset and offset of inspiratory activity earlier than the FL signal, and would therefore be better for the improvement of patient-ventilator synchrony.

JTD Keywords: Estimation, Ventilation, Anesthesia, Dogs, Manuals, Power harmonic filters

Anatomical and/or functional reentries have been proposed as one of the main mechanism of perpetuation of cardiac fibrillation processes. However, technical limitations have difficult the characterization of those reentries and are hampering the development of effective anti-arrhythmic treatments. The goal of this study is to present a novel technology to map with high resolution the center of fibrillation drivers in order to characterize the mechanisms of reentry. Cell cultures of human cardiac-like cells differentiated from pluripotent stem cells were analyzed with a novel microscopic optical mapping system. The pharmacological response to verapamil administration of each type of reentry was analyzed. In all analyzed cell cultures, a reentry was identified as the mechanism of maintenance of the arrhythmia. Interestingly, the administration of verapamil produced opposite effects on activation rate depending on the mechanisms of reentry (i.e. anatomical or functional). Microscopic optical mapping of reentries allows the identification of perpetuation mechanisms which has been demonstrated to be linked with different pharmacological response.

JTD Keywords: Stem cells, Rotors, Microscopy, Optical filters, Calcium, Optical microscopy, Biomedical optical imaging

Obstructive Sleep Apnea (OSA) is a sleep disorder that affects mainly the adult and elderly population. Due to the high percentage of patients who remain undiagnosed and untreated because of limitations of current diagnosis methods, the management of OSA is an important social, scientific and economic problem that will be difficult to be assumed by health systems. On the other hand, smartphone platforms (mHealth systems) are being considered as an innovative solution, thanks to the integration of the essential sensors to obtain clinically relevant parameters in the same device or in combination with wireless wearable devices.

JTD Keywords: Sleep apnea, Microphones, Monitoring, Sensors, Accelerometers, Biomedical monitoring, Band-pass filters

Cardiorespiratory Phase Synchronization (CRPS) is a manifestation of coupling between cardiac and respiratory systems complementary to Respiratory Sinus Arrhythmia. In this work, we investigated CRPS during wake and sleep stages in Polysomnographic (PSG) recordings of 30 subjects suspected from Obstructive Sleep Apnea (OSA). The population was classified into three severity groups according to the Apnea Hypopnea Index (AHI): G1 (AHI<;15), G2 (15<;=AHI<;30) and G3 (AHI>30). The synchrogram between single lead ECG and respiratory abdominal band signals from PSG was computed with the Hilbert transform technique. The different phase locking ratios (PLR) m:n were monitored throughout the night. Ratio 4:1 was the most frequent and it became more dominant as OSA severity increased. CRPS was characterized by the percentage of synchronized time (%Sync) and the average duration of synchronized epochs (AvDurSync) using three different thresholds. Globally, we observed that %Sync significantly decreased and AvDurSync slightly increased with OSA severity. A high synchronization threshold enhanced these population differences. %Sync was significantly higher in NREM than in REM sleep in G2 and G3 groups. Population differences observed during sleep did not translate to the initial wake state. Reduced CRPS could be an early marker of OSA severity during sleep, but further studies are needed to determine whether CRPS is also present during wakefulness.

JTD Keywords: Band-pass filters, Electrocardiography, Heart beat, Sleep apnea, Sociology, Statistics, Synchronization

Extracting clinical information from one single measurement represents a step forward in the assessment of the respiratory muscle function. This attracting idea entails the reduction of the instrumentation and fosters to develop new medical integrated technologies. We present the use of the fixed sample entropy (fSampEn) as a more direct method to non-invasively derive the breathing activity from the diaphragm electromyographic (EMGdi) signal, and thus to extract the respiratory rate, an important vital sign which is cumbersome and time-consuming to be measured by clinicians. fSampEn is a method to evaluate the EMGdi activity that is less sensitive to the cardiac activity (ECG) and its application has proven to be useful to evaluate the load of the respiratory muscles. The behavior of the proposed method was tested in signals from two subjects that performed an inspiratory load protocol, which consists of increments in the inspiratory mouth pressure (P_mouth). Two respiratory signals were derived and compared to the P_mouth signal: the ECG-derived respiration (EDR) signal from the lead-I configuration, and the EMG-derived respiration (EMGDR) signal by applying the fSampEn method over the EMGdi signal. The similitude and the lag between signals were calculated through the cross-correlation between each derived respiratory signal and the P_mouth. The EMGDR signal showed higher correlation and lower lag values (≥ 0.91 and ≤ 0.70 s, respectively) than the EDR signal (≥ 0.83 and ≤0.99 s, respectively). Additionally, the respiratory rate was estimated with the P_mouth, EDR and EMGDR signals showing very similar values. The results from this preliminary work suggest that the fSampEn method can be used to derive the respiration waveform from the respiratory muscle electrical activity.

JTD Keywords: Band-pass filters, Electrocardiography, Electromyography, Entropy, Mouth, Muscles, Protocols

Non-invasive evaluation of respiratory activity is an area of increasing research interest, resulting in the appearance of new monitoring techniques, ones of these being based on the analysis of the diaphragm mechanomyographic (MMGdi) signal. The MMGdi signal can be decomposed into two parts: (1) a high frequency activity corresponding to lateral vibration of respiratory muscles, and (2) a low frequency activity related to excursion of the thoracic cage. The purpose of this study was to apply the empirical mode decomposition (EMD) method to obtain the low frequency of MMGdi signal and selecting the intrinsic mode functions related to the respiratory movement. With this intention, MMGdi signals were acquired from a healthy subject, during an incremental load respiratory test, by means of two capacitive accelerometers located at left and right sides of rib cage. Subsequently, both signals were combined to obtain a new signal which contains the contribution of both sides of thoracic cage. Respiratory rate (RR) measured from the mechanical activity (RR_Mmg) was compared with that measured from inspiratory pressure signal (RR_P). Results showed a Pearson's correlation coefficient (r = 0.87) and a good agreement (mean bias = -0.21 with lower and upper limits of -2.33 and 1.89 breaths per minute, respectively) between RR_mmg and RR_P measurements. In conclusion, this study suggests that RR can be estimated using EMD for extracting respiratory movement from low mechanical activity, during an inspiratory test protocol.

JTD Keywords: Accelerometers, Band-pass filters, Biomedical measurement, Empirical mode decomposition, Estimation, IP networks, Muscles

The study of the mechanomyographic signal of the diaphragm muscle (MMGdi) is a promising technique in order to evaluate the respiratory muscles effort. The relationship between amplitude and frequency parameters of this signal with the respiratory effort performed during respiration is of great interest for researchers and physicians due to its diagnostic potentials. However, MMGdi signals are frequently contaminated by a cardiac vibration or mechanocardiographic (MCG) signal. An adaptive noise cancellation (ANC) can be used to reduce the MCG interference in the recorded MMGdi activity. In this paper, it is evaluated the proposed ANC scheme by means of a synthetic MMGdi signal with a controlled MCG interference. The Pearson's correlation coefficient (PCC) between both root mean square (RMS) and mean frequency (fm) of the synthetic MMGdi signal are considerably reduced with the presence of cardiac vibration noise (from 0.95 to 0.87, and from 0.97 to 0.76, respectively). With the ANC algorithm proposed the effect of the MCG noise on the amplitude and frequency of MMG parameters is reduced considerably (PCC of 0.93 and 0.97 for the RMS and fm, respectively). The ANC method proposed in this work is an interesting technique to attenuate the cardiac interference in respiratory MMG signals. Further investigation should be carried out to evaluate the performance of the ANC algorithm in real MMGdi signals.

JTD Keywords: Adaptive filters, Frequency modulation, Interference, Muscles, Noise cancellation, Vibrations, Cardiology, Medical signal processing, Muscle, Signal denoising, ANC algorithm, MCG interference, Pearson correlation coefficient, Adaptive noise cancellation, Cardiac vibration interference, Cardiac vibration noise, Diaphragm muscle, Mechanocardiographic signal, Mechanomyographic signals, Respiratory muscles effort

Sleep in patients with Sleep Apnea-Hypopnea Syndrome (SAHS) is frequently interrupted with arousals. Increased amounts of arousals result in shortening total sleep time and repeated sleep-arousal change can result in sleep fragmentation. According to the American Sleep Disorders Association (ASDA) an arousal is a marker of sleep disruption representing a detrimental and harmful feature for sleep. The nature of arousals and its role on the regulation of the sleep process raises controversy and has sparked the debate in the last years. In this work, we analyzed and compared the EEG spectral content of respiratory and spontaneous arousals on a database of 45 SAHS subjects. A total of 3980 arousals (1996 respiratory and 1984 spontaneous) were analyzed. The results showed no differences between the spectral content of the two kinds of arousals. Our findings raise doubt as to whether these two kinds of arousals are truly triggered by different organic mechanisms. Furthermore, they may also challenge the current beliefs regarding the underestimation of the importance of spontaneous arousals and their contribution to sleep fragmentation in patients suffering from SAHS.

JTD Keywords: Adaptive filters, Correlation, Databases, Electroencephalography, Hospitals, Sleep apnea, Electroencephalography, Medical signal processing, Pneumodynamics, Sleep, EEG spectral content, Organic mechanism, Respiratory, Sleep apnea hypopnea syndrome, Sleep fragmentation, Spectral content, Spontaneous arousal

In the trend towards miniaturization, a detector module containing multiple IR sensor channels is being built and characterized. In its final form it contains thermopiles, narrow band filters and Fresnel lenses. An important feature of such module is the assembly by flip-chip of the IR filters on top of the thermopiles. The performance of the filter-thermopile ensemble has been assessed by physical simulation and experiments and it has been optimized by the use of an empirically validated model. It has been found that integration of filters (or lenses) too close to the IR detector may lead to degraded performance due to thermal coupling. The impact and extent of this degradation has been thoroughly explored, being the main parameter the distance between the IR sensor and the filter. To avoid such detrimental effects a possibility is to set the device in vacuum conditions, obtaining an improved output response and avoiding the influence of the filters. Another way is to increase the solder joint height. Beyond a certain height, the filter is considered to be isolated from the thermopile.

JTD Keywords: Assembly, Infrared sensor, Infrared filter, Fresnel lenses, FEM simulation, Optimization