Layered Elastomeric Fibrous Scaffolds: An In-Silico Study of the

Aug 10, 2017 - Fiber diameter and alignment are able to modulate cellular morphology, that is, projected cell area, aspect ratio, and long axis length...
0 downloads 0 Views 6MB Size
Subscriber access provided by UNIVERSITY OF ADELAIDE LIBRARIES

Article

Layered elastomeric fibrous scaffolds: An in-silico study of the achievable range of mechanical behaviors James Carleton, G. J. Rodin, and Michael S Sacks ACS Biomater. Sci. Eng., Just Accepted Manuscript • DOI: 10.1021/acsbiomaterials.7b00308 • Publication Date (Web): 10 Aug 2017 Downloaded from http://pubs.acs.org on August 16, 2017

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

ACS Biomaterials Science & Engineering is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 52

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Biomaterials Science & Engineering

Layered Elastomeric Fibrous Scaffolds: An In-Silico Study of the Achievable Range of Mechanical Behaviors James B. Carleton,† Gregory J. Rodin,†,‡ and Michael S. Sacks∗,†,‡ Center for Cardiovascular Simulation, Institute for Computational Engineering and Sciences and the Department of Biomedical Engineering, The University of Texas at Austin, 201 East 24th Street, Austin TX 78712, and Department of Aerospace Engineering and Engineering Mechanics, 210 East 24th Street, Austin, TX 78712 E-mail: [email protected]



To whom correspondence should be addressed Center for Cardiovascular Simulation, Institute for Computational Engineering and Sciences and the Department of Biomedical Engineering, The University of Texas at Austin, 201 East 24th Street, Austin TX 78712 ‡ Department of Aerospace Engineering and Engineering Mechanics, 210 East 24th Street, Austin, TX 78712 †

1

ACS Paragon Plus Environment

ACS Biomaterials Science & Engineering

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ABSTRACT Our goal herein is to understand the mechanisms underlying soft tissue and scaffold behaviors by developing a physically-based micromechanical model as a means to connect the macroscopic behaviors to the underlying microstructural phenomena. Due to its well-documented capacity for generating elastomeric fibrous materials with a wide range of realizable architectures, the electrospun scaffold was used as the exemplar biomaterial. Fibrous network geometries based on a random walk algorithm were first generated to form the basis for subsequent micromechanical simulations. A basic understanding of randomly oriented fibrous network phenomena was then developed, and subsequently expanded on using networks with aligned fibers. Simulation results were then compared with experimental observations of electrospun scaffolds to evaluate the validity of the simulations. The effects of fiber alignment, tortuosity, and material properties on macroscopic mechanical behavior of the material have been presented both individually and in combination. We have seen that all three aspects of the scaffold network can have significant effects on the macroscopic behavior for different load cases. Overall, accurate representation of detailed fibrous network geometry permitted a greater understanding of the complex mechanisms underlying the macroscopic behavior unique to these biomaterials. Insights gained from such simulations can significantly aid the process of designing scaffold network geometries that result in engineered tissues that function as well as or better than the native tissues they are intended to replace.

Keywords scaffolds, micromechanics, simulations, tissue engineering

2

ACS Paragon Plus Environment

Page 2 of 52

Page 3 of 52

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Biomaterials Science & Engineering

INTRODUCTION The term "tissue engineering" was coined in 1988 by the great biomechanical engineering scientist Y.C. Fung. 1 In that early work, Fung underscored the importance of ...the application of principles and methods of engineering and life sciences toward a fundamental understanding of structure-function relationships in normal and pathologic mammalian tissues and the development of biological substitutes to restore, maintain, or improve tissue function. It is thus imperative that fundamental structure-function understanding guides the reproduction of native tissue if it is to emulate its native counterpart successfully. 2,3 Accurate simulations of the mechanical behavior of fibrous networks can be used to inform the design and manufacture of engineered tissues, especially with future fabrication technologies that enable control over important microstructural features. Moreover, while the current approach to manufacturing high quality scaffolds and functional engineered tissues remains somewhat of an art, current and future scaffold production would benefit enormously from the specification of target network geometries that produce desirable mechanical behavior. While the mechanical behavior of both hard and soft tissues can be quite complex, modeling of soft tissues is generally regarded as a more challenging task, due to major effects associated with large strains and significant microstructural changes. Traditionally, the mechanical behavior of soft tissues has been modeled using pseudo-hyperelastic phenomenological models. 4–8 But, like any natural or engineered biomaterial, the complex macroscopic mechanical behavior of soft tissues is a result of multi-scale interactions of the constituent phases. Those include collagen, elastin, nervous and muscular fibers, and, on finer scales, cells and extracellular matrix components such as glycosaminoglycans and proteoglycans. Thus relating the macroscopic mechanical behavior of soft tissues to their microstructure is an important step toward the development of a comprehensive multi-scale modeling framework.

3

ACS Paragon Plus Environment

ACS Biomaterials Science & Engineering

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The fibrous microstructure of many natural soft tissues is often layered, especially for membrane structures (Figure 1-a). Native tissue structures can also include continuous graduations in fiber type, density, and organization within layers. In most native soft tissue models networks primarily act as very long and slender, mechanically non-interacting fibers. For example, we have recently developed a functional gradient long-fiber model (FGM) 9 that incorporated the transmural composition and fiber structure of the aortic valve (AV) leaflet. The FGM model predicted large stress variations both between and within the leaflet layers, suggesting that the continually varying structure of the AV leaflet has an important purpose with regard to valve function. We have also shown that explicit fiber-fiber interactions can be accounted for in related long-fiber models. 10,11 In attempting to mimic such structures, microstructures of engineered soft tissue scaffolds are often layered. As in native tissues, heterogeneity and anisotropy may significantly affect the macroscopic and microscopic mechanical behavior and functionality of synthetic biomaterials. However, in contrast to native tissues, the deformation of each fiber at the local scale is influenced by the positions and deformations of neighboring fibers through interconnections in synthetic biomaterials. 12 Thus understanding of relations between biological functionality and microstructure cannot be attained without characterizing relations between the microstructure and mechanical properties. Previously reported simulations have included three-dimensional assemblies of straight, unconnected fibers as the basis of a multiscale modeling approach in which the microscale geometry is regarded as a unit cell. 13 Volume averaging is used to relate the microscopic forces and displacements to the macroscopic continuum constitutive equation. The larger-scale structure is discretized with finite elements and the constitutive equations are evaluated at each integration point. The predicted macroscopic behavior agreed reasonably well with experimental measurements, but was unable to capture the effects of fiber straightening or interactions due to fiber-to-fiber intersections. Clearly, a more detailed understanding of these multiscale phenomena can only be gained using physically-realistic representations of the actual microstructures.

4

ACS Paragon Plus Environment

Page 4 of 52

Page 5 of 52

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Biomaterials Science & Engineering

Simulations based on directly capturing microstructural details have elucidated the relationship between scaffold fiber deformations and the aspect ratios of cells filling the pore spaces. 14,15 Those results also reveal that the highly loaded fibers in the network form long, straight effective fibers that are reminiscent of the structure of soft tissues. While these simulations provide an accurate geometric representation of actual scaffolds, the generated geometries are limited to specific scaffolds that have been imaged. Thus extending this approach to statistical studies and other microstructures, especially in 3-D, is somewhat problematic. Clearly, better geometric representations are needed for more accurate macroscopic predictions and for evaluation of quantities of interest related to the geometry of the evolving microstructure. Moreover, the ability to parametrically investigate the ranges of possible microstructures that elucidate what can be potentially achievable mechanically at the macro level remains an unmet goal in biomaterials design. It is clear that existing simulation approaches to fibrous networks lack robust geometric representations necessary for more reliable predictions of the macroscopic response and evolving microstructure. Such predictions would enable parametric studies of fibrous microstructures with the objective of identifying optimal microstructures that realize desired macroscopic properties. This mostly unmet objective is central to biomaterials design. Our goal herein is thus to understand the mechanisms underlying scaffold behaviors. For this purpose, physically-based structural models are preferable to phenomenological models since they offer a means to connect the macroscopic behaviors to the underlying microstructural phenomena. Structural models have indicated that the highly nonlinear and anisotropic macroscopic behavior of soft tissues can be explained by alignment and straightening of undulated fibers at the microscale, in addition to material nonlinear effects. 7,16–20 These predictions are based on the assumption of affine (homogeneous) fiber deformations, as supported by recent work on planar collagenous tissues. 21 Due to its widely-documented capacity of generating elastomeric, fibrous materials with a wide range of realizable, controlled architectures, 22–28 the electrospun scaffold was used as the exemplar biomaterial. Fibrous network

5

ACS Paragon Plus Environment

ACS Biomaterials Science & Engineering

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 ACS Paragon Plus Environment

Page 6 of 52

Page 7 of 52

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Biomaterials Science & Engineering

METHODS Approach. In the present study, we first describe the geometric and mechanical models for generalized, planar fibrous networks. These models were based on a random walk algorithm (Figure 1-b) 29 and a mechanical model, which is a variant of Cosserat’s rod theory.Âă We then present simulation results aimed at quantifying the macroscopic response and evolving fibrous microstructure for both initially isotropic and anisotropic networks.Âă Simulation results were then compared with experimental observations of electrospun scaffolds to evaluate the validity of the simulations and are used to explore the effects of different microstructural features on macroscopic behavior. Network simulations were seen to be a useful tool for supplementing experiments when attempting to understand phenomena that are difficult to measure. Finally, we discuss various implications of our simulation results on characterization and design of biomaterials.

Source Material and Geometric Model. In this section, we describe a simulation methodology for quantifying the mechanical response of fibrous networks intended for modeling electrospun polymeric scaffolds, 30,31 whose structural and mechanical properties have been thoroughly characterized previously. 26,28,30,32 In previous studies, cytocompatible poly (ester urethane) urea (PEUU) was synthesized from polycaprolactone diol and 1,4-diisocyanatobutane with subsequent chain extension by putrescine. By syringe pump into a stainless-steel capillary suspended 13-cm vertically over a 4.5 cm diameter aluminum mandrel 5-wt% PEUU solution in hexaïňĆuoroisopropanol was fed at 1.0 mL/h. PEUU was charged with +12 kV and the aluminum target with -7 kV using high voltage generators. Aligned PEUU ïňĄbers were formed by electrospinning onto the target rotating at speeds ranging from 0.0 to 13.8 m/s. Structural and biaxial mechanical characterization were then performed. 26,28,30,32 7

ACS Paragon Plus Environment

ACS Biomaterials Science & Engineering

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 8 of 52

Constructing an accurate representation of the resulting scaffold network geometry is essential to performing high-fidelity simulations that provide insights into the microstructural mechanisms underlying complex macroscopic mechanical behavior. Electrospun scaffolds have been observed to be layered structures. 30,31 Therefore, they are modeled as 3-D fibrous networks formed by stacked layers. The fibrous geometry of each layer is idealized as a collection of random walks of the total length Lt confined to a square LB × LB box. 29 These random walks are characterized by the line density ξ :=

Lt A

and the probability density

function f (ϕ) for the fiber orientation angle ϕ, or simply the orientation density function (ODF). For the purposes of mechanical simulations, it is imperative to rely on periodic random walks. That is, if a walk intersects an edge of the confinement box, it is forced to reenter the box on the opposite side; for details we refer to our previous work. 29 This previously reported algorithm 29 allows one to match a target fiber density ξ t and a target ODF f t (ϕ), including those determined experimentally for the scaffolds. Since the algorithm mimics the actual electro-spinning process, it is capable of capturing many geometric and topological features of the scaffold (Figure 1-b). In addition to the prescribed targets, ξ t and f t (ϕ), the algorithm requires the following input parameters 1. nf - number of random walks per layer. 2. ns - number of segments per random walk. 3. sˆ - segment length normalized by ξ −1 . ˆ B - box edge length normalized by ξ −1 . 4. L 5. α - the maximum orientation change between two sequential segments. These parameters are constrained by nf ns sˆ = 1, ˆ 2B L

8

ACS Paragon Plus Environment

(1)

Page 9 of 52

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Biomaterials Science & Engineering

so that only three of these parameters must be specified. Small values of α result in nearly straight fibers and large α values result in highly tortuous fibers. This geometric feature may be quantified by defining the tortuosity of a fiber chain, which consists of all random walk segments between two adjacent intersections, as Lc τ := ¯ , Lc

(2)

¯ c is the length of the chain’s end-to-end vector. where Lc is the arc length of the chain, and L Although the average tortuosity of the network also depends on ξ, the parameter α allows for tortuosity control. Random walks become fibers once they are assigned a diameter D. In this assignment, we neglect out-of-plane O(D) deflections of the fibers occurring at fiber intersections. Thus each layer is assumed to be planar and of thickness D. With this provision, the 3-D network is naturally formed by stacking the planar layers. The stacking results in forming additional, inter-layer, fiber intersections. The fiber volume fraction φ for the individual layers and the entire 3-D network is φ :=

π ξD. 4

(3)

Mechanical Model and Numerical Implemention. Following the standard practice in continuum mechanics, we model fibers as rods. That is, we introduce simplifications allowing us to describe the fiber deformation as a system of equations defined on the fiber centerline. There is a long history of rod theories going back to Bernoulli and Euler. Modern continuum rod theories, accounting for large strains and rotations, and nonlinear constitutive equations, are classified as either asymptotic or direct approaches. 33 In asymptotic approaches, the rod equilibrium and constitutive equations are derived from the underlying 3-D continuum theory, by exploiting the rod slenderness. In direct approaches, the governing equations are formulated without any reference to 3-D

9

ACS Paragon Plus Environment

ACS Biomaterials Science & Engineering

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

theories. The result of both approaches is a set of governing equations expressed in terms of the position of a point on the centerline and the shape and orientation of the cross section. Theories using a set of vectors to describe the orientation of the cross section are classified as Cosserat rod theories, and the vectors are referred to as the directors. Rod theories of interest to this work result in nonlinear governing equations. Consequently numerical implementations of rod theories are nontrivial, and usually involve tradeoffs among accuracy, efficiency, and robustness. For example, the ABAQUS Euler-Bernoulli beam element uses an asymptotic approach, which includes large 3-D rotations and strains. 34 While in theory this approach is ideal for our simulations, ABAQUS is not an ideal tool for our purposes. A serial implementation may be adequate for purely 2-D fiber networks 35 or 3-D networks with straight fiber geometry between intersections, 36 but this limitation is exceedingly restrictive for large-scale simulations. For single layer simulations considered here, problems with hundreds of thousands of degrees of freedom are common, partly because multiple fiber segments are required to capture the curved fiber geometry. 29 We also needed to perform multilayer simulations 37 to assess the sufficiency of single layer simulations. Multilayer simulations of complete heart valve leaflets could easily require hundreds of millions of degrees of freedom. For a problem with 2.5 million degrees of freedom, the peak parallel speedup (serial time over parallel time) for our code was 6.35 on 64 processors (Figure 2).

Figure 2: Parallel speedup plots for problems with 300,000, 600,000, and 2,500,000 degrees of freedom (dof).

10

ACS Paragon Plus Environment

Page 10 of 52

Page 11 of 52

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Biomaterials Science & Engineering

Our assumptions and approach closely follow those adopted for the ABAQUS EulerBernoulli beam element, but there are some differences, which we highlight here. A complete description of our approach is provided in our previous work. 37 In our approach, the cross section is defined by a position vector and two orthonormal director vectors. We assume that the rod deforms such that plane sections remain plane, no deformations occur within the cross section, and the cross sections remain normal to the centerline. In contrast to the ABAQUS approach, in which the primary kinematical variables are the displacement and rotation vectors at each node, we relied on the position and director vectors the primary kinematical variables. Our approach is computationally less efficient, as it involves more equations and unknowns, but the governing equations are much simpler. We also differ from ABAQUS in the way we enforce the kinematic constraint that the cross section remains normal to the centerline. In the ABAQUS Euler-Bernoulli beam element, this constraint is enforced by direct substitution into the virtual work statement. As a result of this substitution, second derivatives of the position vector with respect to arc length appear in the governing equations. In our approach, we delay applying this constraint and apply it only after performing the finite element discretization, as has been done by others. 38 One advantage of this approach is that only first derivatives of the position vector appear in the governing equations. Another consequence of this choice is that extra terms due to shear deformation appear in our version, which do not appear in the virtual work statement of the ABAQUS Euler-Bernoulli beam element. Terms due to shear deformation do appear in the virtual work statement of the ABAQUS Timoshenko beam element, but we seek to model an Euler-Bernoulli beam, and the applied constraints remove the virtual strain energy due to these terms. One other minor difference is that ABAQUS uses the Green Lagrange strain in the virtual work term due to axial deformation, whereas we choose to use the stretch instead. As in our previous study, 14 an incompressible Yeoh material 39 is used for modeling the axial response, whereas a linear elastic material model used for fiber bending, torsion, and shear. There are many ways of modelling fiber intersections. We simply assumed that at all

11

ACS Paragon Plus Environment

ACS Biomaterials Science & Engineering

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 ACS Paragon Plus Environment

Page 12 of 52

Page 13 of 52

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Biomaterials Science & Engineering

details are provided in our previous work. 37 The finite element method with linear shape functions and a reduced quadrature scheme was used to reformulate the problem as a system of nonlinear algebraic equations, which were solved numerically using Newton’s method. The fiber intersection constraints and periodic boundary conditions were enforced using translational and rotational penalties. The parallel finite element code was developed using PETCSc. 40–42

Model Calibration. In order to make useful observations about material mechanical behavior, one must conduct simulations on sufficiently large networks, which accurately represent the material response. Such networks are referred to as Representative Volume Elements (RVEs). For 3-D networks, ˆ B and (ii) the two sizes are important: (i) the minimum normalized in-plane box size L minimum number of layers. The process of identifying the RVEs is purely empirical. The box size is increased (Figure 3-b) until the error in a quantity of interest falls below an acceptable magnitude. Once the quantity of interest and the corresponding error tolerance have been specified, parametric studies are conducted to determine the minimum sizes. Here we summarize key findings of those studies; for details we refer to our previous work. 37 • The quantity of interest chosen was the strain energy density induced by the macroscopic deformation gradient  

F= 

1.50

0

0

1.25



 . 

(4)

• The relative error tolerance was set at 3%. That is the relative difference in the strain energy density between two RVEs could not exceed 3%. ˆ B ≥ 100. • For single-layered RVEs L

13

ACS Paragon Plus Environment

ACS Biomaterials Science & Engineering

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

• For in-plane loadings considered in this paper, the difference between multi-layer and single-layer networks was minimal, and therefore it is sufficient to consider single-layer RVEs. The last point is somewhat surprising and deserves an explanation. To this end, consider a single layer characterized by Ni fiber intersections. Let us embed this layer into a 3-D multi-layer network. Then the number of fiber intersection for the layer increases roughly by a factor of three, as fibers within the layer can now have intersections with fibers of the neighboring layers. Thus a characteristic spacing between fiber intersections in the isolated layer is roughly three times longer than that in the embedded layer. This spacing reduction significantly increases the bending stiffness but not the axial stiffness of the individual fiber segments connecting neighboring intersections. The ratio of intersection spacing to fiber diameter determines whether the network strain energy will be dominated by stretch or bending deformations in the fibers. 43,44 As a matter of fact, one can argue that the axial stiffness of the isolated and embedded layers must be roughly the same because the embedment does not change the total length of the fibers within the layer. Thus the only scenario under which the overall stiffnesses of the isolated and embedded layers are equal is when the strain energy due to bending is negligible in comparison to the strain energy due to axial deformation. Indeed, this statement held true not only for simulations used for calibrating the model but in entire set of simulations conducted in this study.

RESULTS Macroscopic Material Behavior In this subsection, we present simulation results showing the macroscopic mechanical response of initially isotropic networks. All results were obtained using the data compiled in Table (1). We present results using the first Piola-Kirchhoff macroscopic stress versus the macroscopic stretch curves obtained for several stretch-controlled and stress-controlled 14

ACS Paragon Plus Environment

Page 14 of 52

Page 15 of 52

ACS Biomaterials Science & Engineering

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 ACS Paragon Plus Environment

ACS Biomaterials Science & Engineering

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

kinematics, including straightening and rotation effects, is most easily studied using data for strain-controlled loadings. The nearly linear nature of the stress-stress curves for these loadings suggests that fibers that are initially nearly aligned with the applied load, and thus undergo small rotations, quickly rotate and form a dominant load bearing subnetwork. Other fibers that are not aligned with the applied load rotate more gradually and carry very little load. The nonlinearity of the stress-strain curves under stress-controlled loading may be understood as a non-proportional traversal of the stretch space. This observation is an agreement with tissue 45 and needled non-woven scaffold models, 46 in which any significant network stretching under uniaxial stress-controlled conditions is preceded by lateral contraction of the network. This phenomenon can be explained by the rotation of fibers toward the applied load direction, as described in the next section. In the 3:1 and 4:1 stress controlled loading cases, the initial negative slope of the y-direction stress vs. stretch curves can be viewed as a typical Poisson’s ratio effect observed in a uniaxial tension test, in which axial extensions correspond to lateral contractions. As more load is applied and more fibers rotate toward the x direction, the effective y-direction stiffness decreases, and the network becomes less capable of resisting the load in that direction. Consequently, the stress vs. stretch curve slope becomes positive and the network begins to expand in the y direction. This and other effects associated with microstructural changes are best described in terms of the evolving ODF, as it is done in the next subsection.

Evolution of Microscopic Quantities of Interest. To elucidate some of the more subtle aspects of the macroscopic mechanical material behavior, we examined quantities of interest pertaining to the underlying microstructure. Just as we could not load the network in a state of pure uniaxial stress due to fiber instabilities, pure shear strain loading proved equally difficult to achieve. Instead, we loaded the network in a state of pure shear strain with some equibiaxial strain superimposed. Starting with the

16

ACS Paragon Plus Environment

Page 16 of 52

Page 17 of 52

ACS Biomaterials Science & Engineering

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 ACS Paragon Plus Environment

ACS Biomaterials Science & Engineering

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 ACS Paragon Plus Environment

Page 18 of 52

Page 19 of 52

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Biomaterials Science & Engineering

fiber stretch gives a sense of evolving network mechanics, and shows how the applied load is distributed among the fibers. The average stretch is clearly maximum along the principal loading direction and minimum along the normal direction (Figure 6). Also, the simulation results are in very close agreement with the affine predictedions. The chain tortuosity in the unloaded configuration is roughly uniform for all orientation angles (Figure 6). In the loaded state, the mean tortuosity of chains along the principal load direction is almost exactly equal to one, as these chains become almost completely straight. Along the normal direction, the tortuosity is large, as these chains compress and just begin to buckle.

Comparison with mechanical responses. To compare the simulation and experimental results for actual scaffold configurations, fiber networks were generated to match the properties of electrospun Poly (ester urethane) urea (ES-PEUU) scaffolds that have been previously tested 30 and analyzed. 14,32 Importantly, experimental trends in the relationships between fiber material properties, network geometry, and macroscopic mechanical behavior are comparable with simulation results. This suggests that we have captured the key aspects of the scaffold mechanics. Three scaffolds with different degrees of fiber alignment were studied. In the electrospinning process, greater mandrel speeds produce more highly aligned scaffolds. 28 The fiber diameter D and scaffold fiber volume fraction φ are also affected by mandrel speed, so scaffolds with different degrees of alignment have different fiber diameters and volume fractions. 30 The experimentally measured 30 orientation density functions for the three scaffolds show the trend of increasing geometric anisotropy with increasing mandrel speed (Figure 7a). The scaffold stress vs. stretch plots reported previously 14,30 are for equibiaxial stress loading. We see that mechanical anisotropy and nonlinearity tend to increase with increasing fiber alignment (Figure 7-b). Single-layer networks (Table 2) were generated with target geometric configurations matching the experimental data. As explained in Section , a single layer is sufficient to represent the 3-D scaffold. 19

ACS Paragon Plus Environment

ACS Biomaterials Science & Engineering

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 ACS Paragon Plus Environment

Page 20 of 52

Page 21 of 52

ACS Biomaterials Science & Engineering

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 ACS Paragon Plus Environment

ACS Biomaterials Science & Engineering

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 22 of 52

behavior. This effect was explored by generating networks with both isotropic and aligned fiber ODFs (Table 3 and Figure 9-a,b). The alignment direction is π2 , so the ODF can be characterized by computing the dispersion about π2 , which is defined as 47 Z π π π D( ) = [1 − cos(ϕ − )]f (ϕ) dϕ 2 2 0

For a purely isotropic network (f (ϕ) = π1 ), D( π2 ) = 1 − aligned to

π , 2

D( π2 ) = 0. The dispersions about

π 2

2 π

(6)

≈ 0.363. For a network perfectly

of the generated isotropic and aligned

networks are 0.367 and 0.203 respectively. Table 3: Target geometry, random walk parameters, and fiber material properties for fiber alignment simulations D(µm) 0.33

φt 0.18

ˆB L 100

sˆ 2.0

nf 100

α π 10

Ef (MPa) 7.5

The macroscopic stress versus stretch curves for both equibiaxial stress and 3:1 biaxial stress loading (Figure 9-c,d) show the significant differences between the isotropic and aligned networks. For the aligned network, there are more fibers along the vertical (y) direction, so under equibiaxial stress, the slope of the stress vs. stretch curve is greater for this direction. Fibers rotate toward the horizontal (x) direction, and the stiffening effect along this direction allows the network to traverse a non-proportional load path in strain space. The resulting stress vs. stretch curve is nonlinear. These effects are even more pronounced for nonequibiaxial load case. The aligned network geometry causes the stress versus stretch curves to become significantly nonlinear, and the contraction in the y direction becomes significant. Another key finding was that fiber tortuosity has a moderate effect on the nonlinearity of the macroscopic stress vs. stretch curves for strain controlled loading. This effect was explored by generating isotropic networks with three different levels of tortuosity, which was controlled by varying the maximum random walk angle (Table 4 and Figure 10-a). The 22

ACS Paragon Plus Environment

Page 23 of 52

ACS Biomaterials Science & Engineering

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 ACS Paragon Plus Environment

ACS Biomaterials Science & Engineering

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Another important finding was that fiber material properties have a direct effect on the resulting network properties, so large differences between fiber materials translate into large differences between the identical networks. Two fiber material models were compared: a linear axial force versus displacement model and an incompressible Yeoh model. 39 The Yeoh model constants were taken from extant AFM measurements (c1 = 1.25M P a, c2 = 4.0M P a, c3 = 0.0M P a). 14,32 For the linear model, Young’s modulus was chosen to be equal to the initial slope of the Yeoh model force versus stretch curve. It is clear that that curves are quite different (Figure 11-a) in the interval of typical fiber stretches. The resulting network stress vs. stretch curves under equibiaxial strain loading (Figure 11-b) look similar to the fiber axial force versus stretch curves.

24

ACS Paragon Plus Environment

Page 24 of 52

Page 25 of 52

ACS Biomaterials Science & Engineering

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 ACS Paragon Plus Environment

ACS Biomaterials Science & Engineering

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 ACS Paragon Plus Environment

Page 26 of 52

Page 27 of 52

ACS Biomaterials Science & Engineering

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 ACS Paragon Plus Environment

ACS Biomaterials Science & Engineering

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

network of extracellular matrix (collagen and elastin) fibers. These soft tissue networks can have a wide variety of geometric structures, ranging from highly isotropic to highly aligned (Figure 1-a) along the circumferential direction in heart valve layers 9,52? or helical in chordae tendinae 53 and arteries. 54 In this work, we introduced a method for modeling materials that have a layered, fibrous network microstructure. A primary application is in the field of tissue engineering, where high fidelity models of electrospun scaffolds are needed to better understand how the network geometry affects the mechanical and biological function of the tissues that are grown on the scaffolds. Scaffold geometry has a strong influence on the tissue’s macroscopic mechanical behavior, cell proliferation and attachment, nutrient and waste flows, and extracellular matrix (ECM) generation. This work was focused on characterizing scaffold network geometry and elucidating the impact of geometry on macroscale mechanics.The effects of fiber alignment, tortuosity, and material properties on macroscopic mechanical behavior of the material have been presented both individually and in combination. We have seen that all three aspects of the scaffold network can have significant effects on the bulk behavior for different load cases. Understanding these effects is essential to the process of designing electrospun scaffolds that have optimal material properties for functional engineered tissues. The ultimate goal is to design better engineered tissues that function mechanically and biologically as well as or better than the native tissue. The random walk algorithm 37 used to generate the scaffold geometry mimics the electrospinning process. The 3D, layered geometry that it generates is a reasonable representation of the actual scaffold microstructure. While the generated geometry does not accurately capture fiber undulations through the thickness or interlayer penetrations, it does capture interlayer intersections as well as many in-plane geometric features of each layer, such as fiber curvature, orientation distribution, and intersections, and the fidelity of this geometric representation represents a significant improvement over current models. The network fibers are modeled as Cosserat rods, in which director vectors are used to keep track of the

28

ACS Paragon Plus Environment

Page 28 of 52

Page 29 of 52

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Biomaterials Science & Engineering

orientation of the fiber cross section. This approach allows large, 3D rotations and axial stretches to be handled accurately. Proper treatment of these deformations is important since the large macroscopic strains experienced by tissues during normal function result in large rotations and stretches of the network fibers at the microscale. The individual fibers are connected using constraints at the points of intersection and at the periodic boundary. Periodic boundary conditions are applied to complete the specification of the network boundary value problem. The finite element method is used to reformulate boundary value problem as a system of nonlinear algebraic equations that is solved using Newton’s method. These equations form the basis of the nonlinear, parallel finite element code that was used to perform the mechanical simulations. The ability to run the simulations in parallel allowed problems involving sufficiently large networks with sufficiently intricate representations of fibrous geometry details to be solved in a reasonable time period. A more complete understanding of the mechanisms underlying the macroscopic behavior of scaffolds can be gained using micromechanical simulations than can be gained from analytical models alone. Insight into the relationship between microstructural evolution and mechanical behavior was gained by monitoring microscopic quantities of interest (QOIs). We first developed an understanding of basic fibrous network phenomena for initially isotropic network geometries. The effects of periodic box size, number of samples, and number of layers on the uncertainty in geometric and mechanical QOIs was quantified. Predictions were found to be reasonably accurate for moderate box and sample sizes. For the in-plane loadings considered in this study, multilayer scaffolds were found to be accurately modeled by a single network layer for the purpose of capturing the macroscopic stress vs. stretch response, despite the fact that the monolayer contains one third the number of fiber-to-fiber contacts as the real scaffold. This finding suggests that the fiber interactions are secondary to the effects of fiber volume fraction and ODF. This conclusion is also supported by the finding that the bending energy contribution to the total strain energy in the network is negligible.

29

ACS Paragon Plus Environment

ACS Biomaterials Science & Engineering

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Observations of the loaded network geometry and analyses all support the conclusions that affine fiber kinematics can be used to predict both the orientations and stretches of the fibers. Most of the networks of interest to the field of tissue engineering have aligned fibers in order to duplicate the properties of native tissues, and simulations of aligned fibrous networks based on realistic geometries from tissue engineering applications are used to build on the isotropic network results. Thus, all network geometries and material properties presented here are designed to be good representations of typical electrospun Poly (ester urethane) urea (ES-PEUU) scaffolds, and the investigations presented here are directed toward the ultimate goal of designing scaffolds having optimal properties for functional tissue engineering. Despite the inability of analytical models to capture the details of local deformations, it was found that the assumption of affine kinematics is a reasonable description of average fiber deformations and can be used to predict the evolution of the fiber ODF and average stretch at each orientation with great accuracy. Fibers always rotated toward the direction of the principal eigenvector of the left stretch tensor and always straighten and bear the most load along this direction. We underscore that the transition to affinity as observed in native tissues 21 is a result of the transition to long fiber-like behavior observed in the present simulations. We believe this is a key result, as it will help with the applications to developing simulated native tissue like behaviors. The effects of fiber alignment, tortuosity, and material model, both individually and in combination, on the macroscopic stretch vs. stretch behavior were explored. For moderately tortuous networks of linear fibers, increasing fiber alignment significantly increases mechanical anisotropy, and significantly increases nonlinearity under stress-controlled loading. For isotropic networks of linear fibers, increasing tortuosity moderately increases nonlinearity at small strains. A linear fiber was compared with a fiber modeled as an incompressible Yeoh hyperelastic material. Large differences in the nonlinearity of fiber response translate into large differences in the network response. These simulations allow the effects of geometry and fiber material properties on macroscale behavior to be studied systematically in a way

30

ACS Paragon Plus Environment

Page 30 of 52

Page 31 of 52

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Biomaterials Science & Engineering

that is virtually impossible using current experimental techniques. The effects of both high alignment and high tortuosity were then investigated together by constructing an aligned network with high tortuosity (Figure 12-a). For aligned networks under equibiaxial stress using a linear fiber material model, the effect of tortuosity is significant. Using a Yeoh material model, subjected to equibiaxial stress, we observed that increasing tortuosity significantly increases material nonlinearity, and the added degree of nonlinearity due to tortuosity is comparable to that of using the Yeoh model. The effects of fiber alignment, tortuosity, and material properties on macroscopic mechanical behavior of the material were presented both individually and in combination. We saw that all three aspects of the scaffold network can have significant effects on the bulk behavior for different load cases. Understanding these effects is essential to the process of designing electrospun scaffolds that have optimal material properties for functional engineered tissues.

Implications and broader considerations. Many of the reasons for poor implant performance or failure of prosthetic biomaterials and engineered tissues remain ill-defined, often being a combination of inadequate or mismatched mechanical properties and biological complexities. Despite the multitude of challenges, many early engineered tissue approaches have shown promising results. For example, heart valve prostheses derived from bovine pericardium or porcine aortic valves have long been used to enhance survival and improve the quality of life of patients presenting with a variety of valvular diseases. 55 Similarly, engineered dermal grafts have successfully been used clinically to treat severe burns or wounds that would otherwise be unable to close and heal properly. 56 Decellularized extracellular matrix scaffolds have also shown successes in regenerating organized tissue after severe musculoskeletal tissue loss or injury. 57 In addition to providing invaluable educational experience to guide future efforts, this incremental progress moves the field ever closer towards the ultimate goal of developing technologies for safer and more efficacious tissue repairs and replacements. 26 31

ACS Paragon Plus Environment

ACS Biomaterials Science & Engineering

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Healthy native tissues undergo intricate, multi-scale modes of deformation which work synergistically with biochemical stimuli to drive physiological responses. In order to mimic native tissue structure and organization, scaffolds characteristic lengths on a scale comparable to those observed in nature need to be produced and controllable techniques to achieve those are necessary. Electrospun scaffolds have been employed extensively in tissue engineering as they are able to replicate some of the multi-scale properties that are important characteristic features of native tissues (e.g. among many others, annulus fibrosus, 58 native cartilage, 59 and anistropic cardiac tissues. 60 Electrospun constructs are amenable to modification during and after fabrication to introduce functionality or modify microstructure and mechanical response. 27 Functional groups of peptides can be introduced through surface treatments after fabrication or grafted to the polymer chains prior to solvent processing. Fibers can be patterned or aligned to encourage contact guidance of seeded cells and to produce tunable tensile properties and structural and mechanical anisotropy. However, while electrospinning permits fabrication of biodegradable matrices that resemble the scale, architecture, and mechanical behavior of the native ECM, 30 this morphology also results in pore sizes that are generally smaller (< 5µm) and more irregular than those produced by other non-fibrous production methods. Thus, achieving high cellular density and infiltration remains challenging and several approaches have been proposed to overcome this limitation such as the alteration of manufacturing parameters (e.g. fiber diameter and packing), inclusion of native ECM proteins, labile porogens, or sacrifial fiber populations, and concurrent electrospining of scaffold fibers and electrospraying of cell populations suspended in culture media. 61,62 Fiber diameter and alignment are able to modulate cellular morphology, i.e. projected cell area, aspect ratio, and long axis length. 22 Sub-micron fiber diameters induce diminished cell adhesion and spreading, possibly due ineffective focal adhesions (which can be larger than 1µm). In contrast, increased fiber alignment promoted increased cell spreading. Adhesion and spreading of anchorage-dependent cells is a prerequisite for cell viability and

32

ACS Paragon Plus Environment

Page 32 of 52

Page 33 of 52

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Biomaterials Science & Engineering

proliferation 63 as cells perceive and react to their mechanical environment through adhesion and cellular deformations induced by and into the surrounding mechanical environment. These complex interactions likely play a critical role in the mechanotransduction of proteins required for cell viability and proliferation, and have direct implications towards the development of engineered tissues. 15 Extracellular matrix micro-structure, mechanical properties and cell-ECM interactions play an important role in diverse biological processes, e.g. pathological remodeling of cancer development, 64 or in the vasculature, in atherogenesis and atherosclerosis. 65 Substrate stiffness is well-known to play a substantial role in directing cell differentiation and cell general behavior. 66 Strains imposed to the cells at the cellular and tissue scales (either in vivo or in vitro engineered tissue systems) are propagated into the cell acting through the cytoskeleton and lead to altered stress levels on the nucleus, induce nuclear shape change, and ultimately affect cell function and causing modifications in gene expression and protein synthesis. 67–69 Experimental observations of tissues grown in bioreactors have shown that tissue composition, morphology, and mechanical properties can all be altered by controlling static and dynamic compression and hydrodynamic conditions. 70 Macro-scale coupled multi-physics simulations of tissues grown in the bioreactor environment have also confirmed these effects on tissue growth, as well as the effects of bending loading and nutrient and oxygen concentrations in the fluid phase. 71 The micro-scale scaffold simulations presented herein could in principal be extended to include the coupled mechanical, chemical, and biological phenomena present in bioreactors or even in vivo. The geometry, stiffness, and deformations of the scaffold fibers all provide mechanical cues to the attached cells. These cues can be used to control stem cell differentiation and matrix biosynthesis, 26,72 but the micro-scale mechanisms are still poorly understood. The role of trans-membrane proteins in the adhesion of cells to matrix and scaffold fibers, the transmission of stress to cell nuclei through the cytoskeleton, the chemical, structural, hydrodynamic effects on protein syntheses, and the migration of cells through the scaffold are all phenomena that can be studied in detail using micro-scale

33

ACS Paragon Plus Environment

ACS Biomaterials Science & Engineering

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

models. These models could be combined with macro-scale models of tissue function to create complete multiscale, multiphysics models of the complex phenomena that occur during engineering tissue growth. Alternatively, it would be possible to perform direct numerical simulations of features that extend to the full tissue length scale. Networks containing many layers and reaching the thickness of a heart valve leaflet could be loaded in flexure, torsion, or out of plane shear. The importance of scaffold-cell mechanical coupling becomes readily apparent in the rational design of engineered tissues. Moreover, the unique micromechanics of various scaffolds induce different cell deformation responses, which could correlate to substantial changes in cell proliferation and function. To summarize, the understanding of how microstructure influences cell behavior is important for the following 1. Defining fabrication variables and microstructural features to obtain better performing tissue surrogates that are able to recapitulate native tissues behavior at multiple lengthscales (cellular, tissue and organ), not only structurally and mechanically (in particular, their anisotropy), but also biochemically and biologically. 2. Establishing two-way coupling strategies and multi-scale models to study not only the effects of cells and ECM surrounding the scaffold in the micro- and macro-mechanics of the engineered tissue composite, but also the propagation of macro-scale mechanical stimuli to the micro-scale environment to trigger biological responses that distinguish implant evolution either into healthy integration or pathological rejection/failure. 3. Developing better design tools and methodologies to predict time-evolving scaffold properties and structure of scaffolds (which are desirably biodegradable), not only to tailor the scaffold with sufficient properties to function acutely, but also stress-shielding and smooth transition of load from the degrading scaffold to the developing de novo tissues.

34

ACS Paragon Plus Environment

Page 34 of 52

Page 35 of 52

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Biomaterials Science & Engineering

Future work. Overall, the simulations presented here represent an important advancement in the field of tissue engineering. The accurate representation of detailed fibrous network geometry permits a greater understanding of the complex mechanisms underlying the macroscopic behavior unique to these biomaterials. Furthermore, simulations of scaffold fiber networks form the basis for understanding how interactions with cellular and ECM phases contribute to the growth, remodeling, and ultimate mechanical and biological behavior of the entire engineered tissue. Insights gained from such simulations can significantly aid the process of designing scaffold network geometries that result in engineered tissues that function as well as or better than the native tissues they are intended to replace. Looking forward, the current work may be expanded and enhanced in many ways. The Cosserat rod theory used in this work to model the individual fibers includes the effects of large rotations and large stretches of the fibers, as well as the effects of bending and torsion. 29 However, scaffold networks of practical interest have been found to be axial dominated, with negligible contributions from bending and torsion to the total strain energy density. While fiber bending stiffness is needed to ensure a nonsingular stiffness matrix, it is possible that formulations based simplified rod theories may be used without affecting the accuracy of the quantities of interest. For example, for the terms associated with bending in the expression for the virtual strain energy density, the assumption of small strains could be invoked while still accounting for large fiber rotations and stretches. Since the equations in the current formulation are quite complex, the use of simplified rod theories that still capture the important phenomena would be justified. The basic parameters, such as volume fraction and ODF, that characterize the geometry could be viewed as state variables that evolve as the scaffold network is loaded and the fibers deform. Simulations could be used to help develop equations that describe the evolution of these variables. Many phenomenological 73 and structural 7,16,18 constitutive models of soft tissues have been developed to describe and explain the complex behavior of soft tissues. 35

ACS Paragon Plus Environment

ACS Biomaterials Science & Engineering

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Mesoscale Structural Constitutive Model (MSSCM) 16,19 is known to work well for soft tissues partly because it is a good representation of the microstructure and functioning of these materials . Although electrospun scaffolds are physically very different from soft tissues, the MSSCM has been shown to predict the bulk mechanical behavior of these materials with good fidelity. 30 Simulations could be used to explore the correspondence between soft tissue models and fibrous networks in an attempt to help explain why soft tissue models work as well as they do in predicting the macroscopic mechanical behavior of electrospun scaffolds. In order to better explain how the MSSCM model predictions relate to the simulation results, it would be helpful to consider the correspondence between the parameters of the MSSCM and the parameters that describe the real network geometry. Making these connections would also be helpful to the process of scaffold design. Given MSSCM model parameters that result in a known, desired macroscale behavior, it would be useful to relate those parameters to a corresponding network fiber geometry and material that has the same desired behavior. This network description would provide a target geometry and material for scaffold manufacturers. Another logical step would be to include cellular and ECM phases in the model. Tissue growth and remodeling is known to be strongly influenced by applied stress fields. Ideal engineered heart valve tissues have a dense, organized ECM phase made of aligned collagen and elastin fibers. Replicating the cyclic loading conditions found in vivo may be important in producing these ideal structures. Experimental observations of tissues grown in bioreactors have shown that tissue composition, morphology, and mechanical properties can all be altered by controlling static and dynamic compression and hydrodynamic conditions. 74 Macroscale coupled multiphysics simulations of tissues grown in the bioreactor environment have also confirmed these effects on tissue growth, as well as the effects of bending loading and nutrient and oxygen concentrations in the fluid phase. 71 The microscale scaffold simulations presented herein could in principal be extended to include the coupled mechanical, chemical, and biological phenomena present in bioreactors or even in vivo. The geometry, stiffness, and deformations of the scaffold fibers all provide

36

ACS Paragon Plus Environment

Page 36 of 52

Page 37 of 52

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Biomaterials Science & Engineering

mechanical cues to the attached cells. These cues can be used to control stem cell differentiation and matrix biosynthesis, 62,72 but the microscale mechanisms are still poorly understood. The role of transmembrane proteins in the adhesion of cells to matrix and scaffold fibers, the transmission of stress to cell nuclei through the cytoskeleton, the chemical, structural, hydrodynamic effects on protein syntheses, and the migration of cells through the scaffold are all phenomena that can be studied in detail using microscale simulations. These simulations could be combined with macroscale simulations that treat the tissue as a continuum to create complete multiscale, multiphysics models of the complex phenomena that occur during tissue growth. Alternatively, it would be possible to perform direct numerical simulations of features that extend to the full tissue length scale. Networks containing many layers and reaching the thickness of a heart valve leaflet could be loaded in flexure, torsion, or out of plane shear. The methods introduced may be applied to materials other than electrospun scaffolds. Hydrogel and needled nonwoven tissue scaffolds also have a fibrous network microstructure, and this work could be even more broadly applicable to other materials, such as the extracellular matrix of soft tissues, the actin cytoskeleton, some engineering composite materials, and many textiles. Finally, formal methods could be employed to solve the inverse problem of finding the optimal microstructure that produces materials with desired properties. As biomaterial manufacturing technologies improve, it may one day be possible to produce tissues repeatably that are superior to those produced by nature. Network mechanics simulations will undoubtedly play an important role in understanding and designing the intricate microstructural details of such advanced materials.

37

ACS Paragon Plus Environment

ACS Biomaterials Science & Engineering

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Acknowledgement The authors gratefully acknowledge NIH grant HL068816, and Dr. Joao S. Soares for careful review of the manuscript.

38

ACS Paragon Plus Environment

Page 38 of 52

Page 39 of 52

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Biomaterials Science & Engineering

References References (1) Woo, S. L. Y.; Seguchi, Y. Tissue Engineering - 1989 ; ASME BED: New York, 1989; Vol. 14. (2) Butler, D. L.; Goldstein, S. A.; Guilak, F. J Biomech Eng 2000, 122, 570–5. (3) Guilak, F.; Butler, D. L.; Goldstein, S. A.; Baaijens, F. P. J Biomech 2014, 47, 1933–40. (4) Sacks, M. J Elasticity 2000, 61, 199–246. (5) Criscione, J.; Sacks, M.; Hunter, W. Journal of Biomechanical Engineering 2003, 125, 94–99. (6) Holzapfel, G. A.; Ogden, R. Philosophical Transactions. Series A, Mathematical, Physical, and Engineering Sciences 2009, 367, 3445–3475. (7) Gasser, T. C.; Ogden, R. W.; Holzapfel, G. A. J R Soc Interface 2006, 3, 15–35. (8) Holzapfel, G. A.; Gasser, T. C. Journal of Elasticity 2000, 61, 1–48. (9) Rego, B. V.; Sacks, M. S. Journal of biomechanics 2017, 54, 88–95. (10) Avazmohammadi, R.; Hill, M. R.; Simon, M. A.; Zhang, W.; Sacks, M. S. Biomechanics and modeling in mechanobiology 2017, 16, 561–581. (11) Sacks, M. S.; Zhang, W.; Wognum, S. Interface Focus 2016, 6, 20150090. (12) Preethi, L. C.; Victor, H. B. Journal of Biomechanical Engineering 2007, 129, 137–147. (13) Stylianopoulos, T.; Barocas, V. H. Comput Meth Appl Mech Eng 2007, 196, 2981–2990. (14) D’Amore, A.; Amoroso, N.; Gottardi, R.; Hobson, C.; Carruthers, C.; Watkins, S.; Wagner, W. R.; Sacks, M. S. Journal of the mechanical behavior of biomedical materials 2014, 39, 146–161. 39

ACS Paragon Plus Environment

ACS Biomaterials Science & Engineering

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(15) Stella, J. A.; Liao, J.; Hong, Y.; David Merryman, W.; Wagner, W. R.; Sacks, M. S. Biomaterials 2008, 29, 3228–36. (16) Lanir, Y. J Biomech 1983, 16, 1–12. (17) Lanir, Y. J Biomech Eng 1983, 105, 374–80. (18) Sacks, M. S. J Biomech Eng 2003, 125, 280–7. (19) Fan, R.; Sacks, M. S. Journal of Biomechanics 2014, 47, 2043–2054. (20) Sacks, M. S.; Sun, W. Annu Rev Biomed Eng 2003, 5, 251–84. (21) Lee, C. H.; Zhang, W.; Liao, J.; Carruthers, C. A.; Sacks, J. I.; Sacks, M. S. Biophys J 2015, 108, 2074–87. (22) Bashur, C. A.; Dahlgren, L. A.; Goldstein, A. S. Biomaterials 2006, 27, 5681–5688. (23) Nerurkar, N. L.; Elliott, D. M.; Mauck, R. L. Journal of Orthopaedic Research 2007, 25, 1018–1028. (24) Stylianopoulos, A.; Bashur, C.; Goldstein, A.; Guelcher, S.; Barocas, V. Journal of Mechanical Behavior of Biomedical Materials 2008, 1, 326–335. (25) Stella, J. A.; Wagner, W. R.; Sacks, M. S. Journal of Biomedical Materials Research Part A 2010, 93, 1032–42. (26) Stella, J. A.; D’Amore, A.; Wagner, W. R.; Sacks, M. S. Acta Biomaterialia 2010, 6, 2365–2381. (27) Amoroso, N. J.; D’Amore, A.; Hong, Y.; Rivera, C. P.; Sacks, M. S.; Wagner, W. R. Acta biomaterialia 2012, 8, 4268–4277. (28) Amoroso, N. J.; D’Amore, A.; Hong, Y.; Wagner, W. R.; Sacks, M. S. Advanced Materials 2011, 23, 106–111. 40

ACS Paragon Plus Environment

Page 40 of 52

Page 41 of 52

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Biomaterials Science & Engineering

(29) Carleton, J. B.; D’Amore, A.; Feaver, K. R.; Rodin, G. J.; Sacks, M. S. Acta biomaterialia 2015, 12, 93–101. (30) Courtney, T.; Sacks, M. S.; Stankus, J.; Guan, J.; Wagner, W. R. Biomaterials 2006, 27, 3631–8. (31) Stankus, J. J.; Guan, J.; Wagner, W. R. J Biomed Mater Res 2004, 70A, 603–14. (32) D’Amore, A.; Stella, J. A.; Wagner, W. R.; Sacks, M. S. Biomaterials 2010, 31, 5345– 54. (33) Villaggio, P. Mathematical Models for Elastic Structures; Cambridge University Press, 1997. (34) ABAQUS

6.12

Theory

Manual.

Dassault

Systèmes,

http://abaqus.software.polimi.it/v6.12/books/stm/default.htm. (35) Abhilash, A.; Baker, B. M.; Trappmann, B.; Chen, C. S.; Shenoy, V. B. Biophysical Journal 2014, 107, 1829–1840. (36) Ban, E.; Barocas, V. H.; Shepard, M. S.; Picu, R. C. Journal of the Mechanics and Physics of Solids 2016, 87, 38–50. (37) Carleton, J. B. Microscale Modeling of Layered Fibrous Networks with Applications to Biomaterials for Tissue Engineering. Ph.D. thesis, The University of Texas at Austin, 2015. (38) MacNeal, R. H. Composites and Structures 1978, 8, 175–183. (39) Yeoh, O. Rubber Chemistry and Technology 1993, 66, 754–771. (40) Balay, S.; Gropp, W. D.; McInnes, L. C.; Smith, B. F. Efficient Management of Parallelism in Object Oriented Numerical Software Libraries. Modern Software Tools in Scientific Computing. 1997; pp 163–202. 41

ACS Paragon Plus Environment

ACS Biomaterials Science & Engineering

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(41) Balay, S.; Brown, J.; Buschelman, K.; Eijkhout, V.; Gropp, W. D.; Kaushik, D.; Knepley, M. G.; McInnes, L. C.; Smith, B. F.; Zhang, H. PETSc Users Manual ; 2012. (42) Balay, S.; Brown, J.; Buschelman, K.; Eijkhout, V.; Gropp, W. D.; Kaushik, D.; Knepley, M. G.; McInnes, L. C.; Smith, B. F.; Zhang, H. PETSc Web page. 2012; mcs.anl.gov/petsc. (43) Head, D.; Levine, A.; MacKintosh, F. Physical Review E 2003, 68, 1–15. (44) Broedersz, C.; MacKintosh, F. Reviews of Modern Physics 2014, 86, 995–1036. (45) Billiar, K. L.; Sacks, M. S. J Biomech Eng 2000, 122, 327–35. (46) Soares, J. S.; Zhang, W.; Sacks, M. S. Acta biomaterialia 2017, 51, 220–236. (47) Mardia, K. V.; Jupp, P. E. Directional Statistics; John Wiley & Sons, LTD, 2000. (48) Onck, P. R.; Koeman, T.; van Dillen, T.; van der Giessen, E. Phys Rev Lett 2005, 95, 178102. (49) Kabla, A.; Mahadevan, L. Journal of the Royal Society Interface 2007, 4, 99–106. (50) Ban, E.; Barocas, V. H.; Shephard, M. S.; Picu, C. R. Journal of Applied Mechanics 2016, 83, 1–7. (51) Vacanti, C.; Mikos, A. Tissue Engineering 1995, 1, 1–2. (52) Vesely, I.; Noseworthy, R. Journal of Biomechanics 1992, 25, 101–113. (53) Freed, A. D.; Doehring, T. C. J Biomech Eng 2005, 127, 587–93. (54) Holzapfel, G. A.; Ogden, R. W. Proceedings of the Royal Society a-Mathematical Physical and Engineering Sciences 2010, 466, 1551–1596. (55) Soares, J. S.; Feaver, K. R.; Zhang, W.; Kamensky, D.; Aggarwal, A.; Sacks, M. S. Cardiovascular Engineering and Technology 2016, 7, 309–351. 42

ACS Paragon Plus Environment

Page 42 of 52

Page 43 of 52

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Biomaterials Science & Engineering

(56) Kaefer, M.; Curran, M.; Treves, S. T.; Bauer, S.; Hendren, W. H.; Peters, C. A.; Atala, A.; Diamond, D.; Retik, A. Pediatrics 2000, 105, 800–4. (57) Gilbert, T. W.; Stewart-Akers, A. M.; Simmons-Byrd, A.; Badylak, S. F. J Bone Joint Surg Am 2007, 89, 621–30. (58) Nerurkar, N. L.; Baker, B. M.; Sen, S.; Wible, E. E.; Elliott, D. M.; Mauck, R. L. Nature materials 2009, 8, 986–992. (59) Moutos, F. T.; Freed, L. E.; Guilak, F. Nat Mater 2007, 6, 162–7. (60) Engelmayr, G. C.; Cheng, M.; Bettinger, C. J.; Borenstein, J. T.; Langer, R.; Freed, L. E. Nat Mater 2008, 7, 1003–1010. (61) Stankus, J. J.; Guan, J.; Fujimoto, K.; Wagner, W. R. Biomaterials 2006, 27, 735–44. (62) A, D. A.; Soares, J. S.; Stella, J. A.; Zhang, W.; Amoroso, N. J.; Mayer, J. E., Jr.; Wagner, W. R.; Sacks, M. S. J Mech Behav Biomed Mater 2016, 62, 619–35. (63) Chen, C. S.; Mrksich, M.; Huang, S.; Whitesides, G. M.; Ingber, D. E. Biotechnology progress 1998, 14, 356–363. (64) Huang, S.; Ingber, D. E. Cancer cell 2005, 8, 175–176. (65) Nerem, R. M. Journal of biomechanical engineering 1992, 114, 274–282. (66) Engler, A. J.; Sen, S.; Sweeney, H. L.; Discher, D. E. Cell 2006, 126, 677–689. (67) Thomas, C. H.; Collier, J. H.; Sfeir, C. S.; Healy, K. E. Proc Natl Acad Sci U S A 2002, 99, 1972–7. (68) Guilak, F. J Biomech 1995, 28, 1529–1541. (69) Huang, H. Y.; Liao, J.; Sacks, M. S. J Biomech Eng 2007, 129, 880–89.

43

ACS Paragon Plus Environment

ACS Biomaterials Science & Engineering

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(70) Vunjak-Novakovic, G.; Meinel, L.; Altman, G.; Kaplan, D. Orthod Craniofac Res 2005, 8, 209–18. (71) Soares, J. S.; Sacks, M. S. Biomech Model Mechanobiol 2016, 15, 293–316. (72) Buschmann, M.; Gluzband, Y.; Grodzinsky, A.; Hunziker, E. Journal of Cell Science 1995, 108, 1497–1508. (73) Fung, Y. C. Biomechanics: Mechanical Properties of Living Tissues, 2nd ed.; Springer Verlag: New York, 1993. (74) Vunjak-Novakovic, G.; Martin, I.; Obradovic, B.; Treppo, S.; Grodzinsky, A.; Langer, R.; Freed, L. Journal of Orthopaedic Research 1999, 17, 130–138.

44

ACS Paragon Plus Environment

Page 44 of 52

Page 45 of 52

ACS Biomaterials Science & Engineering

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 ACS Paragon Plus Environment

ACS Biomaterials Science & Engineering

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

7

Page 46 of 52

References

References [1] A, D. A., Soares, J. S., Stella, J. A., Zhang, W., Amoroso, N. J., Mayer, Jr., J. E., Wagner, W. R., and Sacks, M. S. Large strain stimulation promotes extracellular matrix production and stiffness in an elastomeric scaffold model. J Mech Behav Biomed Mater 62 (Sept. 2016), 619–35. [2] Abhilash, A., Baker, B. M., Trappmann, B., Chen, C. S., and Shenoy, V. B. Remodeling of fibrous extracellular matrices by contractile cells: Predictions from discrete fiber network simulations. Biophysical Journal 107 (2014), 1829–1840. [3] Amoroso, N. J., D’Amore, A., Hong, Y., Rivera, C. P., Sacks, M. S., and Wagner, W. R. Microstructural manipulation of electrospun scaffolds for specific bending stiffness for heart valve tissue engineering. Acta biomaterialia 8 (Dec. 2012), 4268–4277. [4] Amoroso, N. J., D’Amore, A., Hong, Y., Wagner, W. R., and Sacks, M. S. Elastomeric electrospun polyurethane scaffolds: the interrelationship between fabrication conditions, fiber topology, and mechanical properties. Advanced Materials 23 (Jan. 2011), 106–111. [5] Avazmohammadi, R., Hill, M. R., Simon, M. A., Zhang, W., and Sacks, M. S. A novel constitutive model for passive right ventricular myocardium: evidence for myofibercollagen fiber mechanical coupling. Biomechanics and modeling in mechanobiology 16 (Apr. 2017), 561–581. [6] Balay, S., Brown, J., Buschelman, K., Eijkhout, V., Gropp, W. D., Kaushik, D., Knepley, M. G., McInnes, L. C., Smith, B. F., and Zhang, H. PETSc Users Manual. Tech. Rep. ANL-95/11 - Revision 3.3, Argonne National Laboratory, 2012. [7] Balay, S., Brown, J., Buschelman, K., Eijkhout, V., Gropp, W. D., Kaushik, D., Knepley, M. G., McInnes, L. C., Smith, B. F., and Zhang, H. PETSc web page, 2012. mcs.anl.gov/petsc. [8] Balay, S., Gropp, W. D., McInnes, L. C., and Smith, B. F. Efficient management of parallelism in object oriented numerical software libraries. In Modern Software Tools in Scientific Computing (1997), E. Arge, A. Bruaset, and H. Langtangen, Eds., Birkh¨auser Press, pp. 163–202.

ACS Paragon Plus Environment

31

Page 47 of 52

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Biomaterials Science & Engineering

[9] Ban, E., Barocas, V. H., Shepard, M. S., and Picu, R. C. Softening in random networks of non-identical beams. Journal of the Mechanics and Physics of Solids 87 (2016), 38–50. [10] Ban, E., Barocas, V. H., Shephard, M. S., and Picu, C. R. Effect of fiber crimp on the elasticity of random fiber networks with and without embedding matrices. Journal of Applied Mechanics 83, 041008 (2016), 1–7. [11] Bashur, C. A., Dahlgren, L. A., and Goldstein, A. S. Effect of fiber diameter and orientation on fibroblast morphology and proliferation on electrospun poly(d,l-lactic-co-glycolic acid) meshes. Biomaterials 27, 33 (2006), 5681–5688. [12] Billiar, K. L., and Sacks, M. S.

Biaxial mechanical properties of the native and

glutaraldehyde-treated aortic valve cusp: Part II–A structural constitutive model. J Biomech Eng 122, 4 (Aug. 2000), 327–35. [13] Broedersz, C., and MacKintosh, F. Modeling semiflexible polymer networks. Reviews of Modern Physics 86 (2014), 995–1036. [14] Buschmann, M., Gluzband, Y., Grodzinsky, A., and Hunziker, E. Mechanical compression modulates matrix biosynthesis in chondrocyte agarose culture. Journal of Cell Science 108, 4 (April 1995), 1497–1508. [15] Butler, D. L., Goldstein, S. A., and Guilak, F. Functional tissue engineering: the role of biomechanics. J Biomech Eng 122, 6 (Dec. 2000), 570–5. [16] Carleton, J. B. Microscale Modeling of Layered Fibrous Networks with Applications to Biomaterials for Tissue Engineering. PhD thesis, The University of Texas at Austin, August 2015. [17] Carleton, J. B., D’Amore, A., Feaver, K. R., Rodin, G. J., and Sacks, M. S. Geometric characterization and simulation of planar layered elastomeric fibrous biomaterials. Acta biomaterialia 12 (Jan. 2015), 93–101. [18] Chen, C. S., Mrksich, M., Huang, S., Whitesides, G. M., and Ingber, D. E. Micropatterned surfaces for control of cell shape, position, and function. Biotechnology progress 14 (1998), 356–363. [19] Courtney, T., Sacks, M. S., Stankus, J., Guan, J., and Wagner, W. R. Design and analysis of tissue engineering scaffolds that mimic soft tissue mechanical anisotropy. Biomaterials 27, 19 (July 2006), 3631–8. ACS Paragon Plus Environment

32

ACS Biomaterials Science & Engineering

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 48 of 52

[20] Criscione, J., Sacks, M., and Hunter, W. Experimentally tractable, pseudo-elastic constitutive law for biomembranes: I. theory. Journal of Biomechanical Engineering 125, 1 (February 2003), 94–99. [21] D’Amore, A., Amoroso, N., Gottardi, R., Hobson, C., Carruthers, C., Watkins, S., Wagner, W. R., and Sacks, M. S. From single fiber to macro-level mechanics: A structural finite-element model for elastomeric fibrous biomaterials. Journal of the mechanical behavior of biomedical materials 39 (Nov. 2014), 146–161. [22] D’Amore, A., Stella, J. A., Wagner, W. R., and Sacks, M. S. Characterization of the complete fiber network topology of planar fibrous tissues and scaffolds. Biomaterials 31, 20 (July 2010), 5345–54. [23] Dassault

`mes. Syste

ABAQUS

6.12

Theory

Manual.

http://abaqus.software.polimi.it/v6.12/books/stm/default.htm. [24] Engelmayr, G. C., Cheng, M., Bettinger, C. J., Borenstein, J. T., Langer, R., and Freed, L. E. Accordion-like honeycombs for tissue engineering of cardiac anisotropy. Nat Mater 7, 12 (2008), 1003–1010. [25] Engler, A. J., Sen, S., Sweeney, H. L., and Discher, D. E. Matrix elasticity directs stem cell lineage specification. Cell 126 (Aug. 2006), 677–689. [26] Fan, R., and Sacks, M. S. Simulation of planar soft tissues using a structural constitutive model: Finite element implementation and validation. Journal of Biomechanics 47 (2014), 2043–2054. [27] Freed, A. D., and Doehring, T. C. Elastic model for crimped collagen fibrils. J Biomech Eng 127, 4 (Aug. 2005), 587–93. [28] Fung, Y. C. Biomechanics: Mechanical Properties of Living Tissues, 2nd ed. Springer Verlag, New York, 1993. [29] Gasser, T. C., Ogden, R. W., and Holzapfel, G. A. Hyperelastic modelling of arterial layers with distributed collagen fibre orientations. J R Soc Interface 3, 6 (Feb. 2006), 15–35. [30] Gilbert, T. W., Stewart-Akers, A. M., Simmons-Byrd, A., and Badylak, S. F. Degradation and remodeling of small intestinal submucosa in canine Achilles tendon repair. J Bone Joint Surg Am 89, 3 (Mar. 2007), 621–30. [31] Guilak, F. Compression-induced changes in the shape and volume of the chondrocyte nucleus. J Biomech 28, 12 (dec 1995), 1529–1541. ACS Paragon Plus Environment

33

Page 49 of 52

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Biomaterials Science & Engineering

[32] Guilak, F., Butler, D. L., Goldstein, S. A., and Baaijens, F. P. Biomechanics and mechanobiology in functional tissue engineering. J Biomech 47, 9 (June 2014), 1933–40. [33] Head, D., Levine, A., and MacKintosh, F. Distinct regimes of elastic response and deformation modes of cross-linked cytoskeletal and semiflexible polymer networks. Physical Review E 68, 061907 (December 2003), 1–15. [34] Holzapfel, G. A., and Gasser, T. C. A new constitutive framework for arterial wall mechanics and a comparative study of material models. Journal of Elasticity 61 (2000), 1–48. [35] Holzapfel, G. A., and Ogden, R. Constitutive modelling of passive myocardium: A structurally based framework for material characterization. Philosophical Transactions. Series A, Mathematical, Physical, and Engineering Sciences 367, 1902 (September 2009), 3445–3475. [36] Holzapfel, G. A., and Ogden, R. W. Constitutive modelling of arteries. Proceedings of the Royal Society a-Mathematical Physical and Engineering Sciences 466, 2118 (June 2010), 1551–1596. [37] Huang, H. Y., Liao, J., and Sacks, M. S. In-situ deformation of the aortic valve interstitial cell nucleus under diastolic loading. J Biomech Eng 129, 6 (Dec. 2007), 880–89. [38] Huang, S., and Ingber, D. E. Cell tension, matrix mechanics, and cancer development. Cancer cell 8 (Sept. 2005), 175–176. [39] Kabla, A., and Mahadevan, L. Nonlinear mechanics of soft fibrous networks. Journal of the Royal Society Interface 4, 12 (2007), 99–106. [40] Kaefer, M., Curran, M., Treves, S. T., Bauer, S., Hendren, W. H., Peters, C. A., Atala, A., Diamond, D., and Retik, A. Sibling vesicoureteral reflux in multiple gestation births. Pediatrics 105, 4 Pt 1 (Apr. 2000), 800–4. [41] Lanir, Y. Constitutive equations for fibrous connective tissues. J Biomech 16 (1983), 1–12. [42] Lanir, Y. Constitutive equations for the lung tissue. J Biomech Eng 105, 4 (Nov. 1983), 374–80. [43] Lee, C. H., Zhang, W., Liao, J., Carruthers, C. A., Sacks, J. I., and Sacks, M. S. On the presence of affine fibril and fiber kinematics in the mitral valve anterior leaflet. Biophys J 108, 8 (Apr. 2015), 2074–87. [44] MacNeal, R. H. A simple quadrilateral shell element. Composites and Structures 8 (1978), 175–183. ACS Paragon Plus Environment

34

ACS Biomaterials Science & Engineering

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 50 of 52

[45] Mardia, K. V., and Jupp, P. E. Directional Statistics. John Wiley & Sons, LTD, 2000. [46] Moutos, F. T., Freed, L. E., and Guilak, F. A biomimetic three-dimensional woven composite scaffold for functional tissue engineering of cartilage. Nat Mater 6, 2 (Feb. 2007), 162–7. [47] Nerem, R. M. Vascular fluid mechanics, the arterial wall, and atherosclerosis. Journal of biomechanical engineering 114 (Aug. 1992), 274–282. [48] Nerurkar, N. L., Baker, B. M., Sen, S., Wible, E. E., Elliott, D. M., and Mauck, R. L. Nanofibrous biologic laminates replicate the form and function of the annulus fibrosus. Nature materials 8 (Dec. 2009), 986–992. [49] Nerurkar, N. L., Elliott, D. M., and Mauck, R. L. Mechanics of oriented electrospun nanofibrous scaffolds for annulus fibrosus tissue engineering. Journal of Orthopaedic Research 25, 8 (August 2007), 1018–1028. [50] Onck, P. R., Koeman, T., van Dillen, T., and van der Giessen, E. Alternative explanation of stiffening in cross-linked semiflexible networks. Phys Rev Lett 95, 17 (Oct. 2005), 178102. [51] Preethi, L. C., and Victor, H. B. Deterministic Material-Based Averaging Theory Model of Collagen Gel Micromechanics. Journal of Biomechanical Engineering 129, 2 (2007), 137–147. [52] Rego, B. V., and Sacks, M. S. A functionally graded material model for the transmural stress distribution of the aortic valve leaflet. Journal of biomechanics 54 (Mar. 2017), 88–95. [53] Sacks, M. Biaxial mechanical evaluation of planar biological materials. J Elasticity 61 (2000), 199–246. [54] Sacks, M. S. Incorporation of experimentally-derived fiber orientation into a structural constitutive model for planar collagenous tissues. J Biomech Eng 125, 2 (Apr. 2003), 280–7. [55] Sacks, M. S., and Sun, W. Multiaxial Mechanical Behavior of Biological Materials. Annu Rev Biomed Eng 5 (Apr. 2003), 251–84. [56] Sacks, M. S., Zhang, W., and Wognum, S. A novel fibre-ensemble level constitutive model for exogenous cross-linked collagenous tissues. Interface Focus 6, 1 (feb 2016), 20150090. [57] Soares, J. S., Feaver, K. R., Zhang, W., Kamensky, D., Aggarwal, A., and Sacks, M. S. Biomechanical behavior of bioprosthetic heart valve heterograft tissues: Characterization, simulation, and performance. Cardiovascular Engineering and Technology 7, 4 (Dec 2016), 309–351. ACS Paragon Plus Environment

35

Page 51 of 52

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Biomaterials Science & Engineering

[58] Soares, J. S., and Sacks, M. S. A triphasic constrained mixture model of engineered tissue formation under in vitro dynamic mechanical conditioning. Biomech Model Mechanobiol 15, 2 (Apr. 2016), 293–316. [59] Soares, J. S., Zhang, W., and Sacks, M. S. A mathematical model for the determination of forming tissue moduli in needled-nonwoven scaffolds. Acta biomaterialia 51 (Mar. 2017), 220–236. [60] Stankus, J. J., Guan, J., Fujimoto, K., and Wagner, W. R. Microintegrating smooth muscle cells into a biodegradable, elastomeric fiber matrix. Biomaterials 27, 5 (2006), 735–44. [61] Stankus, J. J., Guan, J., and Wagner, W. R. Fabrication of biodegradable elastomeric scaffolds with sub-micron morphologies. J Biomed Mater Res 70A, 4 (Sept. 2004), 603–14. [62] Stella, J. A., D’Amore, A., Wagner, W. R., and Sacks, M. S. On the Biomechanical Function of Scaffolds for Engineering Load-Bearing Soft Tissues. Acta Biomaterialia 6, 7 (2010), 2365–2381. [63] Stella, J. A., Liao, J., Hong, Y., David Merryman, W., Wagner, W. R., and Sacks, M. S. Tissue-to-cellular level deformation coupling in cell micro-integrated elastomeric scaffolds. Biomaterials 29, 22 (Aug. 2008), 3228–36. [64] Stella, J. A., Wagner, W. R., and Sacks, M. S. Scale-dependent fiber kinematics of elastomeric electrospun scaffolds for soft tissue engineering. Journal of Biomedical Materials Research Part A 93, 3 (June 2010), 1032–42. [65] Stylianopoulos, A., Bashur, C., Goldstein, A., Guelcher, S., and Barocas, V. Computational predictions of the tensile properties of electrospun fibre meshes: Effect of fibre diameter and fibre orientation. Journal of Mechanical Behavior of Biomedical Materials 1 (2008), 326–335. [66] Stylianopoulos, T., and Barocas, V. H. Volume-averaging theory for the study of the mechanics of collagen networks. Comput Meth Appl Mech Eng 196, 31-32 (2007), 2981–2990. [67] Thomas, C. H., Collier, J. H., Sfeir, C. S., and Healy, K. E. Engineering gene expression and protein synthesis by modulation of nuclear shape. Proc Natl Acad Sci U S A 99, 4 (Feb. 2002), 1972–7. [68] Vacanti, C., and Mikos, A. Letter from the Editors. Tissue Engineering 1, 1 (1995), 1–2. [69] Vesely, I., and Noseworthy, R. Micromechanics of the fibrosa and the ventricularis in aortic valve leaflets. Journal of Biomechanics 25, 1 (1992), 101–113. ACS Paragon Plus Environment

36

ACS Biomaterials Science & Engineering

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 ACS Paragon Plus Environment

Page 52 of 52