Assessing the transferability and reproducibility of 3D in vitro liver models from primary human multi-cellular microtissues to cell-line based HepG2 spheroids.

To reduce, replace, and refine in vivo testing, there is increasing emphasis on the development of more physiologically relevant in vitro test systems to improve the reliability of non-animal-based methods for hazard assessment. When developing new approach methodologies, it is important to standardize the protocols and demonstrate the methods can be reproduced by multiple laboratories. The aim of this study was to assess the transferability and reproducibility of two advanced in vitro liver models, the Primary Human multicellular microtissue liver model (PHH) and the 3D HepG2 Spheroid Model, for nanomaterial (NM) and chemical hazard assessment purposes. The PHH model inter-laboratory trial showed strong consistency across the testing sites. All laboratories evaluated cytokine release and cytotoxicity following exposure to titanium dioxide (TiO2) and zinc oxide (ZnO) nanoparticles. No significant difference was observed in cytotoxicity or IL-8 release for the test materials. The data were reproducible with all three laboratories with control readouts within a similar range. The PHH model ZnO induced the greatest cytotoxicity response at 50.0 µg/mL and a dose-dependent increase in IL-8 release. For the 3D HepG2 spheroid model, all test sites were able to construct the model and demonstrated good concordance in IL-8 cytokine release and genotoxicity data. This trial demonstrates the successful transfer of new approach methodologies across multiple laboratories, with good reproducibility for several hazard endpoints.

Next-generation risk assessment of chemicals - Rolling out a human-centric testing strategy to drive 3R implementation: The RISK-HUNT3R project perspective.

In many industrial sectors, there is a need for reliable ways to evaluate the safety of chemicals with methods anchored to human biology and pathology. For this purpose, many animal-free new approach methods (NAMs) have been developed and implemented in various stages of risk assessment. Now it is time to assemble individual NAMs into a comprehensive next-generation risk assessment (NGRA) strategy. The European Horizon 2020 RISK-HUNT3R project (Risk assessment of chemicals integrating human-centric next-generation testing strategies promoting the 3Rs) has been designed to promote a combination of computational toxicology, in vitro toxicology, and systems biology. It is anticipated that this approach will lead to faster and more accurate risk assessment procedures. The RISK-HUNT3R NGRA strategy will be developed to address the implementation of a comprehensive NAM toolbox into the regulatory framework. Critical conceptual approaches of the project include i) the integration of human-relevant data on biotransformation and elimination, ii) the translation of high-content mode-of-action datasets into predictions of adverse outcomes, iii) development of quantitative adverse outcome pathways (qAOPs), and iv) quantification of uncertainties associated with the predictions based on NGRA strategies. Many of the project steps will be used iteratively to generate datasets with sufficient quality and certainty for NGRA. Scientists and regulators will work together on case studies to evaluate practical applicability of NAMs and strategies to combine information therefrom. Here we delineate how the strategy will be deployed to establish an overall NGRA framework for chemicals, pesticides, food additives, and drugs.

Application of high-throughput transcriptomics for mechanism-based biological read-across of short-chain carboxylic acid analogues of valproic acid.

Chemical read-across is commonly evaluated without specific knowledge of the biological mechanisms leading to observed adverse outcomes in vivo. Integrating data that indicate shared modes of action in humans will strengthen read-across cases. Here we studied transcriptomic responses of primary human hepatocytes (PHH) to a large panel of carboxylic acids to include detailed mode-of-action data as a proof-of-concept for read-across in risk assessment. In rodents, some carboxylic acids, including valproic acid (VPA), are known to cause hepatic steatosis, whereas others do not. We investigated transcriptomics responses of PHHs exposed for 24 h to 18 structurally different VPA analogues in a concentration range to determine biological similarity in relation to in vivo steatotic potential. Using a targeted high-throughput screening assay, we assessed the differential expression of ~3,000 genes covering relevant biological pathways. Differentially expressed gene analysis revealed differences in potency of carboxylic acids, and expression patterns were highly similar for structurally similar compounds. Strong clustering occurred for steatosis-positive versus steatosis-negative carboxylic acids. To quantitatively define biological read-across, we combined pathway analysis and weighted gene co-expression network analysis. Active carboxylic acids displayed high similarity in gene network modulation. Importantly, free fatty acid synthesis modulation and stress pathway responses are affected by active car­boxylic acids, providing coherent mechanistic underpinning for our findings. Our work shows that transcriptomic analysis of cultured human hepatocytes can reinforce the prediction of liver injury outcome based on quantitative and mechanistic biological data and support its application in read-across.

Fluorescent tagging of endogenous Heme oxygenase-1 in human induced pluripotent stem cells for high content imaging of oxidative stress in various differentiated lineages.

Tagging of endogenous stress response genes can provide valuable in vitro models for chemical safety assessment. Here, we present the generation and application of a fluorescent human induced pluripotent stem cell (hiPSC) reporter line for Heme oxygenase-1 (HMOX1), which is considered a sensitive and reliable biomarker for the oxidative stress response. CRISPR/Cas9 technology was used to insert an enhanced green fluorescent protein (eGFP) at the C-terminal end of the endogenous HMOX1 gene. Individual clones were selected and extensively characterized to confirm precise editing and retained stem cell properties. Bardoxolone-methyl (CDDO-Me) induced oxidative stress caused similarly increased expression of both the wild-type and eGFP-tagged HMOX1 at the mRNA and protein level. Fluorescently tagged hiPSC-derived proximal tubule-like, hepatocyte-like, cardiomyocyte-like and neuron-like progenies were treated with CDDO-Me (5.62-1000 nM) or diethyl maleate (5.62-1000 µM) for 24 h and 72 h. Multi-lineage oxidative stress responses were assessed through transcriptomics analysis, and HMOX1-eGFP reporter expression was carefully monitored using live-cell confocal imaging. We found that eGFP intensity increased in a dose-dependent manner with dynamics varying amongst lineages and stressors. Point of departure modelling further captured the specific lineage sensitivities towards oxidative stress. We anticipate that the newly developed HMOX1 hiPSC reporter will become a valuable tool in understanding and quantifying critical target organ cell-specific oxidative stress responses induced by (newly developed) chemical entities.

New approach methods (NAMs) supporting read-across: Two neurotoxicity AOP-based IATA case studies

Read-across approaches are considered key in moving away from in vivo animal testing towards addressing data-gaps using new approach methods (NAMs). Ample successful examples are still required to substantiate this strategy. Here we present and discuss the learnings from two OECD IATA endorsed read-across case studies. They involve two classes of pesticides ­ rotenoids and strobilurins ­ each having a defined mode-of-action that is assessed for its neurological hazard by means of an AOP-based testing strategy coupled to toxicokinetic simulations of human tissue concentrations. The endpoint in question is potential mitochondrial respiratory chain mediated neurotoxicity, specifically through inhibition of complex I or III. An AOP linking inhibition of mitochondrial respiratory chain complex I to the degeneration of dopaminergic neurons formed the basis for both cases but was deployed in two different regulatory contexts. The two cases also exemplify several different read-across concepts: analogue versus category approach, consolidated versus putative AOP, positive versus negative prediction (i.e., neurotoxicity versus low potential for neurotoxicity), and structural versus biological similarity. We applied a range of NAMs to explore the toxicodynamic properties of the compounds, e.g., in silico docking as well as in vitro assays and readouts ­ including transcriptomics ­ in various cell systems, all anchored to the relevant AOPs. Interestingly, although some of the data addressing certain elements of the read-across were associated with high uncertainty, their impact on the overall read-across conclusion remained limited. Coupled to the elaborate regulatory review that the two cases underwent, we propose some generic learnings of AOP-based testing strategies supporting read-across.

In Vitro Three-Dimensional Liver Models for Nanomaterial DNA Damage Assessment.

Whilst the liver possesses the ability to repair and restore sections of damaged tissue following acute injury, prolonged exposure to engineered nanomaterials (ENM) may induce repetitive injury leading to chronic liver disease. Screening ENM cytotoxicity using 3D liver models has recently been performed, but a significant challenge has been the application of such in vitro models for evaluating ENM associated genotoxicity; a vital component of regulatory human health risk assessment. This review considers the benefits, limitations, and adaptations of specific in vitro approaches to assess DNA damage in the liver, whilst identifying critical advancements required to support a multitude of biochemical endpoints, focusing on nano(geno)toxicology (e.g., secondary genotoxicity, DNA damage, and repair following prolonged or repeated exposures).

NAM-supported read-across: From case studies to regulatory guidance in safety assessment.

The use of new approach methodologies (NAMs) in support of read-across (RAx) approaches for regulatory purposes is a main goal of the EU-ToxRisk project. To bring this forward, EU-ToxRisk partners convened a workshop in close collaboration with regulatory representatives from key organizations including European regulatory agencies, such as the European Chemicals Agency (ECHA) and the European Food Safety Authority (EFSA), as well as the Scientific Committee on Consumer Safety (SCCS), national agencies from several European countries, Japan, Canada and the USA, as well as the Organisation for Economic Cooperation and Development (OECD). More than a hundred people actively participated in the discussions, bringing together diverse viewpoints across academia, regulators and industry. The discussion was organized starting from five practical cases of RAx applied to specific problems that offered the oppor-tunity to consider real examples. There was general consensus that NAMs can improve confidence in RAx, in particular in defining category boundaries as well as characterizing the similarities/dissimilarities between source and target substances. In addition to describing dynamics, NAMs can be helpful in terms of kinetics and metabolism that may play an important role in the demonstration of similarity or dissimilarity among the members of a category. NAMs were also noted as effective in providing quanti-tative data correlated with traditional no observed adverse effect levels (NOAELs) used in risk assessment, while reducing the uncertainty on the final conclusion. An interesting point of view was the advice on calibrating the number of new tests that should be carefully selected, avoiding the allure of "the more, the better". Unfortunately, yet unsurprisingly, there was no single approach befitting every case, requiring careful analysis delineating the optimal approach. Expert analysis and assessment of each specific case is still an important step in the process.

Modeling and Simulation Tools: From Systems Biology to Systems Medicine.

Modeling is an integral component of modern biology. In this chapter we look into the role of the model, as it pertains to Systems Medicine, and the software that is required to instantiate and run it. We do this by comparing the development, implementation, and characteristics of tools that have been developed to work with two divergent methodologies: Systems Biology and Pharmacometrics. From the Systems Biology perspective we consider the concept of "Software as a Medical Device" and what this may imply for the migration of research-oriented, simulation software into the domain of human health.In our second perspective, we see how in practice hundreds of computational tools already accompany drug discovery and development at every stage of the process. Standardized exchange formats are required to streamline the model exchange between tools, which would minimize translation errors and reduce the required time. With the emergence, almost 15 years ago, of the SBML standard, a large part of the domain of interest is already covered and models can be shared and passed from software to software without recoding them. Until recently the last stage of the process, the pharmacometric analysis used in clinical studies carried out on subject populations, lacked such an exchange medium. We describe a new emerging exchange format in Pharmacometrics which covers the non-linear mixed effects models, the standard statistical model type used in this area. By interfacing these two formats the entire domain can be covered by complementary standards and subsequently the according tools.

Multiplex Eukaryotic Transcription (In)activation: Timing, Bursting and Cycling of a Ratchet Clock Mechanism.

Activation of eukaryotic transcription is an intricate process that relies on a multitude of regulatory proteins forming complexes on chromatin. Chromatin modifications appear to play a guiding role in protein-complex assembly on chromatin. Together, these processes give rise to stochastic, often bursting, transcriptional activity. Here we present a model of eukaryotic transcription that aims to integrate those mechanisms. We use stochastic and ordinary-differential-equation modeling frameworks to examine various possible mechanisms of gene regulation by multiple transcription factors. We find that the assembly of large transcription factor complexes on chromatin via equilibrium-binding mechanisms is highly inefficient and insensitive to concentration changes of single regulatory proteins. An alternative model that lacks these limitations is a cyclic ratchet mechanism. In this mechanism, small protein complexes assemble sequentially on the promoter. Chromatin modifications mark the completion of a protein complex assembly, and sensitize the local chromatin for the assembly of the next protein complex. In this manner, a strict order of protein complex assemblies is attained. Even though the individual assembly steps are highly stochastic in duration, a sequence of them gives rise to a remarkable precision of the transcription cycle duration. This mechanism explains how transcription activation cycles, lasting for tens of minutes, derive from regulatory proteins residing on chromatin for only tens of seconds. Transcriptional bursts are an inherent feature of such transcription activation cycles. Bursting transcription can cause individual cells to remain in synchrony transiently, offering an explanation of transcriptional cycling as observed in cell populations, both on promoter chromatin status and mRNA levels.

Design principles of nuclear receptor signaling: how complex networking improves signal transduction.

The topology of nuclear receptor (NR) signaling is captured in a systems biological graphical notation. This enables us to identify a number of 'design' aspects of the topology of these networks that might appear unnecessarily complex or even functionally paradoxical. In realistic kinetic models of increasing complexity, calculations show how these features correspond to potentially important design principles, e.g.: (i) cytosolic 'nuclear' receptor may shuttle signal molecules to the nucleus, (ii) the active export of NRs may ensure that there is sufficient receptor protein to capture ligand at the cytoplasmic membrane, (iii) a three conveyor belts design dissipating GTP-free energy, greatly aids response, (iv) the active export of importins may prevent sequestration of NRs by importins in the nucleus and (v) the unspecific nature of the nuclear pore may ensure signal-flux robustness. In addition, the models developed are suitable for implementation in specific cases of NR-mediated signaling, to predict individual receptor functions and differential sensitivity toward physiological and pharmacological ligands.

Stochastic and reversible assembly of a multiprotein DNA repair complex ensures accurate target site recognition and efficient repair.

To understand how multiprotein complexes assemble and function on chromatin, we combined quantitative analysis of the mammalian nucleotide excision DNA repair (NER) machinery in living cells with computational modeling. We found that individual NER components exchange within tens of seconds between the bound state in repair complexes and the diffusive state in the nucleoplasm, whereas their net accumulation at repair sites evolves over several hours. Based on these in vivo data, we developed a predictive kinetic model for the assembly and function of repair complexes. DNA repair is orchestrated by the interplay of reversible protein-binding events and progressive enzymatic modifications of the chromatin substrate. We demonstrate that faithful recognition of DNA lesions is time consuming, whereas subsequently, repair complexes form rapidly through random and reversible assembly of NER proteins. Our kinetic analysis of the NER system reveals a fundamental conflict between specificity and efficiency of chromatin-associated protein machineries and shows how a trade off is negotiated through reversibility of protein binding.

Population-level transcription cycles derive from stochastic timing of single-cell transcription.

Eukaryotic transcription is a dynamic process relying on a large number of proteins. By measuring the cycling expression of the pyruvate dehydrogenase kinase 4 gene in human cells, we constructed a detailed stochastic model for single-gene transcription at the molecular level using realistic kinetics for diffusion and protein complex dynamics. We observed that gene induction caused an approximate 60 min periodicity of several transcription related processes: first, the covalent histone modifications and presence of many regulatory proteins at the transcription start site; second, RNA polymerase II activity; third, chromatin loop formation; and fourth, mRNA accumulation. Our model can predict the precise timing of single-gene activity leading to transcriptional cycling on the cell population level when we take into account the sequential and irreversible multistep nature of transcriptional initiation. We propose that the cyclic nature of population gene expression is primarily based on the intrinsic periodicity of the transcription process itself.

Systems biology towards life in silico: mathematics of the control of living cells.

Systems Biology is the science that aims to understand how biological function absent from macromolecules in isolation, arises when they are components of their system. Dedicated to the memory of Reinhart Heinrich, this paper discusses the origin and evolution of the new part of systems biology that relates to metabolic and signal-transduction pathways and extends mathematical biology so as to address postgenomic experimental reality. Various approaches to modeling the dynamics generated by metabolic and signal-transduction pathways are compared. The silicon cell approach aims to describe the intracellular network of interest precisely, by numerically integrating the precise rate equations that characterize the ways macromolecules' interact with each other. The non-equilibrium thermodynamic or 'lin-log' approach approximates the enzyme rate equations in terms of linear functions of the logarithms of the concentrations. Biochemical Systems Analysis approximates in terms of power laws. Importantly all these approaches link system behavior to molecular interaction properties. The latter two do this less precisely but enable analytical solutions. By limiting the questions asked, to optimal flux patterns, or to control of fluxes and concentrations around the (patho)physiological state, Flux Balance Analysis and Metabolic/Hierarchical Control Analysis again enable analytical solutions. Both the silicon cell approach and Metabolic/Hierarchical Control Analysis are able to highlight where and how system function derives from molecular interactions. The latter approach has also discovered a set of fundamental principles underlying the control of biological systems. The new law that relates concentration control to control by time is illustrated for an important signal transduction pathway, i.e. nuclear hormone receptor signaling such as relevant to bone formation. It is envisaged that there is much more Mathematical Biology to be discovered in the area between molecules and Life.

Mathematical modeling of nucleotide excision repair reveals efficiency of sequential assembly strategies.

Nucleotide excision repair (NER) requires the concerted action of many different proteins that assemble at sites of damaged DNA in a sequential fashion. We have constructed a mathematical model delineating hallmarks and general characteristics for NER. We measured the assembly kinetics of the putative damage-recognition factor XPC-HR23B at sites of DNA damage in the nuclei of living cells. These and other in vivo kinetic data allowed us to scrutinize the dynamic behavior of the nucleotide excision repair process in detail. A sequential assembly mechanism appears remarkably advantageous in terms of repair efficiency. Alternative mechanisms for repairosome formation, including random assembly and preassembly, can readily become kinetically unfavorable. Based on the model, new experiments can be defined to gain further insight into this complex process and to critically test model predictions. Our work provides a kinetic framework for NER and rationalizes why many multiprotein processes within the cell nucleus show sequential assembly strategy.

In vivo dynamics of chromatin-associated complex formation in mammalian nucleotide excision repair.

Chromatin is the substrate for many processes in the cell nucleus, including transcription, replication, and various DNA repair systems, all of which require the formation of multiprotein machineries on the chromatin fiber. We have analyzed the kinetics of in vivo assembly of the protein complex that is responsible for nucleotide excision repair (NER) in mammalian cells. Assembly is initiated by UV irradiation of a small area of the cell nucleus, after which the accumulation of GFP-tagged NER proteins in the DNA-damaged area is measured, reflecting the establishment of the dual-incision complex. The dynamic behavior of two NER proteins, ERCC1-XPF and TFIIH, was studied in detail. Results show that the repair complex is assembled with a rate of approximately 30 complexes per second and is not diffusion limited. Furthermore, we provide in vivo evidence that not only binding of TFIIH, but also its helicase activity, is required for the recruitment of ERCC1-XPF. These studies give quantitative insight into the de novo assembly of a chromatin-associated protein complex in living cells.

Xeroderma pigmentosum group A protein loads as a separate factor onto DNA lesions.

Nucleotide excision repair (NER) is the main DNA repair pathway in mammals for removal of UV-induced lesions. NER involves the concerted action of more than 25 polypeptides in a coordinated fashion. The xeroderma pigmentosum group A protein (XPA) has been suggested to function as a central organizer and damage verifier in NER. How XPA reaches DNA lesions and how the protein is distributed in time and space in living cells are unknown. Here we studied XPA in vivo by using a cell line stably expressing physiological levels of functional XPA fused to green fluorescent protein and by applying quantitative fluorescence microscopy. The majority of XPA moves rapidly through the nucleoplasm with a diffusion rate different from those of other NER factors tested, arguing against a preassembled XPA-containing NER complex. DNA damage induced a transient ( approximately 5-min) immobilization of maximally 30% of XPA. Immobilization depends on XPC, indicating that XPA is not the initial lesion recognition protein in vivo. Moreover, loading of replication protein A on NER lesions was not dependent on XPA. Thus, XPA participates in NER by incorporation of free diffusing molecules in XPC-dependent NER-DNA complexes. This study supports a model for a rapid consecutive assembly of free NER factors, and a relatively slow simultaneous disassembly, after repair.

Large-scale chromatin organization and the localization of proteins involved in gene expression in human cells.

Compartmentalization of the interphase nucleus is an important element in the regulation of gene expression. Here we investigated the functional organization of the interphase nucleus of HeLa cells and primary human fibroblasts. The spatial distribution of proteins involved in transcription (TFIIH and RNA polymerase II) and RNA processing and packaging (hnRNP-U) were analyzed in relation to chromosome territories and large-scale chromatin organization. We present evidence that these proteins are present predominantly in the interchromatin space, inside and between chromosome territories, and are largely excluded by domains of condensed chromatin. We show that they are present throughout the active and inactive X-chromosome territories in primary female fibroblasts, indicating that these proteins can freely diffuse throughout the interchromatin compartment in the interphase nucleus. Furthermore, we established that the in vivo spatial distribution of condensed chromatin in the interphase nucleus does not depend on ongoing transcription. Our data support a conceptually simple model for the functional organization of interphase nuclei.