There is a dire need for targeted approaches to improve trauma patient treatment outcome. Trauma is the leading cause of death between the ages of 1-44 in the U.S.; those who survive suffer huge morbidity and are left with permanent disabilities. Trauma-induced coagulopathy (TIC) occurs after severe trauma and shock, is biologically characterized by perturbations to the balance between clotting and fibrinolysis, and is clinically characterized by uncontrolled bleeding and either death or clotting complications in those who survive. The initial traumatic hemorrhage accounts for the majority of all trauma-related deaths, and 50% of the mortalities of critically injured patients who undergo surgery. Targeting coagulation biology and the resuscitation strategy in the first 24 hours of care are critical, since 80% of deaths from hemorrhage occur within this window.
Many aspects of the present disclosure can be better understood with reference to the following drawings. The components in the drawings are not necessarily to scale, emphasis instead being placed upon clearly illustrating the principles of the present disclosure. Moreover, in the drawings, like reference numerals designate corresponding parts throughout the several views.
The present disclosure relates to coagulation treatment systems and methods for computing coagulation factor concentrations that rectify clotting in a trauma patient in order to quantitatively guide trauma patient coagulation factor levels while accounting for protein interactions. Exemplary systems and methods utilize an improved thrombin dynamics model to capture the coagulation process to control, use rapidly measurable concentrations to help predict patient state and individual clotting dynamics, and account for patient-specific effects and limitations when adding coagulation factors to remedy coagulopathy by following a novel ordering in which to tune coagulation factors. Validation of an exemplary system/method show superior performance over clinical practice in attaining normal coagulation factor concentrations and normal clotting profiles simultaneously.
Memory elements 120 include one or more physical memory devices such as, for example, a local memory and one or more storage devices. Local memory refers to random access memory (RAM) or other non-persistent memory device(s) generally used during actual execution of the program code. Storage device may be implemented as a hard disk drive (HDD), solid state drive (SSD), or other persistent data storage device. Computing system 100 may also include one or more cache memories (not shown) that provide temporary storage of at least some program code in order to reduce the number of times program code must be retrieved from storage device during execution.
Stored in the memory 120 are both data and several components that are executable by the processor 110. In particular, stored in the memory 120 and executable by the processor 110 are code for a predictive thrombin dynamics model 140, a control algorithm 145 and code for outputting a predictive outcome from the predictive model such as a treatment plan 150. Also stored in the memory 120 may be a data store 125 and other data. The data store 125 can include an electronic repository or database relevant to predictive model results. In addition, an operating system may be stored in the memory 120 and executable by the processor 110. In an embodiment, predictive model data are stored in the data store 125, such as model parameters.
For example, a predictive model may include a digitally constructed model of a probability of individual clotting dynamics. In this context, the model refers to an electronic digitally stored set of executable instructions and data values, associated with one another, which are capable of receiving and responding to a programmatic or other digital call, invocation, or request for resolution based upon specified input values, to yield one or more stored output values that can serve as the basis of computer-implemented recommendations, output data displays, or machine control, among other things. Persons of skill in the field find it convenient to express models using mathematical equations, but that form of expression does not confine the models disclosed herein to abstract concepts; instead, each model herein has a practical application in a computer in the form of stored executable instructions and data that implement the model using the computer. The model may include a model of a current status and/or a model of predicted events of one or more fields (e.g., predicted Calibrated Automated Thrombogram (CAT) trajectories from coagulation factor concentrations). Model and field data may be stored in data structures in memory, rows in a database table, in flat files or spreadsheets, or other forms of stored digital data.
Input/output (I/O) devices 160 such as a keyboard, a display device, and a pointing device may optionally be coupled to computing system 100. The I/O devices may be coupled to computing system 100 either directly or through intervening I/O controllers. A network adapter may also be coupled to computing system to enable computing system to become coupled to other systems, computer systems, remote printers, and/or remote storage devices through intervening private or public networks. Modems, cable modems, Ethernet cards, and wireless transceivers are examples of different types of network adapter that may be used with computing system 100.
Current trauma treatments involve rules-of-thumb and lab-based resuscitation guidelines. In most centers, a preset ratio of blood products is administered to rapidly control hemorrhage. Although some studies attribute improved outcomes to such resuscitation control, other studies show the opposite, including conflicting data for the prehospital transfusion of fresh frozen plasma (FFP) and red blood cells (RBC), and for different ratios of blood products. A possible reason is the dynamic nature of patient coagulation state; too much or too little of beneficial static interventions may result in poor outcomes because of a targeting mismatch with resuscitation needs at that timepoint. While well-intentioned, blood product transfusion is linked to inflammatory morbidities and side-effects including acute respiratory distress syndrome and multi-organ failure. Despite much research and vast improvements in clinical care, severely injured patients that require massive transfusions still have 30% mortality.
Hence, trauma patients may benefit from a tailored transfusion strategy, or from innovative treatments that include coagulation factor (blood protein) concentrates. Targeting individual coagulation proteases via coagulation factor concentrates has benefit for hematologic diseases, such as hemophilia. Although the kitchen-sink approach of using FFP in TIC has its proponents, targeted coagulation factor therapy may have better outcomes compared to FFP-based treatments. However, some reports on modulating coagulation factors (including factor VII, factor IX, and factor X) have shown limited benefit in individually correcting coagulopathic hemorrhaging. Moreover, coagulation factor levels cannot be increased in isolation. For example, elevated levels of activated protein C (aPC) inhibit hemorrhaging, but are also associated with undesirable outcomes including pneumonia, multi-organ failure, and death. Thus, there are open requirements to: (1) confirm the benefits of modulating coagulation factors; and (2) develop a new quantitative modulation approach that incorporates interactions between coagulation factors.
Existing protocols for such goal-driven trauma treatment use thromboelastometry, a viscoelastic coagulation assay, but this assay is time-consuming at typically about an hour per run. These protocols using thromboelastometry are non-quantitative, rely on clinical intuition and older standard procedures, and only correct for a small number of coagulation factors and/or their interactions. Such protocols are qualitative because traditional statistical analysis and machine learning on a static trauma patient measures like coagulation factor concentrations and are not informative in diagnosing coagulation or treatment outcomes. These protocols are also non-dynamic, meaning that they are unable to make time-course, patient-specific predictions of recovery, and they do not facilitate future intervention automation.
Because patient responses to trauma are complex and dynamic, with risks in both hemorrhagic and thrombotic states, dynamical systems approaches are preferable since they offer the ability to intervene at any timepoint, or even at multiple timepoints, in a patient’s coagulopathic trajectory. This capability can reduce a need for urgent hospital interventions to improve physiological outcomes, given that there exist numerous unknown or unquantifiable priors such as patient arrival time to hospital, injury severity, co-morbidities, and patient genetics. Dynamical systems models can capture coagulation kinetics and physiological trauma measures to improve treatment. They can also differ in how much mechanistic coagulation knowledge is harnessed, or how much stoichiometry has to be included. In accordance with various embodiments, an exemplary coagulation treatment system 100 of the present disclosure provides a dynamic, goal-oriented, model-based, rapid trauma patient treatment strategy, as demonstrated by architecture in
Toward achieving this vision, the present disclosure describes how to quantitatively attain appropriate trauma patient treatment goals by combining appropriate actuators, blood coagulation sensors 110, thrombin dynamics models, and a control algorithm into an automated treatment delivery platform that can be physically implemented at the point-of-care. For example, blood samples can be readily obtained from trauma patients. By using blood coagulation sensors 110 and coagulation assays, coagulation factor concentrations in the blood sample can be quickly quantified. A controller algorithm can then recommend a personalized treatment plan according to a goal-oriented approach, moving the patient along a recovery trajectory toward healing. An exemplary control algorithm recommends coagulation factor concentrations for administering blood products, which act as interventions to modulate patient coagulation process dynamics, whereby this intervention approach is repeated frequently and the treatment is adjusted dynamically.
Coagulation factors or proteins are central to an exemplary control approach in that they are the actuators for trauma patient treatment and their concentrations can be rapidly measured within a few minutes using blood coagulation sensors and coagulation assays. The process of clot formation after injury proceeds according to the biochemical kinetics of the coagulation cascade, as shown in
Thrombin, factor IIa, is the end product of the coagulation cascade, and thrombin generation measures can be leveraged to predict hemostatic potential and transfusion requirements. Such measures can replace conventional coagulation tests like prothrombin time (PT), partial thromboplastin time (PTT), international normalized ratio (INR), and platelet counts, all of which have limitations. Thrombin is a unique protein that functions as both a procoagulant and an anticoagulants. As a procoagulant, thrombin activates platelets, converts fibrinogen into strands of fibrin, effects the cross-linking of fibrin to produce a firm fibrin clot by activating factor XIII, and catalyzes other coagulation-related reactions, like the activation of factors V, VIII, XI, and protein C (PC), which in turn regulate thrombin generations. As an anticoagulant, thrombin binds to thrombomodulin, a receptor protein on the endothelial membrane of a blood vessel, initiating a series of reactions that leads to fibrinolysis.
The Calibrated Automated Thrombogram (CAT) is a coagulation assay that can measure the concentration time-history of thrombin in a plasma sample. However, this CAT assay takes about 45-60 minutes to run, without including plasma sample preparation time. Such delays are far too long to be used at a patient’s bedside to predict. and guide treatment and outcomes. An exemplary thrombin dynamics model of the present disclosure can mathematically predict the concentration time-history of thrombin from patient plasma sample coagulation factor concentrations, and that thereby capture the dynamics of the coagulation system process while simultaneously replacing the CAT assay. Such a model is integrated in embodiments of an exemplary controller of a point-of-care device that enables clinicians to assess, monitor, and alter trauma patient coagulation status. In accordance with the present disclosure, a treatment process that leverages such a model can provide frequent, personalized, and dynamic recommendations based on sample clotting predictions to move a trauma patient’s coagulation state toward a desired recovery trajectory.
Conventional trauma patient therapy does not yet use quantitative coagulation factor concentration guidance, possibly because common static machine learning approaches on typical patient data with coagulation factor concentrations are uninformative. Moreover, initial biomarker and injury measurements are not correlated to treatment received, and so cannot predict resuscitation need and adverse outcomes. We found that the means of trauma patient coagulation factor concentrations do not indicate if a trauma patient is at high risk for mortality within 28 days, or at high risk for massive transfusion or a thrombotic event. Equally important, coagulation factor concentrations are uncorrelated to treatment and resuscitation: trauma patients who receive fresh frozen plasma (FFP), no matter the number of units they receive, show substantial variation in coagulation factor concentration changes over time, potentially due to a lack of characterization of, and inherent variability in, coagulation factor concentrations per FFP unit. Therefore, FFP units are not predictive of increases or decreases in coagulation factor concentrations. This also substantiates why FFP administration has mixed results for treatment, since units may not deliver required coagulation factors or may oversupply unnecessary coagulation factors in different patients at different timepoints.
Nevertheless, a close examination of the changes in coagulation factor concentrations for subgroups of 252 survivor trauma patients based on initial coagulation factor levels shows clear dynamic information over the first 24 hours after hospital admission. These dynamics can be illustrated using a heatmap of changes (Δ) in coagulation factor (CF) concentrations at different time periods in the first 24 hours (0h-6h, 6h-12h, 12h-24h), as shown by
Specifically, coagulation factor concentrations move toward an equilibrium concentration that is representative of homeostasis, where concentrations that start from a low value increase over time, while concentrations that start from a high value decrease over time, as shown in
To test the significance of the observation that coagulation factor concentrations move toward an equilibrium in patients who recover, a hypothesis (p-value) test was performed that contrasted coagulation factor concentration changes in patients who survived to those who died in the first 24 hours. Four groups were defined as follows: patients who died between 6 and 24 hours, their changes in coagulation factor concentrations between (1) 0 and 6 hours (Deceased 0-6 [hr]); and for patients who were alive at the 24 hour mark post hospital admission time (Alive), their changes in coagulation factor concentrations between (2) 0 and 6 hours, (3) 6 and 12 hours, and (4) 12 and 24 hours. Welch’s t-test was performed and p-values were calculated for α = 0.05. The null hypothesis (H0) was that the mean change in a coagulation factor’s concentration is equal for patients who are dead or alive, i.e., µx= µy where µx and µy are the deceased and alive sample means, respectively.
Given that patients who survive the first 24 hours have coagulation factor concentrations that converge to equilibrium values with the change in coagulation factor concentrations moving to zero (either from inadvertent plasma-based modulation of coagulation factor concentrations, or from innate coagulation factor compensation),
In order to administer coagulation factors to personalize trauma patient treatment, predictions of the effect of administering coagulation factors are required. Menezes et al. proposed a third-order linear dynamical systems model to rapidly predict CAT trajectories from quickly measured coagulation factor concentrations. While this model has satisfactory prediction capability, the inventors hypothesize that an embedded constraint limits its prediction accuracy and have investigated whether model improvement was possible without changing model structure, by adding a degree-of-freedom parameter to remove this underlying constraint to the input-output model as shown below:
where K0, K1, K2, Kn, and Kd are five patient-specific model parameters (the prior model used four parameters with its fifth parameter constrained), Y(s) is the predicted output thrombin concentration time-history in the frequency domain, and U(s) is a 5 pM impulse input tissue factor (TF) concentration in the frequency domain. An impulse input is an input signal with a very high magnitude that is applied to a system over a very short time. Theoretically, this magnitude approaches infinity as time goes to zero. In practice, this magnitude is taken to be some finite value, commonly 5 pM in the CAT literature, a value that also has experimental justification. Practically, the CAT is instantiated with 5 pM of TF in the plasma sample, which then rapidly depletes.
In accordance various embodiments of with the present disclosure, the initial PC concentration is included with the initial concentrations of factors II, V, VII, VIII, IX, X, and antithrombin (ATIII), creating new linear regressions for the five parameters via the same greedy method, the matching pursuit algorithm. The important role of PC in the coagulation cascade motivated this modification and it has been found that the newly updated thrombin dynamics model substantially improves CAT predictions. Among the few existing thrombin-prediction models, however, important coagulation factors like PC are typically excluded.
For the following discussions, studies performed for the present disclosure were based on normal and trauma patient data arranged into nine datasets (datasets 1, 2, 3, ..., 9). Normal data was obtained from a set of plasma samples from healthy individuals with their CAT and coagulation factor concentration measurements characterized according to standard laboratory protocols. Trauma patient data came from the Activation of Coagulation and Inflammation in Trauma study (ACIT), a previously described single-center prospective cohort study that followed severely injured trauma patients from emergency department admission through discharge from hospitalization or death.
For a dataset of 60 samples (20 individual healthy donors and 40 trauma patients), stepwise linear regression was applied that consists of sequentially and greedily adding the linear effect of a coagulation factor concentration measurement that most reduces the error of a least-squares fit to all data for each of thrombin dynamics model (1) parameters. The coagulation factor that minimizes this least square error has the greatest contribution to the system dynamics captured by that particular model parameter. The stepwise process was repeated until further linear additions of coagulation factor concentration measurements no longer improved the fit.
Visual comparisons of model improvement are shown in
Next, validation of the new improved thrombin dynamics model is presented in two ways: first with five-fold cross-validation and second on a separate dataset that was not used for training. Five-fold cross-validation bootstraps available data by subdividing it so that 80% is used for training and the remaining 20% is used for validation. The process is iterated five times for five unique divisions (folds) of the original dataset. The mean model output properties of these five iterations for the combined dataset of 20 normal samples and 40 trauma patient samples (datasets 4 and 5) are reported in
Additional model validation was accomplished with a separate validation dataset (dataset 8) that was not used for model training. This validation set started with normal plasma samples that had coagulation factor concentration and CAT measurements, and into which were spiked increasing concentrations of factors II, VIII, and X that were then quantified. An exemplary thrombin dynamics model trained on the separate 60 samples (datasets 4 and 5) can predict the 20 experimental validation CATs (dataset 8) almost perfectly, as shown in
To examine the dynamic modulation effects of coagulation factors, three experimental datasets (datasets 4, 5 and 7) were used, where dataset 7 started with normal plasma samples that had coagulation factor concentration and CAT measurements, and into which were spiked increasing concentrations of factors II, VIII, and X. New coagulation factor concentration and CAT measurements were taken after each spike.
The effects of different initial TF concentrations and coagulation factor concentration spikes on system poles were examined. The poles of a dynamical system are characteristic parameters that determine the system’s stability and output response. These poles can be obtained from a transfer function model of a system by determining the values for which the denominator of the transfer function becomes zero, i.e., we find the poles of a trauma patient’s coagulation system by setting the denominator of model (1) to zero and solving the resultant equation for s.
Each coagulation factor has a unique effect on system dynamical behavior as described by the movement of pole locations, and is often accompanied by nonlinear limitations, as illustrated by the figures of
Surprisingly, for 20 normal plasma samples from different donors, we found that increased initial TF concentration caused substantial system pole movement away from the origin, essentially recapturing trauma patient variability, as shown in
To assess personalized control of trauma patient thrombin dynamics using coagulation factors, a target goal CAT and an associated region inside which any CAT trajectories can be considered normal was determined by calculating the maximum, minimum, and mean of the experimental data at each time point for all normal plasma samples, dataset 4 and by fitting model (1) to this data. To evaluate how well the identified region represented normal, the identified region was validated using five normal samples (dataset 9) that were different from dataset 4, which was used for identification. The CAT profile of these five validation samples was contrasted against the normal region using mean relative error (MRE), the mean of the error at each time point where the profile was not within normal minimum and maximum bounds.
An exemplary control algorithm (also referred to as a “Goal-oriented Coagulation Management (GCM)” algorithm), was developed to recommend a personalized set of coagulation factor concentration changes to move trauma patients onto a recovery path.
Such an exemplary control algorithm harnesses CAT estimates from coagulation factor concentration measurements via thrombin dynamics model (1), and identifies a patient-specific mapping, as illustrated by
Algorithm treatment goals are defined to simultaneously (a) move coagulation factor concentration values toward normal equilibrium values and (b) achieve a normal thrombin (clotting) profile. To attain these treatment goals, the sequence of control algorithm operations was developed by: (i) prioritizing reaching a normal range of coagulation factor concentrations; (ii) ordering how four thrombin profile properties mimic normal clotting; and (iii) investigating the isolated effects of coagulation factors on thrombin profile properties. To satisfy (i), an exemplary control algorithm first corrects concentrations of coagulation factors that have minimal impact on thrombin profile. Thereafter, the predicted CAT is progressively corrected by modulating a coagulation factor concentration according to our new modeled dynamic interactions with the updated concentration checked to be in the normal range at the end of each thrombin profile correction step. For (ii), the control algorithm sets the order in which the algorithm tunes CAT properties as follows: thrombin generation (peak), response time (peak-time), time delay in system response (time delay), and thrombin potential (area under the curve, which is evaluated and compared using the profile tail, called “sTail”). For (iii), the inventors investigated the most impactful individual coagulation factor concentration changes on estimated CAT properties by performing numerical simulations on datasets 4 and 5, as illustrated by
Referring to
The final step of the control algorithm ensures that the recommended CAT estimate is inside the normal region. If not, then the algorithm chooses to manipulate one of two anticoagulant factors, either protein C or ATIII. The choice is made based on a comparison to the area under the normal CAT curve (thrombin potential) in its post-peak stage, based on the differing ways that protein C and ATIII alter the CAT tail. If this normal area is already surpassed by the patient’s updated CAT estimate, then protein C is selected, otherwise protein ATIII is selected. For all of the above modulations, coagulation factors are modulated only to the extent of their predefined normal limits.
For the requisite four “co-” properties that an algorithm is typically scrutinized for, the exemplary control algorithm is convergent, complete, not complex, and correct. First, the control algorithm is guaranteed to converge to a set of personalized coagulation factor concentration recommendations, because the ordered list of a finite number of coagulation factors is systematically manipulated only once through the list. Next, the control algorithm is complete in the sense that if coagulation factor concentration values exist for all eight coagulation factors to generate a simulated CAT trajectory, then the algorithm will output one possible set. Consider that a set of coagulation factor concentrations always exists consisting of the initial coagulation factor concentrations. Indeed, the control algorithm presumes these concentrations at the start before trying to modulate each concentration in turn. Earlier, we showed that the newly improved thrombin dynamics model can accurately predict CAT trajectories from coagulation factor concentrations, those that are measured before algorithm modulation. Thus, completeness is guaranteed. Third, the algorithm’s complexity is linear in the number of coagulation factors n (i.e., it is O(n) in big O notation); there is only one “for” loop in the pseudocode in
Finally, the control algorithm is correct, and its outputs have been validated against clinical outcomes of CAT profile and normalized coagulation factor concentrations for eight trauma patients (dataset 6) who showed methodical recovery toward our normal goal. This eight-patient validation dataset (dataset 6) is different from, and is not a subset of the 40 trauma patients (dataset 5) used for training the thrombin dynamics model (1). An exemplary control algorithm was validated for the first 24 hours post hospital admission as this time period accounts for 80% of hemorrhage fatalities. Intervention periods of 0, 6, 12, and 24 hours were selected for validation and the control algorithm was compared to clinical data because trauma patient data (in dataset 6) were collected at these time points.
Validation efforts show that control algorithm recommendations drive thrombin generation toward a normal region over time for eight trauma patient samples, validating the correctness and performance of the algorithm. Accordingly,
Correspondingly,
Both panels (
In summary, the pressing need for trauma patient precision medicine treatments is well-documented, but the coagulopathy problem is complex, which restricts clinicians to using rules-of-thumb, generalized treatment protocols, and uncharacterized blood products. As a result, patient recovery often fluctuates between hypo-coagulable and hyper-coagulable states, with conditions that are complicated by the side effects from contemporary non-tailored approaches. This leads to high mortality and poor outcomes in even the best trauma centers. Goal-oriented, frequent, dynamic, and patient-specific interventions are believed to be the solution, especially if the administration of coagulation factors (blood proteins) can transfer a patient onto a desirable healing trajectory.
Accordingly, exemplary systems and methods of the present disclosure can provide quantitative guidance to assist clinicians at the point-of-care by dynamically adjusting coagulation factors using patient-specifics indicated by an embedded thrombin dynamics model and recommending coagulation factor concentrations that are only within a predefined normal range. The thrombin dynamics model was improved by incorporating an additional parameter to increase model flexibility and by adding the effects of an eighth coagulation factor, protein C, because of its known coagulation importance. These modifications substantially improved the model’s thrombin dynamics predictions, which have been validated on data not used for model training. The present disclosure has verified in silico that administering coagulation factor concentrations changed the clotting that was described by the newly improved thrombin dynamics model. Coagulation factor levels have been chosen to comply with generally accepted normal limits when modulated in a treatment scheme due to observed saturating behavior for excessive coagulation factor concentration administration.
Algorithm prediction performance has been validated in silico on data not used for training by contrasting against metrics from actual trauma patients who recovered and also progressed toward normal. Validation efforts show that an exemplary control method not only guides clotting predictions closer to normal, but does so while maintaining all coagulation factor concentrations within normal ranges, which has not been the case in conventional practice.
The present disclosure offers a personalized control approach to trauma patient treatment by utilizing clotting system dynamics, characterizing the effects and limitations of coagulation factor actuators, and then articulating a control algorithm to systematically achieve coagulation goals in precision trauma patient resuscitation. Such methods and systems can be utilized in physical treatment devices at the point-of-care and can facilitate the automation of frequent, tailored clinical interventions in near real-time. An iterative approach of the present disclosure permits quicker model updates, greater personalization, and a responsiveness to uncertainties, all of which will improve patient outcomes and aid in precision trauma patient resuscitation.
In various embodiments, algorithm-recommended coagulation factor concentration increases in blood samples may be achieved by accurately adding specific recombinant coagulation factors, and algorithm-recommended coagulation factor concentration decreases may be achieved by accurately diluting samples or by augmenting inhibiting coagulation factors. In addition to trauma treatments, the results of the disclosed systems and methods may also be for various coagulation disorders, such as hemophilia, von Willebrand disease, factor V Leiden, pulmonary embolism, deep vein thrombosis, stroke, and sickle cell disease.
A comprehensive table of all the reagents and resources that were used to conduct experiments, including for coagulation factor measurements and Calibrated Automated Thrombograms, are presented in the Key Resource Table of
Coagulation factor concentrations were measured using the STA Compact Max® benchtop coagulation analyzer (having blood coagulation sensors 110) as percent activity, which is with respect to the normal coagulation factor concentration in a healthy person. A normal range for coagulation factor concentrations is typically 60-140% activity. Plasma samples were removed from -80° C. storage and thawed at room temperature. Reagents were prepared with DiH20 and left to stabilize for 30-60 minutest. Owren-Koller diluent was used for patient samples, STA-Unicalibrator reagent was used to calibrate the system by measuring/defining ranges of new reagent lots (performed monthly), STA-System Control N+P and STA-Coag Control N+ABN were control reagents measured every four hours and eight hours, respectively, and STA-Deficient reagent was used to measure the activity of a coagulation factor, e.g., STA-Deficient V was used for measuring factor V. The test automatically started after loading sample and reagents into the instrument. Given that quality control was repeated every four hours, coagulation factor concentration measurements were performed once for each sample.
For certain studies, plasma sample thrombin expression experimental data was obtained using the ThermoFisher Fluoroskan Microplate Fluorometer with Calibrated Automated Thrombogram software as follows. Plasma samples were removed from -80° C. storage and thawed at room temperature. Reagents were added to 96-well plates: thrombin calibrator reagent was used for the measurement control, and PPP-reagent was used to measure thrombin in normal or trauma samples. Plasma samples were added to plate wells, with three biological replicates, and the plate was loaded into a Fluoroskan Ascent platereader. Following a ten-minute incubation period at 37° C., the test started automatically when the machine dispensed the FluCa reagent, which was pre-loaded. Ultimately, the generated thrombin generation was measured and recorded every 20 seconds. These measurements included three technical replicates.
Correspondingly, model parameters of the newly improved thrombin dynamics model (1) were fit to experimental data using the MATLAB Simulink Design Optimization (SDO) toolbox. The input was defined as an impulse input with the desired magnitude, e.g., 5 pM of TF. The output to fit was the individual CAT profile experimental data. Solver tolerance was set to 1e-9. Starting from an initial parameter guess, the MATLAB SDO toolbox optimized parameter values of a transfer function model by minimizing the least square error between prediction and actual data using a trust region reflective algorithm. Following convergence, the finalized transfer function model parameters for each experimental sample were recorded and the poles computed.
Poles of a transfer function are the values for which the value of the denominator of the transfer function becomes zero. Therefore, to obtain the pole location values, the denominator of the fitted model was set equal to zero and the resultant equation was solved, i.e., solving s3 + K2s2 + K1s +K0 = 0. Since this is a third order system, the solution is a set of three numbers with real and imaginary parts, which can be plotted in the complex plane as shown in
For statistical analysis, the Welch’s t-test was performed for two unpaired samples (deceased, x, versus any one of the alive groups described in the main text, y) using the following equation:
where µx and µy are the deceased and alive sample means, respectively; Sx and Sy are the sample standard deviations; and n and m are the sample sizes of x and y, respectively. Next, p-values were calculated for α = 0.05, i.e., 95% confidence interval, using MATLAB’s ttest2 function. The results were reported by ns: not significant, p > 0.05; ∗ : p ≤ 0.05; ∗∗ : p ≤ 0.01; and ∗∗∗ : p ≤ 0.001.
Computer program code for carrying out operations of the present disclosure may be written in a variety of computer programming languages. The program code may be executed entirely on at least one computing device (or processor), as a stand-alone software package, or it may be executed partly on one computing device and partly on a remote computer. In the latter scenario, the remote computer may be connected directly to the one computing device via a LAN or a WAN (for example, Intranet), or the connection may be made indirectly through an external computer.
It will be understood that each block of the flowchart illustrations and block diagrams and combinations of those blocks can be implemented by computer program instructions and/or means. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, application specific integrated circuit (ASIC), or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowcharts or block diagrams.
It should be emphasized that the above-described embodiments of the present disclosure are merely possible examples of implementations, merely set forth for a clear understanding of the principles of the disclosure. Many variations and modifications may be made to the above-described embodiment(s) without departing substantially from the principles of the disclosure. All such modifications and variations are intended to be included herein within the scope of this disclosure and protected by the following claims.
This application claims priority to co-pending U.S. provisional application entitled, “Automated Platform to Treat Trauma-Induced Coagulopathy with Personalized Coagulation Factor Concentrations,” having serial number 63/324,323, filed Mar. 28, 2022, which is entirely incorporated herein by reference.
Number | Date | Country | |
---|---|---|---|
63324323 | Mar 2022 | US |