A passive diffusion model of fluorescein derivatives in an in vitro human brain microvascular endothelial cell (HBMEC) monolayer
Jamelle M. Simmons^{1} , Yong Woo Lee^{2} , Luke Achenie^{3}
^{1, 2, 3}Department of Biomedical Engineering and Mechanics, Virginia Polytechnic Institute and State University, 325 Stanger Street, Kelly Hall, Blacksburg, VA 24061, USA
^{3}Department of Chemical Engineering, Virginia Polytechnic Institute and State University, Goodwin Hall, 635 Prices Fork Road, Blacksburg, Virginia 24061, USA
^{1}Corresponding author
Journal of Measurements in Engineering, Vol. 6, Issue 3, 2018, p. 127136.
https://doi.org/10.21595/jme.2018.20186
Received 30 August 2018; accepted 29 September 2018; published 30 September 2018
JVE Conferences
Eukaryotic cells have a protective plasma membrane, which restricts the free movement of molecules from the external environment to the internal environment. This study aims to computationally model the transport of fluorescein derivatives across the monolayer of human brain microvascular endothelial cells (HBMEC). The determination of plausible effective diffusion constants (${D}_{eff}$) will allow models to be built that could be useful beyond in vitro experimentation. Fluorescein5isothiocyanate (FITC) modeling produced a ${D}_{eff}$ range of 1E20 to 5E20 cm^{2}/s at a 1 µm cell monolayer thickness and a ${D}_{eff}$ constant near 5E29 cm^{2}/s at a 5 µm cell monolayer thickness. Both fluorescein and sodium fluorescein (NaFl) modeling at the 1 and 5 µm thicknesses did not produce simulations that closely resembled the HBMEC in vitro model. Overall, it is possible that the fluorescent intensity noted with fluorescein and NaFl may be better explained by a mechanism other than passive diffusion. Simulations of FITC diffusion produced a narrow range of ${D}_{eff}$ constants that closely matched the in vitro HBMEC fluorescent intensity.
Keywords: fluorescein5isothiocyanate, fluorescein, sodium fluorescein, passive diffusion, endothelial cells.
1. Introduction
Eukaryotic cells are enclosed in a plasma membrane which separates the external and internal environments that serves as protection for the cells. The lipid bilayer of the plasma membrane is composed of hydrophilic head groups which are directed towards the aqueous exterior and interior portions of cells and hydrophobic tails separating the hydrophilic groups. The plasma membrane is mostly composed of lipids and proteins that help it form a barrier around the cell’s internal environment [13].
Isothiocyanates (ITCs) are small molecules found in vegetables that are capable of binding to amino groups and proteins and passively diffusion through gastrointestinal epithelium and capillary endothelium [4, 5]. Phenethyl isothiocyanate (PEITC) has been shown to be absorbed into human nonsmall lung carcinoma (A549) cells and ITCs bind to thiol groups at a faster rate compared to amino groups under basic conditions [69]. Research on the cellular accumulation of ITCs have found that the uptake was dependent on intracellular levels of glutathione (GSH) and not dependent on the lipophilicity [1012].
Fluorescein5isothiocyanate (FITC, Isomer I) is a fluorescein fluorochrome conjugated to an isothiocyanate (ITC) that covalently binds to proteins and amine groups. When the pH is raised to an alkaline environment, the covalent bonding is stronger than at physiologic pH [1315]. FITC has been shown to incorporate into the lipid layer and bind to membrane proteins of baby hamster kidney fibroblast (BHK) cells [16, 17] and has been shown to label B and T lymphocytes differently as well bind to internal components of these cells [18, 19]. Labeling of cells using FITC has been shown to be more uniform compared to other lipophilic labels such as 1, 1′dioctadecyl3, 3, 3′, 3′tetramethylindocarbocyanine perchlorate (Dil), but the FITC fluorescence decreases more rapidly compared to labels using tetramethylrhodamine (TRITC) under incubation at 37°C [20]. Additional factors that influence the fluorescence of cells are (A) concentration of fluorescent compound, (B) pH, (C) incubation time, (D) temperature, (E) concentration of cells, and (F) concentration of protein in cell culture medium [18].
Sodium fluorescein (NaFl) is a watersoluble dye that has been used to evaluate the permeability of bloodbrain barrier (BBB) models [21] due to its small size. It has also been used in the study of malignant and normal cells for its internal accumulation due to a pH gradient across the external and internal cell environments [22]. Fluorescein has high binding affinity to intracellular proteins of cells after it is hydrolyzed from a starting material, such as fluorescein diacetate (FDA), which is able to be taken into cells [23, 24] but free diffusion of fluorescein has yet to be proven to our knowledge.
To date, many studies have been conducted on multiple cell types to determine the amount of cell surface attachment and cytoplasmic attachment but work has not been done to determine if fluorescent tracers such as FITC, fluorescein, and NaFl can passively diffuse through the cell surface of human brain microvascular endothelial cells (HBMEC). Brain microvascular endothelial cells help make up the BBB which limits passage across it, to (1) active and passive transport, (2) endocytosis, (3) active and passive diffusion, and (4) paracellular transport, due to tight junctions between endothelial cells. Passive diffusion across cells remains a challenge as it is limited by lipophilicity, size, and charge of the cell and the compound diffusing across it [25, 26].
Additional parameters of interest, cell volume and size fluctuations, of cell monolayers can occur over time that will make an in vitro model more of a dynamic system. Studies conducted on MandinDarby canine kidney (MDCK) cells have volume fluctuations of ± 20% and an average cell thickness of 7.1 ± 0.7 µm (mean ± standard deviation) [27]. Confocal fluorescence microscopy has also been used to measure the cell thickness of C3H10T1/2 and V79 cells grown on various substrates. It was reported that C3H10T1/2 cells, during the exponetial growth phase, grown on mylar have an average cell thickness of 2.9 ± 0.6 µm while V79 cells have an average cell thikness of 6.1 ± 1.0 µm grown on the same substrate during the exponential growth phase [28].
The aim of the present study was to create a simple approach that relies on fluorescence to determine an effective diffusion constant (${D}_{eff}$) of FITC, fluorescein, and NaFl transport across human brain microvascular endothelial cells (HBMEC) by passive diffusion. These constants would allow simulations to be run which could help predict the amount of free diffusion into HBMEC monolayers overtime without loss of fluorescent intensity in signal that can be seen with prolonged in vitro incubation.
2. Materials and methods
2.1. Cell culture
Penicillinstreptomycin solution (100x), MEM vitamin solution (100x), endothelial cell growth supplement (ECGS), fetal bovine serum (FBS), and NuSerum culture supplement were all obtained from Corning Life Sciences. Sodium pyruvate solution (100 mM), RPMI 1640 medium containing 2 mM Lglutamine were purchased from GE Healthcare Life Sciences. Nonessential amino acids were obtained from Gibco laboratories. Human brain microvascular endothelial cells (HBMEC) were isolated and purified as previously described [29]. Cells were cultured in RPMI 1640 medium with 10% fetal bovine serum, 10% NuSerum, 30 µg/ml ECGS, 15 U/ml of heparin, 11.0 mg/ml of sodium pyruvate, 100 U/ml of penicillin, 100 µg/ml of streptomycin, nonessential amino acids, and vitamins. HBMEC cultures were incubated at 37°C in a humid environment with 5% CO_{2}.
2.2. In vitro
Fluorescein isothiocyante isomer I (FITC), fluorescein free acid, and fluorescein sodium salt (NaFl) fluorescent tracers were all purchased through SigmaAldrich. The molecular weights for FITC, fluorescein, and NaFl were 389.38, 332.31, and 376.28 respectively. HBMEC were grown to confluence followed by washing the monolayer with hanks balanced salt solution (HBSS, Corning Life Sciences) and cell detachment with 0.05% trypsin containing EDTA (GE Healthcare Life Sciences). The suspended cells were seeded into Falcon^{TM} 48 well plates at a density of 7.0×10^{4} cells per well and grown to confluence in an incubator for 12 days.
Confluent cells were aspirated of media and replaced with 200 µl of 0.1 mg/ml fluorescent tracer solutions in basal RPMI 1640 medium. Wells were incubated in tracer solutions at 3, 6, 9, 12, 15, 18, 21, and 24 hours with a sample size of 6 wells per time point for all three tracers. All wells were washed twice with sterile phosphate buffered saline (1x) to remove excess fluorescent tracer and replaced with phenol free RPMI 1640 medium (Gibco Laboratories) at a volume of 200 µl. Fluorescent intensity of the wells were recorded using a SpectraMax M2^{e} (Molecular Devices, LLC) plate reader with excitation and emission settings of 492 and 518 nm respectively for FITC, 490 and 514 nm for fluorescein, and 460 and 515 nm for NaFl.
The fluorescent intensity was correlated to concentration (mg/ml) using a dilution curve generated for each fluorescent tracer dissolved in phenol redfree RPMI 1640 for various known concentrations. The fluorescent intensity of tracer found in the HBMEC monolayers were converted to concentration using the dilution curves. The total mass in mg for each monolayer was determined by multiplying the volume the monolayer occupied in each well by the calculated concentration. As cell monolayer thicknesses can vary based on cell type a 1 µm and 5 µm thickness were assumed. The 48 well plates have a growth area of 0.75 cm^{2} and accounting for a 1 µm monolayer thickness results in a volume of 2.36E04 ml being occupied by the monolayer while a 5 µm thickness results in a volume of 1.18E03 ml being occupied by the monolayer. The total mass in mg was averaged and standard deviation calculated for each time point for each tracer (Fig. 1).
Fig. 1. a) Total mass in mg averaged for each time point per fluorescent tracer for a 1 µm thickness; b) total mass in mg averaged for each time point per fluorescent tracer for a 5 µm thickness
a)
b)
2.3. Passive diffusion model
The in vitro system was composed of a cylindrical well where uptake of tracer could occur in the $x$, $y$, and $z$directions (Fig. 2) with the HBMEC. Assuming the transport mechanism is primarily through passive diffusion, the following simplifications were made: (1) no reaction terms, (2) no binding, (3) movement of tracer was uniform in all directions, and (4) an overall diffusion constant is appropriate. These assumptions allow for the reduction of the problem from a three dimensional (3D) space to a one dimensional (1D) space where concentration could be determined in the $z$direction. A 1D partial differential equation (PDE, Eq. (1)) was generated with one initial condition (Eq. (2)) and two boundary conditions (Eqs. (34)). Simulations were run in MATLAB 2015b.
Fig. 2. a) In vitro cell culture of HBMEC where the monolayer thickness spans the $z$direction; b) introduction of fluorescent tracers to the HBMEC for incubation. Uptake of the tracer is assumed to be uniform in all directions
a)
b)
The boundary concentration was set at 0.1 mg/ml for all simulations with the monolayer thickness varying between 1 and 5 µm which resulted in a monolayer volume of 2.36E04 and 1.18E03 ml respectively. A tracer concentration of 0.1 mg/ml was chosen for simulations to match the concentration used in the in vitro setup. The monolayer thickness in the $z$direction was divided into 70 evenly spaced pieces for simulations and the time span was restricted to 24 hours. A range of effective diffusion constants were simulated where the resulting concentration (mass/volume) was converted to mass in mg and the area under the curve for each time point was generated for comparison with the HBMEC in vitro data.
The in vitro data was split into a calibration and validation group where time points at hours 3, 12, 15, 21, and 24 were allocated to the calibration set and time points at hours 6, 9, and 18 were allocated to the validation set. A range of effective diffusion constants (${D}_{eff}$, cm^{2}/s) were selected for simulations and simulated time points at hours 3, 12, 15, 21, and 24 were extracted and compared to the in vitro time points at those hours to determine a range of diffusion constants that produced results not statistically significantly different from the in vitro data. The resulting diffusion constants were used to run simulations and extract time points at hours 6, 9, and 18 and those were compared to the in vitro time points at those hours. The percent error was calculated for the validation sets to determine how far the simulated models deviated from the in vitro model.
2.4. Statistical methods
The means of each experimental time point for each fluorescent tracer were compared using an Analysis of Variance (ANOVA) if equal variance could be assumed at all time points per tracer otherwise the nonparametric equivalent of an ANOVA (Welch Test) was used. Posthoc comparisons were conducted if the $p$value of the ANOVA was less than the alpha ($p<$0.05). The GamesHowell posthoc test was used if the nonparametric ANOVA was conducted and a Tukey Honest Significant Differences (HSD) was used if the parametric ANOVA was conducted.
A linear regression comparison of slopes was used to compare the in vitro time points to the simulated time points for the calibration data at hours 3, 12, 15, 18, 21, and 24. A $p$value less than 0.05 indicates the slopes of the in vitro and simulated models are significantly different while a $p$value larger than 0.05 indicates the slopes are not statistically significantly different. All statistical testing was conducted using Minitab 18 software (Minitab, Inc).
3. Results
3.1. In vitro
With the assumption of a 1 µm monolayer thickness, the calculated total average mass (mg) of tracer for FITC ranged from 5.59E09 mg with a standard deviation of 3.85E10 mg after 3 hours of incubation to 4.85E08 ± 2.76E09 mg at 24 hours (Fig. 1(a)). A nonparametric Welch Test was used to compare average mass among all FITC time points at the 95% confidence level resulting in an $F$statistic = 192.36 and a $p$value of 0.00. A GamesHowell posthoc multiple comparison test was used showing the 24 and 21 hour time points to be significantly greater than all earlier time points ($p$value < 0.05). The 18 hour time point was significantly greater than the 3 through 12 hour time points ($p$value < 0.05). The 15 hour time points were significantly greater than the 3 and 6 hour time points ($p$value < 0.05). The 12, 9, and 6 hour time points were significantly greater than the 3 hour time point with $p$values of 0.001, 0.00, and 0.022 respectively. Fluorescein range from 4.90E08 ± 1.94E09 mg at 3 hours to 5.63E08 ± 4.31E09 mg at 24 hours of incubation (Fig. 1(a)). An ANOVA was used to compare all time points among the fluorescein data set resulting in an $F$statistic = 2.19 and a $p$value = 0.056, no additional posthoc statistical testing was conducted. NaFl range from 5.69E08 ± 2.90E09 mg at 3 hours to 8.39E08 ± 1.12E08 mg at 24 hours of incubation (Fig. 1(a)). A nonparametric Welch Test was conducted resulting in an $F$statistic = 6.33 and a $p$value = 0.001. A GamesHowell posthoc test showed only the 24 hour time point being significantly greater than the 3 and 6 hour time points with $p$values of 0.016 and 0.046 respectively. The 9 hour time point was significantly greater than the 3 hour time point with a $p$value of 0.045.
Under the assumption of a 5 µm monolayer thickness, the calculated total average mass in mg of tracer for FITC ranged from 2.80E08 mg with a standard deviation of 1.94E10 mg after 3 hours of incubation to 2.42E07 ± 1.40E08 mg at 24 hours (Fig. 1(b)). Fluorescein range from 2.45E07 ± 1.10E08 mg at 3 hours to 2.82E07 ± 2.16E08 mg at 24 hours of incubation (Fig. 1(b)). NaFl range from 2.85E07 ± 1.42E08 mg at 3 hours to 4.20E07 ± 5.58E08 mg at 24 hours of incubation (Fig. 1(b)). The statistical testing at the 5 µm monolayer thickness returned the same results as those found above with the 1 µm monolayer thickness.
3.2. Linear Regression: FITC Calibration
For a 1 µm monolayer thickness, FITC calibration simulations run using effective diffusion constants (${D}_{eff}$) of 5.00E21, 1.00E20, and 5.00E20 cm^{2}/s (Fig. 3(a)) created a range of data close to the in vitro calibration data set at hours 3, 12, 15, 18, 21, and 24 hours. The slopes of each simulated data set were statistically compared to the in vitro data set using a linear regression approach. There was not a significant difference in slopes comparing the 5.00E21 cm^{2}/s simulation to the in vitro data set with a $t$statistic of 2.179 and a $p$value of 0.066. The 1.00E20 and 5.00E20 cm^{2}/s simulations yielded similar results with $t$statistics of 2.011 and 0.314 respectively with $p$values of 0.084 and 0.763 respectively.
With a 5 µm monolayer thickness, calibration simulations run using effective diffusion constants (${D}_{eff}$) of 5.00E21 and 1.00E20 cm^{2}/s (Fig. 3(b)) resulted in $t$statistics of 1.647 and 4.865 respectively with $p$values of 0.144 and 0.002 respectively when compared to in vitro data.
Fig. 3. a) Simulations run under a 1 µm monolayer thickness assumption, all simulations were found to not be statistically significantly different from the in vitro data; b) simulations run under a 5 µm monolayer thickness assumption, using a simulated ${D}_{eff}$ of 5.00E21 cm^{2}/s was shown to not be statistically significant while a simulated ${D}_{eff}$ of 1.00E20 cm^{2}/s yielded a significant difference at a $p$value of 0.002 for FITC
a)
b)
3.3. Linear regression: fluorescein calibration
For a 1 µm monolayer thickness, fluorescein calibration simulations run using effective diffusion constants (${D}_{eff}$) of 5.00E21 and 1.00E20 cm^{2}/s (Fig. 4(a)) created a range of data close to the in vitro calibration data set at hours 3, 12, 15, 18, 21, and 24 hours. The slopes of each simulated data set were statistically compared to the in vitro data set using a linear regression approach. There was a significant difference in slopes comparing the 5.00E21 and 1.00E20 cm^{2}/s simulations to the in vitro data set with $t$statistic of 5.122 and 3.594 respectively and $p$values of 0.001 and 0.009 respectively.
Fig. 4. a) Simulations run under a 1 µm monolayer thickness assumption, all simulations were found to be statistically significantly different from the in vitro data with $p$values < 0.05; b) simulations run under a 5 µm monolayer thickness assumption were found to be statistically significant ($p$value < 0.05) for both the 1.00E18 and 1.00E17 cm^{2}/s effective diffusion constants with the 5.00E18 cm^{2}/s diffusion constant gave a $p$value of 0.547
a)
b)
With a 5 µm monolayer thickness, calibration simulations run using effective diffusion constants (${D}_{eff}$) of 1.00E18, 5.00E18, and 1.00E17 cm^{2}/s (Fig. 4(b)) resulted in $t$statistics of 31.364, 0.632, and 3.125 respectively and $p$values of 8.65E9, 0.547, and 0.017 respectively.
3.4. Linear regression: NaFl calibration
For a 1 µm monolayer thickness, fluorescein calibration simulations run using effective diffusion constants (${D}_{eff}$) of 5.00E20 and 5.00E19 cm^{2}/s (Fig. 5(a)) created a range of data close to the in vitro calibration data set at hours 3, 12, 15, 18, 21, and 24 hours. The slopes of each simulated data set were statistically compared to the in vitro data set using a linear regression approach. There was a significant difference in slopes comparing the 5.00E20 and 5.00E19 cm^{2}/s simulations to the in vitro data set with $t$statistic of 20.793 and 3.297 respectively and $p$values of 1.49E07 and 0.013 respectively.
With a 5 µm monolayer thickness, calibration simulations run using effective diffusion constants (${D}_{eff}$) of 5.00E19, 5.00E18, and 5.00E17 cm^{2}/s (Fig. 5(b)) resulted in $t$statistics of 20.131, 0.831, and 4.923 respectively and $p$values of 1.87E07, 0.433, and 0.002 respectively.
Fig. 5. a) Simulations run under a 1 µm monolayer thickness assumption, all simulations were found to be statistically significantly different from the in vitro data with $p$values < 0.05; b) simulations run under a 5 µm monolayer thickness assumption were found to be statistically significant ($p$value < 0.05) for both the 5.00E19 and 5.00E17 cm^{2}/s effective diffusion constants with the 5.00E18 cm^{2}/s diffusion constant gave a $p$value of 0.547
a)
b)
3.5. Percent error: FITC, fluorescein, and NaFl validation
When increasing the effective diffusion (${D}_{eff}$) constant of FITC, under a 1 µm assumption from 1E20 to 5E20 cm^{2}/s shows a reduction in percent error of prediction of simulations (Table 1) at time points of 6 and 9 hours, but an increase in percent error of prediction at 18 hours of the validation data set using Eq. (5). When the monolayer thickness is increased to 5 µm and the ${D}_{eff}$ range is 5.00E19 and 5E18 cm^{2}/s the percent error of prediction increased for all time points:
For fluorescein, under a 1 µm assumption, the ${D}_{eff}$ from 1E19 to 5E19 cm^{2}/s shows a reduction in percent error of prediction of simulations (Table 2) at the 6 hour time point and an increase in error rate at the 9 and 18 hour time points (Table 2). Under a monolayer assumption of 5 µm there is a reduction in percent error at all three time points.
Table 1. FITC percent error (%) of simulations compared to in vitro data at 6, 9, and 18 hours
Hours

Mass (mg)^{A}

Mass
(mg)^{B}

(%)^{C}

Mass
(mg)^{B}

(%)^{C}


Monolayer thickness (µm)

Monolayer volume (ml)

In Vitro

${D}_{eff}=$
1E20 cm^{2}/s

Percent error

${D}_{eff}=$
5E20 cm^{2}/s

Percent error


1.0

2.36E04

6

7.03E09

1.21E09

82.74

6.07E09

13.69

9

7.54E09

1.82E09

75.86

9.10E09

20.68


18

9.94E09

3.64E09

63.38

1.82E08

83.09


Monolayer thickness (µm)

Monolayer volume (ml)

In Vitro

${D}_{eff}=$
5E19 cm^{2}/s

Percent error

${D}_{eff}=$
5E18 cm^{2}/s

Percent Error


5.0

1.18E03

6

3.52E08

1.21E08

65.48

1.21E07

245.23

9

3.77E08

1.82E08

51.72

1.82E07

382.82


18

4.98E08

3.64E08

26.85

3.64E07

631.50


A: Average in vitro mass in milligrams; B: Simulated mass in milligrams;
C: Percent error found using Eq. (5)

Table 2. Fluorescein percent error (%) of simulations compared to in vitro data at 6, 9, and 18 hours
Hours

Mass (mg)^{A}

Mass
(mg)^{B}

(%)^{C}

Mass
(mg)^{B}

(%)^{C}


Monolayer thickness (µm)

Monolayer volume (ml)

In Vitro

${D}_{eff}=$
1E19 cm^{2}/s

Percent error

${D}_{eff}=$
5E19 cm^{2}/s

Percent error


1.0

2.36E04

6

5.30E08

1.21E08

77.12

6.07E08

14.41

9

5.33E08

1.82E08

65.85

9.10E08

70.75


18

5.83E08

3.64E08

37.50

1.82E07

212.47


Monolayer thickness (µm)

Monolayer volume (ml)

In Vitro

${D}_{eff}=$
1E18 cm^{2}/s

Percent error

${D}_{eff}=$
5E18 cm^{2}/s

Percent error


5.0

1.18E03

6

2.65E07

2.43E08

90.85

1.21E07

54.24

9

2.67E07

3.64E08

86.35

1.82E07

31.74


18

2.91E07

7.28E08

75.01

3.64E07

24.96


A: Average in vitro mass in milligrams; B: Simulated mass in milligrams;
C: Percent error found using Eq. (5)

Table 3. NaFl percent error (%) of simulations compared to in vitro data at 6, 9, and 18 hours
Hours

Mass (mg)^{A}

Mass
(mg)^{B}

(%)^{C}

Mass
(mg)^{B}

(%)^{C}


Monolayer thickness (µm)

Monolayer volume (ml)

In Vitro

${D}_{eff}=$
5E20 cm^{2}/s

Percent error

${D}_{eff}=$
5E19 cm^{2}/s

Percent error


1.0

2.36E04

6

6.29E08

6.07E09

90.35

6.07E08

3.54

9

6.44E08

9.10E09

85.86

9.10E08

41.39


18

6.51E08

1.82E08

72.04

1.82E07

179.59


Monolayer thickness (µm)

Monolayer volume (ml)

In Vitro

${D}_{eff}=$
5E19 cm^{2}/s

Percent error

${D}_{eff}=$
5E18 cm^{2}/s

Percent error


5.0

1.18E03

6

3.15E07

1.21E08

96.14

1.21E07

61.44

9

3.22E07

1.82E08

94.34

1.82E07

43.44


18

3.26E07

3.64E08

88.82

3.64E07

11.78


A: Average in vitro mass in milligrams; B: Simulated mass in milligrams;
C: Percent error found using Eq. (5)

NaFl, under an assumed 1 µm monolayer thickness, produced a decrease in percent error at the 6 and 9 hour time points for ${D}_{eff}$ from 5E20 to 5E19 cm^{2}/s and an increase in percent error at the 18 hour time point (Table 3). With a 5 µm monolayer thickness there is a decrease in all time points when the effective diffusion constant (${D}_{eff}$) is increased from 5E19 to 5E18 cm^{2}/s.
4. Conclusions
For both fluorescein and NaFl tracers the in vitro data showed high fluorescence compared to FITC time points with few statistically significant increases in tracer uptake over the 24hour period. Additionally, fluorescein and NaFl simulated diffusion models have trouble narrowing in on a suitable effective diffusion constant range due to very little change in the in vitro data. It is probable that the fluorescein and NaFl tracer mass could be better explained by a more appropriate mechanism other than the passive diffusion such as binding to the cell surface or embedding within the lipid bilayer. Only FITC showed an increase in in vitro tracer that could be closely recreated through computational simulation. For FITC, under a 1 µm thickness, effective diffusion constants between the range of 1.00E20 and 5.00E20 cm^{2}/s can serve as a good approximation of in vitro data. At the 5 µm thickness, effective diffusion constants around 5.00E19 cm^{2}/s would serve as a good approximation for the in vitro data.
Fluorescent tracers have a multitude of uses in biomedical sciences to stain the exterior and interior of cells. This paper attempted to determine the potential magnitude of free diffusion of 3 tracers into HBMEC monolayers at different thicknesses of 1 and 5 µm. It was determined that although higher than FITC, both fluorescein alone and NaFl showed very little increase in tracer concentration from 3 hours to 24 hours while FITC continued to increase at statistically significant levels. Under the assumption of diffusion with no reaction and binding terms, fluorescein and NaFl effective diffusion constants were not able to be narrowed to produce simulated curves similar to the in vitro calibration data but a narrow range can be determined for FITC. This suggests that the retention of fluorescein and NaFl to the HBMEC monolayer may be better explained by another mechanism while FITC uptake make be reasonably approximated using a diffusion model. It is also important to note that in vitro results achieved are highly dependent on a multitude of experimental factors such as pH, temperature, cell type, and concentration levels.
Acknowledgements
We would like to thank the Virginia Tech Initiative for Maximizing Student Development (IMSD) program and the Virginia Tech Center for Autism Research (VTCAR) for partial support of this research project.
References
 Lombard J. Once upon a time the cell membrane: 175 years of cell boundary research. Biology Direct, Vol. 9, 2014, p. 32. [Publisher]
 Owen S. M., Gaus K., Magee A. I., Cebecauer M. Dynamic organization of lymphocyte plasma membrane: lessons from advanced imaging methods. Immunology, Vol. 131, 2010, p. 18. [Search CrossRef]
 Stone M. B., Shelby S. A., Veatch S. L. Superresolution microscopy: shedding light on the cellular plasma membrane. Chemical Reviews, Vol. 117, 2017, p. 74577477. [Publisher]
 Thornalley P. J. Isothiocyanates: mechanism of cancer chemopreventive action. Anticancer Drugs, Vol. 13, 2002, p. 331338. [Publisher]
 Wu X., Zhou Q., Xu K. Are isothiocyanates potential anticancer drugs? Acta Pharmacologica Sinica, Vol. 30, 2009, p. 501512. [Search CrossRef]
 Cross J. V., Foss F. W., Rady J. M., Macdonald T. L., Templeton D. J. The isothiocyanate class of bioactive nutrients covalently inhibit the MEKKI protein kinase. BMC Cancer, Vol. 7, 2007, p. 183. [Publisher]
 Hanschen F. S., Brüggemann N., Brodehl A., Mewis I., Schreiner M., Rohn S., Kroh L. W. Characterization of products from the reaction of glucosinolatederived isothiocyanates with cysteine and lysine derivatives formed in either model systems or broccoli sprouts. Journal of Agricultural and Food Chemistry, Vol. 60, 2012, p. 77357745. [Publisher]
 Ioannides C., Konsue N. A principal mechanism for the cancer chemopreventive activity of phenethyl isothiocyanate is modulation of carcinogen metabolism. Drug Metabolism Reviews, Vol. 47, 2015, p. 356373. [Publisher]
 Mi L., Wang X., Govind S., Hood B. L., Veenstra T. D., Conrads T. P., Saha D. T., Goldman R., Chung F.L. The role of protein binding in induction of apoptosis by phenethyl isothiocyanate and sulforaphane in human nonsmall lung cancer cells. Cancer Research, Vol. 67, 2007, p. 64096416. [Publisher]
 Ji Y. A., Morris M. E. Membrane transport of dietary phenyl isothiocyanate by ABCG2 (breast cancer resistance protein). Molecular Pharmaceutics, Vol. 2, 2005, p. 414419. [Publisher]
 Zhang Y. Molecular mechanism of rapid cellular accumulation of anticarcinogenic isothiocyanates. Carcinogenesis, Vol. 22, 2001, p. 425431. [Publisher]
 Zhang Y., Talalay P. Mechanism of differential potencies of isothiocyanates as inducers of anticarcinogenic phase 2 enzymes. Cancer Research, Vol. 58, 1998, p. 46324639. [Search CrossRef]
 Filatov A. V., Krotov G. I., Zgoda V. G., Volkov Y. Fluorescent immunoprecipitation analysis of cell surface proteins: a methodology compatible with massspectrometry. Journal of Immunological Methods, Vol. 319, 2007, p. 2133. [Publisher]
 Maeda H., Ishida N., Kawauchi H., Tuzimura K. Reaction of fluoresceinisothiocyanate with proteins and amino acids. The Journal of Biochemistry, Vol. 65, 1969, p. 777783. [Publisher]
 Miyata S., Morita S. A new method for visualization of endothelial cells and extravascular leakage in adult mouse brain using fluorescein isothiocyanate. Journal of Neuroscience Methods, Vol. 202, 2011, p. 916. [Publisher]
 Aplin J. D., Hughes R. C. Modification of fibroblast surface amines alters receptormediated cell spreading on proteincoated substrata but not adsorptive endocytosis. Journal of Cell Science, Vol. 49, 1981, p. 283297. [Search CrossRef]
 Butters T. D., Devalia V., Aplin J. D., Hughes R.C. Inhibition of fibronectinmediated adhesion of hamster fibroblasts to substratum: effects of tunicamycin and some cell surface modifying reagents. Journal of Cell Science, Vol. 44, 1980, p. 3358. [Search CrossRef]
 Butcher E. C., Weissman I. L. Direct fluorescent labeling of cells with fluorescein or rhodamine isothiocyanate. I. Technical aspects. Journal of Immunological Methods, Vol. 37, 1980, p. 97108. [Search CrossRef]
 Jerry L. M., Sullivan A. K. The lymphocyte plasma membrane: locus of control in the immune response. In Vitro, Vol. 12, 1976, p. 236259. [Publisher]
 Spotl L., Sarti A., Dierich M. P., Most J. Cell membrane labeling with fluorescent dyes for the demonstration of cytokineinduced fusion between monocytes and tumor cells. Cytometry, Vol. 21, 1995, p. 160169. [Publisher]
 Kaya M., Ahishali B. Assessment of Permeability in Barrier Type of Endothelium in Brain Using Tracers: Evans Blue, Sodium Fluorescein, and Horseradish Peroxidase. Permeability Barrier. Methods in Molecular Biology (Methods and Protocols), Humana Press, Vol 763, 2011, p. 369382. [Search CrossRef]
 Braginskaja O. V., Lazarev V. V., Pershina I. N., Petrov K. V., Rubin L. B., Tinkhonova O. V. Sodium fluorescein accumulation in cultured cells. General Physiology and Biophysics, Vol. 12, 1993, p. 453464. [Search CrossRef]
 Breeuwer P., Drocourt J. L., Bunschoten N., Zwietering M. H., Rombouts F. M., Abee T. Characterization of uptake and hydrolysis of fluorescein diacetate and carboxyfluorescein diacetate by intracellular esterases in Saccharomyces cerevisiae, which result in accumulation of fluorescent product. Applied and Environmental Microbiology, Vol. 61, 1995, p. 16141619. [Search CrossRef]
 Meisingset K. K., Steen H. B. Intracellular binding of fluorescein in lymphocytes. Cytometry, Vol. 1, 1981, p. 272278. [Publisher]
 Cooper G. M. The Cell. Sinauer Associates, 2000. [Search CrossRef]
 Lodish H., Berk A., Zipursky S. L., Matsudaira P., Baltimore D., Darnell J. Molecular Cell Biology. W. H. Freeman, 2000. [Search CrossRef]
 Zehnder S., Suaris M., Bellaire M. M., Angelini T. E. Cell volume fluctuations in mdck monolayers. Biophysical Journal, Vol. 108, 2015, p. 247250. [Publisher]
 Bettega D., Calzolari P., Doglia S. M., Dulio B., Tallone L., Villa A. M. Technical report: cell thickness measurements by confocal fluorescence microscopy on C3H10T1/2 and V79 cells. International Journal of Radiation Biology, Vol. 74, 1998, p. 397403. [Search CrossRef]
 Stins M. F., Gilles F., Kim K. S. Selective expression of adhesion molecules on human brain microvascular endothelial cells. Journal of Neuroimmunology, Vol. 76, 1997, p. 8190. [Publisher]