Culture of 3D bioprinted bone constructs requires an increased fluid dynamic stimulation.

In vitro flow-induced mechanical stimulation of developing bone tissue constructs has been shown to favor mineral deposition in scaffolds seeded with cells directly exposed to the fluid flow. However, the effect of fluid dynamic parameters, such as shear stress (SS), within 3D bioprinted constructs is still unclear. Thus, this study aimed at correlating the SS levels and the mineral deposition in 3D bioprinted constructs, evaluating the possible dampening effect of the hydrogel. Human mesenchymal stem cells (hMSCs) were embedded in 3D bioprinted porous structures made of alginate and gelatin. 3D bioprinted constructs were cultured in an osteogenic medium assessing the influence of different flow rates (0, 0.7 and 7 ml/min) on calcium and collagen deposition through histology, and bone volume (BV) through micro-computed tomography. Uniform distribution of calcium and collagen was observed in all groups. Nevertheless, BV significantly increased in perfused groups as compared to static control, ranging from 0.35±0.28 mm3, 11.90±8.74 mm3 and 25.81±5.02 mm3 at week 3 to 2.28±0.78 mm3, 22.55±2.45 mm3 and 46.05±5.95 mm3 at week 6 in static, 0.7 and 7 ml/min groups, respectively. SS values on construct fibers in the range 10-100 mPa in 7 ml/min samples were twice as high as those in 0.7 ml/min samples showing the same trend of BV. The obtained results suggest that it is necessary to enhance the flow-induced mechanical stimulation of cell-embedding hydrogels to increase the amount of mineral deposited by hMSCs, compared to what is generally reported for the development of in vitro bone constructs. STATEMENT OF SIGNIFICANCE: In this study, we evaluated for the first time how the hydrogel structure dampens the effect of flow-induced mechanical stimulation during the culture of 3D bioprinted bone tissue constructs. By combining computational and experimental techniques we demonstrated that those shear stress thresholds generally considered for culturing cells seeded on scaffold surface, are no longer applicable when cells are embedded in 3D bioprinted constructs. Significantly, more bone volume was formed in constructs exposed to shear stress values generally considered as detrimental than in constructs exposed shear stress values generally considered as beneficial after 3 weeks and 6 weeks of dynamic culture using a perfusion bioreactor.

Haematocrit heterogeneity in blood flows past microfluidic models of oxygenating fibre bundles.

Blood oxygenators act as an extracorporeal artificial lung during certain types of cardiac surgery and intensive care therapies. Inside these devices, blood is forced to flow across an oxygenating bundle, encountering interstitial gaps comparable to those typical of the microvasculature. Despite the well-known effects of such length scales on haemorheology and red blood cell (RBC) behavior, these are generally overlooked in oxygenator modeling and design; it is persistently assumed that RBCs are homogeneously distributed throughout the oxygenating bundle, independently of their microstructure arrangement or main flow directions. The goal of this study is to provide preliminary experimental evidence of heterogeneous RBC distributions inside oxygenating fibre bundles. To this end, a number of microchannels were manufactured inspired by actual oxygenating devices, considering simplified versions of their microstructure. These comprise staggered arrays of micro pillars, which were perfused with RBC suspensions, with feed haematocrit (Ht) and velocities relevant for clinical use. The microchannels were imaged using a microscope and high-speed camera to accurately capture cell distribution. The imaged blood flows revealed the non-uniform nature of RBC distributions in the arrays, characterized by local Ht gradients particularly around the O2 sources inside the bundle. These heterogeneous distributions should be accounted for during oxygenator design, as RBC concentration plays a key role in O2 transport and, ultimately, overall device performance.

µ-Particle tracking velocimetry and computational fluid dynamics study of cell seeding within a 3D porous scaffold.

Cell seeding of 3D scaffolds is a critical step in tissue engineering since the final tissue properties are related to the initial cell distribution and density within the scaffold unit. Perfusion systems can transport cells to the scaffold however; the fact that cells flow inside the scaffold pores does not guarantee cell deposition onto the scaffold substrate and cell attachment. The aim of this study was to investigate how fluid flow conditions modulate cell motion and deposition during perfusion. For such a purpose, a multiphase-based computational fluid dynamics (CFD) model was developed in conjunction with particle tracking velocimetry experiments (PTV) which for the first time were applied to observe cell seeding inside a 3D scaffold. CFD and PTV results showed the strong effect of gravity for lower flow rates leading to cell sedimentation and poor transport of cells to the scaffold. Higher flow rates help overcome the effect of gravity so more cells travelling inside the scaffold were found. Nonetheless, fluid flow drags cells along the fluid streamlines without intercepting the scaffold substrate. As a consequence, if cells do not deposit into the scaffold substrate, cell adhesion cannot occur. Therefore, cell-scaffold interception should be promoted and the present computational model which predicts the effect of gravity and fluid drag on cells trajectories could serve to optimise bioreactors and enhance cell seeding efficiency.

2D µ-Particle Image Velocimetry and Computational Fluid Dynamics Study Within a 3D Porous Scaffold.

Transport properties of 3D scaffolds under fluid flow are critical for tissue development. Computational fluid dynamics (CFD) models can resolve 3D flows and nutrient concentrations in bioreactors at the scaffold-pore scale with high resolution. However, CFD models can be formulated based on assumptions and simplifications. µ-Particle image velocimetry (PIV) measurements should be performed to improve the reliability and predictive power of such models. Nevertheless, measuring fluid flow velocities within 3D scaffolds is challenging. The aim of this study was to develop a µPIV approach to allow the extraction of velocity fields from a 3D additive manufacturing scaffold using a conventional 2D µPIV system. The µ-computed tomography scaffold geometry was included in a CFD model where perfusion conditions were simulated. Good agreement was found between velocity profiles from measurements and computational results. Maximum velocities were found at the centre of the pore using both techniques with a difference of 12% which was expected according to the accuracy of the µPIV system. However, significant differences in terms of velocity magnitude were found near scaffold substrate due to scaffold brightness which affected the µPIV measurements. As a result, the limitations of the µPIV system only permits a partial validation of the CFD model. Nevertheless, the combination of both techniques allowed a detailed description of velocity maps within a 3D scaffold which is crucial to determine the optimal cell and nutrient transport properties.

A novel flow chamber for biodegradable alloy assessment in physiologically realistic environments.

In order to better understand the in vivo corrosion of biodegradable alloys, it is necessary to replicate the physiological environment as closely as possible. In this study, a novel flow chamber system is developed that allows the investigation of biodegradable alloy corrosion in a simulated physiological environment. The system is designed to reproduce flow conditions encountered in coronary arteries using a parallel plate setup and to allow the culturing of cells. Computational fluid dynamics and analytical methods are used as part of the design process to ensure that suitable flow conditions are maintained in the test region. The system is used to investigate the corrosion behavior of AZ31 alloy foils of different thickness, in test media with and without proteins and in static and dynamic solutions. It is observed that pulsatile flows, similar to those in the coronary arteries, significantly increase corrosion rates and lead to a different corrosion surface morphologies relative to static immersion tests.

Simulation of oxygen transfer in stented arteries and correlation with in-stent restenosis.

Computational models are used to study the combined effect of biomechanical and biochemical factors on coronary in-stent restenosis, which is a postoperative remodeling and regrowth pathology of the stented arteries. More precisely, we address numerical simulations, on the basis of Navier-Stokes and mass transport equations, to study the role of perturbed wall shear stresses and reduced oxygen concentration in a geometrical model reconstructed from a real porcine artery treated with a stent. Joining in vivo and in silico tools of investigation has multiple benefits in this case. On one hand, the geometry of the arterial wall and of the stent closely correspond to a real implanted configuration. On the other hand, the inspection of histological tissue samples informs us on the location and intensity of in-stent restenosis. As a result, we are able to correlate geometrical factors, such as the axial variation of the artery diameter and its curvature; the numerical quantification of biochemical stimuli, such as wall shear stresses; and the availability of oxygen to the inner layers of the artery, with the appearance of in-stent restenosis. This study shows that the perturbation of the vessel curvature could induce hemodynamic conditions that stimulate undesired arterial remodeling.

Virtual surgeries in patients with congenital heart disease: a multi-scale modelling test case.

The objective of this work is to perform a virtual planning of surgical repairs in patients with congenital heart diseases--to test the predictive capability of a closed-loop multi-scale model. As a first step, we reproduced the pre-operative state of a specific patient with a univentricular circulation and a bidirectional cavopulmonary anastomosis (BCPA), starting from the patient's clinical data. Namely, by adopting a closed-loop multi-scale approach, the boundary conditions at the inlet and outlet sections of the three-dimensional model were automatically calculated by a lumped parameter network. Successively, we simulated three alternative surgical designs of the total cavopulmonary connection (TCPC). In particular, a T-junction of the venae cavae to the pulmonary arteries (T-TCPC), a design with an offset between the venae cavae (O-TCPC) and a Y-graft design (Y-TCPC) were compared. A multi-scale closed-loop model consisting of a lumped parameter network representing the whole circulation and a patient-specific three-dimensional finite volume model of the BCPA with detailed pulmonary anatomy was built. The three TCPC alternatives were investigated in terms of energetics and haemodynamics. Effects of exercise were also investigated. Results showed that the pre-operative caval flows should not be used as boundary conditions in post-operative simulations owing to changes in the flow waveforms post-operatively. The multi-scale approach is a possible solution to overcome this incongruence. Power losses of the Y-TCPC were lower than all other TCPC models both at rest and under exercise conditions and it distributed the inferior vena cava flow evenly to both lungs. Further work is needed to correlate results from these simulations with clinical outcomes.

Computational models to predict stenosis growth in carotid arteries: which is the role of boundary conditions?

This work addresses the problem of prescribing proper boundary conditions at the artificial boundaries that separate the vascular district from the remaining part of the circulatory system. A multiscale (MS) approach is used where the Navier-Stokes equations for the district of interest are coupled to a non-linear system of ordinary differential equations which describe the circulatory system. This technique is applied to three 3D models of a carotid bifurcation with increasing stenosis resembling three phases of a plaque growth. The results of the MS simulations are compared to those obtained by two stand-alone models. The MS shows a great flexibility in numerically predicting the haemodynamic changes due to the presence of a stenosis. Nonetheless, the results are not significantly different from a stand-alone approach where flows derived by the MS without stenosis are imposed. This is a consequence of the dominant role played by the outside districts with respect to the stenosis resistance.

Computational evaluation of oxygen and shear stress distributions in 3D perfusion culture systems: macro-scale and micro-structured models.

We present a combined macro-scale/micro-scale computational approach to quantify oxygen transport and flow-mediated shear stress to human chondrocytes cultured in three-dimensional scaffolds in a perfusion bioreactor system. A macro-scale model was developed to assess the influence of the bioreactor design and to identify the proper boundary conditions for the micro-scale model. The micro-scale model based on a micro-computed tomography (microCT) reconstruction of a poly(ethylene glycol terephthalate)/poly(butylene terephthalate) (PEGT/PBT) foam scaffold, was developed to assess the influence of the scaffold micro-architecture on local shear stress and oxygen levels within the scaffold pores. Experiments were performed to derive specific oxygen consumption rates for constructs perfused under flow rates of 0.3 and 0.03 ml min(-1). While macro-scale and micro-scale models predicted similar average oxygen levels at different depths within the scaffold, microCT models revealed small local oxygen variations within the scaffold micro-architecture. The combined macro-scale/micro-scale approach indicated that 0.3 ml min(-1), which subjected 95% of the cells to less than 6.3 mPa shear, would maintain the oxygen supply throughout the scaffold above anoxic levels (>1%), with 99.5% of the scaffold supplied with 8-2% O(2). Alternatively, at 0.03 ml min(-1), the macro-scale model predicted 6% of the cells would be supplied with 0.5-1% O(2), although this region of cells was confined to the periphery of the scaffold. Together with local variations predicted by the micro-scale model, the simulations underline that in the current model system, reducing the flow below 0.03 ml min(-1) would likely have dire consequences on cell viability to pronounced regions within the engineered construct. The presented approach provides a sensitive tool to aid efficient bioreactor optimization and scaffold design.

A new bioreactor for the controlled application of complex mechanical stimuli for cartilage tissue engineering.

Mechanical stimuli have been shown to enhance chondrogenesis on both animal and human chondrocytes cultured in vitro. Different mechanical stimuli act simultaneously in vivo in cartilage tissue and their effects have been extensively studied in vitro, although often in a separated manner. A new bioreactor is described where different mechanical stimuli, i.e. shear stress and hydrostatic pressure, can be combined in different ways to study the mechanobiology of tissue engineered cartilage. Shear stress is imposed on cells by forcing the culture medium through the scaffolds, whereas a high hydrostatic pressure up to 15 MPa is generated by pressurizing the culture medium. Fluid-dynamic experimental tests have been performed and successful validation of the bioreactor has been carried out by dynamic culture of tissue-engineered cartilage constructs. The bioreactor system allows the investigation of the combined effects of different mechanical stimuli on the development of engineered cartilage, as well as other possible three-dimensional tissue-engineered constructs.

Chondrocyte response to high regimens of cyclic hydrostatic pressure in 3-dimensional engineered constructs.

PURPOSE: Despite widespread use of 3-dimensional (3D) micro-porous scaffolds to promote their potential application in cartilage tissue engineering, only a few studies have examined the response to hydrostatic pressure of engineered constructs. A high cyclic pressurization, currently believed to be the predominant mechanical signal perceived by cells in articular cartilage, was used here to stimulate bovine articular chondrocytes cultured in a synthetic 3D porous scaffold (DegraPol). METHODS: Construct cultivation lasted 3 days with applied pressurization cycles of amplitude 10 MPa, frequency 0.33 Hz, and stimulation sessions of 4 hours/day. RESULTS: At 3 days of culture, with respect to pre-culture conditions, the viability of the pressurized constructs did not vary, whereas it underwent a 16% drop in the unpressurized controls. Synthesis of alfa-actin was 34% lower in all cultured constructs. Synthesis of collagen II/collagen I did not vary in pressurized constructs, was 76% lower in unpressurized controls, and was around 230% higher in pressurized constructs with respect to unpressurized controls. Chondrocytes showed a phenotypic spherical morphology at time zero and at 3 days of pressurized culture. CONCLUSIONS: Although the passage from 2D expansion to 3D geometry was effective to guide cell differentiation, only mechanical conditioning enabled the maintenance and further cell differentiation toward a mature chondrocytic phenotype.

Modelling drug elution from stents: effects of reversible binding in the vascular wall and degradable polymeric matrix.

Today the most popular approach for the prevention of the restenosis consists in the use of the drug eluting stents. The stent acts as a source of drug, from a coating or from a reservoir, which is transported into and through the artery wall. In this study, the behaviour of a model of a hydrophilic drug (heparin) released from a coronary stent into the arterial wall is investigated. The presence of the specific binding site action is modelled using a reversible chemical reaction that explains the prolonged presence of drug in the vascular tissue. An axi-symmetric model of a single stent strut is considered. First an advection-diffusion problem is solved using the finite element method. Then a simplified model with diffusion only in the arterial wall is compared with: (i) a model including the presence of reversible binding sites in the vascular wall and (ii) a model featuring a drug reservoir made of a degradable polymeric matrix. The results show that the inclusion of a reversible binding for the drug leads to delayed release curves and that the polymer erosion affects the drug release showing a quicker elution of the drug from the stent.

Toward optimal hemodynamics: computer modeling of the Fontan circuit.

The construction of efficient designs with minimal energy losses is especially important for cavopulmonary connections. The science of computational fluid dynamics has been increasingly used to study the hemodynamic performance of surgical operations. Three-dimensional computer models can be accurately constructed of typical cavopulmonary connections used in clinical practice based on anatomic data derived from magnetic resonance scans, angiocardiograms, and echocardiograms. Using these methods, the hydraulic performance of the hemi-Fontan, bidirectional Glenn, and a variety of types of completion Fontan operations can be evaluated and compared. This methodology has resulted in improved understanding and design of these surgical operations.

Use of rapid prototyping models in the planning of percutaneous pulmonary valved stent implantation.

Percutaneous replacement of the pulmonary valve is a recently developed interventional technique which involves the implantation of a valved stent in the pulmonary trunk. It relies upon careful consideration of patient anatomy for both stent design and detailed procedure planning. Medical imaging data in the form of two-dimensional scans and three-dimensional interactive graphics offer only limited support for these tasks. The paper reports the results of an experimental investigation on the use of arterial models built by rapid prototyping techniques. An analysis of clinical needs has helped to specify proper requirements for such model properties as cost, strength, accuracy, elastic compliance, and optical transparency. Two different process chains, based on the fused deposition modelling technique and on the vacuum casting of thermoset resins in rubber moulds, have been tested for prototype fabrication. The use of anatomical models has allowed the cardiologist's confidence in patient selection, prosthesis fabrication, and final implantation to be significantly improved.

MRI-based multiscale models for the hemodynamic and structural evaluation of surgically reconstructed aortic arches.

The surgical reconstruction of the aortic arch is necessary in pediatric patients suffering from different types of congenital heart malformations, in particular, coarctation of the aorta. Among the reconstruction techniques used in surgical practice end-to-end anastomosis (E/E), Gore-tex graft interposition (GGI) and Gore-tex patch graft aortoplasty (GPGA) are compared in this study with a control model, employing a computational fluid-structure-interaction scheme. This study analyzes the impact of introducing synthetic materials on aortic hemodynamics and wall mechanics. Three-dimensional (3D) geometries of a porcine aortic arch were derived from magnetic resonance imaging (MRI) images. Inlet conditions were derived from MRI velocimetry. A multiscale approach was used for the imposition of outlet conditions, wherein a lumped parameter net provided an active afterload. Evidence was found that ring-like repairs increased blood velocity, whereas GPGA limited it. Vortex presence was greater and longer lasting in GGI. The highest power losses corresponded to GPGA. GGI had an intermediate effect, while E/E dissipated only slightly more than the control case. Wall stresses peak in a longitudinal strip on the subject's left side of the vessel, particularly in the frontal area. There was a concentration of stress at the suture lines. All surgical techniques performed equally well in restoring physiological pressures.

Expansion and drug elution model of a coronary stent.

The present study illustrates a possible methodology to investigate drug elution from an expanded coronary stent. Models based on finite element method have been built including the presence of the atherosclerotic plaque, the artery and the coronary stent. These models take into account the mechanical effects of the stent expansion as well as the effect of drug transport from the expanded stent into the arterial wall. Results allow to quantify the stress field in the vascular wall, the tissue prolapse within the stent struts, as well as the drug concentration at any location and time inside the arterial wall, together with several related quantities as the drug dose and the drug residence times.

Multiscale modelling as a tool to prescribe realistic boundary conditions for the study of surgical procedures.

This work was motivated by the problems of analysing detailed 3D models of vascular districts with complex anatomy. It suggests an approach to prescribing realistic boundary conditions to use in order to obtain information on local as well as global haemodynamics. A method was developed which simultaneously solves Navier-Stokes equations for local information and a non-linear system of ordinary differential equations for global information. This is based on the principle that an anatomically detailed 3D model of a cardiovascular district can be achieved by using the finite element method. In turn the finite element method requires a specific boundary condition set. The approach outlined in this work is to include the system of ordinary differential equations in the boundary condition set. Such a multiscale approach was first applied to two controls: (i) a 3D model of a straight tube in a simple hydraulic network and (ii) a 3D model of a straight coronary vessel in a lumped-parameter model of the cardiovascular system. The results obtained are very close to the solutions available for the pipe geometry. This paper also presents preliminary results from the application of the methodology to a particular haemodynamic problem: namely the fluid dynamics of a systemic-to-pulmonary shunt in paediatric cardiac surgery.

Effects of pulmonary afterload on the hemodynamics after the hemi-Fontan procedure.

A computational fluid dynamics study based on the application of the finite volume method has been performed to investigate the effects of the pulmonary afterload on the hemodynamics after the hemi-Fontan procedure. This operation is generally used as part of a series of staged procedures to treat complex congenital malformations of the heart. It consists of re-directing the superior vena caval flow from the right atrium into the pulmonary arteries, by-passing the right ventricle while excluding the inferior caval flow from the lungs. To reproduce correctly the pulmonary afterload conditions, a simplified lumped-parameter mechanical model of the pulmonary circulation has been developed and linked to the finite volume solver. In addition, the effect of a stenosis in the left pulmonary artery was also examined. In this paper the adopted methodology is presented, together with some of the preliminary results. The model has been used to simulate the local fluid dynamics for different values of the pulmonary arteriolar resistance and lung resistances, allowing a quantitative evaluation of the dissipated energy and the flow distribution into the lungs. The results show that both flow distribution into the lungs and energy dissipation after the hemi-Fontan procedure are only minimally affected by the pulmonary arteriolar resistance.

In vitro steady-flow analysis of systemic-to-pulmonary shunt haemodynamics.

A modified Blalock-Taussig shunt is a connection created between the systemic and pulmonary arterial circulations to improve pulmonary perfusion in children with congenital heart diseases. Survival of these patients is critically dependent on blood flow distribution between the pulmonary and systemic circulations which in turn depends upon the flow resistance of the shunt. Previously, we investigated the pressure-flow relationship in rigid shunts with a computational approach. to estimate the pulmonary blood flow rate on the basis of the in vivo measured pressure drop. The present study aims at evaluating, in vitro how the anastomotic distensibility and restrictions due to suture presence affect the shunt pressure-flow relationship. Two actual Gore-Tex shunts (3 and 4 mm diameters) were sutured to compliant conduits by a surgeon and tested at different steady flow rates (0.25-11 min(-1)) and pulmonary pressures (3-34 mmHg). Corresponding computational models were also created to investigate the role of the anastomotic restrictions due to sutures. In vitro experiments showed that pulmonary artery pressure affects the pressure-flow relationship of the anastomoses. particularly at the distal site. However, this occurrence scarcely influences the total shunt pressure drop. Comparisons between in vitro and computational models without anastomotic restrictions show that the latter underestimates the in vitro pressure drops at any flow rate. The addition of the anastomotic restrictions (31 and 47% of the original area of 3 and 4 mm shunts, respectively) to the computational models reduces the gap, especially at high shunt flow rate and high pulmonary pressure.