Soft robotic design strategy applied to mechanically assisting inspiration.
As depicted in the schematic in Fig. 1a, when the diaphragm contracts, the arclength of the diaphragm shortens, and the entire sheet of the diaphragm moves downwards, acting as a pump. The thoracic cavity volume increases and pressure decreases, ultimately driving respiration.
Our strategy aims to harness the contractile function of pneumatic artificial muscles (PAMs) to mimic and augment the native contraction of the diaphragm. We opt for a McKibben type PAM—a classical soft actuator type with a simple fabrication process and high force generation13 that is capable of mimicking and augmenting biological systems.7–9 The McKibben actuators used in this work were capable of generating up to 50 N of contractile force under 20 psi pressurization (see Supplementary Information). At their simplest, McKibben actuators are composed of an expandable weaved mesh surrounding a bladder connected to an airline (Fig. 1b) (see Methods). When the bladder is pressurized, the mesh expands radially and drives linear contraction (Fig. 1c). Conceptually, we harness the linear contraction of these PAMs by placing them superior to the native diaphragm so that the relaxed PAM conforms to the native curvature of the diaphragm (Fig. 1d). Mimicking the native diaphragm, we anchor the ends of the PAMs to the ribs (see Methods). With pressurization, the length of the PAM shortens, the arclength shortens, and the PAM mechanically pushes the diaphragm downwards (shown in situ in Fig. S1).
In contrast to the dielectric artificial diaphragm described by Bashkin, et al.,10 our diaphragm assist system uses a set of two linear PAMs, leaves the native diaphragm intact, and has a low profile presence (deflated: 5 mL volume, inflated: 17 mL volume). To test this concept in a live porcine model, we surgically implanted a pair of McKibben actuators in an anterior-to-posterior direction lateral to the heart. The actuator placement is visualized in a 3D rendering in Fig. 1e. Fluoroscopy of the diaphragm was taken throughout the experiments. The lateral cross-sectional view from the fluoroscopy shows the realization of our soft robotic strategy in an in vivo pig model (Fig. 1f,g).
Augmenting tidal volume and peak inspiratory flow in vivo.
To evaluate the ability of our diaphragm assist system to augment respiratory function, the animals were instrumented to collect physiological data, including respiratory flows, volumes, and pressures within the respiratory system (Fig. S2). The pressurization of the soft robotic actuators was controlled via a custom-built control system; the actuation pressure data was input into the same high-resolution data acquisition system as the physiological data (see Methods).
Ventilation is key to driving CO2 exchange, so we first examine the flow and volume waveforms as metrics of ventilatory function. Flow is measured by a spirometer. Peak inspiratory flow can be used as a clinical metric of inspiratory function12, which yields a direct measurement of the effect of the diaphragm assist system. Integrating the flow with respect to time yields a volume waveform over time. The volume of each breath (tidal volume) and its rate (minute ventilation) are the most relevant parameters of directly measuring ventilation. Pressures within the respiratory system, such as pleural and abdominal pressures, reveal information about the respiratory biomechanics that physically drive ventilation and are discussed later in this work.
To start each study, the animal was anesthetized appropriately with isoflurane and placed on mechanical ventilation. Isoflurane induces a respiratory depression with decreased tidal volumes and increased respiratory rate that ultimately combine to a reduced minute ventilation14. The respiratory depression secondary to the isoflurane is used as our baseline animal model of respiratory insufficiency due to hypoventilation. Each subject has a reduced but non-zero respiratory drive and response to CO2. Mechanical ventilation is used to support the animal throughout the implantation surgery. Within each subject, we introduce a series of respiratory challenges, collecting data during periods of unassisted ventilation (in which any spontaneous respiration is due to the native respiratory drive) and during periods of actuator assisted ventilation. Mechanical ventilation is used to restore and maintain a state of normoventilation after and between respiratory challenges. To investigate the effect of the diaphragm assist system, a representative respiratory challenge was chosen per subject. The phrenic nerve is intact for all data shown in Fig. 2.
In a vignette from the best responding subject (Fig. 2a), we show that the assist system has the direct capacity to augment the peak inspiratory flow from 0.18 L/s to 0.59 L/s and the tidal volume from 55 mL to 161 mL. When the assist is resumed after a short period of unassisted respiration, the augmentation effect of the actuation on the flow and volume waveforms is reestablished nearly immediately over the course of 2 breaths.
An example of a full respiratory challenge is shown in Fig. 2b. During the unassisted ventilation at the start of the challenge, the subject models a state of hypoventilation. During this period, the tidal volumes and flows have a slight increase over time, indicating the baseline respiratory drive is responding to the increasing CO2 status due to the unassisted low minute ventilation (0.9 L/min). When assist is switched on (as indicated by the actuator pressure waveform and the black arrow), there is a clear jump in the peak inspiratory flow (+0.20 L/s, 95%CI +0.19 L/s to +0.22 L/s), tidal volumes (63mL, 95% CI 58 mL to 68 mL), and minute ventilation (0.9 L/min to 3.1 L/min). The actuators cycle between a pressurized and unpressurized state for 10 minutes. At the end of the respiratory challenge when the respiratory effort has reached a steady state, the assist is switched off and we see that the respiratory effort drops slightly (peak inspiratory flow: -0.09 L/s, 95% CI -0.08 to -0.10; tidal volume: -10 mL, 95% CI -7 to -13) but much less than the jump seen at the start of the respiratory challenge.
The respiratory drive is a slow but dynamic factor underlying all of the respiratory physiology data. As seen in the first 200 s of Fig. 2b, the respiratory drive visibly increases as the low minute ventilation leads to CO2 buildup. This response to CO2 is dynamic and varies between subjects based on each animal’s response to isoflurane. By examining the breaths immediately before and after these transition points (off-to-on and on-to-off), we can examine the direct effect of the diaphragm assist system in terms of augmenting volume and peak inspiratory flow while minimizing the influence of the changing baseline.
This analysis was conducted for one representative respiratory challenge per subject. We see a spectrum of responsiveness to the diaphragm assist system across 5 subjects (Fig. 2c,d,e). The subjects are ordered from largest change in tidal volume at the start of the challenge to the smallest (best responder to worst responder according to Fig. 2d). We find the diaphragm assist system generates much larger respiratory augmentations at the beginning of a trial—when mechanical ventilation support has just been removed, minute ventilation drops suddenly, and the animal’s CO2 state rises rapidly—than at the end of the respiratory challenge when the respiratory baseline is relatively more stabilized (Fig. 2c,d,e).
Subject A was much more responsive to the assist system than any other subject. In terms of tidal volume, 4 of the 5 subjects show an augmentation of >30 mL per breath at the beginning, whereas only 1 of the subjects shows substantial augmentation to the tidal volume at the end. Of the 4 less responsive subjects (B,C,D,E), 3 of them show a mild response at the end while in the worst responder (E), the actuation overall decreased the ventilation metrics (Fig. 2c-e). The subject with the weakest response had the highest baseline weight-normalized minute ventilation at the beginning of the trial (Fig. 2e) compared to other subjects.
Body weight normalized minute ventilation is used to compare these results to normal physiology. Minute ventilation is a metric of the ventilation rate, taking into account both tidal volume and the respiratory rate. In a normal, conscious pig, the expected body weight normalized minute ventilation is 198 mL/min/kg ± 41 mL/min/kg with a range of 104 mL/min/kg to 262 mL/min/kg15. Actuator assisted ventilation allowed all 5 subjects to reach the lower range of normal physiology, and 2 of the subjects even achieved a minute ventilation corresponding to one standard deviation below the normal mean (Fig. 2e). However, we note that this minute ventilation is achieved with low tidal volumes and high respiratory rates, which results in a lower alveolar ventilation than the same minute ventilation achieved with high tidal volumes and low respiratory rates.
Synchronizing with the underlying respiratory effort.
Like with standard mechanical ventilation16,17, patient-ventilator synchrony in our system is critical to the ability to augment respiration. Asynchronous ventilation can destructively interfere with the underlying respiratory effort, leading to worse ventilation with assistance than without.
In order to synchronize the actuation of our assist system with the subject’s underlying respiratory effort, we built a control system (Fig. 3a,b) that can actuate based on the respiratory flow rate. The system uses the spirometry flow sensor as the source data. The flow data is read into our data acquisition system. The associated data analysis software allows a user-set threshold voltage; this threshold voltage was manually titrated during every respiratory trial to achieve qualitatively good synchronization. When the flow rate passes this set threshold, a digital pulse is triggered and sent to the microcontroller in our control box. The microcontroller triggers a pre-set actuation pressure waveform of one cycle of pressurization and depressurization in the electropneumatic regulator, filling and emptying the PAMs with pressurized air (further details in Methods).
Our control system can implement both a set, rhythmic control scheme independent of the native respiratory effort or a dynamic control scheme synchronized with the underlying respiratory effort. Due to the phase and frequency mismatch between the independent actuation and the underlying respiratory effort, the mixed interference of the actuator and the underlying respiratory effort can be seen in both the flow and volume waveform (Fig. 3c). Contrastingly, the well synchronized actuation reveals much more homogenous flow and volume waveforms. (Fig. 3d).
Within each subject, we compare the tidal volumes and peak inspiratory flows in one representative challenge of independent actuation with one representative challenge of synchronized actuation (details in methods). We find that synchronized actuation consistently produces much less variance in the tidal volumes (Fig. 3e,f). Although in some subjects—such as subject A—independent actuation achieved a few higher maximum tidal volumes, the independent actuation also achieved lower minimum tidal volumes across all subjects due to the misalignment of actuations with the underlying respiratory effort leading to destructive interference or due to actuation with no underlying breath—representing a breath that is solely actuator driven.
Factors in optimizing synchronization.
As seen by the mixed interference in Fig. 3c, the alignment of the actuation with the underlying respiratory effort will critically determine the constructive versus destructive nature of the interference. In respiratory challenges that had an independent actuation scheme or a poorly synchronized actuation scheme, we found the datasets that provide a natural variation in the timing of the actuation in relationship to the underlying respiratory effort.
Because mechanical respiratory failure exists as a continuous spectrum of loss of function, we looked at the implications of synchronization in different levels of baseline respiratory effort. As seen in Fig. 2, there is variance in the underlying respiratory function between subjects. To simulate a controlled change in the underlying respiratory function within the same subject, we severed the phrenic nerve in some subjects, simulating diaphragm paralysis in combination with the respiratory depression due to the isoflurane (see Methods). Fig. 4 depicts the analysis of aligning the actuator synchronization to the underlying respiratory effort for two respiratory challenges within subject B: (1) the subject with preserved diaphragm function (Fig. 4, left) and (2) the subject with a severed phrenic nerve (Fig. 4, right).
To optimize for maximum inspiratory augmentation, we investigate the relationship of the timing of different waveform features to the resulting tidal volume and peak inspiratory flow of each breath. The high frequency sampling of our data acquisition system (1000 Hz) allows for millisecond temporal resolution. Custom software was written to analyze the actuation pressure, flow, and volume data.
We identify the breath bounds as determined by the local minima in the volume waveform (the locations of V0), and then finds the time distance between identified waveform features for each individual breath (further details in Methods). Waveform features analyzed include the start of an actuation waveform (P0), peak inspiratory flow (Fpk), the start of inspiration (V0), the start of expiration (Vpk), and others (Fig. 4a,b, and Fig. S3).
The distances between features act as different metrics of alignment and elucidate what factors are important to consider in optimizing synchronization. There are many different features and feature distances that can be analyzed. Fig. 4c-f shows the time relationship of the start of expiration to the actuation pressure (Vpk-P0), but other metrics are shown in Fig. S3.
We examine the influence of these time metrics on tidal volume and peak inspiratory flow. We find the most important predictor variables are time metrics related to the start of expiration (Vpk). With diaphragm function preserved, there is a weak linear relationship between Vpk-P0 and the peak inspiratory flow (R2 = 0.31, p<0.001) (Fig. 4c), and no significant relationship to the tidal volume (R2=0.04, p=0.001) (Fig. 4e). However, when the diaphragm function is removed by severing the phrenic nerve, a clear linear relationship emerges between Vpk-P0 and tidal volume (R2 = 0.84, p<0.001) (Fig. 4f) and a weaker relationship with peak inspiratory flow (R2 = 0.30, p<0.001) (Fig. 4d).
Notably, we do not find these relationships when using the timing between the start of actuation and the start of inspiration (P0 - V0) as a metric. There is no linear relationship between P0-V0 and the peak inspiratory flow or tidal volume for both the cases with and without diaphragm function (Fig. S4)
Comparing respiratory biomechanics.
To compare the respiratory biomechanics of different modes of respiration and ventilation, pleural pressure (Ppl), abdominal pressure (Pab), and transdiaphragmatic pressure (Pdi; Pdi = Pab-Ppl) waveforms are analyzed. Transdiaphragmatic pressure is a metric of diaphragm function6,18,19. Pleural pressure and abdominal pressure are approximated by a sensor mounted on a balloon catheter placed in the esophagus and stomach, respectively. As these sensors approximate Ppl and Pab, the measurements are interpreted as relative measurements and not absolute measurements (see Methods for information about instrumentation and normalization). When analyzing relative pressure waveforms, the most informative metric is the maximum change in pressure per each breath.
In Fig. 5a-c, we show that across subjects (subject C was not instrumented for pressure measurements, and is therefore not shown), actuator assisted ventilation more closely matches the respiratory biomechanics of spontaneous respiration than in mechanical ventilation. Mechanical ventilation pushes air into the lungs, increasing pleural pressure with inspiration, whereas both actuator assisted ventilation and spontaneous respiration generate a negative pleural pressure to drive airflow. As the diaphragm is passive in mechanical ventilation, we see a negligible change in the abdominal pressure, whereas the caudal movement of the diaphragm in both actuator assisted ventilation and spontaneous respiration increases abdominal pressure.
In the representative waveforms from subject A (Fig. 5d-f), the case of highest responsiveness as seen in Fig. 3b-f, the actuator assisted ventilation not only more closely resemble that of spontaneous respiration, but also augments all of the pressure waveforms. Actuator assisted ventilation generates more negative changes in pleural pressure, greater increases in abdominal pressure, and ultimately greater increases in transdiaphragmatic pressure per breath.
A graphical technique used to measure work of breathing (WOB) is the Campbell diagram, referencing pleural pressure with lung volume. Using the pressure and volume data from subject A, we generate the pressure-volume (PV) loops of a Campbell diagram (Fig. 5g). Work of breathing is calculated from this PV loop as the internal area between the inspiratory edge of the loop and the passive chest wall compliance derived from the mechanical ventilation PV data. Normal WOB is 0.35-0.7 J/L.12,20,21 During attenuated spontaneous breathing, the subject’s WOB is 0.10 J/L. During actuator assisted ventilation, the assist system shares the WOB and increases the total average WOB to 0.17 J/L, a 66% increase.