Software Tools For Cell Culture-Related 3D Printed Structures
Software Tools For Cell Culture-Related 3D Printed Structures
Software Tools For Cell Culture-Related 3D Printed Structures
* aczirok@gmail.com
Abstract
a1111111111 Three-dimensional (3D) printing technology allowed fast and cheap prototype fabrication in
a1111111111
a1111111111 numerous segments of industry and it also became an increasingly versatile experimental
a1111111111 platform in life sciences. Yet, general purpose software tools to control printer hardware are
a1111111111 often suboptimal for bioprinting applications. Here we report a package of open source soft-
ware tools that we developed specifically to meet bioprinting requirements: Machine move-
ments can be (i) precisely specified using high level programming languages, and (ii) easily
distributed across a batch of tissue culture dishes. To demonstrate the utility of the reported
OPEN ACCESS technique, we present custom fabricated, biocompatible 3D-printed plastic structures that
Citation: Gulyas M, Csiszer M, Mehes E, Czirok A can control cell spreading area or medium volume, and exhibit excellent optical properties
(2018) Software tools for cell culture-related 3D even at 50 ul sample volumes. We expect our software tools to be helpful not only to manu-
printed structures. PLoS ONE 13(9): e0203203. facture customized in vitro experimental chambers, but for applications involving printing
https://doi.org/10.1371/journal.pone.0203203
cells and extracellular matrices as well.
Editor: Marco Livesu, Consiglio Nazionale delle
Ricerche, ITALY
resulting tunels can be populated with endothelial cells and perfused with blood under high-
pressure pulsatile flow [8]. Using a novel coaxial printing method a functional blood vessel can
be fabricated, where the lumen is filled with a water-soluble material, while the bio-ink for the
wall is a composite of ECM proteins and endothelial progenitor cells [9]. The structure of the
extracellular environment can be also shaped by 3D printing technology. As an important
example, follicle explants can be cultured in a suitable gelatin mesh structure, and can be used
as a functional ovarian implant in surgically sterilized mice [10]. Such implanted, follicle-
seeded scaffolds become highly vascularized and ovarian function was fully restored: pups are
born through natural mating and thrive through maternal lactation.
Affordable 3D-printing also allows the development of specialized devices for in vitro cell
technology. As an example, 3D printed inserts can be used to grow and stimulate neurons
within geometrically constrained compartments [11]. By fabricating structures in culture
dishes, one can control cell spreading, hence local cell density, or by restricting medium vol-
ume, decrease the amount of necessary reagents. 3D-printing technology also allows a simple
in-lab fabrication of channels and reservoirs in cell culture dishes—on a cruder scale than lito-
graphy-based microfluidic chambers, but without requiring specialized equipment and at a
fraction of the cost [12, 13]. While 3D-printers, with small modifications, are capable bioengi-
neering tools, these applications also require software tools that have distinct requirements
from general-purpose slicer applications that convert a 3D solid object into a sequence of
machine movements, conventionally encoded in g-code [14]. Most importantly, bioengineer-
ing applications like cell printing or fabricating cell-scale environments often require well-
defined, specialized motion patterns and an ability to reproduce it in parallel targets, like an
array of culture dishes. Prescribing each machine movement by manually editing the low-level
g-code sequence is laborious, error-prone and time-consuming. For this reason we developed
a software package which can represent machine movements or g-code elements by simple
functions of a high level programming language such as python or C#. We also created a
graphic user interface, PetriPrinter, which distributes the programmatically defined printer
movements into several culture dishes organized in a grid pattern.
exposure to various amounts of polylactic acid (PLA), cell monolayers were fixed with 10% tri-
chloroacetic acid and stained for 15 min with sulforhodamine B (SRB). Cells were repeatedly
washed with 1% (v/v) acetic acid to remove excess dye. The protein-bound dye was dissolved
in 10 mM Tris and the optical density of the sample was determined at 570 nm using a micro-
plate reader (EL800, BioTec Instruments, Winooski, VT, USA).
2.6 Sterilization
3D printing in Petri dishes took place in open environment. After the printing process dishes
were sterilized with 70% (v/v) ethanol for 10 min then dried for 2 h under a sterile hood.
Fig 1. Overview of the various software tools and their applications. Geometrical objects can be encoded as python or C# files, as general use Standard Triangle
Language (.stl) files, or as a sequence of machine specific commands (g-code). Internally, extrusion and movement commands are represented as a binarized
gCodeCollection objects (.ptf) that can be edited and distributed into a number of parallel Petri dishes.
https://doi.org/10.1371/journal.pone.0203203.g001
distributes g-code objects into several culture dishes for batch printing. All parts of this soft-
ware package require at least a Windows XP system and a Microsoft.NET Framework 4 pack-
age, and are licensed under the Apache License 2.0.
3.1.1 gCodeAPI.NET. Native computer numerical control (CNC) programming language
(g-code) is a very limited programming language which does not contain variables or logical
operators and it is inconvenient to read by humans. This API encapsulates g-code commands
into functions of high-level programming languages such as C# or python (see code examples
in S1 Appendix). Besides the availability of the entire repertoire (variables, loops, flow control)
of these languages, gCodeAPI.NET adds easily interpretable commands such as line, arc, rela-
tive and absolute position, speed, extrude or temperature. Documentation of the command
library as well as code tutorials are available at the http://petriprinter.elte.hu website and on
github (http://github.com/gulyasmarton/PetriPrinter). As demonstrated in the S1 Appendix,
one can define a simple object like a cylinder surface with a few lines of code.
3.1.2 gCode Editor. Generic purpose CAD programs allow the design of complex struc-
tures in a more convenient way than building it up from the elements of our g-code API. Slicer
programs, which generate g-codes from a 3D solid model, do not offer a detailed control over
the tool movements. When printing delicate structures, however, some extruder head move-
ments may pull thin, undesired PLA threads and deposit them on the tissue culture surface.
Close proximity of the hot extruder head and the tissue culture surface during head travel
can also deteriorate the surface. Thus, there is a need to re-route such movements away from
regions where cells are expected to adhere. Finding such problematic machine movements in
the raw g-code, however, is quite difficult.
To make such tasks easier, our gCode Editor combines features of printing path viewers
and text editors: One can navigate among the g-code commands while a parser function trans-
lates the code to human readable format and the Path Viewer draws the printing path. Using
this visual information and GUI tools, one can identify and replace the problematic head
movements: travel paths can be augmented by inserting new control points and existing con-
trol points can be relocated by mouse clicks. The output of the editor is an API object which is
compatible with the PetriPrinter program.
3.1.3 PetriPrinter. The PetriPrinter application is a g-code generator program which han-
dles model files and aligns them into suitable positions to be printed into culture dishes. Using
a graphic user interface one can assign an arbitrary combination of distinct objects to the avail-
able positions. The software provides collision-safe and optimized entry, exit and movement
between the dishes.
Simple PLA-printable dish positioning objects like grids and rings are distributed as exam-
ples with the program package. The application can be adjusted to various hardware platforms
by specifying settings such as printing temperature, start height, and row/column distances.
When using dish holders, determining the heights of culture surfaces is crucial to obtain
reliable prints. According to our measurments, the average print-by-print variation of the cul-
ture dish surfaces is 0.2 mm (15 wells, 8 independent runs), larger than the height of a single
layer. Thus, we augmented the printer with a microswitch probe to measure the surface height
of each dish (see below), and utilize this information by the PetriPrinter software.
3.1.4 Validation. To test if the gCodeAPI commands are translated into correct g-Codes,
we translated test objects utilizing the basic elements of the API such as lines and ellipses. We
extracted the coordinates from the generated g-Codes, and determined the distance between
the points and the corresponding mathematical curve. For each object we found the distance
between the g-Code points and the mathematical curve to fall wihin the numerical precision
limit (10−8 mm).
Fig 2. Adaptation of an ultimaker original printer for culture dish printing. A: Stage accepting an array of 35 mm cell culture dishes. B: Modified printhead. Yellow
area: modified shroud to increase clearance to the stage and improve airflow. Blue area: a heater block that fits into the tissue culture dish. Red area: holder for a
removable microswitch sensor, used to map the height of each dish. C: Printhead with microswitch sensor installed (orange). Another microswitch sensor, mounted
on the stage (dark red) is used to calculate the height difference between the removable sensor and the nozzle.
https://doi.org/10.1371/journal.pone.0203203.g002
curvature of the medium interface. If wells have an appropriate conical shape, the equilibrium
contact angle between the wall and the medium is achieved by a flat medium-air interface.
To determine the suitable slope angle of the conical surface, we made a series of models
with angles ranging from vertical to outward leaning cones with an aperture angle of 100˚. As
Fig 3 demonstrates, for DMEM within a PLA container the phase contrast quality gradually
improves until the aperture of the cone reaches 80˚, i.e. the slope of the walls is 40 ˚. This
improvement reflects both a decreasing curvature of the medium surface and an increasing
distance of the medium-wall contact area.
3.2.3 PLA structures in tissue culture dishes. We demonstrate the utility of the approach
by some simple 3D printed structures in 35 mm culture dishes that are useful in various exper-
imental studies. Such dishes are convenient for time-lapse studies as they offer better phase
contrast images than 24 or 96 well plates.
Dish 1 (Fig 4A) allows to maintain seven independent cultures within the same 35 mm
dish, without sharing media and without deterioration of optical quality. Each well consists of
a 3 mm tall cylindrical wall, which expands into a conical surface of 1 mm high with the previ-
ously determined apex angle of 100˚. Dish 2 (Fig 4B) contains a small ring to limit the spread-
ing of cells and a higher, larger diameter outer wall to set the volume of medium within the
same 35 mm dish. Limiting cell spreading to a certain area can be used to locally set cell den-
sity. In particular, high cell densities can be maintained with a high medium volume/cell ratio,
enabling long term observation of cultures [18]. Dish 3 (Fig 4C and 4D) contains 1 mm2 mini-
wells that can be used to isolate cell clones, possibly obtained with a single cell sorter [19]. A
microscopic image of the single layer PLA wall indicates the resolution limits obtainable with
such technology: the 0.4 mm wide extrusion nozzle yields 120 μm tall and 920±90 μm wide
PLA barriers on the tissue culture plastic. The mean area between the barriers was determined
to be 1.82 mm2, while the standard deviation of individual cell areas was found to be 1.6.
3.2.4 Biocompatibility of PLA structures. To test the acute (short-term) cytotoxicity of
fusion-deposited PLA structures, we exposed two well-established cell lines, 3T3 fibroblasts
and p31 mesothelioma tumor cells, to various amounts of PLA. After two days in culture, the
total protein content of the cells, as a viability score, was determined by the SRB assay. As Fig 5
indicates, the viability of the cells was not affected by the presence of PLA.
Long term cytotoxicity was assayed using a primary hippocampal neuron culture (Fig 6).
Cells were seeded onto Poly-L-lysine-coated surfaces at a density of 900 cells/mm2. The
Fig 3. Optimization of conical wells for phase contrast optics. A: A series of wells were printed with various wall
slopes or cone aperture angles (left column). 3T3 cells within the wells were imaged using an automated phase contrast
microscope. The tiled mosaic fields are shown in the middle column (scale bar: 500 μm). Red rectangles mark
individual micrographs, that are shown in the right column (scale bar: 100 μm). The best optical image corresponds to
a DMEM-PLA contact angle of 40˚. B: Quantitative analysis of image contrast. For each well we assigned a contrast
quality value, and data is pooled from n = 3 parallel sets of wells. Asterisks denote significant (p < 0.05) difference from
the image obtained in cylindrical wells, error bars indicate standard deviation.
https://doi.org/10.1371/journal.pone.0203203.g003
Fig 4. Various 3D printed culture dish configurations. A: Seven-well dish with a sloped wall to provide optimal optical quality. B: A double ring dish to control cell
spreading and medium volume. C: A grid with 38 rectangular cells to store manually sorted cells. D: Phase contrast micrograph of a single PLA layer adhering to a
tissue culture substrate, part of the cell sorter grid shown in panel C. Air bubbles indicate small areas with imperfect contact between the PLA layer and the underlying
substrate. Scale bar: 500 μm.
https://doi.org/10.1371/journal.pone.0203203.g004
Fig 5. Cytotoxity assay for PLA exposure. Cells were grown in the presence of PLA surfaces for two days after which
their protein content was determined by an SRB assay. Data are shown as average of n = 4 parallel experiments for two
cell ines (A: 3T3 cells, B: p31 cells) and the effect of treatment is expressed relative to untreated controls. Error bars
represent standard deviation.
https://doi.org/10.1371/journal.pone.0203203.g005
Fig 6. Long-term biocompatibility assay with primary hippocampal neuron culture. A: A triple ring dish used for
cell culture. B: To achieve high cell density, cell spreading is restricted within the rings (diameter: 6 mm, height: 1.8
mm). A tiled image of phase contrast micrographs demonstrate the optical quality of the setup. Scalebar: 1000 μm. The
red rectangle marks a region shown in higher magnification in panel C. D: Composite fluorescent image to evaluate
cell death within the same region shown in panel C. Live cells are labelled by calcein (green), while dead cells are
labelled by ethidium (red) and marked by arrows.
https://doi.org/10.1371/journal.pone.0203203.g006
culture dish contained three PLA rings to limit cell spreading and to maintain high medium/
cell ratio. These conditions are favorable for long-term cell culture. After 7 days in vitro, cell
viability was determined using a fluorescent live/dead staining kit, which indicated a cell death
rate of less than 2.2.
4 Conclusions
3D printing is becoming wide-spread in life sciences, either to provide cost-effective alterna-
tive to existing experimental devices, or to fabricate highly specialized structures that are
specifically designed to meet certain experimental requirements. To effectively transfer 3D
printing know-how, some uniformization—a common language—is required. Our approach
can represent a variety of objects by sufficient detail and precision, yet in a conceptionally
clearer form than a sequence of machine-specific low-level movement commands. We hope
that this platform will thus facilitate both the distribution as well as the development of new,
experiment-specific devices or in vitro cell technologies.
Supporting information
S1 Appendix. Example codes using the gCodeAPI. A cylinder surface is defined in C# (A) or
python (B) language.
(PDF)
Acknowledgments
We are grateful to Tamás Garay for his help with the SRB cytotoxicity assay, Krisztián Tárnok
and Katalin Schlett for their gift of the hippocampal neuron culture used in this study.
Author Contributions
Investigation: Marton Gulyas, Elod Mehes.
Methodology: Miklos Csiszer.
Software: Marton Gulyas.
Supervision: Andras Czirok.
Writing – original draft: Marton Gulyas, Andras Czirok.
Writing – review & editing: Elod Mehes, Andras Czirok.
References
1. Shi P, Tan YSE, Yeong WY, Li HY, Laude A. A bilayer photoreceptor-retinal tissue model with gradient
cell density design: A study of microvalve-based bioprinting. J Tissue Eng Regen Med. 2018; 12
(5):1297–1306. https://doi.org/10.1002/term.2661 PMID: 29510003
2. Moldovan NI. Progress in scaffold-free bioprinting for cardiovascular medicine. J Cell Mol Med. 2018;
22(6):2964–2969. https://doi.org/10.1111/jcmm.13598 PMID: 29536627
3. Lopa S, Mondadori C, Mainardi VL, Talò G, Costantini M, Candrian C, et al. Translational Application of
Microfluidics and Bioprinting for Stem Cell-Based Cartilage Repair. Stem Cells Int. 2018;
2018:6594841. https://doi.org/10.1155/2018/6594841 PMID: 29535776
4. Ng WL, Lee JM, Yeong WY, Naing MW. Microvalve-based bioprinting—process, bio-inks and applica-
tions. Biomater Sci. 2017; 5(4):632–647. https://doi.org/10.1039/c6bm00861e PMID: 28198902
5. Jammalamadaka U, Tappa K. Recent Advances in Biomaterials for 3D Printing and Tissue Engineering.
J Funct Biomater. 2018; 9(1).
6. Xia Z, Jin S, Ye K. Tissue and Organ 3D Bioprinting. SLAS Technol. 2018; 23(4):301–314. https://doi.
org/10.1177/2472630318760515 PMID: 29474789
7. Kolesky DB, Homan KA, Skylar-Scott MA, Lewis JA. Three-dimensional bioprinting of thick vascular-
ized tissues. Proc Natl Acad Sci USA. 2016; 113(12):3179–84. https://doi.org/10.1073/pnas.
1521342113 PMID: 26951646
8. Miller JS, Stevens KR, Yang MT, Baker BM, Nguyen DHT, Cohen DM, et al. Rapid casting of patterned
vascular networks for perfusable engineered three-dimensional tissues. Nat Mater. 2012; 11(9):768–
74. https://doi.org/10.1038/nmat3357 PMID: 22751181
9. Gao G, Lee JH, Jang J, Lee DH, Kong JS, Kim BS, et al. Tissue Engineered Bio-Blood-Vessels Con-
structed Using a Tissue-Specific Bioink and 3D Coaxial Cell Printing Technique: A Novel Therapy for
Ischemic Disease. Advanced Functional Materials. 2017; 27(33):1700798. https://doi.org/10.1002/
adfm.201700798
10. Laronda MM, Rutz AL, Xiao S, Whelan KA, Duncan FE, Roth EW, et al. A bioprosthetic ovary created
using 3D printed microporous scaffolds restores ovarian function in sterilized mice. Nat Commun. 2017;
8:15261. https://doi.org/10.1038/ncomms15261 PMID: 28509899
11. Wardyn JD, Sanderson C, Swan LE, Stagi M. Low cost production of 3D-printed devices and electrosti-
mulation chambers for the culture of primary neurons. J Neurosci Methods. 2015; 251:17–23. https://
doi.org/10.1016/j.jneumeth.2015.05.001 PMID: 25962333
12. Rusling JF. Developing Microfluidic Sensing Devices Using 3D Printing. ACS Sens. 2018; 3(3):522–
526. https://doi.org/10.1021/acssensors.8b00079 PMID: 29490458
13. Nie J, Gao Q, Qiu JJ, Sun M, Liu A, Shao L, et al. 3D printed Lego®-like modular microfluidic devices
based on capillary driving. Biofabrication. 2018; 10(3):035001. https://doi.org/10.1088/1758-5090/
aaadd3 PMID: 29417931
14. Oberg E, Jones FD, Horton HL, Ryffel HH, Green RE, McCauley CJ, editors. Machinery’s handbook.
New York, USA: Industrial Press, Inc.; 1996.
15. Czöndör K, Ellwanger K, Fuchs YF, Lutz S, Gulyás M, Mansuy IM, et al. Protein kinase D controls the
integrity of Golgi apparatus and the maintenance of dendritic arborization in hippocampal neurons. Mol
Biol Cell. 2009; 20(7):2108–20. https://doi.org/10.1091/mbc.E08-09-0957 PMID: 19211839
16. Orellana EA, Kasinski AL. Sulforhodamine B (SRB) Assay in Cell Culture to Investigate Cell Prolifera-
tion. Bio Protoc. 2016; 6(21). https://doi.org/10.21769/BioProtoc.1984 PMID: 28573164
17. Schindelin J, Arganda-Carreras I, Frise E, Kaynig V, Longair M, Pietzsch T, et al. Fiji: an open-source
platform for biological-image analysis. Nat Methods. 2012; 9(7):676–82. https://doi.org/10.1038/nmeth.
2019 PMID: 22743772
18. Tarnoki-Zach J, Isai DG, Mehes E, Paku S, Neufeld Z, Hegedus B, et al. Myosin-II dependent cell con-
tractility contributes to spontaneous nodule formation of mesothelioma cells. arXiv:150608090;.
19. Ungai-Salánki R, Gerecsei T, Fürjes P, Orgovan N, Sándor N, Holczer E, et al. Automated single cell
isolation from suspension with computer vision. Sci Rep. 2016; 6:20375. https://doi.org/10.1038/
srep20375 PMID: 26856740