File Name : figs1.tif Caption : figure s1: aso uptake determined based on measured a488 fluorescence intensity from different pores of the scaffold that were treated with 2 µm galnac-aso (a, b) and non-targeted aso (c, d) at 1 h time point, respectively. fov – field of view. scale bar = 50 μm. File Name : figs2.tif Caption : figure s2: cellular phenotypes observed in the liver-on-a-chip scaffold. cells with circular and cuboidal morphologies were observed in the pores of a liver-on-a-chip scaffold. cellular and sub-cellular morphological features were visible from the nad(p)h autofluorescence signal, which was captured in the second detection channel of the custom-built multiphoton imaging system. multiphoton images of microtissue formations obtained from a depth of 30 μm inside the pores of a vehicle control (pbs only) liver-on-a-chip scaffold obtained at 60x magnification. arrow points toward the cuboidal hepatocyte and * denotes circular cells. scale bar: 25 μm. File Name : figs3.tif Caption : figure s3: a representative image obtained from a liver-on-a-chip scaffold that was treated with a488-aso-galnac 1 h post-treatment. a) multiphoton image showing the fluorescence intensity detected in the a488 channel. b) corresponding annotated image showing the manually segmented cuboidal hepatocytes. scale bar = 50 μm. File Name : figs4.tif Caption : figure s4: cellpose algorithm and manual exclusion of the central pore region (to avoid new top layer cells appearing towards the pore center with the increased depth of focus) was used to generate masks quantification of signal loss with the sample depth. the depth information was transformed into a color-coded image using the maximum intensity projection (fiji imagej) along z axis was used to flatten images for both a) 1pf and b) 2pf datasets. b, e) pore center-excluding masks were manually created (dashed white circle). c, f) nuclei masks were generated with a nuclei model trained with >250 nuclei with 1pf and 2pf data, respectively. scale bars = 50 µm. g, h) masks were used to capture and average the fluorescence intensity of each nuclear object per z plane (i.e., 1 dot = 1 nucleus per plane). the intensity scales are relative to the maximum intensity in each 3d dataset. File Name : figs5.tif Caption : figure s5: sample orientation in lsfm. File Name : figs6.tif Caption : figure s6: validation of a) cona and b) hoechst pfa-fixed in untreated hepg2 spheroids. scale bar = 20 µm File Name : figs7.tif Caption : figure s7: a) cona and b) hoechst, c) aso distribution in a488-aso-treated (2 µm, 1 h) hepg2 spheroids. staining was conducted after fixation. no permeabilization was done. scale bar = 20 µm