Beta Phase: Square45 is currently in beta testing. Expect some features or content to be incomplete or missing.
45

Wound Healing Theory

Field: General Surgery

Sequence of Expressions

Let Ci(t,x)C_i(t, \textbf{x}) be the concentration of cytokine ii at time tt and position x\textbf{x}. The spatio-temporal evolution of the cytokine concentration is governed by the reaction-diffusion equation:\n∂Ci∂t=Di∇2Ci+Ri(C)−kiCi+Si(C,Platelets)\frac{\partial C_i}{\partial t} = D_i \nabla^2 C_i + R_i(\mathbf{C}) - k_i C_i + S_i(\mathbf{C}, \text{Platelets})\nwhere DiD_i is the diffusion coefficient, Ri(C)R_i(\mathbf{C}) represents the source/sink terms from cellular activity (e.g., macrophage release), kik_i is the degradation rate, and SiS_i models the signaling dependence on platelet activation and leukocyte presence.
Define the concentrations of key factors: F=(F1,F2,Fibrinogen,Fibrin)\mathbf{F} = (F_1, F_2, \text{Fibrinogen}, \text{Fibrin}). The kinetics of the coagulation cascade are modeled by the following system of coupled ordinary differential equations (ODE), focusing on the conversion of fibrinogen (Fg\text{Fg}) to fibrin (F\text{F}):\nd[F]dt=kII[Prothrombin][FactorXa][Ca2+]−kdeg[F][Plasmin]\frac{d[\text{F}]}{dt} = k_{II} [\text{Prothrombin}] [\text{Factor} \text{Xa}] [\text{Ca}^{2+}] - k_{deg} [\text{F}] [\text{Plasmin}]\nd[Fg]dt=−kI[FactorII][FactorVII][Ca2+][Fg]+Source\frac{d[\text{Fg}]}{dt} = -k_{I} [\text{Factor} \text{II}] [\text{Factor} \text{VII}] [\text{Ca}^{2+}] [\text{Fg}] + \text{Source}
Let M(t,x)M(t, \textbf{x}) be the concentration of the provisional matrix (collagen, fibronectin) at time tt and position x\textbf{x}. The rate of matrix deposition is modeled by a reaction-diffusion equation incorporating cell-mediated synthesis and degradation:\n∂M∂t=DM∇2M+ksynth[Fibroblast]Mmax−MMmax−kdegM[MMP]+Rangi\frac{\partial M}{\partial t} = D_M \nabla^2 M + k_{synth} [\text{Fibroblast}] \frac{M_{max} - M}{M_{max}} - k_{deg} M [\text{MMP}] + R_{angi} \nwhere ksynthk_{synth} is the synthesis rate, MmaxM_{max} is the saturation capacity, and RangiR_{angi} accounts for matrix support provided by newly formed vessels.
Consider the concentration of Vascular Endothelial Growth Factor (VEGF), CVEGFC_{VEGF}, and the vessel density, ρv\rho_v. The dynamics are governed by a reaction-diffusion system:\n∂CVEGF∂t=DVEGF∇2CVEGF−kbindCVEGF[Endothelial]+SVEGF(ρv)\frac{\partial C_{VEGF}}{\partial t} = D_{VEGF} \nabla^2 C_{VEGF} - k_{bind} C_{VEGF} [\text{Endothelial}] + S_{VEGF}(\rho_v) \n∂ρv∂t=kangiCVEGFρv(1−ρvρmax)−kpruneρv\frac{\partial \rho_v}{\partial t} = k_{angi} C_{VEGF} \rho_v (1 - \frac{\rho_v}{\rho_{max}}) - k_{prune} \rho_v \nwhere kangik_{angi} models the growth rate stimulated by VEGF, and kprunek_{prune} accounts for vessel regression.
Let v(t,x)v(t, \textbf{x}) be the velocity of the migrating epithelial sheet boundary. The migration is modeled using a chemotaxis framework, where the velocity is proportional to the gradient of a chemoattractant CchemC_{chem} (e.g., growth factors):\nv=dxdt=β∇Cchem+χ1+e−k(Cchem−Cthreshold)n\mathbf{v} = \frac{d\textbf{x}}{dt} = \beta \nabla C_{chem} + \frac{\chi}{1 + e^{-k(C_{chem} - C_{threshold})}} \textbf{n} \nwhere β\beta is the chemotactic coefficient, χ\chi is the random motility component, and n\textbf{n} is the normal vector to the wound edge.
Define the concentration of active MMPs, [MMP][MMP], and the concentration of the substrate, [ECM][ECM]. The balance is described by the rate of degradation and the regulatory feedback loop:\nd[ECM]dt=ksynth[Collagen]−kMMP[MMP][ECM]+kremodel(1−[ECM][ECM]eq)\frac{d[ECM]}{dt} = k_{synth} [\text{Collagen}] - k_{MMP} [MMP] [ECM] + k_{remodel} (1 - \frac{[ECM]}{[ECM]_{eq}}) \nd[MMP]dt=kind[Cytokines]−kinact[MMP][TIMP]−kclear[MMP]\frac{d[MMP]}{dt} = k_{ind} [\text{Cytokines}] - k_{inact} [MMP] [\text{TIMP}] - k_{clear} [MMP] \nwhere [TIMP][TIMP] represents Tissue Inhibitors of Metalloproteinases, establishing the critical [MMP]/[TIMP][MMP]/[TIMP] ratio.
Let CIIIC_{III} and CIC_{I} be the concentrations of Type III and Type I collagen, respectively. The synthesis rate is modeled as a transition from Type III to Type I, followed by cross-linking (XL\text{XL}). \nd[CIII]dt=ksynth,III[Fibroblast]−ktrans[CIII]+kdeg,III[CIII]\frac{d[C_{III}]}{dt} = k_{synth, III} [\text{Fibroblast}] - k_{trans} [C_{III}] + k_{deg, III} [C_{III}]\nd[CI]dt=ktrans[CIII]−kdeg,I[CI]+kXL[CI]2[Lysyl Oxidase]\frac{d[C_{I}]}{dt} = k_{trans} [C_{III}] - k_{deg, I} [C_{I}] + k_{XL} [C_{I}]^2 [\text{Lysyl Oxidase}]\nwhere kXLk_{XL} is the cross-linking rate, which is dependent on the concentration of enzymatic cross-linkers.
Consider the wound area A(t)A(t) and the contractile stress σc\sigma_c. The rate of change of the wound area is governed by the contractile force exerted by myofibroblasts (Myo\text{Myo}) acting on the wound edges:\ndAdt=−2τ0EALdudt\frac{dA}{dt} = - \frac{2\tau_0}{E} \frac{A}{L} \frac{d\textbf{u}}{dt} \nwhere τ0\tau_0 is the maximum contractile stress, EE is the tissue elasticity modulus, LL is the wound length, and dudt\frac{d\textbf{u}}{dt} is the rate of displacement of the wound edges, which is proportional to the density of active myofibroblasts.
Let M1M_1 and M2M_2 be the concentrations of M1 (inflammatory) and M2 (reparative) macrophages. Their dynamic balance is modeled by a system of coupled ODEs dependent on local cytokine concentrations, CIL−4C_{IL-4} and CIFN−γC_{IFN-\gamma}:\nd[M1]dt=krec[Neutrophils]−kM1→M2[M1]CIL−4+kM2→M1[M2]CIFN−γ\frac{d[M_1]}{dt} = k_{rec} [\text{Neutrophils}] - k_{M1 \to M2} [M_1] C_{IL-4} + k_{M2 \to M1} [M_2] C_{IFN-\gamma} \nd[M2]dt=kM1→M2[M1]CIL−4−krec[Neutrophils]+kM2→M1[M2]CIFN−γ\frac{d[M_2]}{dt} = k_{M1 \to M2} [M_1] C_{IL-4} - k_{rec} [\text{Neutrophils}] + k_{M2 \to M1} [M_2] C_{IFN-\gamma} \n(Note: The signs are simplified to show the dependence on the polarizing cytokines.)
Define the collagen deposition rate RcollR_{coll} and the degradation rate RdegR_{deg}. Pathological scarring is characterized by a positive feedback loop where excessive synthesis overwhelms degradation, leading to a runaway accumulation of matrix MM. \nd[M]dt=Rcoll−Rdeg=ksynth[Myo][TGF-β]−kMMP[MMP][M]\frac{d[M]}{dt} = R_{coll} - R_{deg} = k_{synth} [\text{Myo}] [\text{TGF-}\beta] - k_{MMP} [MMP] [M]\nPathological condition: ksynth[Myo][TGF-β]>kMMP[MMP][M]k_{synth} [\text{Myo}] [\text{TGF-}\beta] > k_{MMP} [MMP] [M] and a failure in the negative feedback mechanism that normally downregulates [TGF-β][\text{TGF-}\beta] and [Myo][\text{Myo}] after wound closure.