Review

Integrating biomarkers for hemostatic disorders into computational models of blood clot formation: A systematic review

  • Received: 04 September 2024 Revised: 07 November 2024 Accepted: 20 November 2024 Published: 04 December 2024
  • In the pursuit of personalized medicine, there is a growing demand for computational models with parameters that are easily obtainable to accelerate the development of potential solutions. Blood tests, owing to their affordability, accessibility, and routine use in healthcare, offer valuable biomarkers for assessing hemostatic balance in thrombotic and bleeding disorders. Incorporating these biomarkers into computational models of blood coagulation is crucial for creating patient-specific models, which allow for the analysis of the influence of these biomarkers on clot formation. This systematic review aims to examine how clinically relevant biomarkers are integrated into computational models of blood clot formation, thereby advancing discussions on integration methodologies, identifying current gaps, and recommending future research directions. A systematic review was conducted following the PRISMA protocol, focusing on ten clinically significant biomarkers associated with hemostatic disorders: D-dimer, fibrinogen, Von Willebrand factor, factor Ⅷ, P-selectin, prothrombin time (PT), activated partial thromboplastin time (APTT), antithrombin Ⅲ, protein C, and protein S. By utilizing this set of biomarkers, this review underscores their integration into computational models and emphasizes their integration in the context of venous thromboembolism and hemophilia. Eligibility criteria included mathematical models of thrombin generation, blood clotting, or fibrin formation under flow, incorporating at least one of these biomarkers. A total of 53 articles were included in this review. Results indicate that commonly used biomarkers such as D-dimer, PT, and APTT are rarely and superficially integrated into computational blood coagulation models. Additionally, the kinetic parameters governing the dynamics of blood clot formation demonstrated significant variability across studies, with discrepancies of up to 1, 000-fold. This review highlights a critical gap in the availability of computational models based on phenomenological or first-principles approaches that effectively incorporate affordable and routinely used clinical test results for predicting blood coagulation. This hinders the development of practical tools for clinical application, as current mathematical models often fail to consider precise, patient-specific values. This limitation is especially pronounced in patients with conditions such as hemophilia, protein C and S deficiencies, or antithrombin deficiency. Addressing these challenges by developing patient-specific models that account for kinetic variability is crucial for advancing personalized medicine in the field of hemostasis.

    Citation: Mohamad Al Bannoud, Tiago Dias Martins, Silmara Aparecida de Lima Montalvão, Joyce Maria Annichino-Bizzacchi, Rubens Maciel Filho, Maria Regina Wolf Maciel. Integrating biomarkers for hemostatic disorders into computational models of blood clot formation: A systematic review[J]. Mathematical Biosciences and Engineering, 2024, 21(12): 7707-7739. doi: 10.3934/mbe.2024339

    Related Papers:

    [1] Shan-Yang Lin . Salmon calcitonin: conformational changes and stabilizer effects. AIMS Biophysics, 2015, 2(4): 695-723. doi: 10.3934/biophy.2015.4.695
    [2] Joshua Holcomb, Nicholas Spellmon, Yingxue Zhang, Maysaa Doughan, Chunying Li, Zhe Yang . Protein crystallization: Eluding the bottleneck of X-ray crystallography. AIMS Biophysics, 2017, 4(4): 557-575. doi: 10.3934/biophy.2017.4.557
    [3] Satoru Ito, Kishio Furuya, Masahiro Sokabe, Yoshinori Hasegawa . Cellular ATP release in the lung and airway. AIMS Biophysics, 2016, 3(4): 571-584. doi: 10.3934/biophy.2016.4.571
    [4] Nily Dan . Bilayer degradation in reactive environments. AIMS Biophysics, 2017, 4(1): 33-42. doi: 10.3934/biophy.2017.1.33
    [5] Timothy Jan Bergmann, Giorgia Brambilla Pisoni, Maurizio Molinari . Quality control mechanisms of protein biogenesis: proteostasis dies hard. AIMS Biophysics, 2016, 3(4): 456-478. doi: 10.3934/biophy.2016.4.456
    [6] Alaa AL-Rahman Gamal, El-Sayed Mahmoud El-Sayed, Tarek El-Hamoly, Heba Kahil . Development and bioevaluation of controlled release 5-aminoisoquinoline nanocomposite: a synergistic anticancer activity against human colon cancer. AIMS Biophysics, 2022, 9(1): 21-41. doi: 10.3934/biophy.2022003
    [7] James C.L. Chow . Computer method and modeling: Medical biophysics applications in cancer therapy, medical imaging and drug delivery. AIMS Biophysics, 2021, 8(3): 233-235. doi: 10.3934/biophy.2021017
    [8] Nikolay Kalaydzhiev, Elena Zlatareva, Dessislava Bogdanova, Svetozar Stoichev, Avgustina Danailova . Changes in biophysical properties and behavior of aging human erythrocytes treated with natural polyelectrolytes. AIMS Biophysics, 2025, 12(1): 14-28. doi: 10.3934/biophy.2025002
    [9] Jonathan Trauth, Johannes Scheffer, Sophia Hasenjäger, Christof Taxis . Strategies to investigate protein turnover with fluorescent protein reporters in eukaryotic organisms. AIMS Biophysics, 2020, 7(2): 90-118. doi: 10.3934/biophy.2020008
    [10] Richard C Petersen . Free-radicals and advanced chemistries involved in cell membrane organization influence oxygen diffusion and pathology treatment. AIMS Biophysics, 2017, 4(2): 240-283. doi: 10.3934/biophy.2017.2.240
  • In the pursuit of personalized medicine, there is a growing demand for computational models with parameters that are easily obtainable to accelerate the development of potential solutions. Blood tests, owing to their affordability, accessibility, and routine use in healthcare, offer valuable biomarkers for assessing hemostatic balance in thrombotic and bleeding disorders. Incorporating these biomarkers into computational models of blood coagulation is crucial for creating patient-specific models, which allow for the analysis of the influence of these biomarkers on clot formation. This systematic review aims to examine how clinically relevant biomarkers are integrated into computational models of blood clot formation, thereby advancing discussions on integration methodologies, identifying current gaps, and recommending future research directions. A systematic review was conducted following the PRISMA protocol, focusing on ten clinically significant biomarkers associated with hemostatic disorders: D-dimer, fibrinogen, Von Willebrand factor, factor Ⅷ, P-selectin, prothrombin time (PT), activated partial thromboplastin time (APTT), antithrombin Ⅲ, protein C, and protein S. By utilizing this set of biomarkers, this review underscores their integration into computational models and emphasizes their integration in the context of venous thromboembolism and hemophilia. Eligibility criteria included mathematical models of thrombin generation, blood clotting, or fibrin formation under flow, incorporating at least one of these biomarkers. A total of 53 articles were included in this review. Results indicate that commonly used biomarkers such as D-dimer, PT, and APTT are rarely and superficially integrated into computational blood coagulation models. Additionally, the kinetic parameters governing the dynamics of blood clot formation demonstrated significant variability across studies, with discrepancies of up to 1, 000-fold. This review highlights a critical gap in the availability of computational models based on phenomenological or first-principles approaches that effectively incorporate affordable and routinely used clinical test results for predicting blood coagulation. This hinders the development of practical tools for clinical application, as current mathematical models often fail to consider precise, patient-specific values. This limitation is especially pronounced in patients with conditions such as hemophilia, protein C and S deficiencies, or antithrombin deficiency. Addressing these challenges by developing patient-specific models that account for kinetic variability is crucial for advancing personalized medicine in the field of hemostasis.



    Antibodies and proteins are crucial modalities of pharmaceutically active ingredients in drug development. Compared with small molecule drugs, antibody and protein drugs have a larger molecular size and substantial surface area. This allows selective binding to target molecules and, in turn, provides high specificity in reducing off-target adverse effects [1],[2]. Antibodies can target protein–protein interactions, whereas small molecules cannot. When the molecular mechanisms of a disease and the target molecules involved in its pathogenesis are identified, protein drugs can directly replace deficient or abnormal target proteins, and antibody drugs can interfere with the target molecules or their interactions [3]. Comparing pharmaceutical clinical trials conducted between 2001 and 2009 and 2015 and 2023, the percentage of small molecule drugs declined from 66% to 47%, while antibodies increased (from 11% to 20%) and proteins/peptides maintained their share (from 9.4% to 9.2%) [1]. Antibodies and proteins/peptides reportedly showed higher success rates in drug clinical development (11% and 5.3%, respectively) compared with small molecule drugs (4.1%) in 2015–2023 [4].

    The large molecular size of antibodies and proteins limits penetration across cellular membranes. The low pH in the stomach, enzymes (proteases and peptidases) in the gastrointestinal tract, and first-pass metabolism by cytochrome P450 in the liver and enterocytes cause the denaturation of protein drugs [5]. This very low oral bioavailability restricts administration routes to injections, such as intravenous (IV), subcutaneous (SC), and intramuscular (IM). However, IV/SC/IM dosing is subject to a number of potential risks, including needle phobia, pain, low patient adherence, phlebitis, and tissue necrosis and infection [6],[7]. To reduce these risks, sustained-release drug delivery systems (DDSs) can allow a decrease in the frequency of injections. In particular, IV infusion carries a high economic burden in medical care and medication costs. These include longer administration times, high costs to patients due to health personnel and clinical occupancy costs, lost workdays, and patient and caregiver travel expenses [8]. Among the IV/SC/IM route, SC dosing is beneficial for patient comfort and compliance and provides time and money savings [9]. SC dosing does not require premedication and requires only a few minutes compared to the several hours required for IV infusion [9]. SC formulations therefore reduce healthcare professional workload and, in turn, decrease hospital costs [10],[11]. Further, patients prefer SC to IV. Also, SC dosing might be safer and have greater tolerance than the IV route in many cases and seems to provide a better health-related quality of life [10],[11]. IM dosing carries the risk of nerve damage and improper drug deposition; it can induce pain, swelling, or irritation at the injection site and tissue damage and has limited injection volumes [9]. These comparisons have made SC dosing an attractive alternative to the IV, IM, and oral routes [9]. Recently, SC formulations have gained attention as secondary formulations after IV products. For instance, atezolizumab, an immunoglobulin (IgG) 1 kappa monoclonal antibody (mAb), which binds to programmed death ligand 1 (PD-L1), was first marketed as an IV product and subsequently as a new SC formulation [12]. Even SC formulations for trastuzumab, rituximab, and daratumumab were secondarily developed after IV products [13],[14]. Sustained-release SC formulations are considered useful technologies that increase the value of marketed drugs as new formulations in product lifecycle management. Moreover, given that some antibody fragments and proteins have a short half-life (t1/2) in blood [15],[16], SC sustained-release DDSs are critical to enabling sufficient efficacy and increasing the probability of successful new drug development.

    Table 1.  Formulation and biophysical parameters that can affect a drug's pharmacokinetics and adverse events.

    No. Parameter

    Formulation parameters 1 Drug loading amount
    2 Drug diffusion in DDS
    3 Drug interactions with polymers
    4 Design of microneedle
    Biophysical parameters 5 In vivo behavior of polymers (swelling, aggregation, dissolution, accumulation)
    6 Drug interactions with extracellular matrix and hypodermis diffusion
    7 Drug unfolding
    8 Hypodermis concentration of drug/inherent proteins/cells
    9 Absorption through lymphatics or blood vessels
    10 Degradation of drugs/polymers by proteases/phagocytic cells
    11 FcRn-mediated recycling
    12 FcγR-mediated endocytosis
    13 Target-mediated drug clearance
    14 Skin thickness/viscoelastic property and individual/site differences
    15 Skin recovery

    DDS: drug delivery system; FcRn: neonatal fragment; FcγR: Fc gamma receptor

     | Show Table
    DownLoad: CSV

    Sustained-release DDS products of small molecule and peptide drugs have been widely established for clinical use. Examples include rods, microspheres, and in situ gel formulations of poly(lactic-co-glycolic acid) (PLGA), as well as Risperdal Consta®, Zoladex®, Lupron Depot®, and Eligard® [17][20]. IV dosing of sustained-release DDS is hampered by the disadvantages of the initial burst drug release by dilution with a large volume of blood and faster systematic elimination through the reticuloendothelial system (liver, spleen, etc.) [21]. IM administration also carries the risks mentioned above. For these reasons, many DDS products that aim for the sustained release of small molecules are administered through the SC route. In contrast to many DDS products of small molecule/peptide drugs, the number of sustained-release DDS for antibody/protein drugs used in clinical practice is limited; examples include Nutropin Depot® [biweekly SC PLGA microsphere of human growth hormone (hGH)] [22] and the INFUSE® bone graft [collagen sponge of recombinant human bone morphogenetic protein-2 (rhBMP-2) in spinal-fusion procedures] [23]. The most challenging points for sustained-release DDS of antibody/protein drugs are drug release control and the maintenance of bioactivity. Among these, the hydrophilicity and high water solubility of antibodies/proteins frequently cause burst drug release. When hydrophobic polymers are used to suppress drug release, their dissolution requires organic solvents, but these in turn affect the conformation and bioactivity of the antibodies/proteins [24]. Most small molecules can rapidly diffuse into the blood without degradation after SC dosing. In contrast, the movement and fate of the DDS and antibody/protein drugs are affected by many formulations and biophysical parameters at the dosing site (Table 1). This review summarizes two kinds of SC dosing of DDS: subcutaneous injection and intradermal microneedle (MN). Several embodiments of sustained-release DDS for antibody/protein drugs are described, including PLGA microspheres for hGH, hydrogel particles for interferon-α2b (IFNα-2b), amphiphilic copolymer hydrogel for bevacizumab, conjugated hydrogel microspheres for anti-tumor necrosis factor alpha (TNFα) single-chain antibody fragment (scFv), and dissolving/hydrogel MN for bevacizumab. The effects of the formulation and biophysical parameters (Table 1) on the in vivo performance of DDS, such as pharmacokinetics (PK) and adverse events, are also discussed. Information and insights from this paper will lead to improvements in the speed and precision of SC sustained-release DDS, which will in turn expand their application to broader antibody/protein drugs and deliver high value to many patients.

    The in vivo fate of antibody/protein drugs after SC injection is affected by various biophysical parameters (Figure 1; Table 1). The first is the drug diffusion into the interstitial space of the hypodermis from the injection site. Tissue interstitial fluid accounts for 70% of the total volume of the hypodermis, while the extracellular matrix (ECM) is the main component of the other 30% [25],[26]. The fluid is an ultrafiltrate of plasma composed of 97% water and 3% proteins, with a flow on the order of 0.1–2 µm/s [26]. The ECM includes proteoglycans, collagen, elastin, laminin, and fibrin. Due to the large size of antibody/protein drugs and the mesh size and high viscosity of the ECM, drug–ECM interactions significantly reduce the diffusion of the drugs to 18%–93% of their free solution values [26] (Table 1, no. 6). Proteins with a larger size [or molecular weight (MW)] are reported to remain for a longer time at the SC injection site [27]. The time required for a 50% loss of protein from the dosing site was 6.8 h for bevacizumab (MW 149 kDa), 2.85 h for bovine serum albumin (BSA) (MW 66 kDa), 1.57 h for ovalbumin (MW 44.3 kDa), and 0.31 h for vascular epithelial growth factor Cys156Ser (VEGF-C156S) (MW 23 kDa) [28]. Proteoglycan molecules contain highly negatively charged glycosaminoglycan chains, which interact with antibody/protein drugs [29]. Antibodies with high positive charges and hydrophobic interactions reduce absorption rates and bioavailability after SC injection [30] (Table 1, no. 6).

    Figure 1.  Subcutaneous injection and structure of skin, blood, and lymph vessels.

    On the other hand, excipients with small molecular sizes, such as histidine, sucrose, and arginine, are frequently used in normal formulations of protein/antibody drugs and diffuse away from the dosing site much faster than antibody/protein drugs [31]. This changes the pH and salt concentration around the drugs after injection, which may in turn cause the precipitation of drug aggregates [32]. These aggregates are retained for a longer time (10–20 h) at the site due to their large size [33][35] (Table 1, no. 6). Moreover, proteases such as matrix metalloproteinases in the hypodermis might degrade antibody/protein drugs (Table 1, no. 10). Degradants of PEGylated erythropoietin (PEG; polyethylene glycol) were observed after incubation with rat subcutaneous tissue homogenate [36]; in contrast, degradation was limited after incubation of bevacizumab, BSA, and ovalbumin with the homogenate [28].

    The length of IgG is around 10 nm, while blood capillaries have a pore size of 5–12 nm, and the optimal molecular size for lymphatic uptake is 10–100 nm [37],[38] (Figure 1). The lymphatic absorption of molecules with different MW after SC dosing to sheep was as follows [37],[39]; 4%, 21%, 17%, 39%, 60%, 62%, 96%, and 99% for 5-fluoro-2′-deoxyuridine (MW 246 kDa), inulin (MW 5.2 kDa), insulin (MW 5.8 kDa), cytochrome c (MW 12.3 kDa), human recombinant IFNα-2a (MW 19.0 kDa), hGH (MW 22 kDa), recombinant human erythropoietin (MW 30 kDa), and darbepoetin α (MW 37 kDa), respectively. This linear relationship (Figure 2) suggests that molecules with MW > 16.5 kDa are absorbed mainly by the lymphatic system (Table 1, no. 9).

    Figure 2.  Association of protein molecular weight with the percentage of lymphatic absorption after subcutaneous dosing to sheep [37],[39].

    Once absorbed, antibody/protein drugs are transported to the collecting lymphatics and pass through lymph nodes [29]. Antigen-presenting cells (APCs) in the lymph nodes constantly sample the surrounding medium by micropinocytosis to degrade proteins [27] (Table 1, no. 10). Although antibodies can bind to neonatal fragment crystallizable receptor (FcRn) to escape lysosomal degradation, when binding is saturated at the administered dose, a significant fraction of antibody is catabolized in APC lysosomes (Table 1, no. 11). In previous work, after incubation of proteins with cell suspensions or homogenate of lymph nodes, some portion of PEGylated erythropoietin and bevacizumab was probably degraded, whereas BSA and ovalbumin were not [28],[36]. Although these data may have been affected by a reduction in APC activity or measurement accuracy, catabolism in the lymph nodes appears to vary with the specific antibody/protein involved.

    Growth hormone (GH or somatotropin) is a single-chain protein containing 191 amino acids (MW 22 kDa) in a four-bundle α-helix protein with no β-sheet structure that is internally cross-linked by two disulfide bonds with an isoelectric point (pI) near 5.3 [40]. GH is secreted episodically from somatotroph cells in the anterior pituitary and stimulates growth and cell reproduction. For patients with GH deficiency, recombinant human GH (rhGH) is the primary treatment. In the 1980s, GH was administered by daily SC injection [41]. However, the plasma t1/2 of rhGH is very short, namely 0.36 h after IV dosing and 3.8 h after SC dosing [42],[43]. The once-monthly product Nutropin Depot® was approved by the United States Food and Drug Administration (FDA) to provide a more stable blood GH concentration over time and reduce dosing frequency [44],[45]. The product consists of PLGA microspheres manufactured using prolease technology, which releases GH gradually.

    Figure 3.  Chemical structures of polymers used in drug delivery systems for sustained release of antibodies and proteins.

    PLGA is a biodegradable synthetic polymer (Figure 3a). PLGA microspheres are one of the most common sustained-release DDSs for small molecule and peptide drugs, used in Lupron Depot® (leuprolide acetate), Risperidal Consta® (risperidone), and Zilretta® (triamcinolone acetoamide), among others [46][48]. The hydrophobic polymer is dissolved in a volatile organic solvent, commonly dichloromethane, and small molecule or peptide drugs are encapsulated in the microspheres by the oil-in-water (O/W) or water-in-oil-in-water (W/O/W) solvent evaporation/extraction method. The large interface between the aqueous and organic phases in the emulsification causes the adsorption and denaturation of proteins [49] (Table 1, no. 7). In O/W emulsification of lysozyme aqueous solution with a PLGA dichloromethane solution, a significant amount (20%–37%) of lysozyme was aggregated at the aqueous/organic interface, and protein recovery was found to be irreproducible [49]. Prolease technology aims to avoid the destruction of the tertiary structures of GH and increase encapsulation efficiency using a cryogenic and nonaqueous preparation process. It consists of the following steps:  (1) freeze-drying, in which a solution of the drug and excipients is sprayed into liquid nitrogen; (2) dispersal of the drug particles into PLGA organic solution and sonication to reduce particle size; (3) atomization of the suspension into liquid nitrogen to produce frozen microspheres; (4) extraction of the solvent with ethanol; and (5) filtration, vacuum-drying, and sieving to remove large particles [45]. The solvents used in this process (dichloromethane and ethanol) do not dissolve rhGH, resulting in a high encapsulation efficiency of more than 98%.

    The microspheres have a median volume diameter of 50 µm and contain 15 wt% of rhGH [42],[45] (Table 1, no. 1). Although this drug loading ratio is not high compared with successful examples of PLGA microspheres for peptides/proteins [50], an amount sufficient for microsphere medication containing 18 mg of rhGH can be dispersed in 1 mL of diluent for dosing [51]. The large size of the hydrophobic microspheres can slow polymer degradation, prevent first burst drug release, and suppress diffusion/release rates of the hydrophilic rhGH (Table 1, no. 2). Several assays have confirmed the maintenance of rhGH integrity after processing and release, including a cell proliferation assay that measures rhGH bioactivity. On the other hand, released drugs showed significant aggregation one day after release in a buffer. To stabilize rhGH from denaturation and reduce water solubility to ensure slower release (Table 1, no. 2), a complex of rhGH with zinc (at a molar ratio of 6: 1) was developed. rhGH solubility sharply decreased at a zinc:rhGH molar ratio above 2. The zinc binds at histidine (His)-18 and His-21 on helix 1 and glutamine (Glu)-174 on helix 4 to form a dimer of rhGH [40]. The precipitated rhGH zinc complex maintained the same secondary structure. Conformation of tyrosine and phenylalanine residues was somewhat more rigid, and the angle of the hydrogen bond between aspartic acid (Asp)-169 and tryptophan (Trp)-86 was altered; however, the tertiary structural changes were fully reversible [40], and complete refolding was confirmed in another study. The zinc complex completely prevented the aggregation of proteins released from the microspheres for 28 days [45] (Table 1, no. 7). The zinc complexation also reduced rhGH solubility to minimize the initial burst release. Moreover, sonication of the spray freeze-dried drug particles decreased their size from 10 to 2 µm, reducing the initial release by 25%–33%. PLGA lactide:glycolide monomer ratio (50: 50), carboxyl polymer terminus, and 1% content of zinc carbonate (insoluble salt of a release modifier) were optimized to ensure the most constant serum levels of rhGH and insulin-like growth factor-I (IGF-I) after SC dosing to monkeys. Bioactive hGH stimulated production of IGF-I in vivo [45].

    After these improvements to the formulation, manufacturability was also improved, and safety and efficacy were confirmed in preclinical and clinical studies. Nutropin Depot® was approved in 1999. In clinical data from 22 children with GH deficiency [52], a single SC injection of the product (0.75 and 1.5 mg/kg GH dose) provided a serum GH level above baseline for 13–15 days (Table 2). Serum IGF-I levels were also maintained above baseline for 16–20 days [52]. PK modeling estimated that the relative bioavailability of Nutropin Depot® was around 50% (vs. daily rhGH formulations; Table 2). No progressive increase in peak or trough levels of GH were observed after repeated dosing, indicating no accumulation at the dosing site or body [52].

    Table 2.  Pharmacokinetic parameters of growth hormone in humans.

    Formulation Dosing route Animal Dose (mg/kg) Cmax (ng/mL) t1/2 (h) TC > 1 µg/mL (days) Relative bioavailability % (vs. IV) References

    Humatrope® IV human 0.02 415 0.36 - [43]
    Humatrope® SC human 0.1 63.3 3.8 ~ 1.3 81 [43]
    Nutropin Depot® SC human 0.75 54 - 13 ~ 50 [52]
    Nutropin Depot® SC human 1.5 80 - 15 ~ 50 [52]

    IV: intravenous dosing; SC: subcutaneous dosing; Cmax: maximum concentration; t1/2: half-life; TC > 1 µg/mL: time when growth hormone concentration is greater than 1 µg/mL.

     | Show Table
    DownLoad: CSV

    In general, PLGA microspheres are gradually hydrolyzed in the hypodermis after SC dosing [53] (Table 1, no. 5). The biodegradation steps include hydrolytic chain cleavage on their surface, pore diffusion, erosion, and autocatalysis by carboxylic end groups [53],[54]. The initial PK profile can be interpreted to mean that rhGH on the microsphere surface was rapidly dissociated, and the initial burst drug release by this surface degradation lasted approximately 48 h [52] (Table 1, no. 5). The long-term sustained GH level was probably caused by subsequent disintegration of the microspheres providing sustained drug release [52]. rhGH is a small anionic inherent protein (pI = 5.3, MW 22 kDa), consequently less likely to interact with the ECM. Its diffusion in the hypodermis was probably rapid (as described in Section 2, Table 1, no. 6). After being released, the rhGH might adsorb on the microspheres to cause unfolding (Table 1, no. 7). Adsorption of rhGH on PLGA nanoparticles was reported in a pH 7.2 buffer [55]. The pI of rhGH is 5.3, and PLGA has anionic charges; accordingly, this adsorption appears to be primarily attributable to hydrophobic interaction. Loss of the tertiary structure of adsorbed rhGH has been suggested, and entropic gain for the unfolding appears to be a primary driving force for adsorption [55]. In the physiological environment, it is unknown whether this adsorption actually occurs or whether the unfolding is reversible or not. It is possible that this adsorption and unfolding might contribute to decreased or slower GH release and low bioavailability.

    The lymphatic absorption of rhGH was 62% [37] (Table 1, no. 9), suggesting the possibility that APCs might degrade the rhGH in the lymph nodes (Table 1, no. 10). After intratracheal dosing of fluorescein isothiocyanate (FITC)-labeled hGH to rats [56], rapid and significant uptake of hGH into alveolar macrophages was observed [56],[57]. The large microparticle size of this product requires a large-diameter 21 G needle (outer diameter 0.813 mm, inner diameter 0.495 mm) for injection [46]. The main adverse events are related to the injection site, including nodules (60% of injections), erythema (54%), and postinjection pain (36%) [42] (Table 1, no. 15). Generally speaking, dosing of sustained-release DDS and PLGA microspheres causes acute and chronic inflammation, granulation tissue formation, foreign body reaction, and fibrous encapsulation in sequence at the administration site [53],[58]. Inflammatory cells respond to the microspheres and the mechanical injury created by injection. Monocytes migrate and differentiate into macrophages to form foreign body giant cells at the microsphere surface. Fibroblasts infiltrate and secrete collagen to form a fibrous capsule surrounding the implant site. Before full degradation, these adverse events are probably caused primarily by reaction to the microspheres and also by injury by the large-diameter needle [42],[46],[58]. The accumulation of macrophages around the microspheres and in the lymph nodes might take up some portion of the released rhGH, which might cause low bioavailability (Table 1, no. 10).

    Human IFNα-2b is a glycoprotein cytokine consisting of 166 amino acids (MW 19.2 kDa) with five α helices, two di-sulfide bridges, and an O-glycosylated threonine (Thr)-106 with a pI around 5.9 [59]. IFNα-2b is produced by leukocytes during viral infection. It has wide biological activity and is used in hepatitis and cancer treatments. Plasma t1/2β of IFNα-2b after IV and SC injection in humans is approximately 1.7 h and 2.9 h, respectively [60] (Table 3). IFNα-2b has a relatively small therapeutic index, and a high dose causes substantial toxicity [61]. To reduce dosing frequency and control drug plasma concentration level, weekly IFNα-2bXL has been developed using the Medusa hydrogel sustained-release DDS.

    Table 3.  Pharmacokinetic parameters of interferon α-2b.

    Formulation Dosing route Animal Dose Cmax (ng/mL) t1/2, T50%AUC, or TSA Relative bioavailability % (vs. IV) References

    Normal solution IV human 5 x 106 IU/m2 188 t1/2α = 0.1 h, t1/2β = 1.7 h [60]
    Normal solution SC human 5 x 106 IU/m2 55.3 t1/2β = 2.9 h > 100 [60]
    Normal solution SC rat 60 µg/kg 1.5 T50%AUC = 1.4 h, TSA = 0.5 day - [62]
    IFNα-2bXL (Medusa gel) SC rat 60 µg/kg 27.9 T50%AUC = 16.8 h, TSA = 2 days 58 [62]
    Normal solution SC dog 60 µg/kg 1.4 T50%AUC = 4.3 h, TSA = 1 day - [62]
    IFNα-2bXL (Medusa gel) SC dog 60 µg/kg 28 T50%AUC = 39 h, TSA = 4 days 43 [62]
    Normal solution SC monkey 60 µg/kg 3.4 T50%AUC = 3.1 h, TSA = 1 day - [62]
    IFNα-2bXL (Medusa gel) SC monkey 60 µg/kg 59.6 T50%AUC = 30 h, TSA = 6 days 53 [62]
    Viraferon® SC human 3 MIU 32 TSA = 1 day - [62]
    IFNα-2bXL (Medusa gel) SC human 18 MIU 23 TSA > 7 days - [62]

    IV: intravenous dosing; SC: subcutaneous dosing; IU: international units; MIU: million international units; Cmax: maximum concentration; t1/2α: half-life during distribution phase; t1/2β: half-life during terminal elimination phase; T50%AUC: time required for 50% of total area under the curve (AUC); TSA: time when interferon-α2b is absorbed; Tmax: time to maximal concentration.

     | Show Table
    DownLoad: CSV

    The Medusa technology uses a polymer consisting of a hydrophilic poly-L-glutamate backbone and grafted hydrophobic α-tocopherols [62] (Figure 3b). The glutamate chains aggregate via the self-assembly of α-tocopherol groups to form hydrogel nanoparticles. The polymer can interact with drugs in the nanoparticles and release them in a sustained manner. SC dosing of the Medusa hydrogel nanoparticles also showed in vivo sustained absorption of IFNα-2b and interleukin-2 (IL-2) in nonclinical and clinical studies [62]. Simple mixing of two aqueous solutions of IFNα-2b and the polymer causes physical binding of the protein to the particles without the use of an organic solvent (Table 1, no. 3). An excess amount of BSA fully displaces the IFNα-2b, and the displaced drug retains the same retention time on HPLC and bioactivity as the native protein [62]. After single SC injection to rats, dogs, and monkeys, IFNα-2bXL prolonged in vivo release by 4, 4, and 6 times compared with normal IFNα-2b solution, respectively (Table 3). The nanoparticles produced a decrease in maximum concentration (Cmax) to around 5% that of the normal solution and showed 43%–58% relative bioavailability to it [62] (Table 3). In a Phase I/II study in 53 patients [63], the absorption time of the immediate formulation of Viraferon® and IFNα-2bXL were 1 day and more than 7 days, respectively. The IFNα-2bXL reduced Cmax to 10% (Table 3) and reduced the rate of adverse events to 3.7%, compared to 5.2% for Viraferon® [62],[63].

    Poly-L-glutamate (MW 10.0–16.9 kDa) is used in this technology [64]; each particle contains 10–15 polymer chains and the product contains IFNα-2b 0.3 mg/mL and 23 mg/mL of the polymer [62]. It was therefore calculated that 14% of particles in the product contained one IFNα-2b molecule and that 86% of particles were empty, without any proteins. Since one IFNα-2b molecule (MW 19.2 kDa) was surrounded by 9–14 polymer molecules (MW 10.0–16.9 kDa), the release of the IFNα-2b is probably quite restricted (Table 1, no. 1). The Medusa hydrogel nanoparticles have a small size of 20–50 nm, and drug diffusion in the particles can be neglected (Table 1, no. 2). Regarding the drug–polymer interaction, the poly-L-glutamate might have high affinity with a conformation of the IFNα-2b (Table 1, no. 3). This Medusa hydrogel particle has also shown in vivo sustained release of different protein drugs of IL-2 and hGH [62],[64]. After SC dosing of IL-2 to rats, dogs, and monkeys, release time was prolonged by about 1.5 times compared with the immediate release product of Proleukin® (Table 4). In a Phase I/II study in renal carcinoma patients, a single SC injection of Medusa particle suspension extended T > 5 ng/mL (time when serum IL-2 concentration is greater than 5 ng/mL) by 2.3 times, and time to maximal concentration (Tmax) by 16 times compared with Proleukin® [62] (Table 4). In SC dosing of hGH to dogs [64], this technology extended T > 5 ng/mL by 7 times and T50%AUC (time required for 50% of total hGH area under the curve, AUC) by 30 times compared with the immediate release solution [62] (Table 5). The effect of this technology on the PK profiles of the three protein drugs was similar. Regarding the biophysical properties of the three proteins, their pIs are different (6.8, 5.9, and 5.3 for IL-2, IFNα-2b, and hGH, respectively); however, they similarly have slight negative charges in the hypodermis at pH 7.4. Poly-L-glutamate has a low pI of 3.9 and exerts a strong negative charge on the particles, raising little concern about ionic interactions or unfolding of the drugs (Table 1, no. 3 and 7). The three proteins have similar sizes (MW 19.2, 15.5, and 22 kDa for IFNα-2b, IL-2, and hGH) and similar tertiary structures. hGH and IL-2 have the same up-up-down-down four-helix bundle motif [65] (Figure 4), and the topology of IFNα-2b resembles it [66]. This motif is classically viewed as an α-helical coiled coil [67], and poly-L-glutamate tracts have been reported to affect the coiled-coil regions and stabilize protein interactions [68]. These protein structures and interactions might contribute to their sustained release from the Medusa particles (Table 1, no. 3).

    Table 4.  Pharmacokinetic parameters of interleukin-2 after dosing of Medusa hydrogel or Proleukin®.

    Formulation Dosing route Animal Dose (IU/m2) Cmax TSA (days) Tmax (h) TC > 5 ng/mL (h) Relative bioavailability % (vs. IV) References

    IL-2XL (Medusa gel) SC rat unreported 1/100 of Proleukin 3 (Proleukin: 2) - - 10 [62]
    IL-2XL (Medusa gel) SC dog unreported 1/13 of Proleukin 5 (Proleukin: 3) - - 40 [62]
    IL-2XL (Medusa gel) SC monkey unreported 1/2 of Proleukin 3–4 (Proleukin: 2–3) - - 100–150 [62]
    Proleukin® SC human 10.6 × 106 ~ 5.3 - 2–4 3 - [62]
    IL-2XL (Medusa gel) SC human 10.6 × 106 ~ 1.6 - 48 7 - [62]

    IV: intravenous dosing; SC: subcutaneous dosing; IU: international units; Cmax: maximum concentration; TSA: time when interferon α-2b is absorbed; Tmax: time to maximal concentration; TC > 5 ng/mL: time when drug concentration is greater than 5 ng/mL.

     | Show Table
    DownLoad: CSV
    Table 5.  Pharmacokinetic parameters of growth hormone in dogs.

    Formulation Dosing route Animal Dose (mg/kg) Cmax (ng/mL) TC > 5 ng/mL (h) T50%AUC (h) References

    Normal solution SC dog 1 582 18 4 [64]
    Medusa gel SC dog 1 25 129 121 [64]

    SC: subcutaneous dosing; Cmax: maximum concentration; TC > 5 ng/mL: time when drug concentration is greater than 5 ng/mL; T50%AUC: time required for 50% of total AUC.

     | Show Table
    DownLoad: CSV
    Figure 4.  Tertiary structures of proteins.

    In hypodermis, the poly-L-glutamates might stably aggregate for an extended period (Table 1, no. 5). Research into the solution-state structures of poly-L-glutamates of 10–20 residues revealed that longer residues increased peptide-peptide hydrogen-bonding contacts and β-strand-like conformation and induced aggregation [69]. Poly-L-glutamates of more than 36 residues long in huntingtin protein are observed in Huntington's disease patients. This aggregation is associated with the faster onset of symptoms [69]. With the Medusa technology, the proportion of grafted tocopherol is 5.1 mol% [64], and the poly-L-glutamate (MW 16.9 kDa) is calculated to contain 130 residues, with an average number of glutamate residues between the tocopherol groups of around 26. Therefore, strong polymer aggregation might contribute to the stable adsorption and sustained release of some protein drugs. Moreover, α-tocopherol may also contribute to the in vivo stability of the polymer (Table 1, no. 5). With regard to α-tocopherol hydrophobic nanodomains, the aggregation number of D-α-tocopheryl polyethylene glycol succinate (TPGS) micelles (number of surfactant molecules in a micelle) is > 100; this is higher than those of other nonionic amphiphiles and is caused by the dense packing of tocopherol moieties within the core [70]. In the Medusa particles, 10–15 α-tocopherol moieties were assembled, probably due to steric hindrance from the polymer length; however, the dense packing of 10 nanodomains/particle probably contributes to the high stability of the particles and the slow water penetration and sustained degradation of the particle [62]. The high stability of the Medusa hydrogel likely helped the sustained drug release in hypodermis.

    Drug release from the hydrogel may be influenced by the hypodermal concentration of inherent BSA (Table 1, no. 8). In an in vitro study, approximately 80% of IFNα-2b was displaced by BSA (MW 66.5 kDa) at [BSA]/[polymer] weight ratio = 10, molar ratio = 2.5 [62]. This suggested that 220 BSA molecules were necessary for the displacement of one IFNα-2b molecule. Comparison of the interstitial albumin concentration (7.36 mg/mL) in human adipose tissue [71] and the polymer concentration of the product (23 mg/mL) also imply a limited in vivo displacement, and this might contribute to sustained release and absorption. After release, IFNα-2b is an anionic and small-sized protein (pI = 5.9, MW 19.2 kDa); as a consequence, it probably did not interact with the ECM, and its hypodermis diffusion was fast (as described in Section 2, Table 1, no. 6). The released IFNα-2b maintained bioactivity, confirming that its folding was correct. The first reason is the mild preparation process for the hydrogel with mixing of the aqueous solutions only. Further, the biocompatible poly-L-glutamates with a negative charge (pI of the polymer = 3.9) probably contribute to the maintenance of IFNα-2b's tertiary structure (Table 1, no. 7). The lymphatic absorption of IFNα-2b is estimated to be approximately 60% (Table 1, no. 9) from the linear correlation (Figure 2), indicating the possibility of drug degradation by APCs in the lymph nodes. Moreover, in the hypodermis, Medusa's polymer is degraded by macrophages [62], and many APCs are probably recruited around the hydrogel. After repeated SC dosing of the hydrogel in rats and monkeys, low-grade inflammation was observed at the injection site [62]. Primarily mononuclear cells, lymphocytes, and/or foamy macrophages were gathered at the cutis and/or subcutis of the site. In the absorption of IFNα-2b, APCs in the hypodermis and the lymph nodes might degrade IFNα-2b to cause sustained absorption and low bioavailability (Table 1, no. 10).

    Bevacizumab is a recombinant humanized anti-VEGF-A monoclonal IgG1 antibody (MW 149 kDa) with pI = 8.3 [72]. Bevacizumab specifically binds to VEGF-A, which interacts with VEGF receptor (VEGFR)-1 and VEGFR-2 and inhibits angiogenesis and tumor growth. It is used in the treatment of various cancers, including metastatic colorectal cancer, non-small cell lung cancer, metastatic renal cell carcinoma, cervical cancer, ovarian/fallopian tube/peritoneal cancer, glioblastoma, and metastatic hepatocellular carcinoma [73]. IV infusion of bevacizumab to patients showed an initial half-life t1/2α of 1.4 days and a terminal half-life t1/2β of 20 days [74] (Table 6). Biweekly infusion of 30–90 min places a high cost and burden on patients and health personnel, and a sustained-release DDS for bevacizumab has been investigated [75][80].

    In one study, a physically cross-linked hydrogel of ABA triblock amphiphilic copolymers was used for delivery [75]. The polymer's A and B blocks were vitamin D-functionalized polycarbonates (V) and PEG (MW 20 kDa), respectively (V-PEG-V) (Figure 3c). The preparation step of the hydrogel involved the addition of bevacizumab aqueous solution to the polymer only. Hydrogel containing 5 wt% polymer allowed quick and easy injection through a 23 G needle and demonstrated a shear-thinning property for recovery of a solid-like state to control drug release after injection. This formulation released the monoclonal antibody (mAb) gradually in a dissolution test. In an efficacy evaluation of a single dose of bevacizumab in tumor-bearing mice, IV and intraperitoneal injections of mAb normal solution showed no efficacy against tumor progression. In contrast, SC injection of the hydrogel remarkably reduced tumor growth [75]. After SC dosing of fluorophore-labeled bevacizumab hydrogel, the mAb was localized to the injection site and gradually decreased with time. In the bevacizumab PK profile after IV dosing of the solution, Cmax was 4.5 µg/mL, t1/2α was 2–6 days, and the concentration decreased to 0.6 µg/mL at 14 days [75] (Table 6). In contrast, the hydrogel SC dosing maintained a constant bevacizumab plasma level (0.5–1.0 µg/mL) for 1–21 days [75] (Table 6). Bioavailability of the SC V-PEG-V hydrogel was estimated to be 50%–70% [75], which was lower than that of the SC bevacizumab solution (> 100%) in another study [81] (Table 6).

    Table 6.  Pharmacokinetic parameters of bevacizumab after dosing in humans and mice.

    Formulation Dosing route Animal Dose (mg/kg) Cmax (µg/mL) t1/2α t1/2β (days) Characteristics for sustained absorption Relative bioavailability % (vs. IV) References

    Normal solution IV human 5 125 1.4 days 20 - [74]
    Avastin® IV mouse 5 ~ 4.5 2–6 days - - 100 [75]
    V-PEG-V hydrogel SC mouse 5 ~ 1.0 - - constant 0.5–1.0 µg/mL for 1–21 days ~ 50–70 [75]
    Normal solution IV mouse 9.3 174 1.2 h 6.8 - 100 [81]
    Normal solution SC mouse 9.3 74.1 6.2 h 6.1 - > 100 [81]

    IV: intravenous dosing; SC: subcutaneous dosing; Cmax: maximum concentration; t1/2α: half-life during distribution phase; t1/2β: half-life during terminal elimination phase.

     | Show Table
    DownLoad: CSV

    Regarding the in vitro bevacizumab release from the V-PEG-V hydrogel (Table 1, no. 2), an initial burst release was not observed, and 65% and 80% of the drug was released in 5 and 9 days, respectively [75]. Hydrogels containing 4, 5, and 6 wt% polymer showed storage moduli (G′) in proportion to concentration and the same release profiles, which suggested that differences in gel stiffness and porosity did not affect drug diffusion in the polymeric network. Since the hydrodynamic radius of mAb (about 7 nm) is smaller than the pore size of the hydrogels (2–20 µm diameter), the network might not act as a diffusion barrier, and the sustained release may be associated with physical interactions between the polymer and mAb (Table 1, no. 3). Another study reported the release profile of tenofovir (MW 476) from the same V-PEG-V hydrogel [82]. Tenofovir is an analog of adenosine monophosphate with an aqueous solubility around 2.5 mg/mL. In the dissolution test, 40% and 80% of tenofovir were released in 0.1 and 2 h, respectively. Although the test conditions were slightly different, the discrepancy in 80% release time between 9 days (bevacizumab) and 2 h (tenofovir) was significant. The in vivo retention of bevacizumab, at about 50% and 25% of the mAb in the dosing site at 2 and 6 days after SC dosing, is similar to the in vitro release profile. Regarding drug solubility, a 100–200 mg/mL mAb aqueous solution without either organic solvents or surfactants is used in several products administered to patients [83]. Accordingly, the water solubility of mAbs is generally high, and solubility likely does not explain the slow release of bevacizumab. In another study, a hydrogel of a similar triblock copolymer gradually released trastuzumab in a dissolution test [84]. In the polymer, vitamin D-functionalized polycarbonates of V-PEG-V were replaced with vitamin E-functionalized polycarbamates. Trastuzumab release showed no initial burst release, and 40% and 60% of the drug was released in 10 and 42 days, respectively [84]. The carbonate and carbamate moieties of the polymers can form hydrogen bonds with amide bonds of the loaded mAb, which might contribute to a sustained release (Table 1, no. 3).

    Regarding the degradation process of the hydrogel (Table 1, no. 5), intact hydrogel remained at 4 days after SC dosing in mice. Gel thickness was noticeably reduced at 1 week, and the hydrogel was completely degraded by 6 weeks [75]. The time course of hydrogel clearance was similar to that of drug retention/clearance at the injection site, which indicated that the in vivo release of bevacizumab was controlled by erosion of the hydrogel [75]. The size of the released bevacizumab is large (MW 149 kDa), and its charge is cationic (pI = 8.3) [72]; accordingly, the drug probably bound with the anionic ECM and slowly diffused in the hypodermis, as described in Section 2 (Table 1, no. 6). In SC dosing of the normal solution to mice, the t1/2 of drug clearance from the dosing site (6.8 h) [28] was consistent with the t1/2α of plasma drug level (6.2 h) [81] (Table 6). This may suggest that the absorption rate of SC-dosed bevacizumab solution is mainly dominated by drug diffusion in the hypodermis. Regarding the folding of the released bevacizumab (Table 1, no. 7), an ELISA test using an IgG-Fc antibody detected the drug at 3 weeks post-dosing. Further, its anti-tumor efficacy lasted 2 months. These findings may support the maintenance of the tertiary structure, albeit no clear evidence was obtained [75]. Regarding the absorption routes (Table 1, no. 9), almost all bevacizumab was probably absorbed via the lymphatic route due to its large size (MW 149 kDa, Figure 2). In another study, 17% of bevacizumab's heavy chain was degraded after 24 hour incubation with lymph node homogenate [28]. Since the V-PEG-V hydrogel gradually released the drug in vivo, bevacizumab concentration in the lymph nodes may be lower after dosing of the gel compared with dosing of the normal drug solution. This lower drug concentration might result in the degradation of a more significant portion of the drug by APCs in the lymph nodes and the lower bioavailability of the V-PEG-V hydrogel (Table 1, no. 8 and 10).

    The same phenomenon in hypodermis might also contribute to low bioavailability (Table 1, no. 8 and 10). Its shear-thinning property enabled the injection of the hydrogel through a thin needle (23 G), which likely, in turn, caused little inflammation (Table 1, no. 14–15). However, SC dosing of the blank hydrogel resulted in infiltration of some inflammatory cells into the hydrogel [75]. CD45-positive cells include macrophages, dendritic cells, and T cells. Although dendritic cells and T cells might be associated with foreign body reactions, macrophages are central mediators and play key roles in the phagocytosis of injected materials [85]. Macrophages express Fc gamma receptor (FcγR) and the highest level of FcRn [86] (Table 1, no. 11–12). Moreover, uptake by APCs was higher for mAb with a more positive charge [87], and bevacizumab (pI = 8.3) is cationic in the hypodermis [72]. All considered, together with the implications of other studies [5],[25],[87],[88], the presence of many macrophages around the hydrogel might reduce the Cmax and bioavailability of bevacizumab in the following manner: at a high concentration of bevacizumab, the free fraction of FcRn in macrophages is decreased. A significant fraction of bevacizumab cannot be recycled by binding to FcRn but rather binds to FcγR and is degraded in lysosomes (Figure 5). If the concentration of bevacizumab in hypodermis after normal solution dosing >> FcγR capacity > the concentration after gel dosing > FcRn capacity, then FcγR-mediated degradation would be negligible in the normal solution but would be more significant in the sustained-release hydrogel (Table 1, no. 8 and 11–12).

    Figure 5.  Recycling and degradation of antibody.

    This lower bioavailability might also be associated with target-mediated clearance (Table 1, no. 13). After mAb binds to the target antigen, the complex may be internalized and degraded in lysosomes [5] (Figure 5). For example, clearance was shown to be important for alemtuzumab (mAb targeting CD52) [89], evolucumab (mAb targeting protein convertase subtilisin kexin type 9: PCSK9) [90], and mAbs targeting high-density and rapidly internalized antigens such as human epidermal growth factor receptor-2 (HER2) [91]. In bevacizumab treatment for patients, bevacizumab serum concentration is approximately 200–400 µg/mL (corresponding to 1–3 µM), and VEGF levels were 0.4–3 pM and 0.3–3 pM in plasma and interstitial muscle of humans, respectively [92]; therefore, the contribution of target-mediated clearance to the overall clearance of bevacizumab was probably insignificant [93]. Similar hypodermis VEGF level was estimated in HCT116 cell-bearing mice, used in the efficacy evaluation of V-PEG-V hydrogel [75], from reported data of plasma VEGF (50 pg/mL, 1.1 pM) in similar mice [94]. On the other hand, a significant effect of the target-mediated clearance on PK has been reported at a low dose of a mAb (2F8) against epidermal growth factor receptor (EGFR) [95]. While a single IV dose of the mAb (2F8) to cynomolgus monkeys showed similar elimination of plasma concentration at 40 and 20 mg/kg, clearance was 3–4 times faster at 2 mg/kg [95]. These findings suggest that the clearance route was negligible at higher mAb plasma levels but became dominant at low plasma levels. mAb (2F8) plasma concentration was 0.2–67 µg/mL at a 2–20 mg/kg dose [95]. Plasma bevacizumab levels were low (approximately 1 µg/mL) after dosing the V-PEG-V hydrogel [75]. These comparisons might suggest that target-mediated clearance plays an important role under a low bevacizumab concentration in plasma and/or hypodermis. Target-mediated clearance of gradually released bevacizumab from the hydrogel might contribute to its low bioavailability (Table 1, no. 8 and 13).

    The scFv is a fragment of mAb [96] with a MW of about 27 kDa. It consists of variable regions of heavy and light chains connected by a peptide linker [97],[98]. Some studies have reported that scFv maintains the binding specificity and potency of the full IgG antibody [97],[99]. DLX105 (ESBA105) is a humanized scFv (246 amino acid, MW 26.3 kDa) for human TNFα [99]. The smaller size of scFv enhances distribution into tissues compared with full-sized mAbs [100] and might enable greater efficacy of scFv especially for the kidney, gastrointestinal tract, pancreas, and skin. On the other hand, scFv lacks FcRn binding and FcRn-mediated recycling and is rapidly eliminated (t1/2 of 0.5–2 h), primarily through renal clearance [98],[101]. Sustained-release DDS is markedly useful in prolonging the t1/2 of scFv and its duration of action and maintaining its primary benefits, namely its great penetration into tissues. In ProLynx β-eliminative DDS, DLX105 was covalently tethered to hydrogel microspheres and gradually released after SC dosing to rats [96].

    ProLynx DDS uses self-cleaving β-eliminative linkers for hydrogel cross-linking and drug-conjugation [96][99]. In the linker, the hydroxide ion catalyzes proton removal, the linker-carbamate bond is rapidly cleaved by β-elimination, and the drug and a substituted alkene are released. Electron-withdrawing modulator groups control the pKa of the acidic carbon-hydrogen bond on the α-carbon as well as the drug release rate. Microspheres (40 µm) were prepared using 4-armed PEG (tetra-PEG) terminated with amines and N-hydroxysuccinimide glutarates and a microfluidic device [102] (Figure 3d). The microspheres were bicyclononyne-derivatized and conjugated with azido-linker-DLX105 [96]. Two different modulators were used for the system: the linker for tetra-PEG cross-linking had a longer cleavage t1/2 compared with a second linker between the drug and microsphere, so the gels degraded following drug release [96]. The gels released the DLX105 with a short amino-propyl remnant on the N-terminus. t1/2 on in vitro release was 11 h at pH 9.4, corresponding to a t1/2 of about 1100 h at pH 7.4 at 37 °C. After IV and SC dosing of DLX105 normal solution to rats, t1/2β was 4.1 and 3.7 h, respectively (Table 7). The ProLynx gels significantly prolonged t1/2β by about 70 times after SC dosing [96] (Table 7).

    Table 7.  Pharmacokinetic parameters of DLX105.

    Formulation Dosing route Animal Dose (nmol/rat) Cmax (µM) t1/2α (h) t1/2β (h) Relative bioavailability % (vs. IV) Reference

    Normal solution IV rat ~ 86 1.2 0.22 4.1 [96]
    Normal solution SC rat ~ 86 0.12 - 3.7 100 [96]
    ProLynx gel SC rat ~ 57 0.0010 - 270 ~ 60 [96]

    IV: intravenous dosing; SC: subcutaneous dosing; Cmax: maximum concentration; t1/2α: half-life during distribution phase; t1/2β: half-life during terminal elimination phase.

     | Show Table
    DownLoad: CSV

    Drug release with the ProLynx DDS is significantly controlled by cleavage of the linkers, and the drug's loading amount, diffusion in DDS, and interactions with polymers consequently make only a minor contribution to the sustained release (Table 1, no. 1–3). The in vivo t1/2β of the DLX105 derivative was 4.1-fold faster than that of in vitro release [96]. From data for tetra-PEG microspheres loaded with exenatide, in vivo cleavage is about 3-fold faster than in vitro [103]. Faster in vivo cleavage was also observed in circulating PEG-conjugates containing these linkers; accordingly, the in vivo rate enhancements were not attributable to the SC environment [104]. Moreover, the rates of linker cleavage were the same in the mouse, rat, and monkey [104],[105]. On the other hand, catalysis of albumin in β-elimination has been reported for several molecules [106][110]. Lysine (Lys) residues of albumin, including Lys-199 in subdomain IIA, are thought to cause a general-base effect and catalysis [109],[110]. Although albumin concentration is lower in adipose tissue (7.36 mg/mL) than in serum (48.9 mg/mL) [71], a study reported that 0.1 mg/mL BSA catalyzed β-elimination, and 2 mg/mL BSA significantly shortened t1/2β of β-elimination from 10 h to 4 min [107]. The study was conducted in vitro at pH 8.8 [107], and the in vivo catalysis of albumin is unclear; however, this might contribute to the shorter t1/2β of DLX105 from ProLynx gels (Table 1, no. 8).

    The in vivo fate of the gel was evaluated by SC injection of the gel containing no drugs to mice [111] (Table 1, no. 5). After four weeks, the gels were retained in the dosing site, and mild infiltration of inflammatory cells was observed [111]. In another study of repeated SC dosing of placebo gels, slight injection site swelling and moderate SC granulomatous inflammation were observed [103]. These results suggested that APCs may have been involved in the polymer's degradation. The pI of DLX105 has not been published, and the possibility of its interaction with ECM cannot be discussed. However, its small size (MW 26.3 kDa) and short t1/2β of the solution on SC dosing (Table 7) suggest a rapid diffusion of DLX105 in hypodermis (Table 1, no. 6). The released DLX105 likely maintained its folding (Table 1, no. 7), which is supported by the finding that DLX105 with a short amino-propyl remnant showed the same bioactivity for TNFα neutralization as free DLX105 [96] and high biocompatibility with the PEG-based hydrogel.

    Lymphatic absorption of DLX105 was estimated to be about 80% by linear correlation (Figure 2) and considering its size (MW 26.3 kDa, Table 1, no. 9). In contrast, another study on SC dosing of anti-HER2 scFv with a similar size (MW 27 kDa) reported 10.8% absorption through the lymphatic pathway [112]. Moreover, the SC bioavailability of DLX105 solution, DLX105 ProLynx gel, and anti-HER2 scFv was 100%, 60% [96], and 22% [112], respectively. The anti-HER2 study reported that the simulation was significantly affected by the catabolism of scFv by macrophages in the hypodermis [113]; therefore, the higher bioavailability of DLX105 may be due to higher stability. Based on this assumption, the lower bioavailability of DLX105 ProLynx gel compared with DLX105 solution might be due to a more significant degradation (Table 1, no. 10). As discussed in the preceding sections, the gel gradually released the scFv, and the resulting low hypodermis scFv concentration might result in the degradation of a larger portion of macrophages [113] (Table 1, no. 8). Abundant accumulation of APCs around the ProLynx gel was also related to more significant drug degradation [111],[114].

    Intradermal MN is another DDS for SC sustained release of mAb and proteins [115],[116]. The skin is a static barrier that protects internal organs. The stratum corneum of the outermost skin layer is the major hindrance to the penetration of macromolecules into the skin, and only lipophilic small molecules (< 500 Da) can passively diffuse through it. MN arrays consist of multiple protrusions of 25–2000 µm in length, which are long enough to penetrate the stratum corneum [115] (Figure 6). MNs transiently create microscopic aqueous channels and thereby allow mAbs and proteins to bypass the barrier and reach the epidermis and dermis for release. Repair of the microchannels occurs soon after MN removal, avoiding long-term skin damage.

    Figure 6.  Microneedle dosing across the stratum corneum.

    MN technologies can be divided into five different categories: solid, coated, hollow, dissolving, and hydrogel [116] (Figure 6). The solid MN is drug-free and is combined with another drug formulation for passive diffusion (patch, solution, cream, or gel). However, the resulting inaccuracy of dosing and the two-step application process are somewhat inconvenient [115]. For coated MNs, drugs are coated on the needle surface; however, drug loading is very low [115] (Figure 6). Hollow MNs are composed of empty cavity needles (5–70 µm wide) that work together with an external auxiliary device (syringe, pump, gas, or electrical assistance) to infuse drug solution through the needles into the skin (Figure 6). The major advantages of hollow MNs are their delivery of a greater volume (10–100 µL/min) of drug solution and precise control at high doses [115],[116]. Among the disadvantages, the bulky device might reduce convenience, and the drug flow may be potentially disturbed by the clogging of needle openings with skin tissue or compression of dense dermal tissue. Further, the liquid drug formulations are at risk of stability issues [115],[116]. Finally, hollow MNs are rigid and brittle, raising the possibility of the breakage of needles in the skin [117].

    Table 8.  Differences between dissolving and hydrogel microneedles.

    Dissolving microneedles Hydrogel microneedles

    Needle tip components Drug and biodegradable polymer Only hydrophilic polymer, the drug is in a reservoir attached to the needle array
    Drug dose Limited by needle volume Not limited, loaded in the reservoir.
    Needle tips after dosing Dissolve Not dissolved, swell with skin interstitial fluid
    Drug release Rapid Relatively slow, requires time for swelling and diffusion
    Polymer deposition in the skin Have risk No risk, needles are completely withdrawn
    Accidental re-piercing No risk, needles dissolved No risk, needles are hydrated and swelled

     | Show Table
    DownLoad: CSV

    In dissolving MNs (D-MNs), the needle tips are composed of drugs and cheap biodegradable polymers and are manufactured by a simple process (Figure 6, Table 8). The tips dissolve after application to the skin, and drug release is controlled by the polymer's dissolution rate. Application is a simple single step, with no risk of accidental re-piercing and no need for sharps disposal [115],[116]. A disadvantage is the fact that drug loading is limited by the needle volume [117]. Further, the effects of polymer deposition on the skin should be considered. Taking D-MN of poly(methylvinylether/maleic acid) copolymer as an example, its administration resulted in the deposition of at least 5 mg/cm2 of polymer in the skin [118]. While this is unlikely to be a concern for vaccinations, most mAb and proteins are given by repeated dosing [115]. Careful selection of polymers and their compositions is therefore important, with particular regard to the accumulation, distribution, and elimination of polymers, as well as drug stability and manufacturability [117].

    For hydrogel MNs (H-MNs), needles are composed of hydrophilic polymers and contain no drugs; rather, a drug reservoir is attached to the needle array (Figure 6, Table 8). After application to the skin, the needles rapidly absorb skin interstitial fluid and swell to form a hydrogel. This fluid penetrates into the reservoir. The drug dissolves in this fluid and then diffuses through the swollen needles into the skin. Control of the cross-linking density of the needle polymers can be used to adjust the degree of swelling and control drug input rates [115],[116]. Fluid uptake of H-MN was 0.9–2.7 µL/h, which is of the same order as hollow MN [119]. The needles are tough and can be completely withdrawn from the skin without deposition. Since swollen needles cannot be re-inserted, there is no risk of re-insertion and no need for sharps disposal [115],[116].

    Regarding the application of MNs, the skin's viscoelastic property causes force distribution, and some parts of the skin resist deformation [120]. A smaller interspace between needles requires a larger insertion force to pierce the skin [121]. Insertion force depends on the skin's resistance to puncture and its biomechanical properties [120]. Young's modulus of skin increases linearly with age. For children and elderly adults, the elasticity modulus was 70 and 60 N/mm2 (MPa), and skin tension was 21 and 17 N/mm2 (MPa), respectively [122]. The mean ultimate skin deformation before bursting was 75% for newborns and 60% for the elderly. These findings indicate that aging makes the skin thinner, stiffer, less tense, and less flexible [122]. On the other hand, the thickness of the stratum corneum differed among body sites [119],[123]. The mean (SD) thickness of the stratum corneum in 71 human volunteers was 18.3 (4.9) µm at the dorsal aspect of the forearm, 11.0 (2.2) µm at the shoulder, and 14.9 (3.4) µm at the buttocks. These findings were positively correlated with pigmentation and were likely attributable to sun exposure [123]. While individual variation was significant, it was much lower than that among body sites. In adults and children, stratum corneum thickness was independent of age and gender [119],[123]. Consideration of these variations in the skin is important in MN design and for dosing instructions aimed at delivering stable efficacy and minimizing the adverse effects of mAb and proteins.

    Among studies of the sustained release of bevacizumab (anti-VEGF-A IgG1, MW 149 kDa, pI = 8.3) [72], a D-MN has been investigated [124] (Figure 6, Table 8). To prepare the D-MN array, an aqueous solution of polyvinyl alcohol (PVA) (Figure 3e) and bevacizumab was filled onto prefabricated silicone molds, from which the arrays were then removed. The in vitro permeation of bevacizumab across dermatomed neonatal porcine skin on a modified Franz cell apparatus was linearly increased over 0–6 h and reached 9% at 24 h [124]. When this D-MN was applied to rats and held in situ for 24 h, the needles completely dissolved by 24 h. In the PK profile, Cmax was 6 h, less than the effective concentration; peak/trough ratio (Cmax/C168 h) was 3.3; and bioavailability was about 0.8% [124] (Table 9). The peak/trough ratio was similar to that with IV dosing (2.7, Table 9), likely due to the fast release of bevacizumab. The PVA is a hydrophilic polymer and the needles are considered to start to dissolve immediately after insertion. This results in the full available payload of bevacizumab being released from the needle tips as a bolus dose [125] (Table 1, no. 5). The low in vitro permeation of D-MN suggests that only bevacizumab contained in the needles is available for delivery. This is supported by the finding that the in vitro permeated drug amount (9%) was consistent with the needle weight ratio (10%) in the total array [124]. In this preparation method, a part of bevacizumab was contained in the baseplate (Figure 6). In the case of small molecule drugs, a proportion of drugs in the baseplate can permeate through the microchannel; in contrast, the movement of bevacizumab was probably prevented by its large molecular size (Table 1, no. 2). To ensure efficient dosing, the preparation method should be improved such that drug loading is limited to the needles only (Figure 6).

    Table 9.  Pharmacokinetic parameters of bevacizumab after intravenous or transdermal microneedle dosing to rats.

    Formulation Dosing route Animal Dose (mg/head) Cmax (µg/mL) Cmax/C168 h Relative bioavailability % (vs. IV) References

    Normal solution IV rat 2.5 75.6 2.7 [124]
    Dissolving microneedle transdermal rat 5 489 3.3 ~ 0.8 [124]
    Hydrogel microneedle transdermal rat 5 81.2 1.4 ~ 0.2 [124]
    Hydrogel microneedle transdermal rat 10 358 1.7 ~ 0.8 [124]

    IV: intravenous dosing; SC: subcutaneous dosing; Cmax: maximum concentration; Cmax/C168 h: ratio of Cmax to concentration at 168 h.

     | Show Table
    DownLoad: CSV

    In the study described in Section 3.1, H-MN was also investigated for the sustained release of bevacizumab [72],[124] (Figure 6, Table 8). An aqueous solution of methyl vinyl ether/maleic anhydride copolymer (Gantrez S-97, Figure 3F, 20 wt%), PEG 10000 (7.5 wt%), and Na2CO3 (3 wt%) was dispensed into the molds and cross-linked. Lyophilized drug reservoirs were separately prepared by dissolving bevacizumab, gelatin, mannitol, NaCl, and sucrose (20 wt%, 10 wt%, 40 wt%, 10 wt%, and 1 wt%, respectively) in water, and the solution was cast into radius cylindrical molds with lyophilization. After lyophilization, the biological activity of bevacizumab was maintained [124]. The in vitro permeation of bevacizumab across porcine skin held in a Franz cell was linearly increased over 0–6 h and reached 24%–28% [124]. After applying the MN to the backs of shaved rats and holding it in situ for 24 h, needles were clearly visible and showed extensive swelling with full dissolution of the reservoirs [124]. Serum concentration after IV dosing of 2.5 mg bevacizumab solution to rats showed a peak/trough ratio (Cmax/C168 h) of 2.7 (Table 9). Dosing of H-MN containing 5 or 10 mg of drugs showed a smaller Cmax, at 48 h, and a peak/trough ratio of 1.4 and 1.7, respectively [124] (Table 9). H-MN showed prolonged Tmax and decreased Cmax compared with D-MN, and the peak/trough ratio was reduced by approximately 50% compared with D-MN and IV dosing, indicating the sustained absorption of bevacizumab. Serum concentration was lower than the effective concentration, and bioavailability was calculated at about 0.2% and 0.8%, respectively.

    The steps for drug release from the H-MN include initial swelling of the array, sufficient diffusion of fluid into the reservoir, dissolution of bevacizumab in the reservoir, and diffusion through the needle into the hypodermis [124] (Table 1, no. 2 and 5). Consistent with these slow steps, H-MN induced sustained in vitro permeation and a slower Tmax and lower Cmax. The cationic bevacizumab (pI = 8.3) probably interacted with the carboxyl group on the hydrogel polymer of Gantrez S-97 (Table 1, no. 2, Figure 3f), which contributed to the in vitro gradual (6 h) and limited permeation (24%–28%) and the in vivo PK profiles (Tmax, Cmax, peak/trough ratio and the low bioavailability) [72],[126]. Since the H-MN was removed at 24 h post-dosing, the Tmax of 24–48 h indicated a delay in systemic absorption and the possibility that bevacizumab was retained within the skin or lymphatic tissues prior to entry into blood [124]. In turn, these might have been due to the binding of cationic bevacizumab (pI = 8.3) with anionic ECM and slow diffusion due to the large size (MW 149 kDa) [72] (Table 1, no. 6). The possibility that the bevacizumab unfolded was low, owing to the mild MN preparation process and use of lyophilized drug reservoirs and hydrogel needles (Table 1, no. 7).

    The lymphatic absorption of bevacizumab (Table 1, no. 9, Figure 2) and its degradation data with lymph node homogenate [28] indicated that APCs in the lymph nodes contributed to degradation (as described in Section 2.3, Table 1, no. 10). Moreover, APCs are recruited around the MN in hypodermis. In another study, MN comprising hyaluronic acid, collagen, and ovalbumin was applied to the skin of mice [127]. Compared with non-treated skin, approximately 10 times more macrophages/monocytes and neutrophils infiltrated into the skin around the MN on the first day after application [127]. As discussed in Section 2.3 above, one possible reason for the sustained absorption and low bioavailability is the degradation of bevacizumab by APCs in hypodermis and/or lymph nodes (Table 1, no. 10). As was also seen with SC V-PEG-V hydrogel dosing, the slow absorption and low bioavailability of H-MN might also be partly attributable to a decrease in FcRn-mediated recycling, increase in FcγR-mediated degradation due to a low drug concentration in hypodermis, and target-mediated clearance (Table 1, no. 8 and 11–13).

    Given that the open pores generated by MN treatment might produce complications such as irritation and infection, consideration of pore closure is an important aspect in the selection and design of DDS [119] (Table 1, no. 4). After application of a MN to rats, the defensive properties of skin were re-established within 3–4 h. Complete pore closure was noted by 15 h, but a delay in closure for up to 72 h was also observed due to the incorporation of occlusive materials [128]. Among animals, rat skin has a similar thickness of the stratum corneum, epidermis, and whole skin to human skin, compared with mouse and pig skin [129]. Rat skin significantly differs from human skin in its intercellular lipid composition of the stratum corneum and in its lower corneocyte surface [130] (Table 1, no. 14). Moreover, rat skin is reported to recover more rapidly than human skin after mild alterations in superficial epidermis; the t1/2 for the first phase of barrier regeneration for rat and human skin was 19 and 40 h, and recovery of transepidermal water loss was 10 and 30 days, respectively [131] (Table 1, no. 15). A clinical study of MNs demonstrated that microchannel repair and resealing were apparent at 8–24 h post-application, whereas disruption triggered by a 25 G hypodermic needle was significantly notable even after 24 h [132]. In 12 clinical trials, adverse events (erythema, swelling, and bleeding) were related to the MN application site. These were mild and generally disappeared a few hours after application [117]. A clinical study examined the repeated application of an H-MN array to the upper arm of human volunteers every day for a period of 5 days. Neither skin reactions nor disruption of the skin barrier function were prolonged [133]. These results suggest that MN research has overcome the barrier of species differences and confirmed its applicability to human medication. In addition to increasing the patch size of the safe MN, the design and manufacturing process will be improved to provide higher loading and increase the efficacy of mAb and proteins.

    Antibodies and proteins are critically important modalities. SC sustained-release DDS has the potential to reduce the dosing frequency of these drugs and contribute to the delivery of high value to patients. This article exemplifies several DDS based around polymeric particles, conjugated hydrogel SC injections, and intradermal dissolving/hydrogel MN with four proteins/antibodies (hGH, IFNα-2b, bevacizumab, and scFv). This review suggests that polymers/excipients and manufacturing processes for formulations should be selected with consideration to interactions with drugs, suppression of drug solubility, and maintenance of drug structure/bioactivity. Other important determinants of safe dosing, sustained release, and drug bioavailability include optimization of the size, strength, morphology, and materials of the DDS, polymer deposition, thickness/viscoelastic properties, and regeneration of skin. In the body, many biophysical parameters affect drug PK and the adverse events of DDS, including the small pore size of blood vessels, negative charge of the ECM, concentration of proteins in the hypodermis, recycling of receptors, infiltration of APCs, and degradation by APCs and protease. Although the contribution of each parameter varies among drugs and types of DDS, consideration of these possible factors will lead to the rapid and reliable development of DDS and many long-acting antibody/protein products for patients.

    The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.



    [1] S. Z. Goldhaber, H. Bounameaux, Pulmonary embolism and deep vein thrombosis, Lancet (London, England), 379 (2012), 1835–1846. https://doi.org/10.1016/S0140-6736(11)61904-1 doi: 10.1016/S0140-6736(11)61904-1
    [2] G. E. Raskob, P. Angchaisuksiri, A. N. Blanco, H. Buller, A. Gallus, B. J. Hunt, et al., Thrombosis: a major contributor to global disease burden, Arterioscler. Thromb. Vasc. Biol., 34 (2014), 2363–2371. https://doi.org/10.1161/ATVBAHA.114.304488 doi: 10.1161/ATVBAHA.114.304488
    [3] D. Voci, U. Fedeli, I. T. Farmakis, L. Hobohm, K. Keller, L. Valerio, et al., Deaths related to pulmonary embolism and cardiovascular events before and during the 2020 COVID-19 pandemic: An epidemiological analysis of data from an Italian high-risk area, Thromb. Res., 212 (2022), 44–50. https://doi.org/10.1016/j.thromres.2022.02.008 doi: 10.1016/j.thromres.2022.02.008
    [4] I. Katsoularis, O. Fonseca-Rodríguez, P. Farrington, H. Jerndal, E. H. Lundevaller, M. Sund, et al., Risks of deep vein thrombosis, pulmonary embolism, and bleeding after covid-19: nationwide self-controlled cases series and matched cohort study, BMJ, 377 (2022). https://doi.org/10.1136/bmj-2021-069590 doi: 10.1136/bmj-2021-069590
    [5] T. N. Nguyen, M. M. Qureshi, P. Klein, H. Yamagami, M. Abdalkader, R. Mikulik, et al., Global impact of the COVID-19 pandemic on cerebral venous thrombosis and mortality, J. Stroke, 24 (2022), 256–265. https://doi.org/10.5853/jos.2022.00752 doi: 10.5853/jos.2022.00752
    [6] E. Berntorp, K. Fischer, D. P. Hart, M. E. Mancuso, D. Stephensen, A. D. Shapiro, et al., Haemophilia, Nat. Rev. Dis. Prim., 7 (2021), 45. https://doi.org/10.1038/s41572-021-00278-x doi: 10.1038/s41572-021-00278-x
    [7] K. G. Link, M. T. Stobb, M. G. Sorrells, M. Bortot, K. Ruegg, M. J. Manco‐Johnson, et al., A mathematical model of coagulation under flow identifies factor V as a modifier of thrombin generation in hemophilia A, J. Thromb. Haemost., 18 (2020), 306–317. https://doi.org/10.1111/jth.14653 doi: 10.1111/jth.14653
    [8] F. W. G. Leebeek, W. Miesbach, Gene therapy for hemophilia: a review on clinical benefit, limitations, and remaining issues, Blood, 138 (2021), 923–931. https://doi.org/10.1182/blood.2019003777 doi: 10.1182/blood.2019003777
    [9] S. S. G. Halfmann, N. Evangelatos, P. Schröder-Bäck, A. Brand, European healthcare systems readiness to shift from 'One-Size Fits All' to personalized medicine, Per. Med., 14 (2017), 63–74. https://doi.org/10.2217/pme-2016-0061 doi: 10.2217/pme-2016-0061
    [10] T. Behl, I. Kaur, A. Sehgal, S. Singh, A. Albarrati, M. Albratty, et al., The road to precision medicine: Eliminating the "One Size Fits All" approach in Alzheimer's disease, Biomed. Pharmacother., 153 (2022), 113337. https://doi.org/10.1016/j.biopha.2022.113337 doi: 10.1016/j.biopha.2022.113337
    [11] N. M. Hamdy, E. B. Basalious, M. G. El-Sisi, M. Nasr, A. M. Kabel, E. S. Nossier, et al., Advancements in current one-size-fits-all therapies compared to future treatment innovations for better improved chemotherapeutic outcomes: a step-toward personalized medicine, Curr. Med. Res. Opin., 40 (2024), 1–19. https://doi.org/10.1080/03007995.2024.2416985 doi: 10.1080/03007995.2024.2416985
    [12] G. Di Minno, E. Tremoli, Tailoring of medical treatment: hemostasis and thrombosis towards precision medicine, Haematologica, 102 (2017), 411–418. https://doi.org/10.3324/haematol.2016.156000 doi: 10.3324/haematol.2016.156000
    [13] D. L. Ornstein, Chapter 41 - Personalized medicine for disorders of hemostasis and thrombosis, in Diagnostic Molecular Pathology (eds. W. B. Coleman and G. J. Tsongalis), Academic Press, (2024), 643–653. https://doi.org/10.1016/B978-0-12-822824-1.00006-7
    [14] S. Nagalla, P. F. Bray, Personalized medicine in thrombosis: back to the future, Blood, 127 (2016), 2665–2671. https://doi.org/10.1182/blood-2015-11-634832 doi: 10.1182/blood-2015-11-634832
    [15] R. J. S. Preston, J. M. O'Sullivan, Personalized approaches to the treatment of hemostatic disorders, Semin. Thromb. Hemostasis, 47 (2021), 117–119. https://doi.org/10.1055/s-0041-1723800 doi: 10.1055/s-0041-1723800
    [16] H. Al‐Samkari, W. Eng, A precision medicine approach to hereditary hemorrhagic telangiectasia and complex vascular anomalies, J. Thromb. Haemost., 20 (2022), 1077–1088. https://doi.org/10.1111/jth.15715 doi: 10.1111/jth.15715
    [17] X. Delavenne, E. Ollier, A. Lienhart, Y. Dargaud, A new paradigm for personalized prophylaxis for patients with severe haemophilia A, Haemophilia, 26 (2020), 228–235. https://doi.org/10.1111/hae.13935 doi: 10.1111/hae.13935
    [18] L. H. Bukkems, L. L. F. G. Valke, W. Barteling, B. A. P. Laros-van Gorkom, N. M. A. Blijlevens, M. H. Cnossen, et al., Combining factor Ⅷ levels and thrombin/plasmin generation: A population pharmacokinetic-pharmacodynamic model for patients with haemophilia A, Br. J. Clin. Pharmacol., 88 (2022), 2757–2768. https://doi.org/10.1111/bcp.15185 doi: 10.1111/bcp.15185
    [19] N. Mackman, W. Bergmeier, G. A. Stouffer, J. I. Weitz, Therapeutic strategies for thrombosis: new targets and approaches, Nat. Rev. Drug. Discov., 19 (2020), 333–352. https://doi.org/10.1038/s41573-020-0061-0 doi: 10.1038/s41573-020-0061-0
    [20] P. S. Wells, R. Ihaddadene, A. Reilly, M. A. Forgie, Diagnosis of venous thromboembolism: 20 years of progress, Ann. Intern. Med., 168 (2018), 131–140. https://doi.org/10.7326/M17-0291 doi: 10.7326/M17-0291
    [21] F. Khan, T. Tritschler, S. R. Kahn, M. A. Rodger, Venous thromboembolism, Lancet (London, England), 398 (2021), 64–77. https://doi.org/10.1016/S0140-6736(20)32658-1 doi: 10.1016/S0140-6736(20)32658-1
    [22] T. D. Martins, S. D. Martins, S. Montalvão, M. Al Bannoud, G. Y. Ottaiano, L. Q. Silva, et al., Combining artificial neural networks and hematological data to diagnose Covid-19 infection in Brazilian population, Neural Comput. Appl., 36 (2024), 4387–4399. https://doi.org/10.1007/s00521-023-09312-3 doi: 10.1007/s00521-023-09312-3
    [23] T. D. Martins, R. Maciel-Filho, S. A. L. Montalvão, G. S. S. Gois, M. Al Bannoud, G. Y. Ottaiano, et al., Predicting mortality of cancer patients using artificial intelligence, patient data and blood tests, Neural Comput. Appl., 36 (2024), 15599–15616. https://doi.org/10.1007/s00521-024-09915-4 doi: 10.1007/s00521-024-09915-4
    [24] F. W. G. Leebeek, New developments in diagnosis and management of acquired hemophilia and acquired von willebrand syndrome, HemaSphere, 5 (2021). https://doi.org/10.1097/HS9.0000000000000586 doi: 10.1097/HS9.0000000000000586
    [25] F. Peyvandi, G. Kenet, I. Pekrul, R. K. Pruthi, P. Ramge, M. Spannagl, Laboratory testing in hemophilia: Impact of factor and non‐factor replacement therapy on coagulation assays, J. Thromb. Haemost., 18 (2020), 1242–1255. https://doi.org/10.1111/jth.14784 doi: 10.1111/jth.14784
    [26] B. Pezeshkpoor, J. Oldenburg, A. Pavlova, Insights into the molecular genetic of hemophilia A and hemophilia B: The relevance of genetic testing in routine clinical practice, Hamostaseologie, 42 (2022), 390–399. https://doi.org/10.1055/a-1945-9429 doi: 10.1055/a-1945-9429
    [27] A. H. Kristoffersen, E. Ajzner, D. Rogic, E. Y. Sozmen, P. Carraro, A. P. Faria, et al., Is D-dimer used according to clinical algorithms in the diagnostic work-up of patients with suspicion of venous thromboembolism? A study in six European countries, Thromb. Res., 142 (2016), 1–7. https://doi.org/10.1016/j.thromres.2016.04.001 doi: 10.1016/j.thromres.2016.04.001
    [28] M. Kafeza, J. Shalhoub, N. Salooja, L. Bingham, K. Spagou, A. H. Davies, A systematic review of clinical prediction scores for deep vein thrombosis, Phlebology, 32 (2017), 516–531. https://doi.org/10.1177/0268355516678729 doi: 10.1177/0268355516678729
    [29] M. T. Greene, A. C. Spyropoulos, V. Chopra, P. J. Grant, S. Kaatz, S. J. Bernstein, et al., Validation of risk assessment models of venous thromboembolism in hospitalized medical patients, Am. J. Med., 129 (2016), 1001.e9–1001.e18. https://doi.org/10.1016/j.amjmed.2016.03.031 doi: 10.1016/j.amjmed.2016.03.031
    [30] P. C. Silveira, I. K. Ip, S. Z. Goldhaber, G. Piazza, C. B. Benson, R. Khorasani, Performance of wells score for deep vein thrombosis in the inpatient setting, JAMA Intern. Med., 175 (2015), 1112–1117. https://doi.org/10.1001/jamainternmed.2015.1687 doi: 10.1001/jamainternmed.2015.1687
    [31] M. Kafeza, J. Shalhoub, N. Salooja, L. Bingham, K. Spagou, A. H. Davies, A systematic review of clinical prediction scores for deep vein thrombosis, Phlebology, 32 (2016), 516–531. https://doi.org/10.1177/0268355516678729 doi: 10.1177/0268355516678729
    [32] I. Nichele, A. Tosetto, Scoring Systems for Estimating the Risk of Recurrent Venous Thromboembolism, Semin. Thromb. Hemost., 43 (2017), 493–499. https://doi.org/10.1055/s-0037-1602662 doi: 10.1055/s-0037-1602662
    [33] A. Muñoz, C. Ay, E. Grilz, S. López, C. Font, V. Pachón, et al., A clinical-genetic risk score for predicting cancer-associated venous thromboembolism: A development and validation study involving two independent prospective cohorts, J. Clin. Oncol., 41 (2023), 2911–2925. https://doi.org/10.1200/JCO.22.00255 doi: 10.1200/JCO.22.00255
    [34] F. Rodeghiero, A. Tosetto, T. Abshire, D. M. Arnold, B. Coller, P. James, et al., ISTH/SSC bleeding assessment tool: a standardized questionnaire and a proposal for a new bleeding score for inherited bleeding disorders, J. Thromb. Haemost., 8 (2010), 2063–2065. https://doi.org/10.1111/j.1538-7836.2010.03975.x doi: 10.1111/j.1538-7836.2010.03975.x
    [35] M. Borhany, N. Fatima, M. Abid, T. Shamsi, M. Othman, Application of the ISTH bleeding score in hemophilia, Transfus. Apher. Sci., 57 (2018), 556–560. https://doi.org/10.1016/j.transci.2018.06.003 doi: 10.1016/j.transci.2018.06.003
    [36] M. Khalifa, M. Albadawy, Artificial intelligence for clinical prediction: Exploring key domains and essential functions, Comput. Methods Programs Biomed. Updat., 5 (2024), 100148. https://doi.org/10.1016/j.cmpbup.2024.100148 doi: 10.1016/j.cmpbup.2024.100148
    [37] T. H. Tan, C. C. Hsu, C. J. Chen, S. L. Hsu, T. L. Liu, H. J. Lin, et al., Predicting outcomes in older ED patients with influenza in real time using a big data-driven and machine learning approach to the hospital information system, BMC Geriatr., 21 (2021), 280. https://doi.org/10.1186/s12877-021-02229-3 doi: 10.1186/s12877-021-02229-3
    [38] C. Guan, F. Ma, S. Chang, J. Zhang, Interpretable machine learning models for predicting venous thromboembolism in the intensive care unit: an analysis based on data from 207 centers, Crit. Care, 27 (2023), 406. https://doi.org/10.1186/s13054-023-04683-4 doi: 10.1186/s13054-023-04683-4
    [39] T. D. Martins, J. M. Annichino-Bizzacchi, A. V. C. Romano, R. Maciel Filho, Artificial neural networks for prediction of recurrent venous thromboembolism, Int. J. Med. Inform., 141 (2020), 104221. https://doi.org/10.1016/j.ijmedinf.2020.104221 doi: 10.1016/j.ijmedinf.2020.104221
    [40] I. Pabinger, C. Ay, Biomarkers and venous thromboembolism, Arter. Thromb. Vasc. Biol., 29 (2009), 332–336. https://doi.org/10.1161/ATVBAHA.108.182188 doi: 10.1161/ATVBAHA.108.182188
    [41] I. Pabinger, J. Thaler, C. Ay, Biomarkers for prediction of venous thromboembolism in cancer, Blood, 122 (2013), 2011–2018. https://doi.org/10.1182/blood-2013-04-460147 doi: 10.1182/blood-2013-04-460147
    [42] F. Galeano-Valle, L. Ordieres-Ortega, C. M. Oblitas, J. del-Toro-Cervera, L. Alvarez-Sala-Walther, P. Demelo-Rodríguez, Inflammatory biomarkers in the short-term prognosis of venous thromboembolism: A narrative review, Int. J. Mol. Sci., 22 (2021), 2627. https://doi.org/10.3390/ijms22052627 doi: 10.3390/ijms22052627
    [43] B. Jacobs, A. Obi, T. Wakefield, Diagnostic biomarkers in venous thromboembolic disease, J. Vasc. Surg. Venous Lymphat. Disord., 4 (2016), 508–517. https://doi.org/10.1016/j.jvsv.2016.02.005 doi: 10.1016/j.jvsv.2016.02.005
    [44] M. A. Bannoud, B. P. Gomes, M. C. de S. P. Abdalla, M. V Freire, K. Andreola, T. D. Martins, et al., Mathematical modeling of drying kinetics of ground Açaí (Euterpe oleracea) kernel using artificial neural networks, Chem. Pap., 78 (2024), 1033–1054. https://doi.org/10.1007/s11696-023-03142-2 doi: 10.1007/s11696-023-03142-2
    [45] M. A. Bannoud, T. D. Martins, B. F. dos Santos, Control of a closed dry grinding circuit with ball mills using predictive control based on neural networks, Digit. Chem. Eng., 5 (2022), 100064. https://doi.org/10.1016/j.dche.2022.100064 doi: 10.1016/j.dche.2022.100064
    [46] J. Berg, K. Nyström, Data-driven discovery of PDEs in complex datasets, J. Comput. Phys., 384 (2019), 239–252. https://doi.org/10.1016/j.jcp.2019.01.036 doi: 10.1016/j.jcp.2019.01.036
    [47] L. Burzawa, L. Li, X. Wang, A. Buganza-Tepole, D. M. Umulis, Acceleration of PDE-based biological simulation through the development of neural network metamodels, Curr. Pathobiol. Rep., 8 (2020), 121–131. https://doi.org/10.1007/s40139-020-00216-8 doi: 10.1007/s40139-020-00216-8
    [48] M. A. Bannoud, C. A. M. da Silva, T. D. Martins, Applications of metaheuristic optimization algorithms in model predictive control for chemical engineering processes: A systematic review, Annu. Rev. Control., 58 (2024), 100973. https://doi.org/10.1016/j.arcontrol.2024.100973 doi: 10.1016/j.arcontrol.2024.100973
    [49] M. A. Bannoud, P. H. N. Ferreira, R. R. de Andrade, C. A. M. da Silva, Control of an integrated first and second-generation continuous alcoholic fermentation process with cell recycling using model predictive control, Chem. Eng. Commun., (2024), 1–24. https://doi.org/10.1080/00986445.2024.2417901 doi: 10.1080/00986445.2024.2417901
    [50] Y. Zhao, C. Li, X. Liu, R. Qian, R. Song, X. Chen, Patient-specific seizure prediction via adder network and supervised contrastive learning, IEEE Trans. Neural Syst. Rehabil. Eng., 30 (2022), 1536–1547. https://doi.org/10.1109/TNSRE.2022.3180155 doi: 10.1109/TNSRE.2022.3180155
    [51] K. Leiderman, S. S. Sindi, D. M. Monroe, A. L. Fogelson, K. B. Neeves, The art and science of building a computational model to understand hemostasis, Semin. Thromb. Hemost., 47 (2021), 129–138. https://doi.org/10.1055/s-0041-1722861 doi: 10.1055/s-0041-1722861
    [52] N. Ratto, A. Bouchnita, P. Chelle, M. Marion, M. Panteleev, D. Nechipurenko, et al., Patient-specific modelling of blood coagulation, Bull. Math. Biol., 83 (2021), 50. https://doi.org/10.1007/s11538-021-00890-8 doi: 10.1007/s11538-021-00890-8
    [53] R. Burghaus, K. Coboeken, T. Gaub, L. Kuepfer, A. Sensse, H.-U. Siegmund, et al., Evaluation of the efficacy and safety of rivaroxaban using a computer model for blood coagulation, PLoS One, 6 (2011), e17626. https://doi.org/10.1371/journal.pone.0017626 doi: 10.1371/journal.pone.0017626
    [54] C. Watson, H. Saaid, V. Vedula, J. C. Cardenas, P. K. Henke, F. Nicoud, et al., Venous thromboembolism: Review of clinical challenges, biology, assessment, treatment, and modeling, Ann. Biomed. Eng., 52 (2024), 467–486. https://doi.org/10.1007/s10439-023-03390-z doi: 10.1007/s10439-023-03390-z
    [55] N. N. Ramli, S. Iberahim, N. H. M. Noor, Z. Zulkafli, T. M. T. M. Shihabuddin, M. H. Din, et al., Haemostasis and inflammatory parameters as potential diagnostic biomarkers for VTE in trauma-immobilized patients, Diagnostics (Basel), 13 (2023), 150. https://doi.org/10.3390/diagnostics13010150 doi: 10.3390/diagnostics13010150
    [56] L. G. R. Ferreira, R. C. Figueiredo, M. das Graças Carvalho, D. R. A. Rios, Thrombin generation assay as a biomarker of cardiovascular outcomes and mortality: A narrative review, Thromb. Res., 220 (2022), 107–115. https://doi.org/10.1016/j.thromres.2022.10.007 doi: 10.1016/j.thromres.2022.10.007
    [57] M. S. Edvardsen, K. Hindberg, E. S. Hansen, V. M. Morelli, T. Ueland, P. Aukrust, et al., Plasma levels of von Willebrand factor and future risk of incident venous thromboembolism, Blood Adv., 5 (2021), 224–232. https://doi.org/10.1182/bloodadvances.2020003135 doi: 10.1182/bloodadvances.2020003135
    [58] L. Anghel, R. Sascău, R. Radu, C. Stătescu, From classical laboratory parameters to novel biomarkers for the diagnosis of venous thrombosis, Int. J. Mol. Sci., 21 (2020), 1920. https://doi.org/10.3390/ijms21061920 doi: 10.3390/ijms21061920
    [59] H. Y. Lim, C. O'Malley, G. Donnan, H. Nandurkar, P. Ho, A review of global coagulation assays — Is there a role in thrombosis risk prediction?, Thromb. Res., 179 (2019), 45–55. https://doi.org/10.1016/j.thromres.2019.04.033 doi: 10.1016/j.thromres.2019.04.033
    [60] H. Hou, Z. Ge, P. Ying, J. Dai, D. Shi, Z. Xu, et al., Biomarkers of deep venous thrombosis, J. Thromb. Thrombolysis., 34 (2012), 335–346. https://doi.org/10.1007/s11239-012-0721-y doi: 10.1007/s11239-012-0721-y
    [61] D. Moher, A. Liberati, J. Tetzlaff, D. G. Altman, Preferred reporting items for systematic reviews and meta-analyses: the PRISMA statement, PLoS Med., 6 (2009), e1000097. https://doi.org/10.1371/journal.pmed.1000097 doi: 10.1371/journal.pmed.1000097
    [62] M. J. Page, J. E. McKenzie, P. M. Bossuyt, I. Boutron, T. C. Hoffmann, C. D. Mulrow, et al., The PRISMA 2020 statement: an updated guideline for reporting systematic reviews, BMJ, 372 (2021), n71. https://doi.org/10.1136/bmj.n71 doi: 10.1136/bmj.n71
    [63] E. N. Sorensen, G. W. Burgreen, W. R. Wagner, J. F. Antaki, Computational simulation of platelet deposition and activation: I. Model development and properties, Ann. Biomed. Eng., 27 (1999), 436–448. https://doi.org/10.1114/1.200 doi: 10.1114/1.200
    [64] K. Boryczko, W. Dzwinel, D. A. Yuen, Modeling fibrin aggregation in blood flow with discrete-particles, Comput. Methods Programs Biomed., 75 (2004), 181–194. https://doi.org/10.1016/j.cmpb.2004.02.001 doi: 10.1016/j.cmpb.2004.02.001
    [65] I. V. Pivkin, P. D. Richardson, G. Karniadakis, Blood flow velocity effects and role of activation delay time on growth and form of platelet thrombi, Proc. Natl. Acad. Sci. USA, 103 (2006), 17164–17169. https://doi.org/10.1073/pnas.0608546103 doi: 10.1073/pnas.0608546103
    [66] Z. Xu, N. Chen, M. M. Kamocka, E. D. Rosen, M. Alber, A multiscale model of thrombus development, J. R. Soc. Interface, 5 (2008), 705–722. https://doi.org/10.1098/rsif.2007.1202 doi: 10.1098/rsif.2007.1202
    [67] Z. Xu, N. Chen, S. C. Shadden, J. E. Marsden, M. M. Kamocka, E. D. Rosen, et al., Study of blood flow impact on growth of thrombi using a multiscale model, Soft. Matter, 5 (2009), 769–779. https://doi.org/10.1039/B812429A doi: 10.1039/B812429A
    [68] Z. Xu, J. Lioi, J. Mu, M. M. Kamocka, X. Liu, D. Z. Chen, et al., A multiscale model of venous thrombus formation with surface-mediated control of blood coagulation cascade, Biophys. J., 98 (2010), 1723–1732. https://doi.org/10.1016/j.bpj.2009.12.4331 doi: 10.1016/j.bpj.2009.12.4331
    [69] A. M. Shibeko, E. S. Lobanova, M. A. Panteleev, F. I. Ataullakhanov, Blood flow controls coagulation onset via the positive feedback of factor Ⅶ activation by factor Xa, BMC Syst. Biol., 4 (2010), 5. https://doi.org/10.1186/1752-0509-4-5 doi: 10.1186/1752-0509-4-5
    [70] S. W. Jordan, E. L. Chaikof, Simulated surface-induced thrombin generation in a flow field, Biophys. J., 101 (2011), 276–286. https://doi.org/10.1016/j.bpj.2011.05.056 doi: 10.1016/j.bpj.2011.05.056
    [71] K. Leiderman, A. L. Fogelson, Grow with the flow: a spatial-temporal model of platelet deposition and blood coagulation under flow, Math. Med. Biol., 28 (2011), 47–84. https://doi.org/10.1093/imammb/dqq005 doi: 10.1093/imammb/dqq005
    [72] A. L. Fogelson, Y. H. Hussain, K. Leiderman, Blood clot formation under flow: The importance of factor XI depends strongly on platelet count, Biophys. J., 102 (2012), 10–18. https://doi.org/10.1016/j.bpj.2011.10.048 doi: 10.1016/j.bpj.2011.10.048
    [73] K. Leiderman, A. L. Fogelson, The influence of hindered transport on the development of platelet thrombi under flow, Bull. Math. Biol., 75 (2013), 1255–1283. https://doi.org/10.1007/s11538-012-9784-3 doi: 10.1007/s11538-012-9784-3
    [74] A. Tosenberger, F. Ataullakhanov, N. Bessonov, M. Panteleev, A. Tokarev, V. Volpert, Modelling of thrombus growth in flow with a DPD-PDE method, J. Theor. Biol., 337 (2013), 30–41. https://doi.org/10.1016/j.jtbi.2013.07.023 doi: 10.1016/j.jtbi.2013.07.023
    [75] A. Sequeira, T. Bodnár, Blood coagulation simulations using a viscoelastic model, Math. Model. Nat. Phenom., 9 (2014), 34–45. https://doi.org/10.1051/mmnp/20149604 doi: 10.1051/mmnp/20149604
    [76] A. Tosenberger, N. Bessonov, V. Volpert, Influence of fibrinogen deficiency on clot formation in flow by hybrid model, Math. Model. Nat. Phenom., 10 (2015), 36–47. https://doi.org/10.1051/mmnp/201510102 doi: 10.1051/mmnp/201510102
    [77] O. S. Rukhlenko, O. A. Dudchenko, K. E. Zlobina, G. T. Guria, Mathematical modeling of intravascular blood coagulation under wall shear stress, PLoS One, 10 (2015), e0134028. https://doi.org/10.1371/journal.pone.0134028 doi: 10.1371/journal.pone.0134028
    [78] J. Pavlova, A. Fasano, J. Janela, A. Sequeira, Numerical validation of a synthetic cell-based model of blood coagulation, J. Theor. Biol., 380 (2015), 367–379. https://doi.org/10.1016/j.jtbi.2015.06.004 doi: 10.1016/j.jtbi.2015.06.004
    [79] A. Piebalgs, X. Y. Xu, Towards a multi-physics modelling framework for thrombolysis under the influence of blood flow, J. R. Soc. Interface, 12 (2015), 20150949. https://doi.org/10.1098/rsif.2015.0949 doi: 10.1098/rsif.2015.0949
    [80] Z. Li, A. Yazdani, A. Tartakovsky, G. E. Karniadakis, Transport dissipative particle dynamics model for mesoscopic advection-diffusion-reaction problems, J. Chem. Phys., 143 (2015), 14101. https://doi.org/10.1063/1.4923254 doi: 10.1063/1.4923254
    [81] A. Bouchnita, K. Bouzaachane, T. Galochkina, P. Kurbatova, P. Nony, V. Volpert, An individualized blood coagulation model to predict INR therapeutic range during warfarin treatment, Math. Model. Nat. Phenom., 11 (2016), 28–44. https://doi.org/10.1051/mmnp/201611603 doi: 10.1051/mmnp/201611603
    [82] J. H. Seo, T. Abd, R. T. George, R. Mittal, A coupled chemo-fluidic computational model for thrombogenesis in infarcted left ventricles, Am. J. Physiol. Heart Circ. Physiol., 310 (2016), H1567–82. https://doi.org/10.1152/ajpheart.00855.2015 doi: 10.1152/ajpheart.00855.2015
    [83] M. N. Ngoepe, Y. Ventikos, Computational modelling of clot development in patient-specific cerebral aneurysm cases, J. Thromb. Haemost., 14 (2016), 262–272. https://doi.org/10.1111/jth.13220 doi: 10.1111/jth.13220
    [84] A. Bouchnita, T. Galochkina, V. Volpert, Influence of antithrombin on the regimes of blood coagulation: Insights from the mathematical model, Acta Biotheor., 64 (2016), 327–342. https://doi.org/10.1007/s10441-016-9291-2 doi: 10.1007/s10441-016-9291-2
    [85] E. V. Dydek, E. L. Chaikof, Simulated thrombin generation in the presence of surface-bound heparin and circulating tissue factor, Ann. Biomed. Eng., 44 (2016), 1072–1084. https://doi.org/10.1007/s10439-015-1377-5 doi: 10.1007/s10439-015-1377-5
    [86] J. Pavlova, A. Fasano, A. Sequeira, Numerical simulations of a reduced model for blood coagulation, Zeitschrift für Angew. Math. und Phys., 67 (2016), 28. https://doi.org/10.1007/s00033-015-0610-2 doi: 10.1007/s00033-015-0610-2
    [87] A. Tosenberger, F. Ataullakhanov, N. Bessonov, M. Panteleev, A. Tokarev, V. Volpert, Modelling of platelet-fibrin clot formation in flow with a DPD-PDE method, J. Math. Biol., 72 (2016), 649–681. https://doi.org/10.1007/s00285-015-0891-2 doi: 10.1007/s00285-015-0891-2
    [88] V. Govindarajan, V. Rakesh, J. Reifman, A. Y. Mitrophanov, Computational study of thrombus formation and clotting factor effects under venous flow conditions, Biophys. J., 110 (2016), 1869–1885. https://doi.org/10.1016/j.bpj.2016.03.010 doi: 10.1016/j.bpj.2016.03.010
    [89] C. Ou, W. Huang, M. M.-F. Yuen, A computational model based on fibrin accumulation for the prediction of stasis thrombosis following flow-diverting treatment in cerebral aneurysms, Med. Biol. Eng. Comput., 55 (2017), 89–99. https://doi.org/10.1007/s11517-016-1501-1 doi: 10.1007/s11517-016-1501-1
    [90] A. Yazdani, H. Li, J. D. Humphrey, G. E. Karniadakis, A general shear-dependent model for thrombus formation, PLoS Comput. Biol., 13 (2017), e1005291. https://doi.org/10.1371/journal.pcbi.1005291 doi: 10.1371/journal.pcbi.1005291
    [91] H. Hosseinzadegan, D. K. Tafti, Prediction of thrombus growth: Effect of stenosis and reynolds number, Cardiovasc. Eng. Technol., 8 (2017), 164–181. https://doi.org/10.1007/s13239-017-0304-3 doi: 10.1007/s13239-017-0304-3
    [92] L. M. Haynes, T. Orfeo, K. G. Mann, S. J. Everse, K. E. Brummel-Ziedins, Probing the dynamics of clot-bound thrombin at venous shear rates, Biophys. J., 112 (2017), 1634–1644. https://doi.org/10.1016/j.bpj.2017.03.002 doi: 10.1016/j.bpj.2017.03.002
    [93] H. Kamada, Y. Imai, M. Nakamura, T. Ishikawa, T. Yamaguchi, Shear-induced platelet aggregation and distribution of thrombogenesis at stenotic vessels, Microcirculation, 24 (2017). https://doi.org/10.1111/micc.12355 doi: 10.1111/micc.12355
    [94] J. D. Horn, D. J. Maitland, J. Hartman, J. M. Ortega, A computational thrombus formation model: application to an idealized two-dimensional aneurysm treated with bare metal coils, Biomech. Model. Mechanobiol., 17 (2018), 1821–1838. https://doi.org/10.1007/s10237-018-1059-y doi: 10.1007/s10237-018-1059-y
    [95] R. Méndez Rojano, S. Mendez, F. Nicoud, Introducing the pro-coagulant contact system in the numerical assessment of device-related thrombosis, Biomech. Model. Mechanobiol., 17 (2018), 815–826. https://doi.org/10.1007/s10237-017-0994-3 doi: 10.1007/s10237-017-0994-3
    [96] B. Gu, A. Piebalgs, Y. Huang, C. Longstaff, A. D. Hughes, R. Chen, et al., Mathematical modelling of intravenous thrombolysis in acute ischaemic stroke: Effects of dose regimens on levels of fibrinolytic proteins and clot lysis time, Pharmaceutics, 11 (2019), 111. https://doi.org/10.3390/pharmaceutics11030111 doi: 10.3390/pharmaceutics11030111
    [97] K. Ayabe, S. Goto, H. Oka, H. Yabushita, M. Nakayama, A. Tomita, et al., Potential different impact of inhibition of thrombin function and thrombin generation rate for the growth of thrombi formed at site of endothelial injury under blood flow condition, Thromb. Res., 179 (2019), 121–127. https://doi.org/10.1016/j.thromres.2019.05.007 doi: 10.1016/j.thromres.2019.05.007
    [98] J. Chen, S. L. Diamond, Reduced model to predict thrombin and fibrin during thrombosis on collagen/tissue factor under venous flow: Roles of γ'-fibrin and factor XIa, PLoS Comput. Biol., 15 (2019), e1007266. https://doi.org/10.1371/journal.pcbi.1007266 doi: 10.1371/journal.pcbi.1007266
    [99] A. Bouchnita, V. Volpert, A multiscale model of platelet-fibrin thrombus growth in the flow, Comput. Fluids, 184 (2019), 10–20. https://doi.org/10.1016/j.compfluid.2019.03.021 doi: 10.1016/j.compfluid.2019.03.021
    [100] H. Hosseinzadegan, D. K. Tafti, A predictive model of thrombus growth in stenosed vessels with dynamic geometries, J. Med. Biol. Eng., 39 (2019), 605–621. https://doi.org/10.1007/s40846-018-0443-5 doi: 10.1007/s40846-018-0443-5
    [101] O. E. Kadri, V. D. Chandran, M. Surblyte, R. S. Voronov, In vivo measurement of blood clot mechanics from computational fluid dynamics based on intravital microscopy images, Comput. Biol. Med., 106 (2019), 1–11. https://doi.org/10.1016/j.compbiomed.2019.01.001 doi: 10.1016/j.compbiomed.2019.01.001
    [102] J. Du, D. Kim, G. Alhawael, D. N. Ku, A. L. Fogelson, Clot permeability, agonist transport, and platelet binding kinetics in arterial thrombosis, Biophys. J., 119 (2020), 2102–2115. https://doi.org/10.1016/j.bpj.2020.08.041 doi: 10.1016/j.bpj.2020.08.041
    [103] W. T. Wu, M. Zhussupbekov, N. Aubry, J. F. Antaki, M. Massoudi, Simulation of thrombosis in a stenotic microchannel: The effects of vWF-enhanced shear activation of platelets, Int. J. Eng. Sci., 147 (2020), 103206. https://doi.org/10.1016/j.ijengsci.2019.103206 doi: 10.1016/j.ijengsci.2019.103206
    [104] A. Bouchnita, K. Terekhov, P. Nony, Y. Vassilevski, V. Volpert, A mathematical model to quantify the effects of platelet count, shear rate, and injury size on the initiation of blood coagulation under venous flow conditions, PLoS One, 15 (2020), e0235392. https://doi.org/10.1371/journal.pone.0235392 doi: 10.1371/journal.pone.0235392
    [105] Z. L. Liu, D. N. Ku, C. K. Aidun, Mechanobiology of shear-induced platelet aggregation leading to occlusive arterial thrombosis: A multiscale in silico analysis, J. Biomech., 120 (2021), 110349. https://doi.org/10.1016/j.jbiomech.2021.110349 doi: 10.1016/j.jbiomech.2021.110349
    [106] V. N. Kaneva, J. L. Dunster, V. Volpert, F. Ataullahanov, M. A. Panteleev, D. Y. Nechipurenko, Modeling thrombus shell: Linking adhesion receptor properties and macroscopic dynamics, Biophys. J., 120 (2021), 334–351. https://doi.org/10.1016/j.bpj.2020.10.049 doi: 10.1016/j.bpj.2020.10.049
    [107] A. Yazdani, Y. Deng, H. Li, E. Javadi, Z. Li, S. Jamali, et al., Integrating blood cell mechanics, platelet adhesive dynamics and coagulation cascade for modelling thrombus formation in normal and diabetic blood, J. R. Soc. Interface., 18 (2021), 20200834. https://doi.org/10.1098/rsif.2020.0834 doi: 10.1098/rsif.2020.0834
    [108] C. Ma, Y. Ren, Q. Zheng, J. Wang, A computational model on cartesian adaptive grid for thrombosis simulation, IEEE Access, 10 (2022), 67694–67702. https://doi.org/10.1109/ACCESS.2022.3184123 doi: 10.1109/ACCESS.2022.3184123
    [109] R. Méndez Rojano, A. Lai, M. Zhussupbekov, G. W. Burgreen, K. Cook, J. F. Antaki, A fibrin enhanced thrombosis model for medical devices operating at low shear regimes or large surface areas, PLoS Comput. Biol., 18 (2022), e1010277. https://doi.org/10.1371/journal.pcbi.1010277 doi: 10.1371/journal.pcbi.1010277
    [110] M. Rezaeimoghaddam, F. N. van de Vosse, Continuum modeling of thrombus formation and growth under different shear rates, J. Biomech., 132 (2022), 110915. https://doi.org/10.1016/j.jbiomech.2021.110915 doi: 10.1016/j.jbiomech.2021.110915
    [111] A. S. Pisaryuk, N. M. Povalyaev, A. V Poletaev, A. M. Shibeko, Systems biology approach for personalized hemostasis correction, J. Pers. Med., 12 (2022). https://doi.org/10.3390/jpm12111903
    [112] M. Zhussupbekov, R. Méndez Rojano, W. T. Wu, J. F. Antaki, von Willebrand factor unfolding mediates platelet deposition in a model of high-shear thrombosis, Biophys. J., 121 (2022), 4033–4047. https://doi.org/10.1016/j.bpj.2022.09.040 doi: 10.1016/j.bpj.2022.09.040
    [113] Y. Wang, J. Luan, K. Luo, T. Zhu, J. Fan, Multi-constituent simulation of thrombosis in aortic dissection, Int. J. Eng. Sci., 184 (2023), 103817. http://dx.doi.org/10.1016/j.ijengsci.2023.103817 doi: 10.1016/j.ijengsci.2023.103817
    [114] R. Petkantchin, A. Rousseau, O. Eker, K. Zouaoui Boudjeltia, F. Raynaud, B. Chopard, A simplified mesoscale 3D model for characterizing fibrinolysis under flow conditions, Sci. Rep., 13 (2023), 13681. https://doi.org/10.1038/s41598-023-40973-1 doi: 10.1038/s41598-023-40973-1
    [115] K. Miyazawa, A. L. Fogelson, K. Leiderman, Inhibition of platelet-surface-bound proteins during coagulation under flow Ⅱ: Antithrombin and heparin, Biophys. J., 122 (2023), 230–240. https://doi.org/10.1016/j.bpj.2022.10.038 doi: 10.1016/j.bpj.2022.10.038
    [116] A. R. Rezaie, S. T. Olson, Calcium enhances heparin catalysis of the antithrombin-factor Xa reaction by promoting the assembly of an intermediate heparin-antithrombin-factor Xa bridging complex: Demonstration by rapid kinetics studies, Biochemistry, 39 (2000), 12083–12090. https://doi.org/10.1021/bi0011126 doi: 10.1021/bi0011126
    [117] M. Anand, K. Rajagopal, K. R. Rajagopal, A model for the formation, growth, and lysis of clots in quiescent plasma: A comparison between the effects of antithrombin Ⅲ deficiency and protein C deficiency, J. Theor. Biol., 253 (2008), 725–738. https://doi.org/10.1016/j.jtbi.2008.04.015 doi: 10.1016/j.jtbi.2008.04.015
    [118] M. J. Griffith, The heparin-enhanced antithrombin Ⅲ/thrombin reaction is saturable with respect to both thrombin and antithrombin Ⅲ, J. Biol. Chem., 257 (1982), 13302–13899.
    [119] M. J. Griffith, Kinetics of the heparin-enhanced antithrombin Ⅲ/thrombin reaction. Evidence for a template model for the mechanism of action of heparin, J. Biol. Chem., 257 (1982), 7360–7365.
    [120] C. M. Danforth, T. Orfeo, K. G. Mann, K. E. Brummel-Ziedins, S. J. Everse, The impact of uncertainty in a blood coagulation model, Math. Med. Biol., 26 (2009), 323–336. https://doi.org/10.1093/imammb/dqp011 doi: 10.1093/imammb/dqp011
    [121] P. P. Naidu, M. Anand, Importance of Ⅷa inactivation in a mathematical model for the formation, growth, and lysis of clots, Math. Model. Nat. Phenom., 9 (2014), 17–33. https://doi.org/10.1051/mmnp/20149603 doi: 10.1051/mmnp/20149603
    [122] F. Saitta, J. Masuri, M. Signorelli, S. Bertini, A. Bisio, D. Fessas, Thermodynamic insights on the effects of low-molecular-weight heparins on antithrombin Ⅲ, Thermochim. Acta, 713 (2022), 179248. https://doi.org/10.1016/j.tca.2022.179248 doi: 10.1016/j.tca.2022.179248
    [123] H. Minakami, M. Morikawa, T. Yamada, T. Yamada, Candidates for the determination of antithrombin activity in pregnant women, J. Perinat. Med., 39 (2011), 369–374. https://doi.org/10.1515/jpm.2011.026 doi: 10.1515/jpm.2011.026
    [124] K. C. Jones, K. G. Mann, A model for the tissue factor pathway to thrombin. Ⅱ: A mathematical simulation, J. Biol. Chem., 269 (1994), 23367–23373.
    [125] M. F. Hockin, K. C. Jones, S. J. Everse, K. G. Mann, A model for the stoichiometric regulation of blood coagulation, J. Biol. Chem., 277 (2002), 18322–18333. https://doi.org/10.1074/jbc.m201173200 doi: 10.1074/jbc.m201173200
    [126] R. Méndez Rojano, S. Mendez, D. Lucor, A. Ranc, M. Giansily-Blaizot, J. F. Schved, et al., Kinetics of the coagulation cascade including the contact activation system: sensitivity analysis and model reduction, Biomech. Model. Mechanobiol., 18 (2019), 1139–1153. https://doi.org/10.1007/s10237-019-01134-4 doi: 10.1007/s10237-019-01134-4
    [127] M. S. Chatterjee, W. S. Denney, H. Jing, S. L. Diamond, Systems biology of coagulation initiation: Kinetics of thrombin generation in resting and activated human blood, PLOS Comput. Biol., 6 (2010), e1000950. https://doi.org/10.1371/journal.pcbi.1000950 doi: 10.1371/journal.pcbi.1000950
    [128] A. Ranc, S. Bru, S. Mendez, M. Giansily-Blaizot, F. Nicoud, R. Méndez Rojano, Critical evaluation of kinetic schemes for coagulation, PLoS One, 18 (2023), e0290531. https://doi.org/10.1371/journal.pone.0290531 doi: 10.1371/journal.pone.0290531
    [129] M. F. Hockin, K. C. Jones, S. J. Everse, K. G. Mann, A model for the stoichiometric regulation of blood coagulation, J. Biol. Chem., 277 (2002), 18322–18333. https://doi.org/10.1074/jbc.m201173200 doi: 10.1074/jbc.m201173200
    [130] D. Luan, M. Zai, J. D. Varner, Computationally derived points of fragility of a human cascade are consistent with current therapeutic strategies, PLoS Comput. Biol., 3 (2007), 1347–1359. https://doi.org/10.1371/journal.pcbi.0030142 doi: 10.1371/journal.pcbi.0030142
    [131] K. G. Link, M. T. Stobb, J. Di Paola, K. B. Neeves, A. L. Fogelson, S. S. Sindi, et al., A local and global sensitivity analysis of a mathematical model of coagulation and platelet deposition under flow, PLoS One, 13 (2018), e0200917. https://doi.org/10.1371/journal.pone.0200917 doi: 10.1371/journal.pone.0200917
  • mbe-21-12-339 supplementary.pdf
  • Reader Comments
  • © 2024 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(1102) PDF downloads(94) Cited by(1)

Figures and Tables

Figures(2)  /  Tables(3)

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog