Search

A major bottleneck in the study of human liver physiology is the provision of stable liver tissue in sufficient quantity. As a result, current approaches to modelling human drug efficacy and toxicity rely heavily on immortalized human and animal cell lines. These models are informative but do possess significant drawbacks. To address the issues presented by those models, researchers have turned to pluripotent stem cells (PSCs). PSCs can be generated from defined genetic backgrounds, are scalable, and capable of differentiation to all the cell types found in the human body, representing an attractive source of somatic cells for in vitro and in vivo endeavours. Although unlimited numbers of somatic cell types can be generated in vitro, their maturation still remains problematic. In order to develop high fidelity PSC-derived liver tissue, it is necessary to better understand the cell microenvironment in vitro including key elements of liver physiology. In vivo a major driver of zonated liver function is the oxygen gradient that exists from periportal to pericentral regions. In this paper, we demonstrate how cell culture conditions for PSC-derived liver sphere systems can be optimised to recapitulate physiologically relevant oxygen gradients by using mathematical modelling. The mathematical model incorporates some often-understated features and mechanisms of traditional spheroid systems such as cell-specific oxygen uptake, media volume, spheroid size, and well dimensions that can lead to a spatially heterogeneous distribution of oxygen. This mathematical modelling approach allows for the calibration and identification of culture conditions required to generate physiologically realistic function within the microtissue through recapitulation of the in vivo microenvironment. [ABSTRACT FROM AUTHOR]

Cancer cells depend on glucose metabolism via glycolysis as a primary energy source, despite the presence of oxygen and fully functioning mitochondria, in order to promote growth, proliferation and longevity. Glycolysis relies upon NAD+ to accept electrons in the glyceraldehyde-3-phosphate dehydrogenase (GAPDH) reaction, linking the redox state of the cytosolic NAD+ pool to glycolytic rate. The free cytosolic NAD+/NADH ratio is involved in over 700 oxidoreductive enzymatic reactions and as such, the NAD+/NADH ratio is regarded as a metabolic readout of overall cellular redox state. Many experimental techniques that monitor or measure total NAD+ and NADH are unable to distinguish between protein-bound and unbound forms. Yet total NAD+/NADH measurements yield little information, since it is the free forms of NAD+ and NADH that determine the kinetic and thermodynamic influence of redox potential on glycolytic rate. Indirect estimations of free NAD+/NADH are based on the lactate/pyruvate (L/P) ratio at chemical equilibrium, but these measurements are often undermined by high lability. To elucidate the sensitivity of the free NAD+/NADH ratio to changes in extracellular substrate, an in silico model of hepatocarcinoma glycolysis was constructed and validated against in vitro data. Model simulations reveal that over experimentally relevant concentrations, changes in extracellular glucose and lactate concentration during routine cancer cell culture can lead to significant deviations in the NAD+/NADH ratio. Based on the principles of chemical equilibrium, the model provides a platform from which experimentally challenging situations may be examined, suggesting that extracellular substrates play an important role in cellular redox and bioenergetic homeostasis. [ABSTRACT FROM AUTHOR]

Paracetamol (acetaminophen (APAP)) is one of the most commonly used analgesics in the United Kingdom and the United States. However, exceeding the maximum recommended dose can cause serious liver injury and even death. Promising APAP toxicity biomarkers are thought to add value to those used currently and clarification of the functional relationships between these biomarkers and liver injury would aid clinical implementation of an improved APAP toxicity identification framework. The framework currently used to define an APAP overdose is highly dependent upon time since ingestion and initial dose; information that is often highly unpredictable. A pharmacokinetic/pharmacodynamic (PK/PD) APAP model has been built in order to understand the relationships between a panel of biomarkers and APAP dose. Visualization and statistical tools have been used to predict initial APAP dose and time since administration. Additionally, logistic regression analysis has been applied to histology data to provide a prediction of the probability of liver injury. [ABSTRACT FROM AUTHOR]

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 HIFp53 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. [ABSTRACT FROM AUTHOR]

[Display omitted] • Developments in drug delivery have prompted interest in chemical permeation via the oral mucosa. • Improvements to in silico modelling could aid research in this area. • In vitro data with mechanistic and finite element modelling drive a new in silico platform. • This in silico approach predicts chemical permeation of oral mucosa. 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. [ABSTRACT FROM AUTHOR]

Redox cycling is an understated mechanism of toxicity associated with a plethora of xenobiotics, responsible for preventing the effective treatment of serious conditions such as malaria and cardiomyopathy. Quinone compounds are notorious redox cyclers, present in drugs such as doxorubicin, which is used to treat a host of human cancers. However, the therapeutic index of doxorubicin is undermined by dose-dependent cardiotoxicity, which may be a function of futile redox cycling. In this study, a doxorubicin-specific in silico quinone redox metabolism model is described. Doxorubicin-GSH adduct formation kinetics are thermodynamically estimated from its reduction potential, while the remainder of the model is parameterised using oxygen consumption rate data, indicative of hydroquinone auto-oxidation. The model is then combined with a comprehensive glutathione metabolism model, facilitating the simulation of quinone redox cycling, and adduct-induced GSH depletion. Simulations suggest that glutathione pools are most sensitive to exposure duration at pharmacologically and supra-pharmacologically relevant doxorubicin concentrations. The model provides an alternative method of investigating and quantifying redox cycling induced oxidative stress, circumventing the experimental difficulties of measuring and tracking radical species. This in silico framework provides a platform from which GSH depletion can be explored as a function of a compound's physicochemical properties. [ABSTRACT FROM AUTHOR]

Abstract Many in vitro liver cell models, such as 2D systems, that are used to assess the hepatotoxic potential of xenobiotics suffer major limitations arising from a lack of preservation of physiological phenotype and metabolic competence. To circumvent some of these limitations there has been increased focus on producing more representative 3D models. Here we have used a novel approach to construct a size-controllable 3D hepatic spheroid model using freshly isolated primary rat hepatocytes (PRH) utilising the liquid-overlay technique whereby PRH spontaneously self-assemble in to 3D microtissues. This system produces viable spheroids with a compact in vivo -like structure for up to 21 days with sustained albumin production for the duration of the culture period. F-actin was seen throughout the spheroid body and P-glycoprotein (P-gp) and multidrug resistance-associated protein 2 (MRP2) transporters had polarised expression on the canalicular membrane of hepatocytes within the spheroids upon formation (day 3). The MRP2 transporter was able to functionally transport 5 μM 5-chloromethylfluorescein diacetate (CMFDA) substrates into these canalicular structures. These PRH spheroids display in vivo characteristics including direct cell-cell contacts, cellular polarisation, 3D cellular morphology, and formation of functional secondary structures throughout the spheroid. Such a well-characterised system could be readily exploited for pre-clinical and non-clinical repeat-dose investigations and could make a significant contribution to replace, reduce and refine the use of animals for applied research. Highlights • Initial cell seeding density allows for size-controllability of PRH spheroids. • PRH spheroids closely resemble in vivo morphology and microarchitecture. • PRH cultured in 3D are an improved system compared to hepatic-derived cell models. • PRH spheroids may be utilised for long-term repeat-dose investigations • Cost-effective methodology for the production of primary cell spheroids [ABSTRACT FROM AUTHOR]

Xenobiotic safety assessment is an area that impacts a multitude of different industry sectors such as medicinal drugs, agrochemicals, industrial chemicals, cosmetics and environmental contaminants. As such there are a number of well-developed in vitro, in vivo and in silico approaches to evaluate their properties and potential impact on the environment and to humans. Additionally, there is the continual investment in multidisciplinary scientists to explore non-animal surrogate technologies to predict specific toxicological outcomes and to improve our understanding of the biological processes regarding the toxic potential of xenobiotics. Here we provide a concise, critical evaluation of a number of in vitro systems utilised to assess the hepatotoxic potential of xenobiotics. [ABSTRACT FROM AUTHOR]

Nanoparticles have the potential to increase the efficacy of anticancer drugs whilst reducing off-target side effects. However, there remain uncertainties regarding the cellular uptake kinetics of nanoparticles which could have implications for nanoparticle design and delivery. Polymersomes are nanoparticle candidates for cancer therapy which encapsulate chemotherapy drugs. Here we develop a mathematical model to simulate the uptake of polymersomes via endocytosis, a process by which polymersomes bind to the cell surface before becoming internalised by the cell where they then break down, releasing their contents which could include chemotherapy drugs. We focus on two in vitro configurations relevant to the testing and development of cancer therapies: a well-mixed culture model and a tumour spheroid setup. Our mathematical model of the well-mixed culture model comprises a set of coupled ordinary differential equations for the unbound and bound polymersomes and associated binding dynamics. Using a singular perturbation analysis we identify an optimal number of ligands on the polymersome surface which maximises internalised polymersomes and thus intracellular chemotherapy drug concentration. In our mathematical model of the spheroid, a multiphase system of partial differential equations is developed to describe the spatial and temporal distribution of bound and unbound polymersomes via advection and diffusion, alongside oxygen, tumour growth, cell proliferation and viability. Consistent with experimental observations, the model predicts the evolution of oxygen gradients leading to a necrotic core. We investigate the impact of two different internalisation functions on spheroid growth, a constant and a bond dependent function. It was found that the constant function yields faster uptake and therefore chemotherapy delivery. We also show how various parameters, such as spheroid permeability, lead to travelling wave or steady-state solutions. [ABSTRACT FROM AUTHOR]