Bioengineered Approaches to Prevent Hypertrophic Scar Contraction

Burn injuries in the United States account for over one million hospital admissions per year, with treatment estimated at four billion dollars. Of severe burn patients, 30-90% will develop hypertrophic scars (HSc). Current burn therapies rely upon the use of bioengineered skin equivalents (BSEs), which assist in wound healing but do not prevent HSc. HSc contraction occurs of 6-18 months and results in the formation of a fixed, inelastic skin deformity, with 60% of cases occurring across a joint. HSc contraction is characterized by abnormally high presence of contractile myofibroblasts which normally apoptose at the completion of the proliferative phase of wound healing. Additionally, clinical observation suggests that the likelihood of HSc is increased in injuries with a prolonged immune response. Given the pathogenesis of HSc, we hypothesize that BSEs should be designed with two key anti-scarring characterizes: (1) 3D architecture and surface chemistry to mitigate the inflammatory microenvironment and decrease myofibroblast transition; and (2) using materials which persist in the wound bed throughout the remodeling phase of repair. We employed electrospinning and 3D printing to generate scaffolds with well-controlled degradation rate, surface coatings, and 3D architecture to explore our hypothesis through four aims.

In the first aim, we evaluate the impact of elastomeric, randomly-oriented biostable polyurethane (PU) scaffold on HSc-related outcomes. In unwounded skin, native collagen is arranged randomly, elastin fibers are abundant, and myofibroblasts are absent. Conversely, in scar contractures, collagen is arranged in linear arrays and elastin fibers are few, while myofibroblast density is high. Randomly oriented collagen fibers native to the uninjured dermis encourage random cell alignment through contact guidance and do not transmit as much force as aligned collagen fibers. However, the linear ECM serves as a system for mechanotransduction between cells in a feed-forward mechanism, which perpetuates ECM remodeling and myofibroblast contraction. The electrospinning process allowed us to create scaffolds with randomly-oriented fibers that promote random collagen deposition and decrease myofibroblast formation. Compared to an in vitro HSc contraction model, fibroblast-seeded PU scaffolds significantly decreased matrix and myofibroblast formation. In a murine HSc model, collagen coated PU (ccPU) scaffolds significantly reduced HSc contraction as compared to untreated control wounds and wounds treated with the clinical standard of care. The data from this study suggest that electrospun ccPU scaffolds meet the requirements to mitigate HSc contraction including: reduction of in vitro HSc related outcomes, diminished scar stiffness, and reduced scar contraction. While clinical dogma suggests treating severe burn patients with rapidly biodegrading skin equivalents, these data suggest that a more long-term scaffold may possess merit in reducing HSc.

In the second aim, we further investigate the impact of scaffold longevity on HSc contraction by studying a degradable, elastomeric, randomly oriented, electrospun micro-fibrous scaffold fabricated from the copolymer poly(l-lactide-co-ε-caprolactone) (PLCL). PLCL scaffolds displayed appropriate elastomeric and tensile characteristics for implantation beneath a human skin graft. In vitro analysis using normal human dermal fibroblasts (NHDF) demonstrated that PLCL scaffolds decreased myofibroblast formation as compared to an in vitro HSc contraction model. Using our murine HSc contraction model, we found that HSc contraction was significantly greater in animals treated with standard of care, Integra, as compared to those treated with collagen coated-PLCL (ccPLCL) scaffolds at d 56 following implantation. Finally, wounds treated with ccPLCL were significantly less stiff than control wounds at d 56 in vivo. Together, these data further solidify our hypothesis that scaffolds which persist throughout the remodeling phase of repair represent a clinically translatable method to prevent HSc contraction.

In the third aim, we attempt to optimize cell-scaffold interactions by employing an anti-inflammatory coating on electrospun PLCL scaffolds. The anti-inflammatory sub-epidermal glycosaminoglycan, hyaluronic acid (HA) was used as a coating material for PLCL scaffolds to encourage a regenerative healing phenotype. To minimize local inflammation, an anti-TNFα monoclonal antibody (mAB) was conjugated to the HA backbone prior to PLCL coating. ELISA analysis confirmed mAB activity following conjugation to HA (HA+mAB), and following adsorption of HA+mAB to the PLCL backbone [(HA+mAB)PLCL]. Alican blue staining demonstrated thorough HA coating of PLCL scaffolds using pressure-driven adsorption. In vitro studies demonstrated that treatment with (HA+mAB)PLCL prevented downstream inflammatory events in mouse macrophages treated with soluble TNFα. In vivo studies using our murine HSc contraction model suggested positive impact of HA coating, which was partiall impeded by the inclusion of the TNFα mAB. Further characterization of the inflammatory microenvironment of our murine model is required prior to conclusions regarding the potential for anti-TNFα therapeutics for HSc. Together, our data demonstrate the development of a complex anti-inflammatory coating for PLCL scaffolds, and the potential impact of altering the ECM coating material on HSc contraction.

In the fourth aim, we investigate how scaffold design, specifically pore dimensions, can influence myofibroblast interactions and subsequent formation of OB-cadherin positive adherens junctions in vitro. We collaborated with Wake Forest University to produce 3D printed (3DP) scaffolds with well-controlled pore sizes we hypothesized that decreasing pore size would mitigate intra-cellular communication via OB-cadherin-positive adherens junctions. PU was 3D printed via pressure extrusion in basket-weave design with feature diameter of ~70 µm and pore sizes of 50, 100, or 150 µm. Tensile elastic moduli of 3DP scaffolds were similar to Integra; however, flexural moduli of 3DP were significantly greater than Integra. 3DP scaffolds demonstrated ~50% porosity. 24 h and 5 d western blot data demonstrated significant increases in OB-cadherin expression in 100 µm pores relative to 50 µm pores, suggesting that pore size may play a role in regulating cell-cell communication. To analyze the impact of pore size in these scaffolds on scarring in vivo, scaffolds were implanted beneath skin graft in a murine HSc model. While flexural stiffness resulted in graft necrosis by d 14, cellular and blood vessel integration into scaffolds was evident, suggesting potential for this design if employed in a less stiff material. In this study, we demonstrate for the first time that pore size alone impacts OB-cadherin protein expression in vitro, suggesting that pore size may play a role on adherens junction formation affiliated with the fibroblast-to-myofibroblast transition. Overall, this work introduces a new bioengineered scaffold design to both study the mechanism behind HSc and prevent the clinical burden of this contractile disease.

Together, these studies inform the field of critical design parameters in scaffold design for the prevention of HSc contraction. We propose that scaffold 3D architectural design, surface chemistry, and longevity can be employed as key design parameters during the development of next generation, low-cost scaffolds to mitigate post-burn hypertrophic scar contraction. The lessening of post-burn scarring and scar contraction would improve clinical practice by reducing medical expenditures, increasing patient survival, and dramatically improving quality of life for millions of patients worldwide.

Méthodes de Bootstrap pour les modèles à facteurs

Cette thèse développe des méthodes bootstrap pour les modèles à facteurs qui sont couram- ment utilisés pour générer des prévisions depuis l'article pionnier de Stock et Watson (2002) sur les indices de diffusion. Ces modèles tolèrent l'inclusion d'un grand nombre de variables macroéconomiques et financières comme prédicteurs, une caractéristique utile pour inclure di- verses informations disponibles aux agents économiques. Ma thèse propose donc des outils éco- nométriques qui améliorent l'inférence dans les modèles à facteurs utilisant des facteurs latents extraits d'un large panel de prédicteurs observés. Il est subdivisé en trois chapitres complémen- taires dont les deux premiers en collaboration avec Sílvia Gonçalves et Benoit Perron. Dans le premier article, nous étudions comment les méthodes bootstrap peuvent être utilisées pour faire de l'inférence dans les modèles de prévision pour un horizon de h périodes dans le futur. Pour ce faire, il examine l'inférence bootstrap dans un contexte de régression augmentée de facteurs où les erreurs pourraient être autocorrélées. Il généralise les résultats de Gonçalves et Perron (2014) et propose puis justifie deux approches basées sur les résidus : le block wild bootstrap et le dependent wild bootstrap. Nos simulations montrent une amélioration des taux de couverture des intervalles de confiance des coefficients estimés en utilisant ces approches comparativement à la théorie asymptotique et au wild bootstrap en présence de corrélation sérielle dans les erreurs de régression. Le deuxième chapitre propose des méthodes bootstrap pour la construction des intervalles de prévision permettant de relâcher l'hypothèse de normalité des innovations. Nous y propo- sons des intervalles de prédiction bootstrap pour une observation h périodes dans le futur et sa moyenne conditionnelle. Nous supposons que ces prévisions sont faites en utilisant un ensemble de facteurs extraits d'un large panel de variables. Parce que nous traitons ces facteurs comme latents, nos prévisions dépendent à la fois des facteurs estimés et les coefficients de régres- sion estimés. Sous des conditions de régularité, Bai et Ng (2006) ont proposé la construction d'intervalles asymptotiques sous l'hypothèse de Gaussianité des innovations. Le bootstrap nous permet de relâcher cette hypothèse et de construire des intervalles de prédiction valides sous des hypothèses plus générales. En outre, même en supposant la Gaussianité, le bootstrap conduit à des intervalles plus précis dans les cas où la dimension transversale est relativement faible car il prend en considération le biais de l'estimateur des moindres carrés ordinaires comme le montre une étude récente de Gonçalves et Perron (2014). Dans le troisième chapitre, nous suggérons des procédures de sélection convergentes pour les regressions augmentées de facteurs en échantillons finis. Nous démontrons premièrement que la méthode de validation croisée usuelle est non-convergente mais que sa généralisation, la validation croisée «leave-d-out» sélectionne le plus petit ensemble de facteurs estimés pour l'espace généré par les vraies facteurs. Le deuxième critère dont nous montrons également la validité généralise l'approximation bootstrap de Shao (1996) pour les regressions augmentées de facteurs. Les simulations montrent une amélioration de la probabilité de sélectionner par- cimonieusement les facteurs estimés comparativement aux méthodes de sélection disponibles. L'application empirique revisite la relation entre les facteurs macroéconomiques et financiers, et l'excès de rendement sur le marché boursier américain. Parmi les facteurs estimés à partir d'un large panel de données macroéconomiques et financières des États Unis, les facteurs fortement correlés aux écarts de taux d'intérêt et les facteurs de Fama-French ont un bon pouvoir prédictif pour les excès de rendement.

State of Connecticut.

Scale 1:125,000: 1 in. equals approx. 2 miles.

Saranac Lake, New York /

Scale 1:25,000; 1 cm. to 250 m.

Saranac Lake, New York /

Scale 1:25,000; 1 cm. to 250 m.

Hydrologic unit map--1974, State of New Jersey.

Scale 1:500,000; 1 in. equals approx. 8 miles.

State of New Jersey /

Scale 1:500,000; 1 in. equals approx. 8 miles.

State of New Jersey /

Scale 1:500,000; 1 in. equals approx. 8 miles.

New Jersey.

Scale 1:500,000; 1 cm. equals 5 kilometers.

Digital revision of USGS 1:500,000-scale New Jersey state base map /

Shows stages and operations undertaken in revising the New Jersey state base map.

Water resources investigations in New Jersey, 1972.

Scales vary.

Transportation map of New Jersey.

Scale ca. 1:538,560; 1 in. represents approx. 8.5 miles.

Transportation map of New Jersey.

Scale ca. 1:538,560; 1 in. represents approx. 8.5 miles.

Burlington County, New Jersey /

UTM and state grids.