Abstract
Introduction
The administration of heparin in large doses to prevent blood clotting within the cardiopulmonary bypass is a defining aspect of cardiac anaesthesiology. These doses frequently exceed 500 units per kilogram in total. Following the separation from the cardiopulmonary bypass (CPB), the excessive heparin in the system must be reversed to avoid ongoing bleeding. Traditionally, a protamine-to-total-heparin ratio of 1:1 has been employed. However, the administration of protamine carries itself inherent risks regarding the coagulation status.
While protamine inhibits heparin by causing dissociation of the heparin-antithrombin complex, it also has the effect of reducing thrombin generation, reducing platelet count and function, and enhancing fibrinolysis.1-3 Accordingly, the surplus of protamine relative to the circulating heparin levels may exert a deleterious impact on the coagulation process. The European Association of Cardiothoracic Surgery (EACTS) and European Association of Cardiothoracic Anaesthesia (EACTA) Guidelines on patient blood management for adult cardiac surgery have determined that a ratio of less than 1:1 should be employed in the majority of cases, citing literature for ratios in the range of 0.6–0.8:1.4,5 Nevertheless, the evidence to support other fixed ratios remains inconclusive.6-8
The pharmacokinetics of heparin are complex and highly dependent on the injected dose. Heparin is eliminated from the bloodstream via two principal mechanisms: rapid renal elimination and uptake by the reticuloendothelial system. The latter occurs mainly in the Kupffer cells of the liver and the sinus endothelial cells of the spleen. 9 The kinetics of heparin elimination are non-linear. In a small study comprising four participants in each group, the half-life of a dose of 200 units per kilogram was approximately 96 minutes, while a dose of 400 units resulted in an approximate half-life of 152 minutes. 10 A search, conducted on PubMed and Google Scholar in December 2024 using the keywords ‘heparin half-life', did not identify any studies of higher doses conducted in humans.
Overview of the Compared Pharmacokinetic Models.
Abbreviation: t1/2, Half-Life.
The objective of this study is to compare the protamine doses calculated with the three presented models in a large sample size, comprising a number of simulated cases with different doses of protamine and different durations of extracorporeal circulation.
Methods
Three models of pharmacokinetics of heparin were compared in our study: the one-compartment model described by Miles et al., 12 and the two-compartment models described by Meesters et al. 11 and in the PRODOSE trial. 13
The high-fidelity open-access VitalDB database of real-world patients published on PhysioNet was employed to provide demographic data for a large cohort of surgical patients.14-16 This database comprises anonymised records of patients who underwent surgical procedures at Seoul National University Hospital between August 2016 and June 2017. The acquisition and subsequent publication of data for further research were approved by the Institutional Review Board of the aforementioned hospital. We selected patients aged 18 years or above undergoing all types of surgeries who required postoperative intensive care unit stay from the database and included them in the further analysis. The ideal body weight for each patient was calculated using the Devine formula. 17
The objective of the study was to compare the hypothetical doses calculated with the three models in two different scenarios. In these simulated scenarios, the typical dosages employed in our institution were utilised. The first scenario (Scenario 1) involved the administration of 500 units of heparin per kilogram of total body weight prior to cardiopulmonary bypass, with an additional 5000 units incorporated into the priming solution of the CPB, and no further heparin administered throughout the course of extracorporeal circulation. The second scenario (Scenario 2) represented a more complex case, with a longer duration of CPB. It entailed the same initial dose and priming solution, but with an additional 5000 units administered as a bolus during the CPB.
Means and Standard Deviations Employed to Generate the Data for the Duration of the Intervals Used in the Calculations.
Abbreviations: tIC, time between initial dose of heparin and initiation of the cardiopulmonary bypass; tCR, time between the initiation of the cardiopulmonary bypass and the reversal of heparin; tC2, time between the initiation of the cardiopulmonary bypass and the administration of the second dose of heparin; t2R, time between the administration of the second dose of heparin and the reversal of heparin; SD, standard deviation.
For each simulated case and scenario, the dose of protamine for heparin reversal was calculated in accordance with the formulas presented in each paper. The calculated doses were rounded up to the nearest ten milligrams. Furthermore, the ratios of calculated doses to initial and total heparin doses were calculated.
Exemplary cases of each scenario were presented on concentration-time plots, thereby facilitating a visual comparison of the decay curves. Furthermore, we plotted a modified formula described by Meesters et al., utilising the t1/2 of 150 minutes instead of 250 minutes stated in the original formula. The t1/2 in the PRODOSE formula is dependent on the initial dose of heparin; therefore, we plotted three curves with doses of 300 units/kg, 400 units/kg and 500 units/kg, respectively.
The primary outcome of the study is the total dose of protamine for heparin reversal, calculated using the cited formulas. Secondary outcome of the study are the ratios of protamine to initial heparin and of protamine to total heparin administered during the case.
Statistical Methods
The statistical analysis of the data was conducted utilising the JASP 0.19 software (JASP Team, Amsterdam, The Netherlands). Descriptive statistics are employed to present the data. A paired Student's t-test is employed to compare data that is normally distributed. Data that is not normally distributed is compared with a Wilcoxon signed-rank test. Normally distributed data are presented as mean ± standard deviation, whereas non-normally distributed data are presented as median [interquartile range]. A P-value of less than 0.05 is assumed to be statistically significant.
Results
Participants and Descriptive Data
Demographic Data of the Selected Patients From the VitalDB Database Along With the Generated Data Regarding the Cardiopulmonary Bypass.
Abbreviations: BMI, body mass index; IBW, ideal body weight; tIC, time between the initial dose of heparin and the initiation of the cardiopulmonary bypass; tCR, time between the initiation of the cardiopulmonary bypass and the reversal of heparin; tC2, time between the initiation of the cardiopulmonary bypass and the administration of the second dose of heparin; t2R, time between the administration of the second dose of heparin and the reversal of heparin.
Main Results and Other Analyses
Comprehensive Results of the Statistical Analysis. All Pairwise Comparisons Were Statistically Significant With P < .001.

Comparison of exemplary heparin decay curves in Scenario 1. After initial bolus of heparin (in this case, 35 000 units), cardiopulmonary bypass was commenced after 10 minutes, with 5000 units heparin in the priming solution. Abbreviations: Meesters 250, original formula of Meesters et al. with t1/2250 minutes; Meesters 150, modified formula of Meesters et al. with t1/2 150 minutes; Miles, formula of Miles et al.; PRODOSE 500, 400, and 300, formula from PRODOSE trial with 500, 400, and 300 units of heparin per kg, respectively; 0.67 and 0.6, two-thirds and six-tenths of the total heparin dose, respectively.].

Comparison of exemplary heparin decay curves in Scenario 2. After initial bolus of heparin (in this case, 35 000 units), cardiopulmonary bypass was commenced after 10 minutes, with 5000 units of heparin in the priming solution; 50 minutes thereafter another 5000 units of heparin were administered; Abbreviations: Meesters 250, original formula of Meesters et al. with t1/2 250 minutes; Meesters 150, modified formula of Meesters et al. with t1/2150 minutes; Miles, formula of Miles et al.; PRODOSE 500, 400, and 300, formula from PRODOSE trial with 500, 400, and 300 units of heparin per kg, respectively; 0.67 and 0.6, two-thirds and six-tenths of the total heparin dose, respectively.].
Discussion
The findings of this study of two modelled scenarios of on-pump surgery on a sample of 1166 patients demonstrate that the protamine doses calculated using the PRODOSE formula are consistently the lowest among all the compared models. In Scenario 1, with a mean total duration between the initial heparin bolus and the administration of protamine of 90 minutes (range 19–166 min), the ratio of protamine to total heparin ranged from 0.46:1 to 0.97:1 across all models. In contrast, in Scenario 2, the mean duration was 160 minutes (range 31–277 min), and the ratio ranged from 0.41:1 to 0.83:1. The ratios were observed to be the highest in cases with the shortest duration.
Despite the disparate assumptions underlying the various study models, the resulting calculations yielded comparable outcomes. In Scenario 1, the discrepancy between the lowest (240 mg in the PRODOSE model) and the highest (270 mg in the Miles’ model) mean calculated protamine dose was 12.5%. In Scenario 2, with 220 mg in the PRODOSE model and 250 mg in the Meesters model, the discrepancy was 13.6%. The issue of whether disparities of such magnitude possess any clinical significance remains unresolved and necessitates investigation in real-life studies.
The half-lives for slow compartment calculated in the Miles’ model and PRODOSE trial (155.2 minutes for 400 units per kilogram ideal body weight and 90.6 minutes for 200 units per kilogram ideal body weight) are in close agreement with the half-times reported by Olsson and Stig in their study.10,12,13 Conversely, Meesters et al. postulated that the impact of non-reversed heparin is more deleterious than that of an excess of protamine. To circumvent the potential underdosing of protamine in cases of intraoperative hypothermia, he employed a prolonged half-time of 250 minutes. 11 Nevertheless, it should be noted that the study also involved cardiopulmonary bypass procedures conducted in patients under normothermic or mild hypothermic conditions, with a core temperature over 34°C. It is important to acknowledge that the model has not been validated for application in cases of moderate or deep hypothermia. In fact, the population of the PRODOSE study had a lower mean core temperature of 33.8 ± 1.5°C. As illustrated in Figures 1 and 2, the decay curve generated by substituting the half-life value in Meesters’ model with 150 minutes exhibits a high degree of similarity to the curve calculated using the PRODOSE model, which involved the administration of 400 units of heparin per kilogram of ideal body weight.
The examined models permitted the calculation of the protamine dose on an individualised basis, thereby achieving the ratios recommended by the current guidelines, which are less than 1:1. 4 Nevertheless, it is essential to take into account the characteristics of the study population in the original trials, as well as the exclusion criteria employed by their authors. The study conducted by Meesters and colleagues excluded patients who had undergone emergency surgery, were underweight or severely overweight, had haematological disorders, or required renal replacement therapy or preoperative heparin. 11 The Miles’ model was not evaluated in patients with preoperative laboratory coagulation abnormalities, surgeries with hypothermic circulatory arrest, surgeries on the aorta other than the aortic root, or implantations of ventricular assist devices. 12 In the PRODOSE study, patients were additionally excluded if they were under the age of 18, had a total body weight exceeding 120 kg, had received adenosine diphosphate receptor antagonists within seven days of surgery, had received an unfractionated heparin infusion or therapeutic low molecular weight heparin within 24 hours before surgery, or were undergoing solid organ transplantation. 13
Hypothermia leads to a deceleration in the rate of heparin elimination.18,19 Cohen et al. demonstrated that at 25°C, there is no notable decline in heparin concentration. 18 Consequently, patients who have undergone surgery in hypothermia may require a greater dose of protamine than would be predicted by the models. Further studies are needed to explore the pharmacokinetics of heparin during hypothermia.
It should also be noted that the calculations are based on a heparin initial dose of 500 units per kilogram. Lower doses result in a shorter half-life in the Miles and PRODOSE models, which would further increase the differences between these models and the Meesters’ model.
In a recent randomised control trial comparing a fixed dose of 250 mg protamine with a 1:1 ratio, researchers found no differences in post-protamine ACT or chest tube output in the first 24 hours postoperatively. 20 The fixed dose chosen by the study authors is close to the median doses calculated using the models analysed, although it is still higher than the calculated dose in some patients.
The efficacy of heparin reversal must be monitored through the use of activated clotting time (ACT) measurements or viscoelastic methods, as some patients across all trials demonstrated the necessity for additional protamine following the initial reversal procedure.
Limitations
The present study aimed to compare protamine dosages administered under three proposed heparin pharmacokinetic models in a large model sample of patients in two clinical scenarios. It should be noted, however, that it did not encompass a comprehensive range of clinical scenarios. Nevertheless, the protamine doses calculated with the PRODOSE model will consistently be the lowest of the three models analysed, due to the presence of a fast compartment (in comparison to the Miles’ model) or a shorter half-time (in comparison to the Meesters’ model). Furthermore, the mathematical models offer only an approximation of the heparin concentration at the time of protamine administration. To confirm the efficacy of the calculated doses, it is essential to employ clinical assessment of oozing and recognised methods to assess completeness of heparin reversal.
Future Perspectives
It is hypothesised that the pharmacokinetic models of heparin, when employed to guide protamine administration, may facilitate individualisation of dosage in accordance with the prevailing European guidelines. Further real-life studies are required to substantiate the safety of this approach, particularly with regard to the incidence of cases where the calculated dose of protamine results in underdosing. In the event of adapting one of the models within an institution, a before-after study may be conducted to assess the real-life impact.
Conclusion
In the simulated cases, the doses of protamine calculated with pharmacokinetic models exhibited discrepancies of up to 13.6%. All the models that were studied yielded a ratio of protamine to heparin less than 1:1. The discrepancies between the models were attributable to the employment of single vs two-compartment models, as well as the utilisation of constant vs dose-dependent half-lives of heparin. While confirmation of complete reversal with alternative methods is still necessary, the models facilitate a more individualised dose selection than the fixed 1:1, 0.8:1, 0.6:1 or other ratios proposed in the literature.
Supplemental Material
Supplemental Material - Comparison of Three Mathematical Models of the Pharmacokinetics of Heparin to Guide the Use of Protamine in a Large Simulated Cohort of Patients
Supplemental Material for Comparison of Three Mathematical Models of the Pharmacokinetics of Heparin to Guide the Use of Protamine in a Large Simulated Cohort of Patients by Stanislaw Vander Zwaag, Imre Kukel, Asen Petrov, and Jens Fassl
Footnotes
Acknowledgements
We acknowledge the use of DeepL Write (DeepL SE, Cologne, Germany), set to the “Academic” style, in the proofreading of the manuscript at the final stage of its preparation.
Authors’ contributions
Declaration of Conflicting Interests
The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Funding
The author(s) received no financial support for the research, authorship, and/or publication of this article.
Ethical Statement
Data Availability Statement
All raw data are made available as Supporting materials.
Supplemental Material
Supplemental material for this article is available online.
Appendix
References
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
