Tuneable hydrogel patterns in pillarless microfluidic devices

Organ-on-chip (OOC) technology has recently emerged as a powerful tool to mimic physiological or pathophysiological conditions through cell culture in microfluidic devices. One of its main goals is bypassing animal testing and encouraging more personalized medicine. The recent incorporation of hydrogels as 3D scaffolds into microfluidic devices has changed biomedical research since they provide a biomimetic extracellular matrix to recreate tissue architectures. However, this technology presents some drawbacks such as the necessity for physical structures as pillars to confine these hydrogels, as well as the difficulty in reaching different shapes and patterns to create convoluted gradients or more realistic biological structures. In addition, pillars can also interfere with the fluid flow, altering the local shear forces and, therefore, modifying the mechanical environment in the OOC model. In this work, we present a methodology based on a plasma surface treatment that allows building cell culture chambers with abutment-free patterns capable of producing precise shear stress distributions. Therefore, pillarless devices with arbitrary geometries are needed to obtain more versatile, reliable, and biomimetic experimental models. Through computational simulation studies, these shear stress changes are demonstrated in different designed and fabricated geometries. To prove the versatility of this new technique, a blood-brain barrier model has been recreated, achieving an uninterrupted endothelial barrier that emulates part of the neurovascular network of the brain. Finally, we developed a new technology that could avoid the limitations mentioned above, allowing the development of biomimetic OOC models with complex and adaptable geometries, with cell-to-cell contact if required, and where fluid flow and shear stress conditions could be controlled.

Experimentally-guided in silico design of engineered heart tissues to improve cardiac electrical function after myocardial infarction

Engineered heart tissues (EHTs) built from human induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs) showed promising results for cardiac function restoration following myocardial infarction. Nevertheless, human iPSC-CMs have longer action potential and lower cell-to-cell coupling than adult-like CMs. These immature electrophysiological properties favor arrhythmias due to the generation of electrophysiological gradients when hiPSC-CMs are injected in the cardiac tissue. Culturing hiPSC-CMs on three-dimensional (3D) scaffolds can promote their maturation and influence their alignment. However, it is still uncertain how on-scaffold culturing influences the overall electrophysiology of the in vitro and implanted EHTs, as it requires expensive and time consuming experimentation. Here, we computationally investigated the impact of the scaffold design on the EHT electrical depolarization and repolarization before and after engraftment on infarcted tissue. We first acquired and processed electrical recordings from in vitro EHTs, which we used to calibrate the modeling and simulation of in silico EHTs to replicate experimental outcomes. Next, we built in silico EHT models for a range of scaffold pore sizes, shapes (square, rectangular, auxetic, hexagonal) and thicknesses. In this setup, we found that scaffolds made of small (0.2 mm2), elongated (30° half-angle) hexagons led to faster EHT activation and better mimicked the cardiac anisotropy. The scaffold thickness had a marginal role on the not engrafted EHT electrophysiology. Moreover, EHT engraftment on infarcted tissue showed that the EHT conductivity should be at least 5% of that in healthy tissue for bidirectional EHT-myocardium electrical propagation. For conductivities above such threshold, the scaffold made of small elongated hexagons led to the lowest activation time (AT) in the coupled EHT-myocardium. If the EHT conductivity was further increased and the hiPSC-CMs were uniformly oriented parallel to the epicardial cells, the total AT and the repolarization time gradient decreased substantially, thus minimizing the likelihood for arrhythmias after EHT transplantation.

The Role of Mechanical Properties and Structure of Type I Collagen Hydrogels on Colorectal Cancer Cell Migration

Mechanical interactions between cells and their microenvironment play an important role in determining cell fate, which is particularly relevant in metastasis, a process where cells invade tissue matrices with different mechanical properties. In vitro, type I collagen hydrogels have been commonly used for modeling the microenvironment due to its ubiquity in the human body. In this work, the combined influence of the stiffness of these hydrogels and their ultrastructure on the migration patterns of HCT-116 and HT-29 spheroids are analyzed. For this, six different types of pure type I collagen hydrogels by changing the collagen concentration and the gelation temperature are prepared. The stiffness of each sample is measured and its ultrastructure is characterized. Cell migration studies are then performed by seeding the spheroids in three different spatial conditions. It is shown that changes in the aforementioned parameters lead to differences in the mechanical stiffness of the matrices as well as the ultrastructure. These differences, in turn, lead to distinct cell migration patterns of HCT-116 and HT-29 spheroids in either of the spatial conditions tested. Based on these results, it is concluded that the stiffness and the ultrastructural organization of the matrix can actively modulate cell migration behavior in colorectal cancer spheroids.

A Mathematical Modelling Study of Chemotactic Dynamics in Cell Cultures: The Impact of Spatio-temporal Heterogeneity

As motivated by studies of cellular motility driven by spatiotemporal chemotactic gradients in microdevices, we develop a framework for constructing approximate analytical solutions for the location, speed and cellular densities for cell chemotaxis waves in heterogeneous fields of chemoattractant from the underlying partial differential equation models. In particular, such chemotactic waves are not in general translationally invariant travelling waves, but possess a spatial variation that evolves in time, and may even oscillate back and forth in time, according to the details of the chemotactic gradients. The analytical framework exploits the observation that unbiased cellular diffusive flux is typically small compared to chemotactic fluxes and is first developed and validated for a range of exemplar scenarios. The framework is subsequently applied to more complex models considering the chemoattractant dynamics under more general settings, potentially including those of relevance for representing pathophysiology scenarios in microdevice studies. In particular, even though solutions cannot be constructed in all cases, a wide variety of scenarios can be considered analytically, firstly providing global insight into the important mechanisms and features of cell motility in complex spatiotemporal fields of chemoattractant. Such analytical solutions also provide a means of rapid evaluation of model predictions, with the prospect of application in computationally demanding investigations relating theoretical models and experimental observation, such as Bayesian parameter estimation.

Modelling cell adaptation using internal variables: Accounting for cell plasticity in continuum mathematical biology

Cellular adaptation is the ability of cells to change in response to different stimuli and environmental conditions. It occurs via phenotypic plasticity, that is, changes in gene expression derived from changes in the physiological environment. This phenomenon is important in many biological processes, in particular in cancer evolution and its treatment. Therefore, it is crucial to understand the mechanisms behind it. Specifically, the emergence of the cancer stem cell phenotype, showing enhanced proliferation and invasion rates, is an essential process in tumour progression. We present a mathematical framework to simulate phenotypic heterogeneity in different cell populations as a result of their interaction with chemical species in their microenvironment, through a continuum model using the well-known concept of internal variables to model cell phenotype. The resulting model, derived from conservation laws, incorporates the relationship between the phenotype and the history of the stimuli to which cells have been subjected, together with the inheritance of that phenotype. To illustrate the model capabilities, it is particularised for glioblastoma adaptation to hypoxia. A parametric analysis is carried out to investigate the impact of each model parameter regulating cellular adaptation, showing that it permits reproducing different trends reported in the scientific literature. The framework can be easily adapted to any particular problem of cell plasticity, with the main limitation of having enough cells to allow working with continuum variables. With appropriate calibration and validation, it could be useful for exploring the underlying processes of cellular adaptation, as well as for proposing favourable/unfavourable conditions or treatments.

A mechanobiological model for tumor spheroid evolution with application to glioblastoma: A continuum multiphysics approach

BACKGROUND

Spheroids are in vitro quasi-spherical structures of cell aggregates, eventually cultured within a hydrogel matrix, that are used, among other applications, as a technological platform to investigate tumor formation and evolution. Several interesting features can be replicated using this methodology, such as cell communication mechanisms, the effect of gradients of nutrients, or the creation of realistic 3D biological structures. The main objective of this work is to link the spheroid evolution with the mechanical activity of cells, coupled with nutrient consumption and the subsequent cell dynamics.

METHOD

We propose a continuum mechanobiological model which accounts for the most relevant phenomena that take place in tumor spheroid evolution under in vitro suspension, namely, nutrient diffusion in the spheroid, kinetics of cellular growth and death, and mechanical interactions among the cells. The model is qualitatively validated, after calibration of the model parameters, versus in vitro experiments of spheroids of different glioblastoma cell lines.

RESULTS

Our model is able to explain in a novel way quite different setups, such as spheroid growth (up to six times the initial configuration for U-87 MG cell line) or shrinking (almost half of the initial configuration for U-251 MG cell line); as the result of the mechanical interplay of cells driven by cellular evolution.

CONCLUSIONS

Glioblastoma tumor spheroid evolution is driven by mechanical interactions of the cell aggregate and the dynamical evolution of the cell population. All this information can be used to further investigate mechanistic effects in the evolution of tumors and their role in cancer disease.

Fabrication of human myocardium using multidimensional modelling of engineered tissues

Biofabrication of human tissues has seen a meteoric growth triggered by recent technical advancements such as human induced pluripotent stem cells (hiPSCs) and additive manufacturing. However, generation of cardiac tissue is still hampered by lack of addequate mechanical properties and crucially by the often unpredictable post-fabrication evolution of biological components. In this study we employ melt electrowriting (MEW) and hiPSC-derived cardiac cells to generate fibre-reinforced human cardiac minitissues. These are thoroughly characterized in order to build computational models and simulations able to predict their post-fabrication evolution. Our results show that MEW-based human minitissues display advanced maturation 28 post-generation, with a significant increase in the expression of cardiac genes such as MYL2, GJA5, SCN5A and the MYH7/MYH6 and MYL2/MYL7 ratios. Human iPSC-cardiomyocytes are significantly more aligned within the MEW-based 3D tissues, as compared to conventional 2D controls, and also display greater expression of CX43. These are also correlated with a more mature functionality in the form of faster conduction velocity. We used these data to develop simulations capable of accurately reproducing the experimental performance. In-depth gauging of the structural disposition (cellular alignment) and intercellular connectivity (CX43) allowed us to develop an improved computational model able to predict the relationship between cardiac cell alignment and functional performance. This study lays down the path for advancing in the development of in silico tools to predict cardiac biofabricated tissue evolution after generation, and maps the route towards more accurate and biomimetic tissue manufacture.

Evidence-Based Selection on the Appropriate FIT Cut-Off Point in CRC Screening Programs in the COVID Pandemic

Background: The COVID pandemic has forced the closure of many colorectal cancer (CRC) screening programs. Resuming these programs is a priority, but fewer colonoscopies may be available. We developed an evidence-based tool for decision-making in CRC screening programs, based on a fecal hemoglobin immunological test (FIT), to optimize the strategy for screening a population for CRC.

Methods: We retrospectively analyzed data collected at a regional CRC screening program between February/2014 and November/2018. We investigated two different scenarios: not modifying vs. modifying the FIT cut-off value. We estimated program outcomes in the two scenarios by evaluating the numbers of cancers and adenomas missed or not diagnosed in due time (delayed).

Results: The current FIT cut-off (20-μg hemoglobin/g feces) led to 6,606 colonoscopies per 100,000 people invited annually. Without modifying this FIT cut-off value, when the optimal number of individuals invited for colonoscopies was reduced by 10–40%, a high number of CRCs and high-risk adenomas (34–135 and 73–288/100.000-people invited, respectively) will be undetected every year. When the FIT cut-off value was increased to where the colonoscopy demand matched the colonoscopy availability, the number of missed lesions per year was remarkably reduced (9–36 and 29–145/100.000 people, respectively). Moreover, the unmodified FIT scenario outcome was improved by prioritizing the selection process based on sex (males) and age, rather than randomly reducing the number invited.

Conclusions: Assuming a mismatch between the availability and demand for annual colonoscopies, increasing the FIT cut-off point was more effective than randomly reducing the number of people invited. Using specific risk factors to prioritize access to colonoscopies should be also considered.

Predicting cell behaviour parameters from glioblastoma on a chip images. A deep learning approach

The broad possibilities offered by microfluidic devices in relation to massive data monitoring and acquisition open the door to the use of deep learning technologies in a very promising field: cell culture monitoring. In this work, we develop a methodology for parameter identification in cell culture from fluorescence images using Convolutional Neural Networks (CNN). We apply this methodology to the in vitro study of glioblastoma (GBM), the most common, aggressive and lethal primary brain tumour. In particular, the aim is to predict the three parameters defining the go or grow GBM behaviour, which is determinant for the tumour prognosis and response to treatment. The data used to train the network are obtained from a mathematical model, previously validated with in vitro experimental results. The resulting CNN provides remarkably accurate predictions (Pearson's ρ > 0.99 for all the parameters). Besides, it proves to be sound, to filter noise and to generalise. After training and validation with synthetic data, we predict the parameters corresponding to a real image of a microfluidic experiment. The obtained results show good performance of the CNN. The proposed technique may set the first steps towards patient-specific tools, able to predict in real-time the tumour evolution for each particular patient, thanks to a combined in vitro-in silico approach.

Force Spectroscopy Imaging and Constriction Assays Reveal the Effects of Graphene Oxide on the Mechanical Properties of Alginate Microcapsules

Microencapsulation of cells in hydrogel-based porous matrices is an approach that has demonstrated great success in regenerative cell therapy. These microcapsules work by concealing the exogenous cells and materials in a robust biomaterial that prevents their recognition by the immune system. A vast number of formulations and additives are continuously being tested to optimize cell viability and mechanical properties of the hydrogel. Determining the effects of new microcapsule additives is a lengthy process that usually requires extensive in vitro and in vivo testing. In this paper, we developed a workflow using nanoindentation (i.e., indentation with a nanoprobe in an atomic force microscope) and a custom-built microfluidic constriction device to characterize the effect of graphene oxide (GO) on three microcapsule formulations. With our workflow, we determined that GO modifies the microcapsule stiffness and surface properties in a formulation-dependent manner. Our results also suggest, for the first time, that GO alters the conformation of the microcapsule hydrogel and its interaction with subsequent coatings. Overall, our workflow can infer the effects of new additives on microcapsule surfaces. Thus, our workflow can contribute to diminishing the time required for the validation of new microcapsule formulations and accelerate their clinical translation.

Mathematical formulation and parametric analysis of in vitro cell models in microfluidic devices: application to different stages of glioblastoma evolution

In silico models and computer simulation are invaluable tools to better understand complex biological processes such as cancer evolution. However, the complexity of the biological environment, with many cell mechanisms in response to changing physical and chemical external stimuli, makes the associated mathematical models highly non-linear and multiparametric. One of the main problems of these models is the determination of the parameters’ values, which are usually fitted for specific conditions, making the conclusions drawn difficult to generalise. We analyse here an important biological problem: the evolution of hypoxia-driven migratory structures in Glioblastoma Multiforme (GBM), the most aggressive and lethal primary brain tumour. We establish a mathematical model considering the interaction of the tumour cells with oxygen concentration in what is called the go or grow paradigm. We reproduce in this work three different experiments, showing the main GBM structures (pseudopalisade and necrotic core formation), only changing the initial and boundary conditions. We prove that it is possible to obtain versatile mathematical tools which, together with a sound parametric analysis, allow to explain complex biological phenomena. We show the utility of this hybrid “biomimetic in vitro-in silico” platform to help to elucidate the mechanisms involved in cancer processes, to better understand the role of the different phenomena, to test new scientific hypotheses and to design new data-driven experiments.

Data-Driven Computational Simulation in Bone Mechanics

The data-driven approach was formally introduced in the field of computational mechanics just a few years ago, but it has gained increasing interest and application as disruptive technology in many other fields of physics and engineering. Although the fundamental bases of the method have been already settled, there are still many challenges to solve, which are often inherently linked to the problem at hand. In this paper, the data-driven methodology is applied to a particular problem in tissue biomechanics, a context where this approach is particularly suitable due to the difficulty in establishing accurate and general constitutive models, due to the intrinsic intra and inter-individual variability of the microstructure and associated mechanical properties of biological tissues. The problem addressed here corresponds to the characterization and mechanical simulation of a piece of cortical bone tissue. Cortical horse bone tissue was mechanically tested using a biaxial machine. The displacement field was obtained by means of digital image correlation and then transformed into strains by approximating the displacement derivatives in the bone virtual geometric image. These results, together with the approximated stress state, assumed as uniform in the small pieces tested, were used as input in the flowchart of the data-driven methodology to solve several numerical examples, which were compared with the corresponding classical model-based fitted solution. From these results, we conclude that the data-driven methodology is a useful tool to directly simulate problems of biomechanical interest without the imposition (model-free) of complex spatial and individually-varying constitutive laws. The presented data-driven approach recovers the natural spatial variation of the solution, resulting from the complex structure of bone tissue, i.e. heterogeneity, microstructural hierarchy and multifactorial architecture, making it possible to add the intrinsic stochasticity of biological tissues into the data set and into the numerical approach.

Force spectroscopy-based simultaneous topographical and mechanical characterization to study polymer-to-polymer interactions in coated alginate microspheres

Cell-laden hydrogel microspheres have shown encouraging outcomes in the fields of drug delivery, tissue engineering or regenerative medicine. Beyond the classical single coating with polycations, many other different coating designs have been reported with the aim of improving mechanical properties and in vivo performance of the microspheres. Among the most common strategies are the inclusion of additional polycation coatings and the covalent bonding of the semi-permeable membranes with biocompatible crosslinkers such as genipin. However, it remains challenging to characterize the effects of the interactions between the polycations and the hydrogel microspheres over time in vitro. Here we use a force spectroscopy-based simultaneous topographical and mechanical characterization to study polymer-to-polymer interactions in alginate microspheres with different coating designs, maintaining the hydrogels in liquid. In addition to classical topography parameters, we explored, for the first time, the evolution of peak/valley features along the z axis via thresholding analysis and the cross-correlation between topography and stiffness profiles with resolution down to tens of nanometers. Thus, we demonstrated the importance of genipin crosslinking to avoid membrane detachment in alginate microspheres with double polycation coatings. Overall, this methodology could improve hydrogel design rationale and expedite in vitro characterization, therefore facilitating clinical translation of hydrogel-based technologies.

Computational Multiscale Solvers for Continuum Approaches

Computational multiscale analyses are currently ubiquitous in science and technology. Different problems of interest-e.g., mechanical, fluid, thermal, or electromagnetic-involving a domain with two or more clearly distinguished spatial or temporal scales, are candidates to be solved by using this technique. Moreover, the predictable capability and potential of multiscale analysis may result in an interesting tool for the development of new concept materials, with desired macroscopic or apparent properties through the design of their microstructure, which is now even more possible with the combination of nanotechnology and additive manufacturing. Indeed, the information in terms of field variables at a finer scale is available by solving its associated localization problem. In this work, a review on the algorithmic treatment of multiscale analyses of several problems with a technological interest is presented. The paper collects both classical and modern techniques of multiscale simulation such as those based on the proper generalized decomposition (PGD) approach. Moreover, an overview of available software for the implementation of such numerical schemes is also carried out. The availability and usefulness of this technique in the design of complex microstructural systems are highlighted along the text. In this review, the fine, and hence the coarse scale, are associated with continuum variables so atomistic approaches and coarse-graining transfer techniques are out of the scope of this paper.

Glioblastoma on a microfluidic chip: Generating pseudopalisades and enhancing aggressiveness through blood vessel obstruction events

Background

Glioblastoma (GBM) is one of the most lethal tumor types. Hypercellular regions, named pseudopalisades, are characteristic in these tumors and have been hypothesized to be waves of migrating glioblastoma cells. These "waves" of cells are thought to be induced by oxygen and nutrient depletion caused by tumor-induced blood vessel occlusion. Although the universal presence of these structures in GBM tumors suggests that they may play an instrumental role in GBM's spread and invasion, the recreation of these structures in vitro has remained challenging.

Methods

Here we present a new microfluidic model of GBM that mimics the dynamics of pseudopalisade formation. To do this, we embedded U-251 MG cells within a collagen hydrogel in a custom-designed microfluidic device. By controlling the medium flow through lateral microchannels, we can mimic and control blood-vessel obstruction events associated with this disease.

Results

Through the use of this new system, we show that nutrient and oxygen starvation triggers a strong migratory process leading to pseudopalisade generation in vitro. These results validate the hypothesis of pseudopalisade formation and show an excellent agreement with a systems-biology model based on a hypoxia-driven phenomenon.

Conclusions

This paper shows the potential of microfluidic devices as advanced artificial systems capable of modeling in vivo nutrient and oxygen gradients during tumor evolution.

Mechanical Stress Redistribution in the First Metatarsal Bone After Autologous Bone Harvesting

BACKGROUND

The first metatarsal bone is a viable source for autologous bone grafting in foot and ankle surgery and may serve as another convenient graft site to correct a flail toe deformity. We aimed to determine how progressive bone removal from the first metatarsal affects the mechanical redistribution of the foot and whether this bone removal increases the risk of fracture.

METHODS

A three-dimensional finite element model developed from computed tomographic images obtained from a healthy man were used to evaluate traction stresses on the first metatarsal bone as a function of applied loads on the talus and Achilles tendon at two phases of the gait cycle (and according to the depth of bone removal).

RESULTS

Simulations indicated that when maximum load was applied to the Achilles tendon, tensile stress increased from 2.049 MPa in the intact foot to 5.941 MPa in the area of maximum bone harvest during the stance phase. Furthermore, as the volume of bone extracted from the first metatarsal increased, there was a redistribution of stress that differed significantly from that of the intact foot.

CONCLUSIONS

Although the maximum stress on the first metatarsal was not significantly affected by increasing the volume of bone harvested, the ankle should be splinted in plantarflexion during the postoperative period to eliminate the stance phase of gait and reduce the risk of metatarsal fracture.

Inhomogeneous Response of Articular Cartilage: A Three-Dimensional Multiphasic Heterogeneous Study

Articular cartilage exhibits complex mechano-electrochemical behaviour due to its anisotropy, inhomogeneity and material non-linearity. In this work, the thickness and radial dependence of cartilage properties are incorporated into a 3D mechano-electrochemical model to explore the relevance of heterogeneity in the behaviour of the tissue. The model considers four essential phenomena: (i) osmotic pressure, (ii) convective and diffusive processes, (iii) chemical expansion and (iv) three-dimensional through-the-thickness heterogeneity of the tissue. The need to consider heterogeneity in computational simulations of cartilage behaviour and in manufacturing biomaterials mimicking this tissue is discussed. To this end, healthy tibial plateaus from pigs were mechanically and biochemically tested in-vitro. Heterogeneous properties were included in the mechano-electrochemical computational model to simulate tissue swelling. The simulation results demonstrated that swelling of the heterogeneous samples was significantly lower than swelling under homogeneous and isotropic conditions. Furthermore, there was a significant reduction in the flux of water and ions in the former samples. In conclusion, the computational model presented here can be considered as a valuable tool for predicting how the variation of cartilage properties affects its behaviour, opening up possibilities for exploring the requirements of cartilage-mimicking biomaterials for tissue engineering. Besides, the model also allows the establishment of behavioural patterns of swelling and of water and ion fluxes in articular cartilage.

Parameter-dependent behavior of articular cartilage: 3D mechano-electrochemical computational model

BACKGROUND AND OBJECTIVE

Changes in mechano-electrochemical properties of articular cartilage play an essential role in the majority of cartilage diseases. Despite of this importance, the specific effect of each parameter into tissue behavior remains still obscure. Parametric computational modeling of cartilage can provide some insights into this matter, specifically the study of mechano-electrochemical properties variation and their correlation with tissue swelling, water and ion fluxes. Thus, the aim of this study is to evaluate the influence of the main mechanical and electrochemical parameters on the determination of articular cartilage behavior by a parametric analysis through a 3D finite element model.

METHODS

For this purpose, a previous 3D mechano-electrochemical model, developed by the same authors, of articular cartilage behavior has been used. Young's modulus, Poisson coefficient, ion diffusivities and ion activity coefficients variations have been analyzed and quantified through monitoring tissue simulated response.

RESULTS

Simulation results show how Young's modulus and Poisson coefficient control tissue behavior rather than electrochemical properties. Meanwhile, ion diffusivity and ion activity coefficients appear to be vital in controlling velocity of incoming and outgoing fluxes.

CONCLUSIONS

This parametric study establishes a basic guide when defining the main properties that are essential to be included into computational modeling of articular cartilage providing a helpful tool in tissue simulations.

Abstract B04: An in vitro model for glioblastoma using microfluidics: Generating pseudopalisades on a chip

Glioblastoma (GBM), also named grade IV astrocytoma, is the most common and lethal malignant primary brain tumor. GBM is characterized by two main histopathological conditions: necrotic foci typically surrounded by areas of high cellularity and microvascular proliferation.

The causes of such densely populated regions, known a “pseudopalisades”, remain poorly understood. Firstly, high cellularity was thought to be due to the rapid proliferation of GBM cells, however recent histological studies have shown that proliferation in pseudopalisading areas is significantly lower than in adjacent regions. Additionally, in pseudopalisades, apoptosis is substantially larger than in neighboring regions. These evidences suggest that pseudopalisades are due to other causes than simply higher proliferation or survival rates. Recent reports have pointed out that, in histological slices from GBM patients, more than 50% of pseudopalisades present clearly a central obstructed blood vessel and in more than 90% microscopic evidences of thrombosis are observed.

As a consequence, oxygen and nutrients supply is compromised around the surrounding zones. It has been proposed that one of the driven forces of glioma aggressiveness is these nutrient and oxygen starvation. According to this hypothesis, GBM cell proliferation and secretion of pro-coagulant signals would cause thrombotic events, leading to hypoxia and nutrient depletion in the microenvironment. As a consequence, cells migrate towards nutrients and oxygen enriched regions, creating the characteristic GBM pseudopalisades. Eventually, these migrating cells would reach other blood vessels and eventually cause the collapse of these vessels, restarting the process and creating an expansive wave of tumor cells across the brain.

This complex process is hardly reproducible using “in vitro” models because the conventional migration assays are unable to mimic the complex microenvironment described. Recently, microfabrication and microfluidic technologies have arisen as interesting alternatives for creating high-performance cell culture systems. In this paper we describe the design, fabrication and biological validation of a microfluidic device using SU-8 technology for three-dimensional GBM cell cultures under thrombotic conditions. The fabricated microdevice possesses a central microchamber to locate the cells embedded within a hydrogel, mimicking the ECM and allowing migration in three dimensions. On both sides of the microchamber, two lateral microchannels are filled with culture medium, allowing nutrient and oxygen diffusion towards the microchamber, mimicking the brain blood vessels.

In this study we have demonstrated the GBM cells (U-251-MG) viability within the microdevice. Moreover, controlling medium flow through lateral microchannels we can mimic the GBM-associated thrombotic pathophysiological conditions. Under these thrombotic circumstances nutrient starvation leads to a chemotactic process and the formation of a migratory front similar to the pseudopalisades observed in vivo and validated with mathematical algorithms. Moreover, our results suggest that the whole process stimulates a more aggressive behavior of GBM cells. In an early stage when nutrients are plenty, GBM cells remain in a non-invasive phenotype. When nutrients are depleted, GBM cells initiate a migration process towards enriched regions, leading to the pseudopalisade formation. When compared with patient´s tumor pathology images, GBM cell invasion behavior follow the same pattern observed in vivo. This novel technique could help us to understand the mechanism of the pseudopalisade formation and to suggest novel therapeutic targets to avoid tumor progression. Besides, these microfluidic devices could represent an extremely useful platform to test new anticancer agents in a preclinical setting that mimics the complex GBM microenvironment.

Citation Format: Jose M. Ayuso, Rosa Monge, Alicia Martínez-González, Guillermo A. Llamazares, Javier Berganzo, Aurelio Hernández-Laín, Jorge Santolaria, Manuel Doblaré, Pilar Sánchez-Gómez, Víctor M. Pérez-García, Ignacio Ochoa, Luis J. Fernández. An in vitro model for glioblastoma using microfluidics: Generating pseudopalisades on a chip. [abstract]. In: Proceedings of the AACR Special Conference: Advances in Brain Cancer Research; May 27-30, 2015; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2015;75(23 Suppl):Abstract nr B04.

Reduction of Dorsal Displacement of the Proximal and Middle Phalanges Using a Neutral or Angled Implant for Joint Arthrodesis to Treat Hammertoe Deformity A Finite Element Study

BACKGROUND

We used finite element analysis to study the mechanical displacements at three planes of the second through fourth hammertoes during the push-off phase of gait using a new neutral or 10° angled memory alloy intramedullary implant (FDA K070598) used for proximal interphalangeal joint arthrodesis.

METHODS

After geometric reconstruction of the foot skeleton from computed tomographic images of a 36-year-old man, an intramedullary implant was positioned in the virtual model at the neutral and 10° angled positions at the proximal interphalangeal joints of the second through fourth hammertoes during the push-off phase of gait. The obtained displacement results in three planes were compared with those derived from the nonsurgical foot model using finite element analysis.

RESULTS

These results support the successful use of either a neutral or angled implant for proximal interphalangeal joint arthrodesis, with the neutral implant yielding slightly better results.

CONCLUSIONS

The neutral implant reduced vertical displacement to a greater extent than did the angled implant. We also highlight the potential risk of iatrogenic curly toe when performing a proximal interphalangeal joint arthrodesis using an angled implant specifically at the fourth toe.

Altered swelling and ion fluxes in articular cartilage as a biomarker in osteoarthritis and joint immobilization: a computational analysis

In healthy cartilage, mechano-electrochemical phenomena act together to maintain tissue homeostasis. Osteoarthritis (OA) and degenerative diseases disrupt this biological equilibrium by causing structural deterioration and subsequent dysfunction of the tissue. Swelling and ion flux alteration as well as abnormal ion distribution are proposed as primary indicators of tissue degradation. In this paper, we present an extension of a previous three-dimensional computational model of the cartilage behaviour developed by the authors to simulate the contribution of the main tissue components in its behaviour. The model considers the mechano-electrochemical events as concurrent phenomena in a three-dimensional environment. This model has been extended here to include the effect of repulsion of negative charges attached to proteoglycans. Moreover, we have studied the fluctuation of these charges owning to proteoglycan variations in healthy and pathological articular cartilage. In this sense, standard patterns of healthy and degraded tissue behaviour can be obtained which could be a helpful diagnostic tool. By introducing measured properties of unhealthy cartilage into the computational model, the severity of tissue degeneration can be predicted avoiding complex tissue extraction and subsequent in vitro analysis. In this work, the model has been applied to monitor and analyse cartilage behaviour at different stages of OA and in both short (four, six and eight weeks) and long-term (11 weeks) fully immobilized joints. Simulation results showed marked differences in the corresponding swelling phenomena, in outgoing cation fluxes and in cation distributions. Furthermore, long-term immobilized patients display similar swelling as well as fluxes and distribution of cations to patients in the early stages of OA, thus, preventive treatments are highly recommended to avoid tissue deterioration.

Structural biology response of a collagen hydrogel synthetic extracellular matrix with embedded human fibroblast: computational and experimental analysis

Adherent cells exert contractile forces which play an important role in the spatial organization of the extracellular matrix (ECM). Due to these forces, the substrate experiments a volume reduction leading to a characteristic shape. ECM contraction is a key process in many biological processes such as embryogenesis, morphogenesis and wound healing. However, little is known about the specific parameters that control this process. With this aim, we present a 3D computational model able to predict the contraction process of a hydrogel matrix due to cell-substrate mechanical interaction. It considers cell-generated forces, substrate deformation, ECM density, cellular migration and proliferation. The model also predicts the cellular spatial distribution and concentration needed to reproduce the contraction process and confirms the minimum value of cellular concentration necessary to initiate the process observed experimentally. The obtained continuum formulation has been implemented in a finite element framework. In parallel, in vitro experiments have been performed to obtain the main model parameters and to validate it. The results demonstrate that cellular forces, migration and proliferation are acting simultaneously to display the ECM contraction.

Evaluation of the stiffnesses of the Achilles tendon and soleus from the apparent stiffness of the triceps surae

The triceps surae plays an important role in the performance of many sports. Although the apparent average mechanical properties of the triceps surae may be a satisfactory parameter for estimating the training level of an athlete, a knowledge of the mechanical properties of the individual constituents of the triceps surae (in particular the Achilles tendon and soleus) permits a more detailed and in-depth control of the effects of training from more physically based parameters. The objective of this work is therefore the estimation of the individual viscoelastic properties (stiffness and viscosity) of soleus and Achilles tendon from the apparent properties of the triceps surae obtained by free vibration techniques. Different procedures have been developed and discussed, showing a high degree of robustness in the predictions. The results obtained for a non-oriented set of subjects present a high level of variability, depending on the training conditions and anthropometric features, although the corresponding average values compare well with data previously reported in the literature, particularly those associated with the tendon stiffness.

Zeolites for the selective adsorption of sulfur hexafluoride

Molecular simulations have been used to investigate at the molecular level the suitability of zeolites with different topology on the adsorption, diffusion and separation of a nitrogen–sulfur hexafluoride mixture containing the latter at low concentration.

Stress at the second metatarsal bone after correction of hammertoe and claw toe deformity: a finite element analysis using an anatomical model

BACKGROUND

We used finite element analysis to evaluate three techniques for the correction of hammertoe and claw toe deformities: flexor digitorum longus tendon transfer (FDLT), flexor digitorum brevis tendon transfer (FDBT), and proximal interphalangeal joint arthrodesis (PIPJA).

METHODS

We performed a finite element analysis of FDLT and FDBT compared with PIPJA of the second toe using multislice computed tomography and 93 tomographic images of the foot obtained in a healthy 36-year-old man.

RESULTS

The PIPJA showed a significantly higher increase in traction and compressive stresses and strain at the medial aspect of the shaft of the second metatarsal bone compared with FDLT or FDBT (P < .01). Mean ± SD compressive stresses increased to -4.35 ± 7.05 MPa compared with the nonsurgical foot (-3.10 ± 4.90 MPa). It can, therefore, be hypothesized that if PIPJA is used to correct the hammertoe and claw toe deformities, it could also increase traction and compressive stresses and strain in the metatarsals during running and other vigorous activities.

CONCLUSIONS

There is a biomechanical advantage to performing FDLT or FDBT instead of PIPJA to surgically treat a hammertoe or claw toe deformity. In addition, tensile strain at the dorsal aspect of the second metatarsal bone when performing PIPJA increases the risk of metatarsalgia or stress fracture in patients at risk.

Computational modelling of multi-cell migration in a multi-signalling substrate

Cell migration is a vital process in many biological phenomena ranging from wound healing to tissue regeneration. Over the past few years, it has been proven that in addition to cell-cell and cell-substrate mechanical interactions (mechanotaxis), cells can be driven by thermal, chemical and/or electrical stimuli. A numerical model was recently presented by the authors to analyse single cell migration in a multi-signalling substrate. That work is here extended to include multi-cell migration due to cell-cell interaction in a multi-signalling substrate under different conditions. This model is based on balancing the forces that act on the cell population in the presence of different guiding cues. Several numerical experiments are presented to illustrate the effect of different stimuli on the trajectory and final location of the cell population within a 3D heterogeneous multi-signalling substrate. Our findings indicate that although multi-cell migration is relatively similar to single cell migration in some aspects, the associated behaviour is very different. For instance, cell-cell interaction may delay single cell migration towards effective cues while increasing the magnitude of the average net cell traction force as well as the local velocity. Besides, the random movement of a cell within a cell population is slightly greater than that of single cell migration. Moreover, higher electrical field strength causes the cell slug to flatten near the cathode. On the other hand, as with single cell migration, the existence of electrotaxis dominates mechanotaxis, moving the cells to the cathode or anode pole located at the free surface. The numerical results here obtained are qualitatively consistent with related experimental works.

Epicardial delivery of collagen patches with adipose-derived stem cells in rat and minipig models of chronic myocardial infarction

Although transplantation of adipose-derived stem cells (ADSC) in chronic myocardial infarction (MI) models is associated with functional improvement, its therapeutic value is limited due to poor long-term cell engraftment and survival. Thus, the objective of this study was to examine whether transplantation of collagen patches seeded with ADSC could enhance cell engraftment and improve cardiac function in models of chronic MI. With that purpose, chronically infarcted Sprague-Dawley rats (n = 58) were divided into four groups and transplanted with media, collagen scaffold (CS), rat ADSC, or CS seeded with rat ADSC (CS-rADSC). Cell engraftment, histological changes, and cardiac function were assessed 4 months after transplantation. In addition, Göttingen minipigs (n = 18) were subjected to MI and then transplanted 2 months later with CS or CS seeded with autologous minipig ADSC (CS-pADSC). Functional and histological assessments were performed 3 months post-transplantation. Transplantation of CS-rADSC was associated with increased cell engraftment, significant improvement in cardiac function, myocardial remodeling, and revascularization. Moreover, transplantation of CS-pADSC in the pre-clinical swine model improved cardiac function and was associated with decreased fibrosis and increased vasculogenesis. In summary, transplantation of CS-ADSC resulted in enhanced cell engraftment and was associated with a significant improvement in cardiac function and myocardial remodeling.

Culture of human bone marrow-derived mesenchymal stem cells on of poly(L-lactic acid) scaffolds: potential application for the tissue engineering of cartilage

PURPOSE

Due to the attractive properties of poly(L-lactic acid) (PLLA) for tissue engineering, the aim was to determine the growth and differentiation capacity of mesenchymal stromal cells (MSCs) in PLLA scaffolds and their potential use in the treatment of cartilage diseases.

METHODS

MSCs were cultured in PLLA films and thin porous membranes to study adherence and proliferation. Permeability and porosity were determined for the different scaffolds employed. The optimal conditions for cell seeding were first determined, as well as cell density and distribution inside the PLLA. Scaffolds were then maintained in expansion or chondrogenic differentiation media for 21 days. Apoptosis, proliferation and chondrogenic differentiation was assessed after 21 days in culture by immunohistochemistry. Mechanical characteristics of scaffolds were determined before and after cell seeding.

RESULTS

MSCs uniformly adhered to PLLA films as well as to porous membranes. Proliferation was detected only in monolayers of pure PLLA, but was no longer detected after 10 days. Mechanical characterization of PLLA scaffolds showed differences in the apparent compression elastic modulus for the two sizes used. After determining high efficiencies of seeding, the production of extracellular matrix (ECM) was determined and contained aggrecan and collagens type I and X. ECM produced by the cells induced a twofold increase in the apparent elastic modulus of the composite.

CONCLUSIONS

Biocompatible PLLA scaffolds have been developed that can be efficiently loaded with MSCs. The scaffold supports chondrogenic differentiation and ECM deposition that improves the mechanics of the scaffold. Although this improvement does not met the expectations of a hyaline-like cartilage ECM, in part due to the lack of a mechanical stimulation, their potential use in the treatment of cartilage pathologies encourages to improve the mechanical component.

3D computational modelling of cell migration: a mechano-chemo-thermo-electrotaxis approach

Single cell migration constitutes a fundamental phenomenon involved in many biological events such as wound healing, cancer development and tissue regeneration. Several experiments have demonstrated that, besides the mechanical driving force (mechanotaxis), cell migration may be also influenced by chemical, thermal and/or electrical cues. In this paper, we present an extension of a previous model of the same authors adding the effects of chemotaxis, thermotaxis and electrotaxis to the initial mechanotaxis model of cell migration, allowing us to predict cell migration behaviour under different conditions and substrate properties. The present model is based on the balance of effective forces during cell motility in the presence of the several guiding cues. This model has been applied to several numerical experiments to demonstrate the effect of the different drivers onto the cell path and final location within a certain three-dimensional substrate with heterogeneous properties. Our findings indicate that the presence of the chemotaxis, thermotaxis and/or electrotaxis reduce, in general, the random component of cell movement, being this reduction more important in the case of electrotaxis that can be considered a dominating signal during cell migration (besides the underlying mechanical effects). These results are qualitatively in agreement with well-known experimental ones.

Preparation and characterization of collagen-based ADSC-carrier sheets for cardiovascular application

The use of scaffolds composed of natural biodegradable matrices represents an attractive strategy to circumvent the lack of cell engraftment, a major limitation of stem cell therapy in cardiovascular diseases. Bovine-derived non-porous collagen scaffolds with different degrees of cross-linking (C0, C2, C5 and C10) were produced and tested for their mechanical behavior, in vitro biocompatibility with adipose-derived stem cells (ADSCs) and tissue adhesion and inflammatory reaction. Uniaxial tensile tests revealed an anisotropic behavior of collagen scaffolds (2×0.5cm) and statistically significant differences in the mechanical behavior between cross-linked and non-cross-linked scaffolds (n=5). In vitro, ADSCs adhered homogenously and showed a similar degree of proliferation on all four types of scaffolds (cells×10(3)cm(-2) at day 7: C0: 94.7±37.1; C2: 91.7±25.6; C5: 88.2±6.8; C10: 72.8±10.7; P=n.s.; n=3). In order to test the in vivo biocompatibility, a chronic myocardial infarction model was performed in rats and 1.2×1.2cm size collagen scaffolds implanted onto the heart 1month post-infarction. Six animals per group were killed 2, 7 and 30days after transplant. Complete and long-lasting adhesion to the heart was only observed with the non-cross-linked scaffolds with almost total degradation 1month post-transplantation. After 7 and 30days post-implantation, the degree of inflammation was significantly lower in the hearts treated with non-cross-linked scaffolds (day 7: C0: 10.2±2.1%; C2: 16.3±2.9%; C5: 15.9±4.8%; C10: 17.4±4.1%; P<0.05 vs. C0; day 30: C0: 1.3±1.3%; C2: 9.4±3.0%; C5: 7.0±2.1%; C10: 9.8±2.5%; P<0.01 vs. C0). In view of the results, the non-cross-linked scaffold (C0) was chosen as an ADSC-carrier sheet and tested in vivo. One week post-implantation, 25.3±7.0% of the cells transplanted were detected in those animals receiving the cell-carrier sheet whereas no cells were found in animals receiving cells alone (n=3 animals/group). We conclude that the biocompatibility and mechanical properties of the non-cross-linked collagen scaffolds make them a useful cell carrier that greatly favors tissue cell engraftment and may be exploited for cell transplantation in models of cardiac disease.

Cell migration and cell-cell interaction in the presence of mechano-chemo-thermotaxis

Although there are several computational models that explain the trajectory that cells take during migration, till now little attention has been paid to the integration of the cell migration in a multi-signaling system. With that aim, a generalized model of cell migration and cell-cell interaction under multisignal environments is presented herein. In this work we investigate the spatio-temporal cell-cell interaction problem induced by mechano-chemo-thermotactic cues. It is assumed that formation of a new focal adhesion generates traction forces proportional to the stresses transmitted by the cell to the extracellular matrix. The cell velocity and polarization direction are calculated based on the equilibrium of the effective forces associated to cell motility. It is also assumed that, in addition to mechanotaxis signals, chemotactic and thermotactic cues control the direction of the resultant traction force. This model enables predicting the trajectory of migrating cells as well as the spatial and temporal distributions of the net traction force and cell velocity. Results indicate that the tendency of the cells is firstly to reach each other and then migrate towards an imaginary equilibrium plane located near the source of the signal. The position of this plane is sensitive to the gradient slope and the corresponding efficient factors. The cells come into contact and separate several times during migration. Adding other cues to the substrate (such as chemotaxis and/or thermotaxis) delays that primary contact. Moreover, in all states, the average local velocity and the net traction force of the cells decrease while the cells approach the cues source. Our findings are qualitatively consistent with experimental observations reported in the related literature.

Response of sheep chondrocytes to changes in substrate stiffness from 2 to 20 Pa: effect of cell passaging

AIM

The influence of culture substrate stiffness (in the kPa range) on chondrocyte behavior has been described. Here we describe the response to variations in substrate stiffness in a soft range (2-20 Pa), as it may play a role in understanding cartilage physiopathology.

METHODS

We developed a system for cell culture in substrates with different elastic moduli using collagen hydrogels and evaluated chondrocytes after 2, 4, and 7 days in monolayer and three-dimensional (3D) cultures. Experiments were performed in normoxia and hypoxia in order to describe the effect of a low oxygen environment on chondrocytes. Finally, we also evaluated if dedifferentiated cells preserve the capacity for mechanosensing.

RESULTS

Chondrocytes showed less proliferating activity when cultured in monolayer in the more compliant substrates. Expression of the cartilage markers Aggrecan (Acan), type II collagen (Col2a1), and Sox9 was upregulated in the less stiff gels (both in monolayer and in 3D culture). Stiffer gels induced an organization of the actin cytoskeleton that correlated with the loss of a chondrocyte phenotype. When cells were cultured in hypoxia, we observed changes in the cellular response that were mediated by HIF-1α. Results in 3D hypoxia cultures were opposite to those found in normoxia, but remained unchanged in monolayer hypoxic experiments. Similar results were found for dedifferentiated cells.

CONCLUSIONS

Chondrocytes respond differently according to the stiffness of the substrate. This response depends greatly on the oxygen environment and on whether the chondrocyte is embedded or grown onto the hydrogel, since mechanosensing capacity was not lost with cell expansion.

Computational modelling and analysis of mechanical conditions on cell locomotion and cell-cell interaction

Between other parameters, cell migration is partially guided by the mechanical properties of its substrate. Although many experimental works have been developed to understand the effect of substrate mechanical properties on cell migration, accurate 3D cell locomotion models have not been presented yet. In this paper, we present a novel 3D model for cells migration. In the presented model, we assume that a cell follows two main processes: in the first process, it senses its interface with the substrate to determine the migration direction and in the second process, it exerts subsequent forces to move. In the presented model, cell traction forces are considered to depend on cell internal deformation during the sensing step. A random protrusion force is also considered which may change cell migration direction and/or speed. The presented model was applied for many cases of migration of the cells. The obtained results show high agreement with the available experimental and numerical data.

Numerical simulation of bone remodelling around dental implants

Crestal bone loss can result in the failure of dental implants and can be caused, by among other factors, the development of non-physiological mechanical conditions. Bone remodelling (BR) is the physiological process through which bone adapts itself to the mechanical environment. A previously published mathematical model of BR is used in this work to study the homogenized structural evolution of peri-implant bone. This model is used to study the influence of the diameter and length of a dental implant of pure titanium on its long-term stability. The temporal evolution of porosity and microstructural damage of the peri-implant bone are the variables analysed in this study. The results show that damage and porosity increase as the implant length decreases and, more pronouncedly, as its diameter decreases. The increase in damage and porosity levels is localized, as many other studies confirm, at the implant neck due to the stress concentration that is created in that area. The main conclusion of this study is that in implants with a diameter equal to or greater than 3 mm the damage is under control and there is no mechanical failure of the peri-implant bone in the long term.

Mechanical behaviour of synthetic surgical meshes: finite element simulation of the herniated abdominal wall

The material properties of meshes used in hernia surgery contribute to the overall mechanical behaviour of the repaired abdominal wall. The mechanical response of a surgical mesh has to be defined since the haphazard orientation of an anisotropic mesh can lead to inconsistent surgical outcomes. This study was designed to characterize the mechanical behaviour of three surgical meshes (Surgipro®, Optilene® and Infinit®) and to describe a mechanical constitutive law that accurately reproduces the experimental results. Finally, through finite element simulation, the behaviour of the abdominal wall was modelled before and after surgical mesh implant. Uniaxial loading of mesh samples in two perpendicular directions revealed the isotropic response of Surgipro® and the anisotropic behaviour of Optilene® and Infinit®. A phenomenological constitutive law was used to reproduce the measured experimental curves. To analyze the mechanical effect of the meshes once implanted in the abdomen, finite element simulation of the healthy and partially herniated repaired rabbit abdominal wall served to reproduce wall behaviour before and after mesh implant. In all cases, maximal displacements were lower and maximal principal stresses higher in the implanted abdomen than the intact wall model. Despite the fact that no mesh showed a behaviour that perfectly matched that of abdominal muscle, the Infinit® mesh was able to best comply with the biomechanics of the abdominal wall.

FSI Analysis of a Human Trachea Before and After Prosthesis Implantation

In this work we analyzed the response of a stenotic trachea after a stent implantation. An endotracheal stent is the common treatment for tracheal diseases such as stenosis, chronic cough, or dispnoea episodes. Medical treatment and surgical techniques are still challenging due to the difficulties in overcoming potential complications after prosthesis implantation. A finite element model of a diseased and stented trachea was developed starting from a patient specific computerized tomography (CT) scan. The tracheal wall was modeled as a fiber reinforced hyperelastic material in which we modeled the anisotropy due to the orientation of the collagen fibers. Deformations of the tracheal cartilage rings and of the muscular membrane, as well as the maximum principal stresses, are analyzed using a fluid solid interaction (FSI) approach. For this reason, as boundary conditions, impedance-based pressure waveforms were computed modeling the nonreconstructed vessels as a binary fractal network. The results showed that the presence of the stent prevents tracheal muscle deflections and indicated a local recirculatory flow on the stent top surface which may play a role in the process of mucous accumulation. The present work gives new insight into clinical procedures, predicting their mechanical consequences. This tool could be used in the future as preoperative planning software to help the thoracic surgeons in deciding the optimal prosthesis type as well as its size and positioning.

Mechanical characterization and numerical simulation of polyether-ether-ketone (PEEK) cranial implants

Cranial implants have experienced a significant evolution in the last decade in different aspects such as materials, method of fixation, and the structure. In addition, patient-specific cranial implants have recently been started to be developed. To achieve this objective, efficient mechanical characterization and numerical modeling of the implant are required to guarantee its functionality on each patient as well as to facilitate further developments. In this work, mechanical characterization and numerical models have been performed for patient-specific Polyaryletherketone (PEEK) scaffold cranial implants. Mechanical characterization has been performed at the scaffold and the whole implant levels under displacement control tests. Two different implant designs for the same patient but with different scaffold structure were experimentally characterized, and finite element models of the implants were developed within the framework of linear elasticity. Two types of finite element models were developed: a detailed finite element model with the actual scaffold geometry, and a solid shell-like model with effective material properties. These effective material properties were obtained by means of the Asymptotic Expansion Homogenization (AEH) theory which accounts for the periodicity of the underlying structure of the material. Experimental results showed a linear response of the material and the implant up to failure, therefore supporting the use of linear elastic models for simulation. Numerical models showed excellent agreement with experiments regarding load-displacement response. Models also showed a very consistent behavior with regard to the location and the value of the maximum principal stress in the implant when subjected to the maximum load of the experiments. The two numerical models were compared. The homogenized model gave results that were very close to those obtained with the detailed model, while reducing the number of degrees of freedom by 90%, and therefore the overall computational burden. The results showed that the models are able to reproduce experimental results conducted on actual implants, offering a valid alternative to be used in the design of customized cranial implants with a scaffold structure.

Experimental study and constitutive modelling of the passive mechanical properties of the porcine carotid artery and its relation to histological analysis: Implications in animal cardiovascular device trials

The present study focusses on the determination, comparison and constitutive modelling of the passive mechanical properties of the swine carotid artery over very long stretches in both proximal and distal regions. Special attention is paid to the histological and mechanical variations of these properties depending on the proximity to the heart. The results can have clinical relevance, especially in the research field of intravascular device design. Before the final clinical trials on humans, research in the vascular area is conducted on animal models, swine being the most common due to the similarities between the human and swine cardiovascular systems as well as the fact that the swine size is suitable for testing devices, in this case endovascular carotid systems. The design of devices usually involves numerical techniques, and an important feature is the appropriate modelling of the mechanical properties of the vessel. Fourteen carotid swine arteries were harvested just after sacrifice and cyclic uniaxial tension tests in longitudinal and circumferential directions were performed for distal and proximal samples. The stress-stretch curves obtained were fitted with a hyperelastic anisotropic model. Stress-free configuration states were also analyzed. Finally, human and swine samples were processed in a histological laboratory and images were used to quantify their microconstituents. The statistical analysis revealed significant differences between the mechanical behavior of proximal and distal locations in the circumferential but not in the longitudinal direction. Circumferential direction samples show clear differences both in residual stretches and tensile curves between the two locations, while the features of longitudinal specimens are independent of the axial position. The statistical analysis provides significant evidence of changes depending on the position of the sample, mainly in elastin and SMC quantification.