Hemodynamic environments from opposing sides of human aortic valve leaflets evoke distinct endothelial phenotypes in vitro.

The regulation of valvular endothelial phenotypes by the hemodynamic environments of the human aortic valve is poorly understood. The nodular lesions of calcific aortic stenosis (CAS) develop predominantly beneath the aortic surface of the valve leaflets in the valvular fibrosa layer. However, the mechanisms of this regional localization remain poorly characterized. In this study, we combine numerical simulation with in vitro experimentation to investigate the hypothesis that the previously documented differences between valve endothelial phenotypes are linked to distinct hemodynamic environments characteristic of these individual anatomical locations. A finite-element model of the aortic valve was created, describing the dynamic motion of the valve cusps and blood in the valve throughout the cardiac cycle. A fluid mesh with high resolution on the fluid boundary was used to allow accurate computation of the wall shear stresses. This model was used to compute two distinct shear stress waveforms, one for the ventricular surface and one for the aortic surface. These waveforms were then applied experimentally to cultured human endothelial cells and the expression of several pathophysiological relevant genes was assessed. Compared to endothelial cells subjected to shear stress waveforms representative of the aortic face, the endothelial cells subjected to the ventricular waveform showed significantly increased expression of the "atheroprotective" transcription factor Kruppel-like factor 2 (KLF2) and the matricellular protein Nephroblastoma overexpressed (NOV), and suppressed expression of chemokine Monocyte-chemotactic protein-1 (MCP-1). Our observations suggest that the difference in shear stress waveforms between the two sides of the aortic valve leaflet may contribute to the documented differential side-specific gene expression, and may be relevant for the development and progression of CAS and the potential role of endothelial mechanotransduction in this disease.

Pulmonary tissue engineering using dual-compartment polymer scaffolds with integrated vascular tree.

OBJECTIVES: The persistent shortage of donor organs for lung transplantation illustrates the need for new strategies in organ replacement therapy. Pulmonary tissue engineering aims at developing viable hybrid tissue for patients with chronic respiratory failure. METHODS: Dual-chamber polymer constructs that mimic the characteristics of the pulmonary air-blood interface were fabricated by microfabrication techniques using the biocompatible polymer polydimethylsiloxane. One compartment ("vascular chamber") was designed as a capillary network to mimic the pulmonary microvasculature. The other compartment ("parenchymal chamber") was designed to permit gas exchange. Immortalized mouse lung epithelium cells (MLE-12) were cultured on the surface of polystyrene microcarrier beads. These beads were subsequently injected into the parenchymal chamber of the dual-chamber microsystems. The vascular compartment was perfused with cell culture medium in a bioreactor and the construct was maintained in culture for 1 week. RESULTS: The microcarriers evenly distributed MLE-12 cells on the parenchymal compartment surface. Confluent cell layers were confirmed by fluorescent and electron microscopy. Adequate proliferation of MLE-12 cells within the construct was monitored via the DNA content. Viability of the cells was maintained over 1 week. Finally, cellular specificity and functional capacity in situ were demonstrated by immunostaining for proSP-B and proSP-C (alveolar epithelium), and by using MLE-12 cells transfected to overexpress green fluorescent protein. CONCLUSION: We conclude that functional hybrid microsystems mimicking the basic building plan of alveolar tissue can be engineered in vitro.

A multiscale computational comparison of the bicuspid and tricuspid aortic valves in relation to calcific aortic stenosis.

Patients with bicuspid aortic valve (BAV) are more likely to develop a calcific aortic stenosis (CAS), as well as a number of other ailments, as compared to their cohorts with normal tricuspid aortic valves (TAV). It is currently unknown whether the increase in risk of CAS is caused by the geometric differences between the tricuspid and bicuspid valves or whether the increase in risk is caused by the same underlying factors that produce the geometric difference. CAS progression is understood to be a multiscale process, mediated at the cell level. In this study, we employ multiscale finite-element simulations of the valves. We isolate the effect of one geometric factor, the number of cusps, in order to explore its effect on multiscale valve mechanics, particularly in relation to CAS. The BAV and TAV are modeled by a set of simulations describing the cell, tissue, and organ length scales. These simulations are linked across the length scales to create a coherent multiscale model. At each scale, the models are three-dimensional, dynamic, and incorporate accurate nonlinear constitutive models of the valve leaflet tissue. We compare results between the TAV and BAV at each length scale. At the cell-scale, our region of interest is the location where calcification develops, near the aortic-facing surface of the leaflet. Our simulations show the observed differences between the tricuspid and bicuspid valves at the organ scale: the bicuspid valve shows greater flexure in the solid phase and stronger jet formation in the fluid phase relative to the tricuspid. At the cell-scale, however, we show that the region of interest is shielded against strain by the wrinkling of the fibrosa. Thus, the cellular deformations are not significantly different between the TAV and BAV in the calcification-prone region. This result supports the assertion that the difference in calcification observed in the BAV versus TAV may be due primarily to factors other than the simple geometric difference between the two valves.

Microfabrication of three-dimensional engineered scaffolds.

One of the principal challenges facing the field of tissue engineering over the past 2 decades has been the requirement for large-scale engineered constructs comprising precisely organized cellular microenvironments. For vital organ assist and replacement devices, microfluidic-based systems such as the microcirculation, biliary, or renal filtration and resorption systems and other functional elements containing multiple cell types must be generated to provide for viable engineered tissues and clinical benefit. Over the last several years, microfabrication technology has emerged as a versatile and powerful approach for generating precisely engineered scaffolds for engineered tissues. Fabrication process tools such as photolithography, etching, molding, and lamination have been established for applications involving a range of biocompatible and biodegradable polymeric scaffolding materials. Computational fluid dynamic designs have been used to generate scaffold designs suitable for microvasculature and a number of organ-specific constructs; these designs have been translated into 3-dimensional scaffolding using microfabrication processes. Here a brief overview of the fundamental microfabrication technologies used for tissue engineering will be presented, along with a summary of progress in a number of applications, including the liver and kidney.

A finite shell element for heart mitral valve leaflet mechanics, with large deformations and 3D constitutive material model.

This paper presents a shell finite element formulation appropriate for simulating the heart valve leaflet mechanics, including three-dimensional (3D) stress and strain effects. A 4-node mixed-interpolation shell is formulated in convected coordinates. This shell model is made capable of handling arbitrary 3D material models by use of an algorithm that satisfies the shell stress assumption at every element integration point. A method for tracking the fiber direction is incorporated. The resulting shell element operates under the same conditions as a standard 4-node shell element with 5 degrees of freedom per node, but extends the modeling capabilities to handle large-deformation and anisotropic behavior.

A combined FEM/genetic algorithm for vascular soft tissue elasticity estimation.

Tissue elasticity reconstruction is a parameter estimation effort combining imaging, elastography, and computational modeling to build maps of soft tissue mechanical properties. One application is in the characterization of atherosclerotic plaques in diseased arteries, wherein the distribution of elastic properties is required for stress analysis and plaque stability assessment. In this paper, a computational scheme is proposed for elasticity reconstruction in soft tissues, combining finite element modeling (FEM) for mechanical analysis of soft tissues and a genetic algorithm (GA) for parameter estimation. With a model reduction of the discrete elasticity values into lumped material regions, namely the plaque constituents, a robust, adaptive strategy can be used to solve inverse elasticity problems involving complex and inhomogeneous solution spaces. An advantage of utilizing a GA is its insistence on global convergence. The algorithm is easily implemented and adaptable to more complex material models and geometries. It is meant to provide either accurate initial guesses of low-resolution elasticity values in a multi-resolution scheme or as a replacement for failing traditional elasticity estimation efforts.

A coarse-grained model for force-induced protein deformation and kinetics.

Force-induced changes in protein conformation are thought to be responsible for certain cellular responses to mechanical force. Changes in conformation subsequently initiate a biochemical response by alterations in, for example, binding affinity to another protein or enzymatic activity. Here, a model of protein extension under external forcing is created inspired by Kramers' theory for reaction rate kinetics in liquids. The protein is assumed to have two distinct conformational states: a relaxed state, C(1), preferred in the absence of external force, and an extended state, C(2), favored under force application. In the context of mechanotransduction, the extended state is a conformation from which the protein can initiate signaling. Appearance and persistence of C(2) are assumed to lead to transduction of the mechanical signal into a chemical one. The protein energy landscape is represented by two harmonic wells of stiffness kappa(1) and kappa(2), whose minima correspond to conformations C(1) and C(2). First passage time t(f) from C(1) to C(2) is determined from the Fokker-Plank equation employing several different approaches found in the literature. These various approaches exhibit significant differences in behavior as force increases. Although the level of applied force and the energy difference between states largely determine equilibrium, the dominant influence on t(f) is the height of the transition state. Distortions in the energy landscape due to force can also have a significant influence, however, exhibiting a weaker force dependence than exponential as previously reported, approaching a nearly constant value at a level of force that depends on the ratio kappa(1)/kappa(2). Two model systems are used to demonstrate the utility of this approach: a short alpha-helix undergoing a transition between two well-defined states and a simple molecular motor.

A large-strain finite element formulation for biological tissues with application to mitral valve leaflet tissue mechanics.

This paper presents a finite element formulation suitable for large-strain modeling of biological tissues and uses this formulation to implement an accurate finite element model for mitral valve leaflet tissue. First, an experimentally derived strain energy function is obtained from literature. This function is implemented in finite elements using the mixed pressure-displacement formulation. A modification is made to aid in maintaining positive definiteness of the stiffness matrix at low strains. The numerical implementation is shown to be accurate in representing the analytical model of material behavior. The mixed formulation is useful for modeling of soft biological tissues in general, and the model presented here is applicable to finite element simulation of mitral valve mechanics.

Endothelialized microvasculature based on a biodegradable elastomer.

Vital organs maintain dense microvasculature to sustain the proper function of their cells. For tissue- engineered organs to function properly, artificial capillary networks must be developed. We have microfabricated capillary networks with a biodegradable and biocompatible elastomer, poly(glycerol sebacate) (PGS). We etched capillary patterns onto silicon wafers by standard micro-electromechanical systems (MEMS) techniques. The resultant silicon wafers served as micromolds for the devices. We bond the patterned PGS film with a flat film to create capillary networks that were perfused with a syringe pump at a physiological flow rate. The devices were endothelialized under flow conditions, and part of the lumens reached confluence within 14 days of culture. This approach may lead to tissue-engineered microvasculature that is critical in vital organs engineering.

Distinct endothelial phenotypes evoked by arterial waveforms derived from atherosclerosis-susceptible and -resistant regions of human vasculature.

Atherosclerotic lesion localization to regions of disturbed flow within certain arterial geometries, in humans and experimental animals, suggests an important role for local hemodynamic forces in atherogenesis. To explore how endothelial cells (EC) acquire functional/dysfunctional phenotypes in response to vascular region-specific flow patterns, we have used an in vitro dynamic flow system to accurately reproduce arterial shear stress waveforms on cultured human EC and have examined the effects on EC gene expression by using a high-throughput transcriptional profiling approach. The flow patterns in the carotid artery bifurcations of several normal human subjects were characterized by using 3D flow analysis based on actual vascular geometries and blood flow profiles. Two prototypic arterial waveforms, "athero-prone" and "athero-protective," were defined as representative of the wall shear stresses in two distinct regions of the carotid artery (carotid sinus and distal internal carotid artery) that are typically "susceptible" or "resistant," respectively, to atherosclerotic lesion development. These two waveforms were applied to cultured EC, and cDNA microarrays were used to analyze the differential patterns of EC gene expression. In addition, the differential effects of athero-prone vs. athero-protective waveforms were further characterized on several parameters of EC structure and function, including actin cytoskeletal organization, expression and localization of junctional proteins, activation of the NF-kappaB transcriptional pathway, and expression of proinflammatory cytokines and adhesion molecules. These global gene expression patterns and functional data reveal a distinct phenotypic modulation in response to the wall shear stresses present in atherosclerosis-susceptible vs. atherosclerosis-resistant human arterial geometries.

Mechanical analysis of atherosclerotic plaques based on optical coherence tomography.

Finite element analysis is a powerful tool for investigating the biomechanics of atherosclerosis and has thereby provided an improved understanding of acute myocardial infarction. Structural analysis of arterial walls is traditionally performed using geometry contours derived from histology. In this paper we demonstrate the first use of a new imaging technique, optical coherence tomography (OCT), as a basis for finite element analysis. There are two primary benefits of OCT relative to histology: 1) imaging is performed without excessive tissue handling, providing a more realistic geometry than histology and avoiding structural artifacts common to histologic processing, and 2) OCT imaging can be performed in vivo, making it possible to study disease progression and the effect of therapeutic treatments in animal models and living patients. Patterns of mechanical stress and strain distributions computed from finite element analysis based on OCT were compared with those from modeling based on "gold standard" histology. Our results indicate that vascular structure and composition determined by OCT provides an adequate basis for investigating the biomechanical factors relevant to atherosclerosis and acute myocardial infarction.

How flexible is alpha-actinin's rod domain?

Alpha-actinin, an actin binding protein, plays a key role in cell migration, cross-links actin filaments in the Z-disk, and is a major component of contractile muscle apparatus. The flexibility of the molecule is critical to its function. The flexibility of various regions of the molecule, including the linker connecting central subunits is studied using constant force steered molecular dynamics simulations. The linker, whose structure has been a subject of debate, is predicted to be semi-flexible. The flexibility of the linker is compared to all possible segments of equal length throughout the molecule. The stretching profile of the molecule at different forces suggests that loops and regions adjacent to the loops are much more rigid than the helices in the protein. Amino acid composition analysis of most flexible and most rigid regions of the molecule reveals that the rigid regions are rich in Ser, Val and Ile whereas the flexible regions are rich in Ala, Leu and Glu.

On the molecular basis for mechanotransduction.

Much is currently known about the signaling pathways that are excited when cells are subjected to a mechanical stimulus, yet we understand little of the process by which the mechanical perturbation is transformed into a biochemical signal. Numerous theories have been proposed, and each has merit. While cells may possess many different ways of responding to stress, the existence of a single unifying principle has much appeal. Here we propose the hypothesis that cells sense mechanical force through changes in protein conformation, leading to altered binding affinities of proteins, ultimately initiating an intracellular signaling cascade or producing changes in the proteins localized to regions of high stress. More generally, this represents an alternative to transmembrane signaling through receptor-ligand interactions providing the cell with a means of reacting to changes in its mechanical, as opposed to biochemical, environment. One example is presented showing how the binding affinity between the focal adhesion targeting domain of focal adhesion kinase and the LD motif of paxillin is influenced by externally applied force.

A three-dimensional viscoelastic model for cell deformation with experimental verification.

A three-dimensional viscoelastic finite element model is developed for cell micromanipulation by magnetocytometry. The model provides a robust tool for analysis of detailed strain/stress fields induced in the cell monolayer produced by forcing one microbead attached atop a single cell or cell monolayer on a basal substrate. Both the membrane/cortex and the cytoskeleton are modeled as Maxwell viscoelastic materials, but the structural effect of the membrane/cortex was found to be negligible on the timescales corresponding to magnetocytometry. Numerical predictions are validated against experiments performed on NIH 3T3 fibroblasts and previous experimental work. The system proved to be linear with respect to cytoskeleton mechanical properties and bead forcing. Stress and strain patterns were highly localized, suggesting that the effects of magnetocytometry are confined to a region extending <10 microm from the bead. Modulation of cell height has little effect on the results, provided the monolayer is >5 micro m thick. NIH 3T3 fibroblasts exhibited a viscoelastic timescale of approximately 1 s and a shear modulus of approximately 1000 Pa.

Dynamic rotational seeding and cell culture system for vascular tube formation.

Optimization of cell seeding and culturing is an important step for the successful tissue engineering of vascular conduits. We evaluated the effectiveness of using a hybridization oven for rotational seeding and culturing of ovine vascular myofibroblasts onto biodegradable polymer scaffolds suitable for replacement of small- and large-diameter blood vessels. Large tubes (12 mm internal diameter and 60 mm length, n = 4) and small tubes (5 mm internal diameter and 20 mm length, n = 4) were made from a combination of polyglycolic acid/poly-4-hydroxybutyrate and coated with collagen solution. Tubes were then placed in culture vessels containing a vascular myofibroblast suspension (10(6) cells/cm(2)) and rotated at 5 rpm in a hybridization oven at 37 degrees C. Light and scanning electron microscopy analyses were performed after 5, 7, and 10 days. Myofibroblasts had formed confluent layers over the outer and inner surfaces of both large and small tubular scaffolds by day 5. Cells had aligned in the direction of flow by day 7. Multiple spindle-shaped cells were observed infiltrating the polymer mesh. Cell density increased between day 5 and day 10. All conduits maintained their tubular shape throughout the experiment. We conclude that dynamic rotational seeding and culturing in a hybridization oven is an easy, effective, and reliable method to deliver and culture vascular myofibroblasts onto tubular polymer scaffolds.

A microfabricated array bioreactor for perfused 3D liver culture.

We describe the design, fabrication, and performance of a bioreactor that enables both morphogenesis of 3D tissue structures under continuous perfusion and repeated in situ observation by light microscopy. Three-dimensional scaffolds were created by deep reactive ion etching of silicon wafers to create an array of channels (through-holes) with cell-adhesive walls. Scaffolds were combined with a cell-retaining filter and support in a reactor housing designed to deliver a continuous perfusate across the top of the array and through the 3D tissue mass in each channel. Reactor dimensions were constructed so that perfusate flow rates meet estimated values of cellular oxygen demands while providing fluid shear stress at or below a physiological range (<2 dyne cm(2)), as determined by comparison of numerical models of reactor fluid flow patterns to literature values of physiological shear stresses. We studied the behavior of primary rat hepatocytes seeded into the reactors and cultured for up to 2 weeks, and found that cells seeded into the channels rearranged extensively to form tissue like structures and remained viable throughout the culture period. We further observed that preaggregation of the cells into spheroidal structures prior to seeding improved the morphogenesis of tissue structure and maintenance of viability. We also demonstrate repeated in situ imaging of tissue structure and function using two-photon microscopy.