Adrian Neagu

Tissue Engineering by Self-Assembly of Cells Printed into Topologically Defined Structures

Understanding the principles of biological self-assembly is indispensable for developing efficient strategies to build living tissues and organs. We exploit the self-organizing capacity of cells and tissues to construct functional living structures of prescribed shape. In our technology, multicellular spheroids (bio-ink particles) are placed into biocompatible environment (bio-paper) by the use of a three-dimensional delivery device (bio-printer). Our approach mimics early morphogenesis and is based on the realization that the genetic control of developmental patterning through self-assembly involves physical mechanisms. Three-dimensional tissue structures are formed through the postprinting fusion of the bio-ink particles, in analogy with early structure-forming processes in the embryo that utilize the apparent liquid-like behavior of tissues composed of motile and adhesive cells. We modeled the process of self-assembly by fusion of bio-ink particles, and employed this novel technology to print extended cellular structures of various shapes. Functionality was tested on cardiac constructs built from embryonic cardiac and endothelial cells. The postprinting self-assembly of bio-ink particles resulted in synchronously beating solid tissue blocks, showing signs of early vascularization, with the endothelial cells organized into vessel-like conduits.

Developmental Biology and Tissue Engineering

Morphogenesis implies the controlled spatial organization of cells that gives rise to tissues and organs in early embryonic development. While morphogenesis is under strict genetic control, the formation of specialized biological structures of specific shape hinges on physical processes. Tissue engineering (TE) aims at reproducing morphogenesis in the laboratory, i.e., in vitro, to fabricate replacement organs for regenerative medicine. The classical approach to generate tissues/organs is by seeding and expanding cells in appropriately shaped biocompatible scaffolds, in the hope that the maturation process will result in the desired structure. To accomplish this goal more naturally and efficiently, we set up and implemented a novel TE method that is based on principles of developmental biology and employs bioprinting, the automated delivery of cellular composites into a three-dimensional (3D) biocompatible environment. The novel technology relies on the concept of tissue liquidity according to which multicellular aggregates composed of adhesive and motile cells behave in analogy with liquids: in particular, they fuse. We emphasize the major role played by tissue fusion in the embryo and explain how the parameters (surface tension, viscosity) that govern tissue fusion can be used both experimentally and theoretically to control and simulate the self-assembly of cellular spheroids into 3D living structures. The experimentally observed postprinting shape evolution of tube- and sheet-like constructs is presented. Computer simulations, based on a liquid model, support the idea that tissue liquidity may provide a mechanism for in vitro organ building.

Fusion of uniluminal vascular spheroids: A model for assembly of blood vessels

We evaluated the self‐assembly properties of uniluminal vascular spheroids having outer layers of vascular smooth muscle cells and a contiguous inner layer of endothelial cells lining a central lumen. We showed that while pairs of uniluminal vascular spheroids suspended in culture medium fused to form a larger diameter spheroidal structure, spheroids in collagen hydrogels formed elongated structures. These findings highlight the potential use of uniluminal vascular spheroids as modules to engineer blood vessels. We also demonstrate that uniluminal vascular spheroid fusion conforms to models describing the coalescence of liquid drops. Furthermore, the fusion of uniluminal vascular spheroids in vitro closely resembled the in vivo process by which the descending aorta forms from the fusion of the paired dorsal aortae during embryonic development. Together, the findings indicate that tissue liquidity underlies uniluminal vascular spheroid fusion and that in vivo anastomosis of blood vessels may involve a similar mechanism. Developmental Dynamics 239:398–406, 2010.

Fluctuations and the Hofmeister Effect

The Hofmeister effect consists in changes of protein solubility triggered by neutral electrolyte cosolutes. Based on the assumption that salts cause stochastic fluctuations of the free energy barrier profiles, a kinetic theory of this phenomenon is proposed. An exponentially correlated noise, of amplitude proportional to the salt concentration, is added to each energy level, and the time-dependence of the mean protein concentration is calculated. It is found that the theory yields the well-known Setschenow equation if the noise correlation time is low in comparison to the aggregation time constant. Experimental data on salting-in agents are well fitted, whereas, in the case of salting-out cosolutes, two independent dichotomic fluctuations are needed to fit the data. This may result from the fact that, in both cases, the low-concentration regime is dominated by salting-in electrostatic contributions, whereas, at high salt concentrations, electron donor/acceptor interactions become important; these have opposite effects. The theory offers a novel way to metricate Hofmeister effects and also leads to thermodynamic quantities, which account for the influence of salts. The formalism may also be applied to describe kinetic phenomena in the presence of cosolutes.

Kinetic Monte Carlo and cellular particle dynamics simulations of multicellular systems.

Computer modeling of multicellular systems has been a valuable tool for interpreting and guiding in vitro experiments relevant to embryonic morphogenesis, tumor growth, angiogenesis and, lately, structure formation following the printing of cell aggregates as bioink particles. Here we formulate two computer simulation methods: (1) a kinetic Monte Carlo (KMC) and (2) a cellular particle dynamics (CPD) method, which are capable of describing and predicting the shape evolution in time of three-dimensional multicellular systems during their biomechanical relaxation. Our work is motivated by the need of developing quantitative methods for optimizing postprinting structure formation in bioprinting-assisted tissue engineering. The KMC and CPD model parameters are determined and calibrated by using an original computational-theoretical-experimental framework applied to the fusion of two spherical cell aggregates. The two methods are used to predict the (1) formation of a toroidal structure through fusion of spherical aggregates and (2) cell sorting within an aggregate formed by two types of cells with different adhesivities.

Relating Biophysical Properties Across Scales

A distinguishing feature of a multicellular living system is that it operates at various scales, from the intracellular to organismal. Genes and molecules set up the conditions for the physical processes to act, in particular to shape the embryo. As development continues the changes brought about by the physical processes lead to changes in gene expression. It is this coordinated interplay between genetic and generic (i.e., physical and chemical) processes that constitutes the modern understanding of early morphogenesis. It is natural to assume that in this multiscale process the smaller defines the larger. In case of biophysical properties, in particular, those at the subcellular level are expected to give rise to those at the tissue level and beyond. Indeed, the physical properties of tissues vary greatly from the liquid to solid. Very little is known at present on how tissue level properties are related to cell and subcellular properties. Modern measurement techniques provide quantitative results at both the intracellular and tissue level, but not on the connection between these. In the present work we outline a framework to address this connection. We specifically concentrate on the morphogenetic process of tissue fusion, by following the coalescence of two contiguous multicellular aggregates. The time evolution of this process can accurately be described by the theory of viscous liquids. We also study fusion by Monte Carlo simulations and a novel Cellular Particle Dynamics (CPD) model, which is similar to the earlier introduced Subcellular Element Model (SEM; Newman, 2005). Using the combination of experiments, theory and modeling we are able to relate the measured tissue level biophysical quantities to subcellular parameters. Our approach has validity beyond the particular morphogenetic process considered here and provides a general way to relate biophysical properties across scales.

Experimental evaluation of apparent tissue surface tension based on the exact solution of the Laplace equation

The notion of apparent tissue surface tension offered a systematic way to interpret certain morphogenetic processes in early development. It also allowed deducing quantitative information on cellular and molecular parameters that is otherwise difficult to obtain. To accurately determine such tensions we combined novel experiments with the exact solution of the Laplace equation for the profile of a liquid drop under the employed experimental conditions and used the exact solution to evaluate data collected on tissues. Our results confirm that tissues composed of adhesive and motile cells indeed can be characterized in terms of well-defined apparent surface tension. Our experimental technique presents a way to measure liquid interfacial tensions under conditions when known methods fail.

COMPUTATIONAL MODELING OF TISSUE SELF-ASSEMBLY

Starting from the beginning of the 20th century, theoretical models of living tissues have evolved along two distinct conceptual lines. The first of these views considers the tissue as a set of discrete, interacting cells, whereas the other treats it as a continuum, and monitors cell densities instead of individual cells1. Here we briefly describe a few of these models. The interested reader can find further details in the cited literature.

Computational modeling of epithelial–mesenchymal transformations

An epithelial-mesenchymal transformation (EMT) involves alterations in cell-cell and cell-matrix adhesion, the detachment of epithelial cells from their neighbors, the degradation of the basal lamina and acquisition of mesenchymal phenotype. Here we present Monte Carlo simulations for a specific EMT in early heart development: the formation of cardiac cushions. Cell rearrangements are described in accordance with Steinbergs differential adhesion hypothesis, which states that cells possess a type-dependent adhesion apparatus and are sufficiently motile to give rise to the tissue conformation with the largest number of strong bonds. We also implement epithelial and mesenchymal cell proliferation, cell type change and extracellular matrix production by mesenchymal cells. Our results show that an EMT is promoted more efficiently by an increase in cell-substrate adhesion than by a decrease in cell-cell adhesion. In addition to cushion tissue formation, the model also accounts for the phenomena of matrix invasion and mesenchymal condensation. We conclude that in order to maintain epithelial integrity during EMT the number of epithelial cells must increase at a controlled rate. Our model predictions are in qualitative agreement with available experimental data.

Spatial distribution of VEGF isoforms and chemotactic signals in the vicinity of a tumor

We propose a mathematical model that describes the formation of gradients of different isoforms of vascular endothelial growth factor (VEGF). VEGF is crucial in the process of tumor-induced angiogenesis, and recent experiments strongly suggest that the molecule is most potent when bound to the extracellular matrix (ECM). Using a system of reaction-diffusion equations, we study diffusion of VEGF, binding of VEGF to the ECM, and cleavage of VEGF from the ECM by matrix metalloproteases (MMPs). We find that spontaneous gradients of matrix-bound VEGF are possible for an isoform that binds weakly to the ECM (i.e. VEGF(165)), but cleavage by MMPs is required to form long-range gradients of isoforms that bind rapidly to the ECM (i.e. VEGF(189)). We also find that gradient strengths and ranges are regulated by MMPs. Finally, we find that VEGF molecules cleaved from the ECM may be distributed in patterns that are not conducive to chemotactic migration toward a tumor, depending on the spatial distribution of MMP molecules. Our model elegantly explains a number of in vivo observations concerning the significance of different VEGF isoforms, points to VEGF(165) as an especially significant therapeutic target and indicator of a tumors angiogenic potential, and enables predictions that are subject to testing with in vitro experiments.