Lung tissue engineering and preservation of alveolar microstructure using a novel casting method

Abstract

We used a rat model to decellularize and seed alveolar cells on a three-dimensional lung scaffold to preserve alveolar microarchitecture. We verified the preservation of terminal respiratory structure by casting and by scanning electron microscopy (SEM) of the casts after decellularization. Whole lungs were obtained from 12 healthy Sprague-Dawley rats, cannulated through the trachea under sterile conditions, and decellularized using a detergent-based method. Casting of both natural and decellularized lungs was performed to verify preservation of the inner microstructure of scaffolds for further cell seeding. Alveolar cell seeding was performed using green fluorescent protein (GFP) lung cells and non-GFP lung cells, and a peristaltic pump. We assessed cell seeding using histological and immunohistochemical staining, and enzymatic evaluation. All cellular components were removed completely from the scaffolds, and histological staining and SEM of casts were used to verify the preservation of tissue structure. Tensile tests verified conservation of biomechanical properties. The hydroxyproline content of decellularized lungs was similar to native lung. Histological and immunohistochemical evaluations showed effective cell seeding on decellularized matrices. Enzymatic measurement of trypsin and alpha 1 antitrypsin suggested the potential functional properties of the regenerated lungs. Casts produced by our method have satisfactory geometrical properties for further cell seeding of lung scaffolds. Preservation of micro-architecture and terminal alveoli that was confirmed by SEM of lung casts increases the probability of an effective cell seeding process.     

The effect of autologous bone marrow stromal cells differentiated on scaffolds for canine tibial bone reconstruction

Bone marrow contains mesenchymal stem cells that form many tissues. Various scaffolds are available for bone reconstruction by tissue engineering. Osteoblastic differentiated bone marrow stromal cells (BMSC) promote osteogenesis on scaffolds and stimulate bone regeneration. We investigated the use of cultured autologous BMSC on different scaffolds for healing defects in tibias of adult male canines. BMSC were isolated from canine humerus bone marrow, differentiated into osteoblasts in culture and loaded onto porous ceramic scaffolds including hydroxyapatite 1, hydroxyapatite gel and calcium phosphate. Osteoblast differentiation was verified by osteonectine and osteocalcine immunocytochemistry. The scaffolds with stromal cells were implanted in the tibial defect. Scaffolds without stromal cells were used as controls. Sections from the defects were processed for histological, ultrastructural, immunohistochemical and histomorphometric analyses to analyze the healing of the defects. BMSC were spread, allowed to proliferate and differentiate to osteoblasts as shown by alizarin red histochemistry, and osteocalcine and osteonectine immunostaining. Scanning electron microscopy showed that BMSC on the scaffolds were more active and adhesive to the calcium phosphate scaffold compared to the others. Macroscopic bone formation was observed in all groups, but scaffolds with stromal cells produced significantly better results. Bone healing occurred earlier and faster with stromal cells on the calcium phosphate scaffold and produced more callus compared to other scaffolds. Tissue healing and osteoblastic marker expression also were better with stromal cells on the scaffolds. Increased trabecula formation, cell density and decreased fibrosis were observed in the calcium phosphate scaffold with stromal cells. Autologous cultured stromal cells on the scaffolds were useful for healing of canine tibial bone defects. The calcium phosphate scaffold was the best for both cell differentiation in vitro and bone regeneration in vivo. It may be possible to improve healing of bone defects in humans using stem cells from bone marrow.     

A computational model of in vitro angiogenesis based on extracellular matrix fibre orientation

Recent interest in the process of vascularisation within the biomedical community has motivated numerous new research efforts focusing on the process of angiogenesis. Although the role of chemical factors during angiogenesis has been well documented, the role of mechanical factors, such as the interaction between angiogenic vessels and the extracellular matrix, remains poorly understood. In vitro methods for studying angiogenesis exist; however, measurements available using such techniques often suffer from limited spatial and temporal resolutions. For this reason, computational models have been extensively employed to investigate various aspects of angiogenesis. This paper outlines the formulation and validation of a simple and robust computational model developed to accurately simulate angiogenesis based on length, branching and orientation morphometrics collected from vascularised tissue constructs. Microvessels were represented as a series of connected line segments. The morphology of the vessels was determined by a linear combination of the collagen fibre orientation, the vessel density gradient and a random walk component. Excellent agreement was observed between computational and experimental morphometric data over time. Computational predictions of microvessel orientation within an anisotropic matrix correlated well with experimental data. The accuracy of this modelling approach makes it a valuable platform for investigating the role of mechanical interactions during angiogenesis.     

A mixed finite element formulation for a non-linear,transversely isotropic material model for the cardiac tissue

In this paper we present a mixed finite element method for modeling the passive properties of the myocardium. The passive properties are described by a non-linear, transversely isotropic, hyperelastic material model, and the myocardium is assumed to be almost incompressible. Single-field, pure displacement-based formulations are known to cause numerical difficulties when applied to incompressible or slightly compressible material cases. This paper presents an alternative approach in the form of a mixed formulation, where a separately interpolated pressure field is introduced as a primary unknown in addition to the displacement field. Moreover, a constraint term is included in the formulation to enforce (almost) incompressibility. Numerical results presented in the paper demonstrate the difficulties related to employing a pure displacement-based method, applying a set of physically relevant material parameter values for the cardiac tissue. The same problems are not experienced for the proposed mixed method. We show that the mixed formulation provides reasonable numerical results for compressible as well as nearly incompressible cases, also in situations of large fiber stretches. There is good agreement between the numerical results and the underlying analytical models.     

Predicting the structural integrity of bone defects repaired using bone graft materials

Bone defects create stress concentrations which can cause fracture under impact or cyclic loading. Defects are often repaired by filling them with a bone graft material; this will reduce the stress concentration, but not completely, because these materials have lower stiffness than bone. The fracture risk decreases over time as the graft material is replaced by living bone. Many new bone graft materials are being developed, using tissue engineering and other techniques, but currently there is no rational way to compare these materials and predict their effectiveness in repairing a given defect. This paper describes, for the first time, a theoretical model which can be used to predict failure by brittle fracture or fatigue, initiating at the defect. Preliminary results are presented, concentrating on the prediction of stress fracture during the crucial post-operative period. It is shown that the likelihood of fracture is strongly influenced by the shape of the defect as well as its size, and also by the level of post-operative exercise. The most important finding is that bone graft materials can be successful in preventing fracture even when their mechanical properties are greatly inferior to those of bone. Future uses of this technique include pre-clinical assessment of bone replacement materials and pre-operative planning in orthopaedic surgery.     

The use of preserved tissue in finite element modelling of fresh ovine vertebral behaviour

The aim of this study was to investigate whether the predicted finite element (FE) stiffness of vertebral bone is altered when using images of preserved rather than fresh tissue to generate specimen-specific FE models. Fresh ovine vertebrae were used to represent embalmed (n = 3) and macerated dry-bone (n = 3) specimens and treated accordingly. Specimens were scanned pre- and post-treatment using micro-computed tomography. A constant threshold level derived from these images was used to calculate the respective bone volume fraction (BV/TV) from which the conversion factor validated for fresh tissue was used to determine material properties that were assigned to corresponding FE models. Results showed a definite change in the BV/TV between the fresh and the preserved bone. However, the changes in the predicted FE stiffness were not generally greater than the variations expected from assignment of loading and boundary conditions. In conclusion, images of preserved tissue can be used to generate FE models that are representative of fresh tissue with a tolerable level of error.     

Determination of material parameters of the two-dimensional Holzapfel–Weizsäcker type model based on uniaxial extension data of arterial walls

Soft tissues are anisotropic materials yet a majority of mechanical property tests have been uniaxial, which often failed to recapitulate the tensile response in other directions. This paper aims to study the feasibility of determining material parameters of anisotropic tissues by uniaxial extension with a minimal loss of anisotropic information. We assumed that by preselecting a certain constitutive model, we could give the constitutive parameters based on uniaxial extension data from orthogonal strip samples. In our study, the Holzapfel–Weizsäcker type strain energy density function (H–W model) was used to determine the material parameters of arterial walls from two fresh donation bodies. The key points we applied were the relationships between strain components in uniaxial tensile tests and the methods of stochastic optimisation. Further numerical experiments were taken. The estimate–effect ratio, defined by the number of data with the precision of estimation less than 0.5% over whole size of data, was calculated to demonstrate the feasibility of our method. The material parameters for Chinese aorta and pulmonary artery were given with the maximum root mean square (RMS) errors 0.042, and the minimal estimate–effect ratio in numerical experiments was 90.79%. Our results suggest that the constitutive parameters of arterial walls can be determined from uniaxial extension data, given the passive mechanical behaviour governed by H–W model. This method may apply to other tissues using different constitutive models.     

Which factors influence the ability of a computational model to predict the in vivo deformation behaviour of skeletal muscle?

Deep tissue injury (DTI) is a severe form of pressure ulcer where tissue damage starts in deep tissues underneath intact skin. Tissue deformation may play an important role in the aetiology, which can be investigated using an experimental–numerical approach. Recently, an animal-specific finite element model has been developed to simulate experiments in which muscle tissue was compressed with an indenter. In this study, the material behaviour and boundary conditions were adapted to improve the agreement between model and experiment and to investigate the influence of these adaptations on the predicted strain distribution. The use of a highly nonlinear material law and including friction between the indenter and the muscle both improved the quality of the model and considerably influenced the estimated strain distribution. With the improved model, the required sample size to detect significant differences between loading conditions can be diminished, which is clearly relevant in experiments involving animals.     

How does muscle stiffness affect the internal deformations within the soft tissue layers of the buttocks under constant loading?

Mechanical loading of soft tissues covering bony prominences can cause skeletal muscle damage, ultimately resulting in a severe pressure ulcer termed deep tissue injury (DTI). Deformation plays an important role in the aetiology of DTI. Therefore, it is essential to minimise internal muscle deformations in subjects at risk of DTI. As an example, spinal cord-injured (SCI) individuals exhibit structural changes leading to a decrease in muscle thickness and stiffness, which subsequently increase the tissue deformations. In the present study, an animal-specific finite element model, where the geometry and boundary conditions were derived from magnetic resonance images, was developed. It was used to investigate the internal deformations in the muscle, fat and skin layers of the porcine buttocks during loading. The model indicated the presence of large deformations in both the muscle and the fat layers, with maximum shear strains up to 0.65 in muscle tissue and 0.63 in fat. Furthermore, a sensitivity analysis showed that the tissue deformations depend considerably on the relative stiffness values of the different tissues. For example, a change in muscle stiffness had a large effect on the muscle deformations. A 50% decrease in stiffness caused an increase in maximum shear strain from 0.65 to 0.99, whereas a 50% increase in stiffness resulted in a decrease in maximum shear strain from 0.65 to 0.49. These results indicate the importance of restoring tissue properties after SCI, with the use of, for example, electrical stimulation, to prevent the development of DTI.     

Detailed finite element modelling of deep needle insertions into a soft tissue phantom using a cohesive approach

Detailed finite element modelling of needle insertions into soft tissue phantoms encounters difficulties of large deformations, high friction, contact loading and material failure. This paper demonstrates the use of cohesive elements in high-resolution finite element models to overcome some of the issues associated with these factors. Experiments are presented enabling extraction of the strain energy release rate during crack formation. Using data from these experiments, cohesive elements are calibrated and then implemented in models for validation of the needle insertion process. Successful modelling enables direct comparison of finite element and experimental force–displacement plots and energy distributions. Regions of crack creation, relaxation, cutting and full penetration are identified. By closing the loop between experiments and detailed finite element modelling, a methodology is established which will enable design modifications of a soft tissue probe that steers through complex mechanical interactions with the surrounding material.