A Dynamic Model of Rescuer Parameters for Optimizing Blood Gas Delivery during Cardiopulmonary Resuscitation.
Cardiopulmonary resuscitation (CPR) involves delivering chest compressions and positive pressure ventilation to cardiac arrest victims to maintain circulatory blood flow and oxygen delivery . Studies have shown that the quality of CPR can impact post-CPR outcomes and survival [2-7]. CPR optimization remains an active topic of resuscitation research because survival rates for post-CPR hospital inpatients remain low . Various approaches have been proposed to modify and optimize the CPR process. Simultaneous ventilation and compression as well as interposed abdominal compression are two examples of CPR modification [9-11].
Another approach to optimizing CPR involves defining and establishing reliable, appropriate CPR performance measures and then modifying CPR parameters to enhance these measures. Brain ischemia is a primary contributor to postarrest morbidity; thus, oxygen delivery and carbon dioxide elimination (end-tidal carbon dioxide [ETCO2]) have been established as important factors for measuring CPR performance [12, 13]. Other performance measures include mean coronary perfusion pressure (CPP), nitric oxide (NO), and balance of systemic and pulmonary perfusion with ventilation [14-17]. The CPR compression to ventilation ratio has been shown to be a CPR parameter that influences postarrest outcomes [18-21]. Another important CPR parameter is ventilation time, which is defined as the pause for ventilation between compressions. Experimental models have been used to assess the effects of ventilation time on coronary perfusion pressure .
In this study, we aimed to define a reliable CPR performance measure, optimize the quality of CPR based on the defined performance measure, and design a dynamically optimized scheme that varies CPR parameters to optimize CPR performance.
2. Materials and Methods
The Institutional Review Board at the Children's Hospital of Philadelphia approved this study.
2.1. Model-Based CPR Optimization. The combination of systemic oxygen delivery and carbon dioxide elimination at the lungs was selected as an objective function for maximization. Delivery of systemic carbon dioxide to the lung for elimination is associated with partial pressure of end-tidal carbon dioxide PETCO2 and consequently has been correlated with CPR outcome . Maximizing the defined objective function will result in improving CPR performance by
(i) ensuring sufficient oxygen delivery to vital organs and (ii) preventing carbon dioxide from accumulating in the body.
2.2. Calculating Blood Gas Delivery. During CPR, the blood gas (oxygen and carbon dioxide) delivery to any organ in the body can be estimated using the following simplified equation :
[D.sub.BG] = [bar.Q][DELTA][C.sub.BG], (1)
where [bar.Q] is the mean blood flow during CPR, [DELTA]C is the change in the blood gas concentration, and subscript BG refers to blood gas. Mean blood flow or [bar.Q] during CPR is
[bar.Q] = [Q.sub.max] x/(T/t)+x, (2)
where T is the average artificial ventilation time, t is the time for one full compression, (1/t) is the compression speed, and x is the compression to ventilation ratio . For systemic oxygen delivery, [Q.sub.max] is the maximum systemic blood flow Qs, and for carbon dioxide delivery to the lungs, [Q.sub.max] is the maximum pulmonary blood flow [Q.sub.p]. On average, we assume that systemic blood flow is equal to pulmonary blood flow; hence, [bar.Q] is equal for both the cases of systemic oxygen delivery and carbon dioxide delivery to the lungs.
To calculate [DELTA]C[O.sub.2] and [DELTA][C.sub.CO2], we need to start with pulmonary gas exchange equations presented in . The presented model assumes a steady-state condition as well as continuous ventilation and perfusion. The fraction of alveolar gas ([f.sub.A]) is equal to end-tidal gas,
[mathematical expression not reproducible], (3)
In the steady-state condition, the oxygen and carbon dioxide balance in the lungs can be expressed using the following equations:
[mathematical expression not reproducible], (4)
where subscript in refers to tracheal inflow, subscript out refers to tracheal outflow, and [mathematical expression not reproducible] are slopes of oxygen and carbon dioxide dissociation curves, respectively. In CPR-related parameters,
[mathematical expression not reproducible], (5)
where R is the average ventilation rate, [v.sub.t] is the tidal volume, and [v.sub.d] is the dead-space volume. Solving for alveolar fractions, we get
[mathematical expression not reproducible]. (6)
Since changes in concentrations are small, they are in the linear range of dissociation curves, or in mathematical terms:
[mathematical expression not reproducible]. (7)
Finally, combining the above equations will result in mathematical expressions for systemic oxygen delivery and carbon dioxide delivery:
[mathematical expression not reproducible], (8)
where [Q.sub.s,max] and [Q.sub.p,max] are maximum systemic and pulmonary blood flow, respectively. In steady-state CPR condition, [Q.sub.s,max] and [Q.sub.p,max] are equal to ensure balanced blood flow to the body and the lungs. The objective function we sought to maximize was the total blood gas delivery, defined as the summation of carbon dioxide delivery to the lung and systemic oxygen delivery. It should be noted that a differentially weighted sum of the two deliveries might also be a possible criterion.
[mathematical expression not reproducible]. (9)
2.3. Optimization. In the next step, CPR parameters were divided into three categories: rescuer dependent, patient dependent, and constant parameters. Rescuer-dependent parameters are parameters that depend on rescuer performance, such as [mathematical expression not reproducible]. Patient-dependent parameters are parameters that vary from patient to patient, such as [Q.sub.max]. Constant parameters are fixed parameters, such as [mathematical expression not reproducible].
In the optimization procedure, we varied the rescuer-dependent parameters for each patient (i.e., with constant patient-dependent parameters) over the entire range of patients. The described parameters and their values are listed in Table 1.
Next, two different optimization schemes were developed: a global optimization scheme and a sequential optimization scheme. Both of the methods used simulated annealing as the optimization technique. Simulated annealing is a computational intelligence technique that aims to solve combinatorial optimization problems to minimize the defined cost function which is characterized by a large number of solutions . The simulated annealing method has been widely applied to problems in the fields of engineering, science, and biomedical research [26-29].
Unlike most of the optimization techniques which could be considered analogous to rapid cooling methods, simulated annealing is analogous to a process of slowly cooling a physical system by providing ample time in order to obtain states with globally minimum energy . As a result, simulated annealing is able to find solutions near a global minimum for very large optimization problems. To define the simulated annealing, we first defined the energy of the system and determined how it related to the temperature of the system. The Helmholtz free energy of a physical system represented by function f at a given temperature T is defined by the following equation:
[E.sub.[empty set]] = [summation over ([empty set])][E.sub.[empty set]]+T[summation over ([empty set])][P.sub.[empty set]]log[P.sub.[empty set]], (10)
where [p.sub.[empty set]] is the probability of occurrence of state [empty set], which could be calculated using the following equation:
[mathematical expression not reproducible], (11)
[PHI] is a set of all possible states of the system. This equation will result in the conclusion that low energy ordered states are strongly favored at low temperatures.
2.3.1. Global Optimization. First, we considered global optimization. In this method, for each patient group, i.e., a fixed Qmax, we sought to maximize the total blood gas delivery that is expressed by Equation (1). We used rescuer-dependent parameters as a set of free parameters of the model for a range of fixed [Q.sub.max] to find the global maximum for total delivery. Instead of a set of fixed values for all of the patients as recommended by American Heart Association CPR guidelines, the proposed global optimization method tried to find the set of parameters for each patient group that resulted in maximum blood gas delivery.
2.3.2. Sequential Optimization. This method involved applying a sequential optimization scheme to actively vary the CPR parameters in order to maximize the defined objective function and thereby improve CPR performance. In this scheme, instead of finding the global maximum of the objective function for each patient group, at each CPR sequence, the rescuer-dependent parameters were varied based on a predefined protocol to find the global maximum of blood gas delivery. We defined a CPR sequence as one cycle of compressions and positive pressure ventilation. During each CPR sequence, a compression to ventilation ratio that maximized the total delivery was found, and then the CPR was performed. Total delivery was then maximized again by varying other rescuer-dependent parameters. These new parameters were used as the starting point of the next CPR sequence, and this procedure was continued until reaching the maximum total delivery for each group.
The steps of the sequential optimization algorithm can hence be summarized as follows:
(1) Start with nominal values.
(2) Optimize the total delivery based on varying x.
(3) Rescuer performs the CPR cycle.
(4) Optimize the remaining rescuer parameters to maximize the oxygen delivery with a fixed x.
(5) Go to Step 2 and repeat until objective function shown in Equation (9) reaches the maximum.
To simulate a wide range of patient groups, we varied Qmax in the range of 700 to 1100 mL blood per minute.
3. Results and Discussion
Results obtained from two different optimization prospective are presented in Figure 1. Results show (1) that the proposed sequential algorithm is more effective in optimizing CPR performance, and (2) increasing [Q.sub.max] will increase total blood gas delivery.
Figure 1 shows that, in all the patient groups, the sequential optimization scheme will result in significant enhancement in the effectiveness of the CPR procedure. The advantage of the sequential optimization technique stemmed from its dynamically adapting effect. Compression to ventilation ratios for three different [Q.sub.max] are presented in Figure 2.
Figure 2 shows the compression to ventilation ratio for three different [Q.sub.max,] (i.e., three different patient groups) during CPR. During initiation of CPR, the rescuer should place more emphasis on chest compressions than on ventilation (high compression to ventilation ratio) and gradually decrease the compression to ventilation ratio and carry out the positive pressure ventilation more frequently.
The results obtained for average ventilation time T for three different [Q.sub.max] are presented in Figure 3. The plot shows that the ventilation time does not change over the time and remains constant for all patient groups.
The results obtained for compression speed 1/t for three different [Q.sub.max] are presented in Figure 4. This figure also shows that there should be interactive changes in the rescuer's performance with progression of the CPR. The plot shows that the compression speed should generally increase as time progresses.
In this study, we present an analysis of dynamic modeling of CPR parameters during CPR to improve CPR performance. While previous research efforts to optimize CPR have focused on finding the best compression to ventilation ratio and keeping the ratio constant during the CPR procedure, we propose a sequential optimization scheme to vary the rescuer-dependent parameters as needed during CPR sequences to optimize CPR performance.
Our results suggest that as CPR progresses, the compression to ventilation ratio should decrease over time to optimize performance. One possible explanation for this observation in our model is that when CPR typically commences, blood oxygen content is sufficient for satisfactory systemic delivery and the carbon dioxide level has only started to rise above normal levels. As CPR progresses, oxygen decreases in the blood and carbon dioxide starts to accumulate, so ventilation becomes increasingly important. Hence, our model suggests that the optimal compression to ventilation ratio should begin higher than is currently practiced (greater than 30 compressions per ventilation) and then decrease over time.
In this study, we followed to a large extent the methods as described in a conference paper . In that earlier study, we used systemic oxygen delivery as a performance measure for optimization . A model developed earlier by Babbs et al.  was used to calculate the oxygen delivery through CPR. The free parameters of the previous model that depended on rescuer performance were ventilation time, compression speed, tidal volume, and fraction of oxygen in the inspired air. Two different optimization problems were carried out. Results obtained from that study showed that the potential of the sequential optimization procedure to enhance the performance of the CPR, but the results of that study only recommended one drop in compression to ventilation ratio during dynamic phase. Promising results from that study encouraged us to explore the subject of dynamic variation of CPR parameters with a broader, more realistic objective function: total blood gas delivery.
There were several limitations to this study. First, we used one published systemic oxygen delivery model to develop analytical expressions for both systemic oxygen delivery and carbon dioxide delivery to the lungs , and we developed and analyzed only two approaches for optimization. Our findings are not necessarily reproducible using different physiological models and analytical approaches. Second, we selected parameter ranges for the optimization process based on expert knowledge and our experiences; others may find different parameter ranges more relevant to their practice population. In the next step, we plan to apply the sequential optimization method on a more complex model of chest that we developed to add systemic blood flow to the objective function for more comprehensive CPR optimization . Lastly, our findings are based on mathematical modeling and have not been validated in vivo, and they should not be used to adjust current CPR practices until such validation has been completed.
In summary, our study illustrates the potential benefit of considering dynamic changes in rescuer-dependent parameters during CPR in order to improve performance. Validation in an animal model and an assessment of generalizability to real data sets are important future steps.
The data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
This research was supported by a grant from the National Institutes of Health (no. 1 R01 NS 72338 01A1).
 J. A. Cooper, J. D. Cooper, and J. M. Cooper, "Cardiopulmonary resuscitation: history, current practice, and future direction," Circulation, vol. 114, no. 25, pp. 2839-2849, 2006.
 R. M. Sutton, M. R. Maltese, D. Niles et al., "Quantitative analysis of chest compression interruptions during in-hospital resuscitation of older children and adolescents," Resuscitation, vol. 80, no. 11, pp. 1259-1263, 2009.
 R. M. Sutton, D. Niles, J. Nysaether et al., "Quantitative analysis of CPR quality during in-hospital resuscitation of older children and adolescents," Pediatrics, vol. 124, no. 2, pp. 494-499, 2009.
 L. J. Morrison, R. W. Neumar, J. L. Zimmerman et al., "Strategies for improving survival after in-hospital cardiac arrest in the United States: 2013 consensus recommendations: a consensus statement from the American Heart Association," Circulation, vol. 127, no. 14, pp. 1538-1563, 2013.
 A. C. Glatz, A. Nishisaki, D. E. Niles et al., "Sternal wall pressure comparable to leaning during CPR impacts intrathoracic pressure and haemodynamics in anaesthetized children during cardiac catheterization," Resuscitation, vol. 84, no. 12, pp. 1674-1679, 2013.
 O. Badaki-Makun, F. Nadel, A. Donoghue et al., "Chest compression quality over time in pediatric resuscitations," Pediatrics, vol. 131, no. 3, pp. e797-e804, 2013.
 R. M. Sutton, S. H. Friess, U. Bhalala et al., "Hemodynamic directed CPR improves short-term survival from asphyxia-associated cardiac arrest," Resuscitation, vol. 84, no. 5, pp. 696-701, 2013.
 E. Jung, C. F. Babbs, S. Lenhart, and V. A. Protopopescu, "Optimal strategy for cardiopulmonary resuscitation with continuous chest compression," Academic Emergency Medicine, vol. 13, no. 7, pp. 715-721, 2006.
 N. Chandra, M. L. Weisfeldt, J. Tsitlik et al., "Augmentation of carotid flow during cardiopulmonary resuscitation by ventilation at high airway pressure simultaneous with chest compression," American Journal of Cardiology, vol. 48, no. 6, pp. 1053-1063, 1981.
 R. C. Koehler, N. Chandra, A. D. Guerci et al., "Augmentation of cerebral perfusion by simultaneous chest compression and lung inflation with abdominal binding after cardiac arrest in dogs," Circulation, vol. 67, no. 2, pp. 266-275, 1983.
 C. F. Babbs, "Interposed abdominal compression CPR: a comprehensive evidence based review," Resuscitation, vol. 59, no. 1, pp. 71-82, 2003.
 C. Sandroni, F. Cavallaro, C. W. Callaway et al., "Predictors of poor neurological outcome in adult comatose survivors of cardiac arrest: a systematic review and meta-analysis. part 2: patients treated with therapeutic hypothermia," Resuscitation, vol. 84, no. 10, pp. 1324-1338, 2013.
 C. Sandroni, F. Cavallaro, C. W. Callaway et al., "Predictors of poor neurological outcome in adult comatose survivors of cardiac arrest: a systematic review and meta-analysis. part 1: patients not treated with therapeutic hypothermia," Resuscitation, vol. 84, no. 10, pp. 1310-1323, 2013.
 K. B. Kern, G. A. Ewy, W. D. Voorhees, C. F. Babbs, and W. A. Tacker, "Myocardial perfusion pressure: a predictor of 24-hour survival during prolonged cardiac arrest in dogs," Resuscitation, vol. 16, no. 4, pp. 241-250, 1988.
 N. M. Tsoukias and S. C. George, "A two-compartment model of pulmonary nitric oxide exchange dynamics," Journal of Applied Physiology, vol. 85, no. 2, pp. 653-666, 1998.
 C. F. Babbs and V. Nadkarni, "Optimizing chest compression to rescue ventilation ratios during one-rescuer CPR by professionals and lay persons: children are not just little adults," Resuscitation, vol. 61, no. 2, pp. 173-181, 2004.
 D. S. Karbing, S. Kjrgaard, S. Andreassen, K. Espersen, and S. E. Rees, "Minimal model quantification of pulmonary gas exchange in intensive care patients," Medical Engineering & Physics, vol. 33, no. 2, pp. 240-248, 2011.
 E. Dorph, L. Wik, T. A. Stremme, M. Eriksen, and P. A. Steen, "Oxygen delivery and return of spontaneous circulation with ventilation:compression ratio 2:30 versus chest compressions only CPR in pigs," Resuscitation, vol. 60, no. 3, pp. 309-318, 2004.
 I. Turner and S. Turner, "Optimum cardiopulmonary resuscitation for basic and advanced life support: a simulation study," Resuscitation, vol. 62, no. 2, pp. 209-217, 2004.
 I. U. Haque, J. P. Udassi, S. Udassi, D. W. Theriaque, J. J. Shuster, and A. L. Zaritsky, "Chest compression quality and rescuer fatigue with increased compression to ventilation ratio during single rescuer pediatric CPR," Resuscitation, vol. 79, no. 1, pp. 82-89, 2008.
 M. Botran, J. Lopez-Herce, J. Urbano, M. Solana, A. Garcia, and A. Carrillo, "Chest compressions versus ventilation plus chest compressions: a randomized trial in a pediatric asphyxial cardiac arrest animal model," Intensive Care Medicine, vol. 37, no. 11, pp. 1873-1880, 2011.
 R. A. Berg, A. B. Sanders, K. B. Kern et al., "Adverse hemodynamic effects of interrupting chest compressions for rescue breathing during cardiopulmonary resuscitation for ventricular fibrillation cardiac arrest," Circulation, vol. 104, no. 20, pp. 2465-2470, 2001.
 K. Lah, M. Krizmaric, and S. Grmec, "The dynamic pattern of end-tidal carbon dioxide during cardiopulmonary resuscitation: difference between asphyxial cardiac arrest and ventricular fibrillation/pulseless ventricular tachycardia cardiac arrest," Critical Care, vol. 15, no. 1, p. R13, 2011.
 A. Jalali, R. A. Berg, V. Nadkarni, and C. Nataraj, "Model based optimization of the cardiopulmonary resuscitation (CPR) procedure," in Proceedings of IEEE Engineering in Medicine and Biology, vol. 2012, pp. 715-718, San Diego, CA, USA, August 2012.
 S. O. Haykin, Neural Networks and Learning Machines, Prentice Hall, Upper Saddle River, NJ, USA, 3rd edition, 2008.
 K. Bryan, P. Cunningham, and N. Bolshakova, "Application of simulated annealing to the biclustering of gene expression data," IEEE Transactions on Information Technology in Biomedicine, vol. 10, no. 3, pp. 519-525, 2006.
 N. Mishra, S. Petrovic, and S. Sundar, "A self-adaptive case-based reasoning system for dose planning in prostate cancer radiotherapy," Medical Physics, vol. 38, no. 12, pp. 6528-6538, 2011.
 T. D. C. Martins, E. D. L. B. de Camargo, R. G. Lima, M. B. P. Amato, and M. D. S. G. Tsuzuki, "Image reconstruction using interval simulated annealing in electrical impedance tomography," IEEE Transactions on Biomedical Engineering, vol. 59, no. 7, pp. 1861-1870, 2012.
 M. Hasegawa, "Evaluation of the physical annealing strategy for simulated annealing: a function-based analysis in the landscape paradigm," Physical Review E: Statistical, Nonlinear, and Soft Matter Physics, vol. 85, no. 2, article 056704, 2012.
 S. P. Brooks and B. J. T. Morgan, "Optimization using simulated annealing," Journal of the Royal Statistical Society. Series D (The Statistician), vol. 44, no. 2, pp. 241-257, 1995.
 A. Jalali, R. A. Berg, V. Nadkarni, and C. Nataraj, "Improving cardiopulmonary resuscitation (CPR) by dynamic variation of CPR parameters," in Proceedings from American Society of Mechanical Engineers System Identification and Therapeutic Control in Bio-Systems Conference, Boston, MA, USA, April 2013.
 C. F. Babbs and K. B. Kern, "Optimum compression to ventilation ratios in CPR under realistic, practical conditions: a physiological and mathematical analysis," Resuscitation, vol. 54, no. 2, pp. 147-157, 2002.
 A. Jalali, A. F. Simpao, V. M. Nadkarni, R. A. Berg, and C. Nataraj, "A novel nonlinear mathematical model of thoracic wall mechanics during cardiopulmonary resuscitation based on a porcine model of cardiac arrest," Journal of Medical Systems, vol. 41, no. 2, p. 20, 2017.
Ali Jalali (iD),(1) Allan F. Simpao (iD),(2) Jorge A. Galvez, (2) Robert A. Berg, (2) Vinay M. Nadkarni, (2) and Chandrasekhar Nataraj (iD)(3)
(1) Health Informatics Core, Johns Hopkins All Children's Hospital, 501 6th Avenue South, St. Petersburg, FL 33701, USA
(2) Department of Anesthesiology & Critical Care Medicine, Children's Hospital of Philadelphia, Perelman School of Medicine at the University of Pennsylvania, 3401 Civic Center Blvd., Philadelphia, PA 19104, USA
(3) Villanova Center for Analytics of Dynamic Systems, Department of Mechanical Engineering, Villanova University, Villanova, PA 19085, USA
Correspondence should be addressed to Ali Jalali; firstname.lastname@example.org
Received 29 August 2018; Accepted 11 November 2018; Published 29 November 2018
Academic Editor: Thomas Desaive
Caption: Figure 1: Results obtained from the two proposed schemes for optimizing total blood gas delivery during cardiopulmonary resuscitation. For systemic oxygen delivery, [Q.sub.max] is the maximum systemic blood flow [Q.sub.s], and for carbon dioxide delivery to the lungs, [Q.sub.max] is the maximum pulmonary blood flow [Q.sub.p].
Caption: Figure 2: Plot of compression to ventilation ratio for three different [Q.sub.max], (i.e., three different patient groups) during cardiopulmonary resuscitation (CPR).
Caption: Figure 3: Plot of average ventilation time for three different [Q.sub.max] values during cardiopulmonary resuscitation (CPR).
Caption: Figure 4: Plot of compression speed over time for three different [Q.sub.max] values during cardiopulmonary resuscitation (CPR).
Table 1: Parameter values and their respective ranges. Parameter Range CPR parameters T [0.132,0.134] [v.sub.t] [600,1000] t [0.008,0.012] [mathematical expression not reproducible] [0.13,0.19] [mathematical expression not reproducible] [0.0002,0.0006] Patient specific parameter [Q.sub.max] [400,1000] Constant parameters [v.sub.d] 150 [mathematical expression not reproducible] 0.8 [mathematical expression not reproducible] 1.5 Parameter Unit CPR parameters T Min [v.sub.t] mL t Min [mathematical expression not reproducible] mLx[O.sub.2]/mL gas [mathematical expression not reproducible] mLx[O.sub.2]/mL gas Patient specific parameter [Q.sub.max] mL/min Constant parameters [v.sub.d] mL [mathematical expression not reproducible] mL/mL [mathematical expression not reproducible] mL/mL The range for each parameter shows the lower and upper limit of variable changes during the optimization. Except for the constant parameters, parameter values change during the CPR delivery.
|Printer friendly Cite/link Email Feedback|
|Title Annotation:||Research Article|
|Author:||Jalali, Ali; Simpao, Allan F.; Galvez, Jorge A.; Berg, Robert A.; Nadkarni, Vinay M.; Nataraj, Chand|
|Publication:||Computational and Mathematical Methods in Medicine|
|Date:||Jan 1, 2018|
|Previous Article:||Does the Temporal Asymmetry of Short-Term Heart Rate Variability Change during Regular Walking? A Pilot Study of Healthy Young Subjects.|
|Next Article:||True Random Number Generation from Bioelectrical and Physical Signals.|