Abstract
Mechanical metamaterials constitute a nascent category of architected structures comprising arranged periodic components with tailored geometrical features. These materials are now being employed as advanced medical implants due to their extraordinary mechanical properties over traditional devices. Nevertheless, to achieve desired tissue integration and regeneration, it is critical to study how the microarchitecture affects interactions between metamaterial scaffolds and living biological tissues. Based on human induced pluripotent stem cell technology and multiphoton lithography, we report the establishment of an in vitro microtissue model to study the integration and remodeling of human mesenchymal tissues on metamaterial scaffolds with different unit geometries. Microtissues showed distinct tissue morphologies and cellular behaviors between architected octet-truss and bowtie structures. Under the active force generated from mesenchymal tissues, the octet-truss and bowtie metamaterial scaffolds demonstrated unique instability phenomena, significantly different from uniform loading using conventional mechanical testing.
Introduction
Mechanical metamaterials are architected structures possessing microscale or nanoscale topological features that are assembled in a specialized order. The material properties of metamaterials depend on the geometry of the unit structure. The rational design of the unit structure and arrangement enables superior material properties at macroscales over natural materials.1,2 For instance, pentamode metamaterials, which are made of face-centered cubic unit cell structures, have high bulk modulus, but very low shear modulus. 3 This enables them to function as liquid structures and inherit controllable wave propagation. Moreover, novel metamaterial designs encompassing hyperelastic behavior have been proposed to employ brittle materials for fabricating devices that require large, but recoverable, deformations.
For applications in biomedical engineering, mechanical metamaterials have been mostly used in design and fabrication of metaimplants.4,5 By combining metamaterials with different degrees of auxeticity in a hip implant design, implant–bone integration and implant longevity were substantially improved. 6 To improve host tissue integration by manipulating the anisotropy of the implant (i.e., the Zener ratio), 7 a combination of auxetic and isotropic structures makes it possible to fabricate the metaimplant possessing the failure resistance and adjustable anisotropy that could match that of specific tissues.
Despite these incipient results of metaimplants, the use of metamaterials is still focused on optimization of implant adaptivity to native conditions. Tissue–scaffold interaction between mechanical metamaterials and biological tissues is still underexplored, while unraveling this response is critical for the purpose of tissue repair and regeneration.
Biomechanical studies are interlaced with the development of in vitro models using architected scaffolds.8–10 With the emergence of multiphoton lithography (MPL), manufacturing complex three-dimensional (3D) scaffolds enabled the design of cell niches to control and regulate cell behaviors and functions. 11 The ultrahigh printing resolution of MPL makes it possible to recapitulate the nanoscopic topography of the extracellular matrix, which is necessary to create physiologically relevant models. For example, using MPL, the architecture of 3D lattice microscaffolds was found to trigger β-catenin activity of breast cancer cells through mechanotransduction pathways, which further promoted their proliferation and invasiveness. 12
In another study, tetrakaidecahedral nanolattices with better compliance were shown to enhance the expression of f-actin and calcium secretion of osteoblast-like cells. 13 Two-dimensional auxetic structures were also found to affect the spreading of mesenchymal stem cells and division of fibroblasts.14,15 Although many studies demonstrated the effect of lattice architecture on cellular behaviors and phenotypes, how the tissue remodeling process would affect the scaffold architecture through biomechanical interactions has not been elucidated yet.
To quantitatively understand the relationship between tissue remodeling and tissue mechanical conditions, the biomechanical models based on finite element analysis (FEA) elucidated the stress distribution within 3D microtissues and predicted their remodeling behaviors.23–26 The accuracy of FEA is dependent on the comprehensive understanding of mechanical properties for both living tissues and biomaterial scaffolds. For example, formation of sarcomeres has been linked with the local stress state of cardiac microtissues in an FEA model to computationally simulate maturation of stem cell-derived cardiac microtissues. Furthermore, this computational simulation was used to guide the disease modeling of arrhythmogenic cardiomyopathy by modulating the tissue geometry. 27 More recently, real-time stress distribution within mouse femurs was traced by a micro-FEA model and then utilized to adapt the mechanical loading and lower the variance across all mice. This helped to prevent overloading of individual femurs, which would create bone fractures. 28
In this work, we developed an in vitro tissue model based on different designs of mechanical metamaterials to study how the metamaterial architecture would affect tissue formation. Two types of thoroughly investigated metamaterials (octet-truss and auxetic bowtie structures) were fabricated by MPL using SZ2080TM, which is an inorganic–organic hybrid photoresist used for high-resolution laser fabrication. Compared with other commonly used photoresist materials, it shows negligible shrinkage during solvent washing and supports the growth of various cell types.29–31
Microtissues were generated by growing mesenchymal stromal cells (MSCs) derived from human induced pluripotent stem cells (hiPSCs) on metamaterial scaffolds. MSCs are widely used for general mechanobiological studies due to their presence in many tissue types and differentiation capability in response to different mechanical environments.16–19 However, primary MSCs obtained from human donors have limited proliferation potency and high variability in cell quality. 20 To overcome these restrictions, hiPSCs have become a promising resource for MSCs because they can be expanded to a large quantity while maintaining the differentiation potential. There are already many protocols available to robustly generate MSCs from hiPSCs using small molecules and defined media.21,22 Consistent production of MSCs from the same hiPSC source ensures model development with high robustness, consistency, and reproducibility. We found that both overall tissue morphology and local cell behaviors were highly dependent on the microarchitecture of metamaterials. Furthermore, the mechanical force generated from microtissues induced unique deformation patterns on metamaterial scaffolds, leading to structural instability due to buckling. This in vitro model provides the avenue to obtain a fundamental understanding of biomechanical interactions between living microtissues and metamaterial scaffolds, which will potentially be utilized to formulate new design principles to guide generation of artificial tissues for various applications.
Materials and Methods
MPL fabrication of metamaterial scaffolds
All of the metamaterial scaffolds were fabricated by the MPL process with an organic–inorganic hybrid resin, SZ2080TM, as the basal material. 32 The experimental apparatus and fabrication conditions have been reported previously. 33 Specifically, the MPL apparatus is equipped with an near infrared laser (FemtoFiber pro, Toptica) with 780 nm wavelength, 100 fs pulse width, and 80 MHz repetition rate. The laser beam was focused using a 100 × microscope objective lens (Plan Apochromat 100 × /1.40 Oil M27; Zeiss). The stage was translated such that the laser beam could polymerize inside the material and fabricate the geometry. Each structure was designed in a CAD file using SolidWorks 2019x64. Through an STL file, it was converted into G-code so that the stages can be translated accordingly.
hiPSC-MSC differentiation
The hiPSC line was obtained from Dr. Conklin's laboratory at the University of California, San Francisco (UCSF). The committee on Human Research at UCSF approved the hiPSC research protocol (#10-02521). hiPSCs were maintained on 6-well plates coated with growth factor-reduced Geltrex (Ca# A1413302; Life Technologies) in Essential 8 (E8) medium (Ca# A1517001; Life Technologies). The protocol to differentiate hiPSC-MSCs has been published previously. 23 Briefly, hiPSCs were treated with Essential 6 (E6) medium (A1516401; Life Technologies) supplemented with 10 ng/mL basic fibroblast growth factor (Ca# 233-FB; R&D Systems), 4 μM SB431542 (Ca# 04-0010-10; Stemgent), and 4 μM CHIR99021 (Ca# 04-2004; Stemgent) for 5 days. Next, differentiated cells were dissociated and replated in a serum-free MSC culture medium (CTS StemPro MSC SFM) (A1033201; Life Technologies) for four more passages to obtain differentiated hiPSC-MSCs.
Generation of mesenchymal microtissues
To generate mesenchymal microtissues, one mechanical metamaterial scaffold was placed into one well of a 6-well plate and sterilized with 70% ethanol for 1 h. After rinsing three times with Dulbecco's phosphate-buffered saline (PBS) (Gibco), the scaffolds were then coated with diluted Matrigel for 1 h. To seed hiPSC-MSCs onto the metamaterial scaffolds, the cells were dissociated, suspended, concentrated with a density of 1.5 × 107 cells/mL, and seeded with a total volume of 20 μL. After cell seeding, scaffolds were incubated at 37°C for 1.5 h to promote initial cell attachment before more cell culture medium was added.
Immunostaining and confocal microscopy
Architected microtissues were fixed with 4% (vol/vol) paraformaldehyde, permeabilized with 0.2% (vol/vol) Triton X-100, and blocked with 2% (wt/vol) bovine serum albumin. Next, microtissues were incubated with phalloidin for 1 h to stain F-actin and then incubated with DAPI for 10 min to stain nuclei. Confocal microscopy (U880; Zeiss) was used to capture z-stacks of the microtissues architected by the metamaterial scaffolds, which can be observed simultaneously due to autofluorescence.
Scanning electron microscopy
The microtissues were fixed in 4% paraformaldehyde overnight with 0.1 M PBS at room temperature. After washing three times with PBS, samples were then dehydrated in a series of concentrations of ethanol (15%, 30%, 50%, 70%, 95%, and 100%) at room temperature for 15 min at each concentration, plus two more 100% ethanol dehydration steps at the end. After dehydration, samples were dried in a vacuum oven for one day. Prepared samples were placed on stubs and sputter-coated with a thin layer of gold with a thickness of 10 nm. Finally, a scanning electron microscope (JSM-IT100LA; JEOL USA, Inc.) was used to image the microtissues architected by metamaterial scaffolds.
Helium ion microscopy
After sample preparation for SEM, the samples were imaged using helium ion microscopy (HIM; Zeiss ORION NanoFab). To avoid charging effects, imaging was conducted using the electron flood gun to negate the accumulating positive charge from a scanning ion beam, permitting imaging of samples without additional sputtering and allowing extremely high-resolution imaging required to observe the nanofeatures of cell membranes.
Finite element analysis
To investigate the effect of cell attachment to the scaffold, FEA simulations were performed using the multiphysics software, ANSYS R18.1. Both bowtie and octet-truss structures were discretized by 3D, 10-node, tetrahedral solid elements. Specifically, the octet-truss structure was discretized with 173018 nodes and 84558 elements, while the bowtie structure was discretized with 32301 nodes and 16648 elements. Mechanical properties were set using the same values that have been reported in the previous work.
31
Based on SEM images, cells within the bowtie structure were distributed along the edges of beam members, while cells within the octet-truss structure occupied the sides of the lattice. For this reason, a quasi-static displacement field was applied at these specific members of both lattices for large nonlinear deformations. For the bowtie structure, distributed displacement was set at the nodes of beam members at the unit edges to apply compression, while the beam members at the top were subjected to distributed displacement that would cause them to bend inward, which was represented by HIM images. For the octet-truss structure, the side beam members were subjected to distributed displacement that would cause compression of the array. Since the structures were fixed on a glass substrate, the bottom beam members were defined as fixed supports. After calculation of nonzero components of the displacement vector at each node, the force vector was obtained by the following equation:
where
Results and Discussion
We employed two types of unit structure designs, namely octet-truss and auxetic bowtie structures, to create metamaterial scaffolds. The octet-truss structure, furnishing ultrastiff and ultralight material properties, has been used as the building block of load-bearing biomaterials for hard tissue implants.34,35 The bowtie structure, on the other hand, is heralded by macroscopically negative Poisson's ratios, fracture resistance, and high-energy absorption performance36,37 (Fig. 1A). Previously reported dimensions of structural units were used in our designs (octet-truss structure: beam thickness: 1.5 μm, and base length: 50 μm; and bowtie structure: beam thickness: 8 μm, and base length: 50 μm),31,38 which ensure that the scaffolds possess the minimum volume possible while simultaneously having the structural integrity to prevent collapse. The bowtie structures had beam members with a larger diameter to prevent collapse during the fabrication process. The 3D metamaterial scaffolds were fabricated using a Ti: sapphire femtosecond laser scanning system (power: 1 mW, and scanning speed: 10 μm/s) to achieve submicron resolution in scaffold printing (Fig. 1B). 31 An organic–inorganic hybrid resin, SZ2080TM, was used as the basal material for fabricating all scaffolds. 32

Demonstration of octet-truss and auxetic metamaterials.
Fabricated scaffolds were imaged using scanning electron microscopy (SEM), showing precise microstructures and great mechanical integrity (Fig. 1C). The length scale of metamaterial structural units was similar to the size of individual biological cells, which would provide a direct association between biomechanical responses of cells and the structural geometry of the metamaterial unit. hiPSC-MSCs showed a spindle-like shape and expression of standard MSC markers, including CD73, CD90, and CD105 (Supplementary Fig. S1). hiPSC-MSCs were then seeded onto Geltrex-coated metamaterial scaffolds with a density of 1.5 × 107 cells/mL. The high cell seeding density resulted in formation of thick 3D tissue on metamaterial scaffolds. We estimated around four hiPSC-MSCs surrounding each structural unit of metamaterials.
On the octet-truss metamaterial scaffolds, we observed efficient penetration of hiPSC-MSCs throughout the 3D scaffolds based on the immunostaining of actin filaments and cell nuclei (Fig. 2A). Although many structural branches were present in the octet-truss structure, cells were able to traverse the struts, invade into pores within the scaffolds, and exhibit an elongated morphology. Compared with standard SEM imaging, HIM enhances the contrast, resolution, and depth of field, which enable visualization of tissue–material interactions at microscopic levels. 39 To maximize attachment to metamaterials with a low material density, hiPSC-MSCs grew and elongated following the longitudinal orientation of beams, which is a common phenomenon for cells on most architected scaffolds (Supplementary Fig. S2A).10,40,41 To further facilitate cell adhesion, thin branches of plasma membranes were extended from cell bodies to attach to adjacent beams (Supplementary Fig. S2B). In addition, micropores were present on cell bodies, which might have been created under high mechanical tension when hiPSC-MSCs spanned over different beams (Supplementary Fig. S2C).

Biomechanical interaction of mesenchymal microtissues on octet-truss metamaterials.
The mechanical load generated by the microtissues led to a large deformation of the entire octet-truss metamaterial scaffold, which manifested as a decrease in its height and the inward buckling of struts at the corners (Fig. 2A). The stability of center posts maintained an isotropic stress distribution across the octet-truss structure to keep the overall integrity of the metamaterial scaffold. To define the deformation modes of metamaterials, FEA simulations were performed using the multiphysics software, ANSYS R18.1. The mechanical properties were set with the same values reported in previous work. 33 We first simulated the mechanical behaviors of the octet-truss structure under uniform mechanical load. The octet-truss structure can sustain a high mechanical force without fracture with a linear relationship between force and displacement (Supplementary Fig. S4), which confirmed the ultrastiff property of this structure.
To recapitulate the deformation characteristics of octet-truss metamaterial scaffolds under the mechanical loading from microtissues, the beam members at the top part of the octet-truss structure were subjected to distributed displacement, which caused compression of the array (Fig. 2B). The bottom part attached on the coverslip, which limited the spatial moving capacity of beam members, was simulated as a fixed support. The force–displacement relationship showed that the octet-truss structure experienced a bifurcation buckling mode evidenced by the positive–negative–positive transition of the slope. 42 This distinct deformation mode indicated a unique way of mechanical loading from living biological tissues on metamaterial scaffolds.
The mesenchymal tissues growing on scaffolds formed tight connections with monolayer tissues grown on the surface of the coverslip, which applied a dominant compressive force transmitted through the main axis of beam members of the metamaterial scaffolds. Since the octet-truss structure is a stretch-dominated metastructure, the compression-induced displacement was limited by the interior hardness of materials. 43 However, the top part of the octet-truss structure had better compliance, showing that most of the beams buckled under the biological force from mesenchymal microtissues. The buckling of individual beams disrupted the integrity of the unit structure and then modulated the load-bearing mode of the octet-truss metamaterial scaffolds. Once the force was applied perpendicular to the beam axis, it was possible to deform the scaffold to a large extent with less mechanical loading.
On the auxetic bowtie metamaterial scaffolds, thicker tissue was formed with a bulky morphology and a clear tissue boundary following the contour of the scaffolds (Fig. 3A). The cells tended to aggregate around the bowtie structures with dense tissue formation instead of extending cell bodies with thin plasma membrane branches that were observed on the octet-truss metamaterial scaffolds. When imaging with HIM, we observed that cells wrapped around the beams and micropores were formed on the cell membrane (Supplementary Fig. S3A, B). The flat edges of bowtie structures were significantly deformed inward by microtissues, leading to shrinkage of their diagonal ribs and closing of pores (Fig. 3A). In addition, the vertical parts flipped 90° and met their neighbors due to the high compression force.

Biomechanical interaction of mesenchymal microtissues on the auxetic bowtie metamaterials.
To simulate the deformation characteristics of bowtie metamaterial scaffolds under mechanical loading from microtissues, the distributed displacement was set at the nodes of beam members at the flat edges, while the top beams were subjected to displacement that would cause them to bend inward (Fig. 3B). The mechanical force of the microtissue caused a snap-through buckling mode on the bowtie structure demonstrated as a plateau on the force–displacement curve. 44
Compared with the octet-truss structure, comprising multidirectional beam members that provide a homogeneous mechanical behavior, the high anisotropy degree of the bowtie structure based on its beam arrangement assured sufficient area for cell attachment without stretching the cell membrane to form thin branches. Due to disruption of structural units, the metamaterial scaffold became more compliant. However, the bowtie structure had less stress relaxation compared with the octet-truss structure since bowtie structures contained fewer buckled beams under the mechanical loading. More importantly, the negative Poisson's ratio of the bowtie structure allowed tissues to aggregate in a natural way, which brought the cells in close contact with strong intercellular connections. Hence, the entire scaffold was dramatically compressed, shown as inward buckling of individual beams at scaffold edges, flipping of vertical flat beams, and closing up of unit pores.
Results described above evidenced the fact that the mechanical load from biological tissues was prone to cause buckling of individual fibers, triggering the instability of metamaterials. However, the metamaterial scaffolds could sustain significant distortion without fracture, which possibly resulted from efficient tissue–scaffold integration as a single living composite. The active remodeling of mesenchymal microtissues made it possible to establish an equilibrium state of the composite under minimal mechanical force and energy potential. In comparison, uniform mechanical loading was passive and applied externally to the metamaterial scaffolds, thus it was not able to create the buckling forms that were generated internally from the living composite. The concept of living tissue–scaffold composite should be taken into consideration for future mechanistic analysis on tissue–material mechanical interaction and evolution of tissue instability.
Meanwhile, optimization of computational models is dependent on the accurate experimental measurement of spatial stiffness of the living tissue constructs. However, conventional atomic force microscopy or nanoindentation is limited for this measurement due to the highly hierarchical textures of 3D tissue constructs. Bioindenter devices are featured by their large displacement range, high-resolution force control, and capability of testing samples in a liquid condition, which are essential for measuring the mechanical properties of soft living tissues.45–47 In future, our microtissue model can be used to establish a correlation of metamaterial geometry with soft tissue mechanics based on bioindentation measurement.
Conclusions
In summary, we have generated the first 3D human microtissue model on metamaterial scaffolds fabricated using MPL technology. hiPSC-MSCs could assemble into 3D microtissues on different metamaterial scaffolds, but tissue morphology and scaffold deformation mode were significantly different between octet-truss and bowtie structures. On the microscopic level, we observed the microtissue remodeling adapting to the complex scaffold geometry. FEA simulations revealed that mesenchymal microtissues induced different deformation modes on the metamaterial scaffolds and caused a large transformation of their geometries. Further quantification of cell surface tension, immunostaining of other cell markers, and simulation of force distribution within the metastructures would benefit a deeper understanding of how unit geometries of metamaterials affect tissue phenotypes.
Data Availability
The authors declare that all data supporting the findings of this study are available. Source data for the figures are available from the corresponding author on reasonable request.
Footnotes
Authors' Contributions
Z.M., C.W., Z.V., and C.P.G. conceived and designed the experiments. C.W. conducted the biological experiments and confocal microscopy. Z.V. performed the fabrication of metamaterial scaffolds. T.W. differentiated and provided the hiPSC-MSCs. S.S. performed the SEM imaging. C.W. and Z.V. analyzed and interpreted the data. C.W., Z.V., C.P.G., and Z.M. wrote the manuscript with discussions and improvements from all authors. Z.M. and C.P.G. supervised the project development and funded the study.
Acknowledgments
The authors thank Dr. Elmina Kambouraki, George Flamourakis, and Dr. Maria Farsari from FORTH for assisting in fabrication of scaffolds. The HIM experiments were conducted at the California Institute for Quantitative Biosciences (QB3).
Author Disclosure Statement
No competing financial interests exist.
Funding Information
This work was supported by the NSF (CBET-1804875, CBET-1804922, CBET-1943798, and SNM-1449305), NIH NICHD (R01HD101130), and a Syracuse University intramural CUSE Grant.
References
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
