Use of a systems model of drug-induced liver injury (DILIsym®) to elucidate the mechanistic differences between acetaminophen and its less-toxic isomer, AMAP, in mice

AbstractAcetaminophen (APAP) has been used as a probe drug to investigate drug-induced liver injury (DILI). In mice, 3′-hydroxyacetanilide (AMAP), a less-toxic isomer of APAP, has also been studied as a negative control. Various mechanisms for the divergence in toxicological response between the two isomers have been proposed. This work utilized a mechanistic, mathematical model of DILI to test the plausibility of four mechanistic hypotheses. Simulation results were compared to an array of measured endpoints in mice treated with APAP or AMAP. The four hypotheses included: (1) quantitative differences in drug metabolism profiles as a result of different affinities for the relevant enzymes; (2) differences in the amount of reactive metabolites produced due to cytochrome P450 (CYP450) inhibition by the AMAP reactive metabolites; (3) differences in the rate of conjugation between the reactive metabolites and proteins; (4) differences in the downstream effects or potencies of the reactive metabolites on vital components within hepatocytes. The simulations did not support hypotheses 3 or 4 as the most likely hypotheses underlying the difference in hepatoxic potential of APAP and AMAP. Rather, the simulations supported hypotheses 1 and 2 (less reactive metabolite produced per mole of AMAP relative to APAP). Within the simulations, the difference in reactive metabolite formation was equally likely to have occurred from differential affinities for the relevant drug metabolism enzymes or from direct CYP450 inhibition by the AMAP reactive metabolite. The demonstrated method of using simulation tools to probe the importance of possible contributors to toxicological observations is generally applicable across species

Investigating the exposure and impacts of chemical UV-filters in coral reef ecosystems : Review and research gap prioritisation

Coral reefs are amongst the world's most productive and biologically diverse ecosystems and in recent decades have experienced an unparalleled decline due to various anthropogenically induced stressors. Ultraviolet (UV) filters found within personal care products such as sunscreen are an emerging chemical pollutant with a growing concern of toxicity to reef organisms. In this study a systematic literature review was conducted to [1] determine the current understanding of spatial distribution and occurrence of UV-filters exposed to the marine environment, [2] synthesise current ecotoxicological thresholds of relevant reef organisms under various UV filter exposures, [3] identify research gaps related to both exposure and toxicity of UV-filters in coral reef ecosystems. With gaps identified, a survey was developed and distributed to experts in the field representing academic, governmental, not-for-profit, and industry researchers in order to prioritise research gaps and inform future research efforts. The survey identified the need for better understanding of the impacts of co-stressors, long-term exposure, mixture and degradation product exposure and realistic environmental conditions. Ultimately, this review, will help guide priority research efforts to understand the risks of UV-filter exposure to coral reef ecosystems. This article is protected by copyright. All rights reserved

Shockwave/Boundary-Layer Interaction Studies Performed in the NASA Langley 20-Inch Mach 6 Air Tunnel

This paper highlights results from a collaborative study performed by The University of Tennessee Space Institute (UTSI) and NASA Langley Research Center on the Shockwave/Boundary-Layer Interaction (SWBLI) generated by a cylindrical protuberance on a flat plate in a Mach 6 flow. The study was performed in the 20-Inch Mach 6 Air Tunnel at NASA Langley Research Center and consisted of two separate entries. In the first entry, simultaneous high-speed schlieren and high-speed pressure-sensitive paint (PSP) imaging which was performed for the first time in the 20-Inch Mach 6 facility at NASA Langley were performed as well as simultaneous high-speed schlieren and oil-flow imaging. In the second entry, the model configuration was modified to increase the size of the interaction region. High-speed schlieren and infrared thermography (IR) surface imaging were performed in this second entry. The goal of these tests was to characterize the SBLI in the presence of a laminar, transitional, and turbulent boundary layer using high-speed optical imaging techniques. AoA = sting angle-of-attack () dcylinder = cylinder diameter (mm) dtrip = cylindrical tripping element diameter (mm) shock = shock stand-off distance (mm) hcylinder = cylinder height (mm) htrip = cylindrical tripping element height (mm) HSS = high-speed schlieren M = freestream Mach number PSP = pressure-sensitive paint Re = freestream unit Reynolds number (m-1) SWBLI = shockwave/boundary-layer interaction plate = model plate angle () Introduction his paper highlights two experimental entries performed in the 20-Inch Mach 6 Air Blowdown Tunnel at NASA Langley Research Center in collaboration with The University of Tennessee Space Institute (UTSI). The purpose of these entries was to characterize the dynamic shockwave/boundary-layer interaction (SWBLI) between a vertical cylinder on a flat plate and laminar, transitional (XSWBLI), and turbulent (SWTBLI) boundary layers with a freestream Mach number of 6 using non-intrusive optical diagnostics. Experiments performed by Murphree et al.1,2 were among the first to specifically characterize XSWBLI induced by a vertical cylinder on a flat plate geometry using several optical measurement techniques. Recent optical studies of XSWBLI phenomenon have been performed by UTSI at Mach 2 in their low-enthalpy blow wind tunnel3-8 and by Texas A&M University and UTSI at Mach numbers of 6 and 7 in their Adjustable Contour Expansion wind tunnel.9 The experiments described in this paper were intended to complement previous studies by expanding the freestream unit Reynolds number range, Re, over which the XSWBLI phenomena has been observed. Additionally these experiments, made possible under NASAs new facility funding model under the Aeronautics Evaluation and Test Capabilities (AETC) project, promoted collaboration between university and NASA researchers. The initial entry in the 20-Inch Mach 6 Air Tunnel at NASA Langley occurred in December of 2016. Originally, testing was to occur in November of 2016 in the 31-Inch Mach 10 Air Tunnel at NASA Langley. This facility was chosen so that the XSWBLI phenomenon could be observed at much higher Mach numbers than had previously been attempted in ground test experiments. The model selected for this experiment, a 10 half-angle wedge with a sharp leading edge (described in detail in section II.B), had previously been used by Danehy et al. [10] for boundary layer transition studies using the nitric oxide planar laser-induced fluorescence (NO PLIF) flow visualization technique. In that work, it was determined that transition could be induced downstream of a single htrip = 1-mm tall, dtrip = 4-mm diameter cylindrical tripping element and that the streamwise location of the transition could be changed for a single Re by changing the model angle-of-attack (AoA) (see Fig. A3 in Ref. [10] for more details). Based on the findings of that work, a decision was made to use the wedge model with the cylindrical tripping element to trip the boundary layer flow ahead of a cylindrical protuberance in order to achieve a XSWBLI. Unfortunately, the 31-Inch Mach 10 facility had been taken offline for repairs in October of 2016 and a decision was made to move the test to the 20-Inch Mach 6 facility. Since the behavior of the boundary layer with the chosen model configuration had not been studied before in that facility and the available test time was limited, the entry was considered to be exploratory and was used to collect spatially-resolved and time-resolved flow and surface visualization data that would be used to inform a second entry. Test techniques included simultaneous high-speed schlieren (HSS) captured at 160 kHz and high-speed pressure sensitive paint captured at 10 kHz as well as oil flow visualization, captured at 750 Hz. The second entry in the 20-Inch Mach 6 facility occurred in June and July of 2017. In this follow-on test, modifications to the wind tunnel model were made based on observations made during the first entry and included removing the cylindrical tripping element, increasing the size of the cylinder used to induce the SWBLI to increase the size of the interaction while simultaneously improving spatial resolution, and using a swept ramp array, similar to that described in Ref. [11], to trip the flow to turbulence. Simultaneous HSS (captured at 140 kHz, 100 kHz, and 40 kHz) and conventional IR thermography (captured at 30 Hz) imaging were performed simultaneously in this follow-on entry. This paper is intended to serve as a summary of the work performed during these two entries, to detail lessons learned from each entry, and to highlight some of the datasets acquired. Details on the experimental setup, model configuration, and techniques used are provided. Papers providing a more rigorous analysis of data acquired during the second entry, including statistical, spectral, and modal decomposition methods, can be found in Refs. [12,13]. An entry examining XSWBLI in the 31-Inch Mach 10 Blowdown Wind Tunnel facility is currently planned for mid-to-late calendar year 2019, pending the success of facility repairs. The work performed and described in this paper and the upcoming entry in the 31-Inch Mach 10 facility at NASA Langley have been made possible by NASAs new facility funding model under the Aeronautics Evaluation and Test Capabilities (AETC) project. Wind Tunnel Facility All experiments discussed in this paper were performed in the 20-Inch Mach 6 Air Tunnel at NASA Langley Research Center. Specific details pertaining to this facility can be found in Refs. [14,15], with only a brief description of the facility provided here. For both entries, the nominal freestream unit Reynolds number was varied between 1.8106 m-1 (0.5106 ft-1) and 26.3106 m-1 (8106 ft-1). The nominal stagnation pressure was varied between 0.21 MPa and 3.33 MPa and the nominal stagnation temperature was varied between 480 K and 520 K to achieve the desired Re condition. For all runs, the nominal freestream Mach number was 6. The nearly square test section is 520.7-mm (20.5-inches) wide by 508-mm (20-inches) high. Two 431.8-mm (17-inch) diameter windows made of Corning 7940, Grade 5F schlieren-quality glass serve as the side walls of the tunnel and provide optical access for the high-speed schlieren measurements. A rectangular window made of the same material as the side windows served as the top wall of the test section and provided optical access for the high-speed PSP and oil flow measurements. For the second entry, this top window was replaced with a Zinc Selenide (ZnSe) window with an anti-reflection coating capable of passing IR wavelengths between 8m and 12m with greater than 98% transmittance. The model was sting supported by a strut attached to a hydraulic system that allows for the model pitch angle to be adjusted between -5 to +55. For the first entry, an initial pitch/pause sweep of the model AoA was performed to observe the resulting SWBLI. Ultimately, however, the sting pitch angle for this entry was fixed at +10.0 so that the angle of the top surface of the wedge relative to the streamwise axis of the tunnel (referred to herein as the plate angle, plate), was plate = 0. For the second entry, plate = 0 and plate = -13.25 were initially tested with the swept ramp array (discussed in the following section) to determine which orientation produced conditions most favorable for XSWBLI to occur based on the heating signatures observed over the top surface of the model in the IR thermography images. Based on these initial tests, plate = -13.25 was set for the remainder of the runs in the second entry. For both entries, any model changes were performed in a housing located beneath the closed test section. Prior to performing a run of the tunnel, the housing was sealed and the tunnel started. Once the appropriate freestream conditions were achieved, the model was injected into the test section using a hydraulic injection system. B. Model Geometry For all runs, a 10 half-angle (20 full-angle) wedge model with a sharp leading edge was used. The model is described in detail in Refs. [10,16]. The top surface of the sharp leading edge of the model extended 47.8 mm from its upstream-most edge to a junction with the upstream edge of a stainless steel top plate that then extended an (a) (c) (b) Fig. 1 (a) Schematic of top surface of wedge model with gas seeding insert, (b) perspective view of the model in the 20-Inch Mach 6 tunnel with centerline pressure orifices on sharp leading edge, and (c) a perspective view of the model with stainless steel (top) and SLA middle insert (bottom) during the first entry. Flow occurs from left to right

Evaluation of an ambulatory system for the quantification of cough frequency in patients with chronic obstructive pulmonary disease

BACKGROUND: To date, methods used to assess cough have been primarily subjective, and only broadly reflect the impact of chronic cough and/or chronic cough therapies on quality of life. Objective assessment of cough has been attempted, but early techniques were neither ambulatory nor feasible for long-term data collection. We evaluated a novel ambulatory cardio-respiratory monitoring system with an integrated unidirectional, contact microphone, and report here the results from a study of patients with COPD who were videotaped in a quasi-controlled environment for 24 continuous hours while wearing the ambulatory system. METHODS: Eight patients with a documented history of COPD with ten or more years of smoking (6 women; age 57.4 ± 11.8 yrs.; percent predicted FEV(1 )49.6 ± 13.7%) who complained of cough were evaluated in a clinical research unit equipped with video and sound recording capabilities. All patients wore the LifeShirt(® )system (LS) while undergoing simultaneous video (with sound) surveillance. Video data were visually inspected and annotated to indicate all cough events. Raw physiologic data records were visually inspected by technicians who remained blinded to the video data. Cough events from LS were analyzed quantitatively with a specialized software algorithm to identify cough. The output of the software algorithm was compared to video records on a per event basis in order to determine the validity of the LS system to detect cough in COPD patients. RESULTS: Video surveillance identified a total of 3,645 coughs, while LS identified 3,363 coughs during the same period. The median cough rate per patient was 21.3 coughs·hr(-1 )(Range: 10.1 cghs·hr(-1 )– 59.9 cghs·hr(-1)). The overall accuracy of the LS system was 99.0%. Overall sensitivity and specificity of LS, when compared to video surveillance, were 0.781 and 0.996, respectively, while positive- and negative-predictive values were 0.846 and 0.994. There was very good agreement between the LS system and video (kappa = 0.807). CONCLUSION: The LS system demonstrated a high level of accuracy and agreement when compared to video surveillance for the measurement of cough in patients with COPD

Sandwich-Cultured Hepatocytes as a Tool to Study Drug Disposition and Drug-Induced Liver Injury

Sandwich-cultured hepatocytes (SCH) are metabolically competent and have proper localization of basolateral and canalicular transporters with functional bile networks. Therefore, this cellular model is a unique tool that can be used to estimate biliary excretion of compounds. SCH have been used widely to assess hepatobiliary disposition of endogenous and exogenous compounds and metabolites. Mechanistic modeling based on SCH data enables estimation of metabolic and transporter-mediated clearances, which can be employed to construct physiologically-based pharmacokinetic models for prediction of drug disposition and drug-drug interactions in humans. In addition to pharmacokinetic studies, SCH also have been employed to study cytotoxicity and perturbation of biological processes by drugs and hepatically-generated metabolites. Human SCH can provide mechanistic insights underlying clinical drug-induced liver injury (DILI). In addition, data generated in SCH can be integrated into systems pharmacology models to predict potential DILI in humans. In this review, applications of SCH in studying hepatobiliary drug disposition and bile acid-mediated DILI are discussed. An example is presented to show how data generated in the SCH model was used to establish a quantitative relationship between intracellular bile acids and cytotoxicity, and how this information was incorporated into a systems pharmacology model for DILI prediction

Hypersonic Laminar Boundary Layer Velocimetry with Discrete Roughness on a Flat Plate

Laminar boundary layer velocity measurements are made on a 10-degree half-angle wedge in a Mach 10 flow. Two types of discrete boundary layer trips were used to perturb the boundary layer gas. The first was a 2-mm tall, 4-mm diameter cylindrical trip. The second was a scaled version of the Orbiter Boundary Layer Transition (BLT) Detailed Test Objective (DTO) trip. Both 1-mm and 2.5-mm tall BLT DTO trips were tested. Additionally, side-view and plan-view axial boundary layer velocity measurements were made in the absence of these tripping devices. The free-stream unit Reynolds numbers tested for the cylindrical trips were 1.7x10(exp 6)/m and 3.3x10(exp 6)/m. The free-stream unit Reynolds number tested for the BLT DTO trips was 1.7x10(exp 6)/m. The angle of attack was kept at approximately 5-degrees for most of the tests resulting in a Mach number of approximately 8.3. These combinations of unit Reynolds numbers and angle of attack resulted in laminar flowfields. To study the precision of the measurement technique, the angle of attack was varied during one run. Nitric-oxide (NO) molecular tagging velocimetry (MTV) was used to obtain averaged axial velocity values and associated uncertainties. These uncertainties are as low as 20 m/s. An interline, progressive scan CCD camera was used to obtain separate images of the initial reference and shifted NO molecules that had been tagged by the laser. The CCD configuration allowed for sub-microsecond sequential acquisition of both images. The maximum planar spatial resolution achieved for the side-view velocity measurements was 0.07-mm in the wall-normal direction by 1.45-mm in the streamwise direction with a spatial depth of 0.5-mm. For the plan-view measurements, the maximum planar spatial resolution in the spanwise and streamwise directions was 0.69-mm by 1.28-mm, respectively, with a spatial depth of 0.5-mm. Temperature sensitive paint (TSP) measurements are provided to compliment the velocity data and to provide further insight into the behavior of the boundary layers. The experiments were performed at the NASA Langley Research Center 31-Inch Mach 10 Air tunnel

Substrate complexity reduces prey consumption in functional response experiments: Implications for extrapolating to the wild

Understanding the density-dependent impacts of an invasive predator is integral for predicting potential consequences for prey populations. Functional response experiments are used to assess the rate of prey consumption and a predator’s ability to search for and consume prey at different resource densities. However, results can be highly context-dependent, limiting their extrapolation to natural ecosystems. Here, we examined how simulated habitat complexity, through the addition of substrate in which prey can escape predation, affects the functional response of invasive European green crabs (Carcinus maenas) foraging on two different bivalve species. Green crabs feeding on varnish clams (Nuttallia obscurata) shifted from a Type II hyperbolic functional response in the absence of substrate to density-independent consumption when prey could bury. Green crabs ate few Japanese littleneck clams (Venerupis philippinarum) under all densities, such that no functional response curve of any type could be produced and their total consumption was always density independent. However, the probability of at least one Japanese littleneck clam being consumed increased significantly with initial clam density and crab claw size across all treatments. At mean crab claw size and compared to trials without substrate, the proportion of varnish clams consumed were 4.2 times smaller when substrate was present, but substrate had a negligible effect (1.2 times) on Japanese littlenecks. The proportion of varnish clams consumed increased with crab claw size and were higher across both substrate conditions than the proportion of Japanese littlenecks consumed; however, the proportion of Japanese littleneck clams consumed increased faster with claw size than that of varnish clams. Our results suggest that including environmental features and variation in prey species can influence the density-dependent foraging described by functional response experiments. Incorporating replicable features of the natural environment into functional response experiments is imperative to make more accurate predictions about the impact of invasive predators on prey populations

Application of a Mechanistic Model to Evaluate Putative Mechanisms of Tolvaptan Drug-Induced Liver Injury and Identify Patient Susceptibility Factors

Tolvaptan is a selective vasopressin V2 receptor antagonist, approved in several countries for the treatment of hyponatremia and autosomal dominant polycystic kidney disease (ADPKD). No liver injury has been observed with tolvaptan treatment in healthy subjects and in non-ADPKD indications, but ADPKD clinical trials showed evidence of drug-induced liver injury (DILI). Although all DILI events resolved, additional monitoring in tolvaptan-treated ADPKD patients is required. In vitro assays identified alterations in bile acid disposition and inhibition of mitochondrial respiration as potential mechanisms underlying tolvaptan hepatotoxicity. This report details the application of DILIsym software to determine whether these mechanisms could account for the liver safety profile of tolvaptan observed in ADPKD clinical trials. DILIsym simulations included physiologically based pharmacokinetic estimates of hepatic exposure for tolvaptan and2 metabolites, and their effects on hepatocyte bile acid transporters and mitochondrial respiration. The frequency of predicted alanine aminotransferase (ALT) elevations, following simulated 90/30 mg split daily dosing, was 7.9% compared with clinical observations of 4.4% in ADPKD trials. Toxicity was multifactorial as inhibition of bile acid transporters and mitochondrial respiration contributed to the simulated DILI. Furthermore, simulation analysis identified both pre-treatment risk factors and on-treatment biomarkers predictive of simulated DILI. The simulations demonstrated that in vivo hepatic exposure to tolvaptan and the DM-4103 metabolite, combined with these 2 mechanisms of toxicity, were sufficient to account for the initiation of tolvaptan-mediated DILI. Identification of putative risk-factors and potential novel biomarkers provided insight for the development of mechanism-based tolvaptan risk-mitigation strategies

Exploring BSEP inhibition-mediated toxicity with a mechanistic model of drug-induced liver injury

Inhibition of the bile salt export pump (BSEP) has been linked to incidence of drug-induced liver injury (DILI), presumably by the accumulation of toxic bile acids in the liver. We have previously constructed and validated a model of bile acid disposition within DILIsym®, a mechanistic model of DILI. In this paper, we use DILIsym® to simulate the DILI response of the hepatotoxic BSEP inhibitors bosentan and CP-724,714 and the non-hepatotoxic BSEP inhibitor telmisartan in humans in order to explore whether we can predict that hepatotoxic BSEP inhibitors can cause bile acid accumulation to reach toxic levels. We also simulate bosentan in rats in order to illuminate potential reasons behind the lack of toxicity in rats compared to the toxicity observed in humans. DILIsym® predicts that bosentan, but not telmisartan, will cause mild hepatocellular ATP decline and serum ALT elevation in a simulated population of humans. The difference in hepatotoxic potential between bosentan and telmisartan is consistent with clinical observations. However, DILIsym® underpredicts the incidence of bosentan toxicity. DILIsym® also predicts that bosentan will not cause toxicity in a simulated population of rats, and that the difference between the response to bosentan in rats and in humans is primarily due to the less toxic bile acid pool in rats. Our simulations also suggest a potential synergistic role for bile acid accumulation and mitochondrial electron transport chain (ETC) inhibition in producing the observed toxicity in CP-724,714, and suggest that CP-724,714 metabolites may also play a role in the observed toxicity. Our work also compares the impact of competitive and noncompetitive BSEP inhibition for CP-724,714 and demonstrates that noncompetitive inhibition leads to much greater bile acid accumulation and potential toxicity. Our research demonstrates the potential for mechanistic modeling to contribute to the understanding of how bile acid transport inhibitors cause DILI

Investigating bile acid-mediated cholestatic drug-induced liver injury using a mechanistic model of multidrug resistance protein 3 (MDR3) inhibition

Inhibition of the canalicular phospholipid floppase multidrug resistance protein 3 (MDR3) has been implicated in cholestatic drug-induced liver injury (DILI), which is clinically characterized by disrupted bile flow and damage to the biliary epithelium. Reduction in phospholipid excretion, as a consequence of MDR3 inhibition, decreases the formation of mixed micelles consisting of bile acids and phospholipids in the bile duct, resulting in a surplus of free bile acids that can damage the bile duct epithelial cells, i.e., cholangiocytes. Cholangiocytes may compensate for biliary increases in bile acid monomers via the cholehepatic shunt pathway or bicarbonate secretion, thereby influencing viability or progression to toxicity. To address the unmet need to predict drug-induced bile duct injury in humans, DILIsym, a quantitative systems toxicology model of DILI, was extended by representing key features of the bile duct, cholangiocyte functionality, bile acid and phospholipid disposition, and cholestatic hepatotoxicity. A virtual, healthy representative subject and population (n = 285) were calibrated and validated utilizing a variety of clinical data. Sensitivity analyses were performed for 1) the cholehepatic shunt pathway, 2) biliary bicarbonate concentrations and 3) modes of MDR3 inhibition. Simulations showed that an increase in shunting may decrease the biliary bile acid burden, but raise the hepatocellular concentrations of bile acids. Elevating the biliary concentration of bicarbonate may decrease bile acid shunting, but increase bile flow rate. In contrast to competitive inhibition, simulations demonstrated that non-competitive and mixed inhibition of MDR3 had a profound impact on phospholipid efflux, elevations in the biliary bile acid-to-phospholipid ratio, cholangiocyte toxicity, and adaptation pathways. The model with its extended bile acid homeostasis representation was furthermore able to predict DILI liability for compounds with previously studied interactions with bile acid transport. The cholestatic liver injury submodel in DILIsym accounts for several processes pertinent to bile duct viability and toxicity and hence, is useful for predictions of MDR3 inhibition-mediated cholestatic DILI in humans