A Mathematical Model for Postimplant Collagen Remodeling in an Autologous Engineered Pulmonary Arterial Conduit

This study was undertaken to develop a mathematical model of the long-term in vivo remodeling processes in postimplanted pulmonary artery (PA) conduits. Experimental results from two extant ovine in vivo studies, wherein polyglycolic-acid (PGA)/poly-L-lactic acid tubular conduits were constructed, cell seeded, incubated for 4 weeks, and then implanted in mature sheep to obtain the remodeling data for up to two years. Explanted conduit analysis included detailed novel structural and mechanical studies. Results in both studies indicated that the in vivo conduits remained dimensionally stable up to 80 weeks, so that the conduits maintained a constant in vivo stress and deformation state. In contrast, continued remodeling of the constituent collagen fiber network as evidenced by an increase in effective tissue uniaxial tangent modulus, which then stabilized by one year postimplant. A mesostructural constitute model was then applied to extant planar biaxial mechanical data and revealed several interesting features, including an initial pronounced increase in effective collagen fiber modulus, paralleled by a simultaneous shift toward longer, more uniformly length-distributed collagen fibers. Thus, while the conduit remained dimensionally stable, its internal collagen fibrous structure and resultant mechanical behaviors underwent continued remodeling that stabilized by one year. A time-evolving structural mixture-based mathematical model specialized for this unique form of tissue remodeling was developed, with a focus on time-evolving collagen fiber stiffness as the driver for tissue-level remodeling. The remodeling model was able to fully reproduce (1) the observed tissue-level increases in stiffness by time-evolving simultaneous increases in collagen fiber modulus and lengths, (2) maintenance of the constant collagen fiber angular dispersion, and (3) stabilization of the remodeling processes at one year. Collagen fiber remodeling geometry was directly verified experimentally by histological analysis of the time-evolving collagen fiber crimp, which matches model predictions very closely. Interestingly, the remodeling model indicated that the basis for tissue homeostasis was maintenance of the collagen fiber ensemble stress for all orientations, and not individual collagen fiber stresses. Unlike other growth and remodeling models that traditionally treat changes in the external boundary conditions (e.g., changes in blood pressure) as the primary input stimuli, the driver herein is changes to the internal constituent collagen fiber themselves due to cellular mediated cross-linking.