Innovative organotypic in vitro models for safety assessment: aligning with regulatory requirements and understanding models of the heart, skin, and liver as paradigms.

The development of improved, innovative models for the detection of toxicity of drugs, chemicals, or chemicals in cosmetics is crucial to efficiently bring new products safely to market in a cost-effective and timely manner. In addition, improvement in models to detect toxicity may reduce the incidence of unexpected post-marketing toxicity and reduce or eliminate the need for animal testing. The safety of novel products of the pharmaceutical, chemical, or cosmetics industry must be assured; therefore, toxicological properties need to be assessed. Accepted methods for gathering the information required by law for approval of substances are often animal methods. To reduce, refine, and replace animal testing, innovative organotypic in vitro models have emerged. Such models appear at different levels of complexity ranging from simpler, self-organized three-dimensional (3D) cell cultures up to more advanced scaffold-based co-cultures consisting of multiple cell types. This review provides an overview of recent developments in the field of toxicity testing with in vitro models for three major organ types: heart, skin, and liver. This review also examines regulatory aspects of such models in Europe and the UK, and summarizes best practices to facilitate the acceptance and appropriate use of advanced in vitro models.

Tight control of hypoxia-inducible factor-α transient dynamics is essential for cell survival in hypoxia.

Intracellular signaling involving hypoxia-inducible factor (HIF) controls the adaptive responses to hypoxia. There is a growing body of evidence demonstrating that intracellular signals encode temporal information. Thus, the dynamics of protein levels, as well as protein quantity and/or localization, impacts on cell fate. We hypothesized that such temporal encoding has a role in HIF signaling and cell fate decisions triggered by hypoxic conditions. Using live cell imaging in a controlled oxygen environment, we observed transient 3-h pulses of HIF-1α and -2α expression under continuous hypoxia. We postulated that the well described prolyl hydroxylase (PHD) oxygen sensors and HIF negative feedback regulators could be the origin of the pulsatile HIF dynamics. We used iterative mathematical modeling and experimental analysis to scrutinize which parameter of the PHD feedback could control HIF timing and we probed for the functional redundancy between the three main PHD proteins. We identified PHD2 as the main PHD responsible for HIF peak duration. We then demonstrated that this has important consequences, because the transient nature of the HIF pulse prevents cell death by avoiding transcription of p53-dependent pro-apoptotic genes. We have further shown the importance of considering HIF dynamics for coupling mathematical models by using a described HIF-p53 mathematical model. Our results indicate that the tight control of HIF transient dynamics has important functional consequences on the cross-talk with key signaling pathways controlling cell survival, which is likely to impact on HIF targeting strategies for hypoxia-associated diseases such as tumor progression and ischemia.

Transport of fibroblast growth factor 2 in the pericellular matrix is controlled by the spatial distribution of its binding sites in heparan sulfate.

The heparan sulfate (HS) chains of proteoglycans are a key regulatory component of the extracellular matrices of animal cells, including the pericellular matrix around the plasma membrane. In these matrices they regulate transport, gradient formation, and effector functions of over 400 proteins central to cell communication. HS from different matrices differs in its selectivity for its protein partners. However, there has been no direct test of how HS in the matrix regulates the transport of its partner proteins. We address this issue by single molecule imaging and tracking in fibroblast pericellular matrix of fibroblast growth factor 2 (FGF2), stoichiometrically labelled with small gold nanoparticles. Transmission electron microscopy and photothermal heterodyne imaging (PHI) show that the spatial distribution of the HS-binding sites for FGF2 in the pericellular matrix is heterogeneous over length scales ranging from 22 nm to several µm. Tracking of individual FGF2 by PHI in the pericellular matrix of living cells demonstrates that they undergo five distinct types of motion. They spend much of their time in confined motion (∼110 nm diameter), but they are not trapped and can escape by simple diffusion, which may be slow, fast, or directed. These substantial translocations (µm) cover distances far greater than the length of a single HS chain. Similar molecular motion persists in fixed cells, where the movement of membrane PGs is impeded. We conclude that FGF2 moves within the pericellular matrix by translocating from one HS-binding site to another. The binding sites on HS chains form non-random, heterogeneous networks. These promote FGF2 confinement or substantial translocation depending on their spatial organisation. We propose that this spatial organisation, coupled to the relative selectivity and the availability of HS-binding sites, determines the transport of FGF2 in matrices. Similar mechanisms are likely to underpin the movement of many other HS-binding effectors.

An advanced in silico model of the oral mucosa reveals the impact of extracellular spaces on chemical permeation.

Accurately predicting the permeation of chemicals through human epithelial tissues is crucial for pharmaceutical therapeutic design and toxicology. Current mathematical models of multi-layered stratified squamous epithelium such as those in the oral cavity use simplistic 'bricks and mortar' geometries that do not fully account for the complex cellular architecture that may affect chemical permeation in these tissues. Here we aimed to develop a new, advanced mechanistic mathematical model of the human epithelium that more accurately represents chemical tissue permeation. Using measurements of cell size and tortuosity from micrograph images of both human oral (buccal) and tissue-engineered buccal mucosa along with mechanistic mathematical modelling, we show that the convoluted geometry of the extracellular spaces within the epithelium significantly impacts chemical permeation. We next developed an advanced histologically and physiologically-relevant in silico model of buccal mucosal chemical permeation using partial differential equations, fitted to chemical permeation from in vitro assay data derived from tissue-engineered buccal mucosal models and chemicals with known physiochemical properties. Our novel in silico model can predict epithelial permeation kinetics for chemicals with different physicochemical properties in the absence or presence of permeability enhancers. This in vitro - in silico approach constitutes a step-change in the modelling of chemical tissue permeation and has the potential to expedite pharmaceutical innovation by improved and more rapid screening of chemical entities whilst reducing the need for in vivo animal experiments.

Diffusion about the mean drift location in a biased random walk.

Random walks are used to model movement in a wide variety of contexts: from the movement of cells undergoing chemotaxis to the migration of animals. In a two-dimensional biased random walk, the diffusion about the mean drift position is entirely dependent on the moments of the angular distribution used to determine the movement direction at each step. Here we consider biased random walks using several different angular distributions and derive expressions for the diffusion coefficients in each direction based on either a fixed or variable movement speed, and we use these to generate a probability density function for the long-time spatial distribution. We demonstrate how diffusion is typically anisotropic around the mean drift position and illustrate these theoretical results using computer simulations. We relate these results to earlier studies of swimming microorganisms and explain how the results can be generalized to other types of animal movement.

Multiscale modelling of drug transport and metabolism in liver spheroids.

In early preclinical drug development, potential candidates are tested in the laboratory using isolated cells. These in vitro experiments traditionally involve cells cultured in a two-dimensional monolayer environment. However, cells cultured in three-dimensional spheroid systems have been shown to more closely resemble the functionality and morphology of cells in vivo. While the increasing usage of hepatic spheroid cultures allows for more relevant experimentation in a more realistic biological environment, the underlying physical processes of drug transport, uptake and metabolism contributing to the spatial distribution of drugs in these spheroids remain poorly understood. The development of a multiscale mathematical modelling framework describing the spatio-temporal dynamics of drugs in multicellular environments enables mechanistic insight into the behaviour of these systems. Here, our analysis of cell membrane permeation and porosity throughout the spheroid reveals the impact of these properties on drug penetration, with maximal disparity between zonal metabolism rates occurring for drugs of intermediate lipophilicity. Our research shows how mathematical models can be used to simulate the activity and transport of drugs in hepatic spheroids and in principle any organoid, with the ultimate aim of better informing experimentalists on how to regulate dosing and culture conditions to more effectively optimize drug delivery.

Modeling the dynamics of hypoxia inducible factor-1α (HIF-1α) within single cells and 3D cell culture systems.

HIF (hypoxia inducible factor) is an oxygen-regulated transcription factor that mediates the intracellular response to hypoxia in human cells. There is increasing evidence that cell signaling pathways encode temporal information, and thus cell fate may be determined by the dynamics of protein levels. We have developed a mathematical model to describe the transient dynamics of the HIF-1α protein measured in single cells subjected to hypoxic shock. The essential characteristics of these data are modeled with a system of differential equations describing the feedback inhibition between HIF-1α and prolyl hydroxylases (PHD) oxygen sensors. Heterogeneity in the single-cell data is accounted through parameter variation in the model. We previously identified the PHD2 isoform as the main PHD sensor responsible for controlling the HIF-1α transient response, and make here testable predictions regarding HIF-1α dynamics subject to repetitive hypoxic pulses. The model is further developed to describe the dynamics of HIF-1α in cells cultured as 3D spheroids, with oxygen dynamics parameterized using experimental measurements of oxygen within spheroids. We show that the dynamics of HIF-1α and transcriptional targets of HIF-1α display a non-monotone response to the oxygen dynamics. Specifically we demonstrate that the dynamic transient behavior of HIF-1α results in differential dynamics in transcriptional targets.