_{2.5}in 2011, with a societal cost of $886 billion, highlighting the importance of modeling emissions at fine spatial scales to prioritize emissions mitigation efforts.

## New Research In

### Physical Sciences

### Social Sciences

#### Featured Portals

#### Articles by Topic

### Biological Sciences

#### Featured Portals

#### Articles by Topic

- Agricultural Sciences
- Anthropology
- Applied Biological Sciences
- Biochemistry
- Biophysics and Computational Biology
- Cell Biology
- Developmental Biology
- Ecology
- Environmental Sciences
- Evolution
- Genetics
- Immunology and Inflammation
- Medical Sciences
- Microbiology
- Neuroscience
- Pharmacology
- Physiology
- Plant Biology
- Population Biology
- Psychological and Cognitive Sciences
- Sustainability Science
- Systems Biology

# Assembly of complex cell microenvironments using geometrically docked hydrogel shapes

Edited by Robert Langer, Massachusetts Institute of Technology, Cambridge, MA, and approved February 5, 2013 (received for review January 10, 2013)

## Abstract

Cellular communities in living tissues act in concert to establish intricate microenvironments, with complexity difficult to recapitulate in vitro. We report a method for docking numerous cellularized hydrogel shapes (100–1,000 µm in size) into hydrogel templates to construct 3D cellular microenvironments. Each shape can be uniquely designed to contain customizable concentrations of cells and molecular species, and can be placed into any spatial configuration, providing extensive compositional and geometric tunability of shape-coded patterns using a highly biocompatible hydrogel material. Using precisely arranged hydrogel shapes, we investigated migratory patterns of human mesenchymal stem cells and endothelial cells. We then developed a finite element gradient model predicting chemotactic directions of cell migration in micropatterned cocultures that were validated by tracking ∼2,500 individual cell trajectories. This simple yet robust hydrogel platform provides a comprehensive approach to the assembly of 3D cell environments.

Biological tissues are composed of cellular “building blocks” that cooperate to provide tissue-specific functions (1⇓⇓⇓⇓–6). Specific cells, molecules, and their geometric assembly establish a biological system, whether it is a vascular network surrounded by parenchymal cells (7), a developing tissue (8⇓⇓–11), or a metastatic tumor (12⇓⇓⇓–16). In vitro culture systems designed to control the 3D presentation of multiple cells and molecular species in a biologically relevant matrix are needed to faithfully recapitulate intricate biological niches (17⇓⇓⇓⇓⇓–23). Previous attempts to assemble complex tissue structures in vitro lacked the specificity and yield to control large numbers of cells and soluble factors (22⇓–24), had limited resolution at the microscale (25, 26), and used matrices without adequate biological function (27⇓⇓⇓⇓–32). Therefore, spatial microenvironmental control of biological systems has been difficult to achieve. Such control is important in many processes, including the migratory formation of vasculature, where gradient patterns dictate growth of vascular sprouts. For example, endothelial cell (EC) stabilization via mesenchymal stem cell (MSC) interactions is known to facilitate the maturation of blood vessels impacting many physiologic systems, from tumors to engineered tissues (33⇓–35). Recent studies clearly showed the importance of microscale gradients for vasculogenesis (36⇓–38).

We report a method by which microsized 3D hydrogels are shape-coded for their biological and physical properties and docked by iterative sedimentation into shape-matching hydrogel templates. Microenvironmental niches were fabricated using gelatin methacrylate (GelMA), a modified native protein with excellent biocompatibility, tunable mechanical properties, and micrometer-scale patterning resolution (39, 40). GelMA was molded into diverse geometric shapes with dimensions of 100–1,000 µm, after encapsulating cells or labeled molecular species of a desired type and concentration. This way, biological niches were shape-coded, and various shapes could be assembled into collective structures with emergent properties.

Shape-coded cell niches were docked into hydrogel templates containing congruent wells of different geometries using gravity-based sedimentation. Over multiple sedimentation events, docked shapes were shear protected and remained in their wells, whereas mismatched shapes were resuspended to settle again. The geometric uniqueness of each shape enabled exclusivity and specificity of docking into congruent wells, allowing simple yet precise fabrication of geometric patterns mimicking native cell microenvironments. These techniques were used to form radially organized 3D patterns to study sprouting and the homing of MSCs to ECs, by using the diffusive signaling patterns to dictate spatial orientation of cell migration. A finite element gradient model was built by accounting for the repelling and attracting forces originating from the docked hydrogels. Large numbers of single cells were tracked using advanced image processing algorithms to validate the dependence of cell migration on the spatial gradient.

## Results

### Fabrication of Hydrogel Shapes and Templates.

Cells and labeled molecular species were encapsulated in GelMA, a naturally occurring peptide-based hydrogel, which was fabricated using elastomeric poly(dimethylsiloxane) (PDMS) molds into a range of nonsuperimposable geometric shapes (Fig. 1 *A–C* and Fig. S1). Each hydrogel shape was coded for the specific types and concentrations of encapsulated cells and cytokine surrogates, and was able to preserve cell viability and biological function over weeks of in vitro culture. Templates containing combinations of negative reliefs of shape geometries were designed to mimic different potential microenvironments (Fig. 1 *D* and *E*) and fabricated using methacrylated glass, which covalently cross-linked with the GelMA to form a substrate-attached hydrogel template (Fig. 1 *D* and *E*).

### Assembly of Hydrogel Shapes.

To assemble microenvironments, shapes in suspension were pipetted onto a hydrogel template and docked into their respective wells by sequential sedimentation–resuspension on an orbital shaker. Shapes, used in excess of their corresponding wells, were first agitated for 1 min and then allowed to sediment for another minute. At the end of each cycle, a fraction of shapes docked into their congruent wells. The undocked shapes, which settled outside of wells, being geometrically excluded or incorrectly oriented, were resuspended and allowed to sediment again, whereas the correctly docked shapes stayed locked in their wells due to shear protection (Fig. 1*F*). This process allowed for shapes to be sorted into their congruent wells with high fidelity (Fig. 1 *G* and *H*).

### Specificity and Yield of Docking.

To maximize specificity of docking, different sets of canonical shapes were designed. Shapes sized at 60% of the surface area of their congruent well sorted with poor specificity, such that only 55% docked into their correct wells (Fig. 1*I*). When the shapes were sized to 85% of the surface area of their congruent wells, the specificity increased to 90% (Fig. 1*J*). The specificity remained unchanged over sequential rounds of shaking, confirming geometric exclusion as the primary determinant of specificity. Shapes correctly sorted with a yield of ∼20% in the first cycle, and an additional 10% in each subsequent cycle, to a maximum yield of 85%, a limit imposed by the maximum specificity. The yields and specificities of iterative sedimentation surpass those achieved with any previous technique. The change in hydrogel composition (percentage GelMA) had significant effect on mechanical properties of the hydrogel without affecting the sorting specificity or yield (Fig. S2). After sorting, another layer of GelMA was applied to permanently lock the sorted shapes into a cell microenvironment composed of a homogenous biocompatible hydrogel matrix.

### Molecular Diffusion from Hydrogel Shapes.

The ability to encapsulate molecular species within GelMA, along with the precision of the sorting technique, provided a method to generate unique spatial gradients. Fluorescently tagged dextran or BSA, serving as surrogates for cytokines of large molecular size, were encapsulated into circular and rectangular shapes, docked to form specific patterns, and were allowed to outwardly diffuse. Over the course of 100 min, red fluorescent dextran emanated from the circle and green fluorescent dextran emanated from the rectangle, providing a spatial gradient field over the millimeter-size scale (Fig. 2 *A–D* and Movies S1 and S2). The evolution of concentration gradients demonstrated that the concentrations of molecular species changed with time in a predictable manner (Fig. 2*E*). Hydrogel composition modulated the diffusivity, with faster transport rates of encapsulated BSA measured for hydrogels made with a lower percentage of GelMA (Fig. S3).

### Engineering Cell Behavior Using Hydrogel Shapes.

We demonstrated that hydrogel shapes can be used to form biomimetic cell environments to modulate the function of encapsulated cells. As an example, neonatal rat cardiomyocytes were encapsulated within rectangular hydrogels of two width/length aspect ratios: 1:3 and 1:10. After 7 d of culture, all shapes exhibited spontaneous macroscopic beating, with high fractions of viable cells (>98%; Fig. 2 *F* and *G*) and remarkable differences in contractile cell function. The 1:3 rectangles contracted at a significantly higher frequency and to only ∼12% of their length, whereas the 1:10 rectangles contracted more slowly and to ∼36% of their total length. These data suggest differential phenotypic maturation of cardiac subunits by simple changes in hydrogel geometry (Fig. 2 *H* and *I*) and capability to mature the cellular subunits before assembly into docked patterns.

### Patterning MSC Migration.

The hydrogel sorting platform was then used to control the coculture of human MSCs and human ECs. MSCs were encapsulated in circular hydrogels, sorted into centrally positioned circular wells, and allowed to migrate out of the hydrogel toward the surrounding rectangles containing either acellular or EC-laden hydrogels (Fig. 2 *J–M* and Fig. S4). Images in both groups were taken at the same time points, with the only variable being the presence of ECs. MSCs by themselves migrated out in a constant symmetric, radial fashion; however, MSCs cultured with EC-containing hydrogels initially migrated radially, curved their trajectories over time, creating long, arcing sprouts homing toward ECs. Confocal microscopy of actin (Fig. 2*N*) provided a closer view of the cytoskeletal architecture of the migratory MSC sprouts. This visually striking, spatially dependent influence on MSC growth suggested a gradient-based communication between the two cell types.

### Finite Element Model of Cell Migration.

To explore gradient effects between the ECs and MSCs, a 3D finite-element model was built using the two hydrogel shape types as gradient sources (Fig. 3*A*). The model accounted for two driving forces: (*i*) repulsive homotypic interaction between MSCs dictating cell migration from a region of high to low cell density (MSC–MSC repulsion), and (*ii*) the chemotactic gradient of EC-secreted factors dictating MSC migration from a region of low to a region of high concentration (MSC–EC attraction). The MSC-only system was driven solely by the repulsive cell gradient, whereas the coculture system combined the two driving forces to form a gradient field akin to a dipole electric field. The gradient field model was used to check the relative independence of specific shape geometry on the subsequent formation of the gradient field. Different shapes—rectangles, circles and squares—all formed similar gradient fields when positioned equally, validating that the mutual exclusivity of shape geometry did not necessitate variable gradient formation (Fig. S5 *A–C*). Changing individual shape geometry had only a small effect compared with changing spatial position of wells in the hydrogel template, which caused grossly different gradient patterns (Fig. S5 *D–F*). When biologically assessed, cell migration in linear configurations matched the model of Fig. S5*E* accordingly (Fig. S6).

### Cell Migration Vectors.

To validate the accuracy of the gradient model of cell migration behavior, cell trajectories were quantified in the time-lapse images by generating cell-associated vectors using a MATLAB-based image segmentation algorithm. Contrast thresholding, image smoothing, Sobel edge detection function, and a watershed transformation were implemented to create segmented image areas for each cell (Fig. S7). Ellipses were inscribed within segmented areas, and the resultant long axis vector was used to represent the cell (Fig. 3 *B* and *C*). The vectors formed a comprehensive dataset (corresponding to ∼2,500 cells), to allow quantitative comparisons to be made for cell migratory speed, angular distribution, and correlation of cell migration to the developed gradient model. Using the vector dataset, an average rate of MSC migration was calculated in the presence or absence of ECs. The MSCs alone migrated with an average speed of 25.9 µm/d, but when in the presence of ECs, they migrated 36.1 µm/d (a 40% increase; *P* = 4.2 × 10^{−5}). The angular distributions of cells were examined using histograms of the cell angles relative to the right horizontal axis (Fig. 3 *F* and *G*). MSCs cultured alone exhibited uniform distributions, whereas MSCs cocultured with ECs clustered at angles correlating to the locations of ECs.

### Correlation of Cell Migration to Finite Element Model.

The measured directions of cell migration were compared with model predictions by superimposing cell vector data onto the model vector fields. Overall, the migrating cells followed the gradient field lines (Fig. 3 *D* and *E*). Along radial diagonal lines between rectangles in the coculture model, the migratory deviations reflected ambiguous cell responses to competing gradients (equal attraction to both rectangles) (Fig. 3 *D* and *E*).

For quantitative comparison, the angles between cell data vectors and model vectors were measured. The angle between the cell data vector and the finite-element model vector ranged from 0° to 90°, 0 corresponding to perfect agreement between the data and the model (Fig. 3 *D* and *E*, *Insets*). Comparison of the MSC-only data with the single-gradient model resulted in an average angle deviation of 21.1°, substantiating that ∼20° deviation is intrinsic to MSC migration. Interestingly, the average angle was 22° and 29° when comparing the coculture cell data to the coculture model and radial model, respectively. A statistically significant decrease (by ∼7°, *P* = 2.23 × 10^{−7}) of the average angle validates that two simple driving forces can establish a gradient model accounting for the interactions between MSCs and ECs.

To spatially resolve where the MSCs deviated from the radial model and migrated toward the ECs, heat maps of cell direction identified the locations where the cells followed the coculture model over the radial model (Fig. 3*I*). The red areas of the heat map indicate the specific regions where the MSCs responded to the paracrine signaling effects of ECs by changing course from a radial migratory pattern, which occurred close to the EC-containing rectangles. By comparison, the heat map of the angle magnitudes (Fig. 3*H*) for the MSC-only data compared with the single-gradient model showed no spatial localization.

## Discussion

This simple assembly technique solves some of the problems of fidelity and stepwise construction of cellular microenvironments, by virtue of precision, versatility, and scalability. Precision was achieved using mutual geometric exclusion, where shape-coded hydrogels could be uniquely tuned to ensure high specificity of spatial docking within the hydrogel template (Fig. 1). Versatility was achieved using combinatorial patterning of multiple shape geometries with customizable amounts of encapsulated cells and fluorescently labeled surrogates of cytokines. Additionally, the hydrogel material, GelMA, is highly biocompatible, supports physiologic cell functions, and has tunable stiffness and diffusivity. Efficiency and scalability were achieved through the use of an iterative docking–resuspension protocol (Fig. 1). Shear protection allowed for rapid conjugate-well filling due to accumulation of successfully docked shapes in the final assembly. Because each docking event was independent, organization occurred in parallel, solving scalability problems, as the building blocks were placed not sequentially, but simultaneously. Thus, both microsized cell niches and large tissue constructs can now be assembled over short time periods.

Cell migration was particularly apt to study diffusive cell–cell communication during vasculogenic network formation, due to the visible spatial organization of vascular cells along the concentration gradient fields (Figs. 2 and 3). Despite potentially complicated signaling networks between MSCs and ECs, the pattern formation was predictable using only two driving forces—MSC–MSC repulsion and MSC–EC attraction. The mechanism of outward cell migration is not fully understood (41, 42) and has been described as a general diffusive process (43, 44). Using the microscale hydrogel system and an image-processing algorithm, we generated a large and reproducible cell vector dataset, to quantitate the rates of migration that increased by ∼40% due to the chemotactic effects of ECs that attracted MSCs. As literature suggests that the PDGF-BB secreted by ECs is a key recruitment signal for MSCs (33, 45⇓⇓⇓–49), our future experiments will investigate the role of PDGF-BB in EC–MSC interactions, among multiple factors that may be at play in microscale hydrogel environments we describe here.

We present a universal fabrication platform for precisely generating complex 3D cell microenvironments over multiple size scales, by patterning any number of encapsulated cells and molecular species, in any configuration, simultaneously. The characterization of paracrine cell–cell interactions was done using advanced image processing, finite element modeling, and complex 3D cell patterning. Two model growth patterns for MSC migration provided the basis for mathematical models defining spatial maps of diffusive cell communication. Due to the pervasive effect of the cell microenvironment on cell function, we envision that this method of subunit-based microenvironmental assembly will be used to explore numerous processes, such as tumor metastasis, stem cell differentiation, and scaling for large, spatially designed tissue engineered constructs.

## Materials and Methods

### Cell Culture.

Human MSCs were isolated from fresh bone marrow aspirates obtained from Cambrex Life Sciences as previously described (38, 50), grown in high-glucose DMEM with 10% (vol/vol) FBS and 1% penicillin/streptomycin (all from Invitrogen) with 0.1 ng/mL basic FGF (Invitrogen), and expanded to passage 8 or lower. Human umbilical vein endothelial cells were isolated from postnatal discarded umbilical veins and fully de-identified following the protocol IRB-AAAC4839 approved by the Columbia University Institutional Review Board, and grown using EGM-2 Media (Lonza) to passage 8 or lower. Rat neonatal cardiomyocytes were harvested from day 1–2 neonatal rats as in our previous studies (51) under protocol AC-AAAB2757 approved by the Columbia University Institutional Animal Care and Use Committee.

### Micromolding of Hydrogel Shapes and Templates.

GelMA was synthesized according to a previously published method (39, 40). PDMS (Fisher) molds for generating the hydrogel shapes and wells were prepared using standard soft lithography. The masks used in SU-8 (Microchem) photolithography to fabricate master molds were generated by (Cad/Art Services). Specifically, molds for circle shapes were 600 µm in diameter × 500 µm deep, rectangles were 300 µm wide × 750 µm long × 300 µm deep, and triangles were equilateral with each side being 600 µm × 300 µm deep. Shapes used in the cardiomyocyte study were 100 µm wide × 300 µm or 1,000 µm long × 100 µm deep. Each shape was generated in ∼400 repeating copies per batch, using a 4 cm × 2.5 cm mold. Molds for the well templates were fabricated using the same soft lithography technique. PDMS molds were created from the master patterns, washed, dried, and plasma cleaned for 60 s.

Hydrogel shapes were fabricated using a 10% (wt/vol) prepolymer solution of GelMA, 1% (wt/vol) 2-hydroxy-1-[4-(2-hydroxyethoxy)phenyl]-2-methyl-1-propanone (Irgacure-2959; BASF), and 1% (vol/vol) 1-vinyl-2-pyrrolidone (Sigma) in PBS (Invitrogen). For each geometric shape, a defined concentration of cells or fluorescently labeled molecules was added to each prepolymerization GelMA solution. The GelMA solution containing a desired concentration of specific cells or molecular species was pipetted onto the PDMS mold, placed on a glass slide excluding excess GelMA, and exposed to UV light for 30 s (Omnicure 100; Lumen Dynamics). Following polymerization, the shapes were harvested from the wells by resuspension and pipetting. To generate conjugate well templates, masks with a specific configuration of multiple geometric shapes were used to generate the PDMS molds. The 10% (wt/vol) prepolymer GelMA solution was placed between the PDMS mold and a glass slide treated with 3-(trimethoxysilyl)propyl methacrylate (Sigma) to permit a covalent linkage of the well hydrogel to the glass substrate upon polymerization. For additional specifications, please see Fig. 1 and Fig. S1.

### Sorting and Docking.

After fabrication, a mixture of geometric shapes needed for a specific experiment was suspended in DMEM at a concentration of ∼100 shapes per mL. A suspension of the three canonical shapes—circle, triangle and rectangle—was pipetted on the hydrogel template containing a specific configuration of geometrically different wells, and allowed to settle for 1 min. The hydrogel template was then exposed to agitation using an orbital shaker (Innova 4000; New Brunswick Scientific) set to 90 rpm, which allowed sufficient agitation to resuspend incorrectly sorted shapes, but was gentle enough to maintain shear protection for the correctly sorted shapes. After 1 min of agitation, the shapes were allowed to settle again for 1 min. This process was repeated for up to 10 cycles, by which time the well templates were occupied by properly docked shapes. Following sorting, excess shapes were washed away from the wells, and an additional layer of un–cross-linked GelMA was UV polymerized for 30 s to permanently lock the sorted shapes.

### Shape Sorting Analysis.

A suspension containing excess amount of hydrogel shapes of any specific type was pipetted onto a template in form of a Cartesian grid of circles, triangles, and rectangles at a density of 92 wells per cm^{2}. For each shape and gel concentration, the docking specificity was calculated as a ratio of correctly sorted shapes and the total number of docked shapes. The docking yield was calculated as a ratio of the number of the docked shapes of a particular type and the number of congruent wells.

### Diffusion Studies.

Circular shapes were suspended in a 100 µg/mL Texas Red-tagged 10-kDa dextran or Texas Red-tagged BSA. Rectangular shapes were suspended in a 100 µg/mL Alexa 488-tagged 10-kDa dextran or Alexa 488-tagged BSA for 8 h to allow the fluorescently tagged molecules to fully diffuse into the gel (all reagents, Invitrogen). Shapes were sorted and assembled into the template wells, and fluorescence images of the patterns were taken every 10 min over 100 min at the same image settings. Using MATLAB, areas with the shapes were masked out, and fluorescence intensity of the remaining images was normalized to the maximum intensity of all images in the experiment. As an additional metric comparing 10% and 20% GelMA, the distance at which fluorescence intensity dropped below 50% was measured for each time point.

### Cell Migration.

MSCs were encapsulated into circular shapes; ECs were encapsulated into rectangular shapes at a concentration of 10^{6} cells/mL. For the migration experiments, the hydrogel template contained patterns with a central circle (600 µm diameter × 750 µm deep) surrounded by four rectangles (300 µm wide × 750 µm long × 750 µm deep) that were spaced 750 µm from the center of the circle. Cells were cultured in EGM-2, and phase contrast images were taken at regular intervals for up to 14 d.

### Numerical Simulation of Diffusion Gradients.

The multiphysics software COMSOL (COMSOL AB) was used to solve the partial differential equation (PDE) for diffusion of bioactive molecular species (encapsulated in hydrogel or released by the cells) over the entire domain [, where *D* is the diffusion coefficient, and *c* is the concentration of molecular species]. We chose the whole template to serve as a simulation domain, because we were interested in cellular migratory behavior in the template and not within the shapes themselves. We assume that the molecular diffusion within the shapes did not affect the gradient within the template. The circumferential boundary conditions of different shapes—positive flux for the growth factor secreting cells (ECs as a sink) and negative flux for the growth factor-consuming cells (MSCs as a source)—were defined using the following equation:, where is the normal vector and is the constant inward flux at the boundary. We further set the initial growth factor concentration at zero throughout the entire domain (), defined the domain boundaries with convective flux (which denotes that the growth factor can freely diffuse in and out of the domain), and held the diffusion coefficient constant over the domain (). To achieve short simulation times with high accuracy, we chose a mesh size of 2 µm in the vicinity of the boundaries and 10 µm in the bulk. We used “Direct (UMFPACK)” solver to solve the PDE with the above boundary conditions and exported the concentration gradient vector map to compare with the experimental data.

### Image Analysis of Cell Migration.

Cell migration images were processed by MATLAB to inscribe vectors quantifying their orientation. Images were masked to remove the area containing shapes, because the shapes gave obvious noncell vectors due to edge effects. The subsequent image was smoothed using a reconstruction algorithm, with a close-open filtering sequence. A serial opening by reconstruction using disk-shaped structuring elements of 2–7 pixels in diameter allowed for discolorations within the cell to be removed without affecting the cell boundaries, effectively smoothing the cell image while preserving clear edges of the cell boundaries. A Gaussian filter further reduced noise in the images by using a rotationally symmetric Gaussian blur for each pixel. Next, the images were thresholded using a high-pass filter to remove low-intensity pixels from the background. A Sobel operator was used to accentuate cell edges, with the gradient of pixel intensity used to highlight the regions with large differences in pixel intensity, corresponding to a cell edge. A watershed transformation allowed the segmentation of each cell into an enclosed shape for further analysis (Fig. S7). Using built-in ImageJ functions, an ellipse was fit to each segmented cell, where the centroid of the ellipse demarcated the location of the cell and the long axis was used as the cell vector.

### Determination of Migration Rate.

Images were quantified for overall migration rates using vectors generated from day 14 images. The total distance from the circle edge to the leading migration edge was measured using ImageJ by connecting vectors head to tail, resulting in a total distance traveled. This analysis was performed for all points around the circle and averaged over time, resulting in an average migration rate for the MSC alone and MSC–EC coculture groups.

### Cell Vector Analysis.

To compare the measured angles of cell vectors with the model predictions, each cell location was correlated to the model by aligning the hydrogel rectangles on bright-field images to the rectangles in the model, and scaling the distances so all of the cell images would be spatially compatible. Each cell vector was compared with the associated model vector and the dot product was used to determine the angle between the two vectors, to be used to generate heat maps and histograms.

### Heat Maps.

Heat maps were generated using a MATLAB script with the grid-fit function to derive the plots from the cell vector angle data. For monoculture heat maps (Fig. 3*G*), the angle between each cell vector and its associated radial vector was determined from the dot vector product, and plotted as a heat map (with red corresponding to 90°, and blue corresponding to 0°). For coculture heat maps (Fig. 3*I*), the angle between the measured and modeled coculture cell vectors and the coculture model was compared with the radial diffusion model. The heat maps highlight in red the regions where the coculture cell vectors showed agreement with the coculture diffusion model, but not the radial model.

### Statistics.

To compare the effect of shape geometry or hydrogel composition on the yield of shape sorting, two-way ANOVA was performed (*n* = 3, *P* < 0.05). Furthermore, to compare efficiency, two-way ANOVA was performed at each round of shaking (*n* = 3, *P* < 0.05). To establish differences between the experimental groups for cardiac contraction studies, the nonparametric, Wilcoxon rank sum test (also known as Mann–Whitney *U* test) was used (*n* = 27–50, α = 0.05, *P* < 0.01). Statistical differences between migration rates for MSC-alone and MSC–EC groups were determined using a two-sample *t* test (*n* > 26, *P* < 0.05). A circular statistics toolbox designed for MATLAB was used to compare the average angle differences between the cell data and models (*n* = 2,500) (52). The Watson–Williams two-sample test (circular analog to the Student *t* test) was used to compare the groups, with the angle binning of <10°. This test assumes an appropriate von Mises distribution or circular analog of the normal distribution (52). With α = 0.05, a final *P* value was <0.01. Normality testing using the Kuiper test (circular invariant analog to Kolmogorov–Smirnov test) demonstrated normal distributions of compared with the von Mises distribution.

For additional details on GelMA synthesis, mechanical testing, and cardiac contraction analysis, please see *SI Materials and Methods*.

## Acknowledgments

This work has been supported by National Institutes of Health Grants HL076485 and EB002520 (to G.V.-N.) and HL112505 (to G.E.).

## Footnotes

- ↵
^{1}To whom correspondence should be addressed. E-mail: gv2131{at}columbia.edu.

Author contributions: G.E. and G.V.-N. designed research; G.E., B.W.L., J.S., and G.L. performed research; G.E., B.W.L., H.P., C.D.C., S.K.S., and G.V.-N. analyzed data; and G.E., B.W.L., H.P., S.K.S., and G.V.-N. wrote the paper.

The authors declare no conflict of interest.

This article is a PNAS Direct Submission.

This article contains supporting information online at www.pnas.org/lookup/suppl/doi:10.1073/pnas.1300569110/-/DCSupplemental.

## References

- ↵
- Turing AM

- ↵
- Kondo S,
- Miura T

- ↵
- ↵
- ↵
- Nelson CM,
- et al.

- ↵
- ↵
- Madden LR,
- et al.

- ↵
- Marolt D,
- et al.

- ↵
- Grayson WL,
- et al.

- ↵
- Radisic M,
- et al.

- ↵
- ↵
- ↵
- ↵
- ↵
- Chen CS,
- Mrksich M,
- Huang S,
- Whitesides GM,
- Ingber DE

- ↵
- ↵
- ↵
- Khademhosseini A,
- Langer R,
- Borenstein J,
- Vacanti JP

- ↵
- ↵
- Derda R,
- et al.

- ↵
- Alsberg E,
- Anderson KW,
- Albeiruti A,
- Rowley JA,
- Mooney DJ

- ↵
- ↵
- ↵
- Du Y,
- Lo E,
- Ali S,
- Khademhosseini A

- ↵
- ↵
- Jakab K,
- Neagu A,
- Mironov V,
- Markwald RR,
- Forgacs G

- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- Kim D-H,
- et al.

- ↵
- ↵
- Armulik A,
- Abramsson A,
- Betsholtz C

- ↵
- ↵
- ↵
- Zervantonakis IK,
- et al.

- ↵
- ↵
- ↵
- ↵
- ↵
- Ridley AJ,
- et al.

- ↵
- ↵
- ↵
- Hirschi KK,
- Rohovsky SA,
- Beck LH,
- Smith SR,
- D’Amore PA

- ↵
- ↵
- ↵
- Cool J,
- DeFalco TJ,
- Capel B

- ↵
- ↵
- ↵
- ↵
- Berens P

## Citation Manager Formats

## Article Classifications

- Biological Sciences
- Applied Biological Sciences

- Physical Sciences
- Engineering