Mathematical Modeling of Perifusion Cell Culture Experiments

In perifusion cell cultures, the culture medium flows continuously through a chamber containing immobilized cells and the effluent is collected at the end. In our main applications, gonadotropin releasing hormone (GnRH) or oxytocin is introduced into the chamber as the input. They stimulate the cells to secrete luteinizing hormone (LH), which is collected in the effluent. To relate the effluent LH concentration to the cellular processes producing it, we develop and analyze a mathematical model consisting of coupled partial differential equations describing the intracellular signaling and the movement of substances in the cell chamber. We analyze three different data sets and give cellular mechanisms that explain the data. Our model indicates that two negative feedback loops, one fast and one slow, are needed to explain the data and we give their biological bases. We demonstrate that different LH outcomes in oxytocin and GnRH stimulations might originate from different receptor dynamics. We analyze the model to understand the influence of parameters, like the rate of the medium flow or the fraction collection time, on the experimental outcomes. We investigate how the rate of binding and dissociation of the input hormone to and from its receptor influence its movement down the chamber. Finally, we formulate and analyze simpler models that allow us to predict the distortion of a square pulse due to hormone-receptor interactions and to estimate parameters using perifusion data. We show that in the limit of high binding and dissociation the square pulse moves as a diffusing Gaussian and in this limit the biological parameters can be estimated.

The Effects of Brominated Flame Retardants on Thyroid Hormone Homeostasis in Human Placenta Tissues and Cell Culture

Polybrominated diphenyl ethers (PBDEs) are a class of brominated flame retardants (BFRs) that have been heavily used in consumer products such as furniture foams, plastics, and textiles since the mid-1970’s. BFRs are added to products in order to meet state flammability standards intended to increase indoor safety in the event of a fire. The three commercial PBDE mixtures, Penta-, Octa-, and DecaBDE, have all been banned in the United States, however, limited use of DecaBDE is still permitted. PBDEs were phased out of production and added to the Stockholm Convention due to concerns over their environmental persistence and toxicity. Human exposure to PBDEs occurs primarily through the inadvertent ingestion of contaminated house dust, as well as though dietary sources. Despite the phase-out and discontinued use of PBDEs, human exposure to this class of chemicals is likely to continue for decades due to the continued use of treated products and existing environmental reservoirs of PBDEs. Extensive research over the years has shown that PBDEs disrupt thyroid hormone (TH) levels and neurodevelopmental endpoints in rodent and fish models. Additionally, there is growing epidemiological evidence linking PBDE exposure in humans to altered TH homeostasis and neurodevelopmental impairments in children. Due to the importance of THs throughout gestation, there is a great need to understand the effects of BFRs on the developing fetus. Specifically, the placenta plays a critical role in the transport, metabolism, and delivery of THs to the fetal compartment during pregnancy and is a likely target for BFR bioaccumulation and endocrine disruption. The central hypothesis of this dissertation research is that BFRs disrupt the activity of TH sulfotransferase (SULT) enzymes, thereby altering TH concentrations in the placenta.

In the first aim of this dissertation research, the concentrations of PBDEs and 2,4,6-TBP were measured in a cohort of 102 placenta tissue samples from an ongoing pregnancy cohort in Durham, NC. Methods were developed for the extraction and analysis of the BFR analytes. It was found that 2,4,6-TBP was significantly correlated with all PBDE analytes, indicating that 2,4,6-TBP may share common product applications with PBDEs or that 2,4,6-TBP is a metabolite of PBDE compounds. Additionally, this was the first study to measure 2,4,6-TBP in human placenta tissues.

In the second aim of this dissertation research, the placenta tissue concentrations of THs, as well as the endogenous activity of deiodinase (DI) and TH SULT enzymes were quantified using the same cohort of 102 placenta tissue samples. Enzyme activity was detected in all samples and this was the first study to measure TH DI and SULT activity in human placenta tissues. Enzyme activities and TH concentrations were compared with BFR concentrations measured in Aim 1. There were few statistically significant associations observed for the combined data, however, upon stratifying the data set based on infant sex, additional significant associations were observed. For example, among males, those with the highest concentrations of BDE-99 in placenta had T3 levels 0.80 times those with the lowest concentration of BDE-99 (95% confidence interval (CI): 0.59, 1.07). Whereas females with the highest concentrations of BDE-99 in placenta had T3 levels 1.50 times those with the lowest concentration of BDE-99 (95% CI: 1.10, 2.04). Additionally, all BFR analyte concentrations were higher in the placenta of males versus females and they were significantly higher for 2,4,6-TBP and BDE-209. 3,3’-T2 SULT activity was significantly higher in female placenta tissues, while type 3 DI activity was significantly higher in male placenta tissues. This research is the first to show sex-specific differences in the bioaccumulation of BFRs in human placenta tissue, as well as differences in TH concentrations and endogenous DI and SULT activity. The underlying mechanisms of these observed sex differences warrant further investigation.

In the third aim of this dissertation research, the effects of BFRs were examined in a human choriocarcinoma placenta cell line, BeWo. Michaelis-Menten parameters and inhibition curves were calculated for 2,4,6-TBP, 3-OH BDE-47, and 6-OH BDE-47. 2,4,6-TBP was shown to be the most potent inhibitor of 3,3’-T2 SULT activity with a calculated IC50 value of 11.6 nM. It was also shown that 2,4,6-TBP and 3-OH BDE-47 exhibit mixed inhibition of 3,3’-T2 sulfation in BeWo cell homogenates. Next, a series of cell culture exposure experiments were performed using 1, 6, 12, and 24 hour exposure durations. Once again, 2,4,6-TBP was shown to be the most potent inhibitor of basal 3,3’-T2 SULT activity by significantly decreasing activity at the high and medium dose (1 M and 0.5 M, respectively) at all measured time points. Interestingly, BDE-99 was also shown to inhibit basal 3,3’-T2 SULT activity in BeWo cells following the 24 hour exposure, despite exhibiting no inhibitory effects in the BeWo cell homogenate experiments. This indicates that BDE-99 must act through a pathway other than direct enzyme inhibition. Following exposures, the TH concentrations in the cell culture growth media and mRNA expression of TH-related genes were also examined. There was no observed effect of BFR treatment on these endpoints. Future work should focus on determining the downstream biological effects of TH SULT disruption in placental cells, as well as the underlying mechanisms of action responsible for reductions in basal TH SULT activity following BFR exposure.

This was one of the first studies to measure BFRs in a cohort of placenta tissue samples from the United States and the first study to measure THs, DI activity, and SULT activity in human placenta tissues. This research provides a novel contribution to our growing understanding of the effects of BFRs on TH homeostasis within the human placenta, and provides further evidence for sex-specific differences within this important organ. Future research should continue to investigate the effects of environmental contaminants on TH homeostasis within the placenta, as this represents the most critical and vulnerable stage of human development.

Scaffold-free, Human Mesenchymal Stem Cell-Based Tissue Engineered Blood Vessels.

Tissue-engineered blood vessels (TEBV) can serve as vascular grafts and may also play an important role in the development of organs-on-a-chip. Most TEBV construction involves scaffolding with biomaterials such as collagen gel or electrospun fibrous mesh. Hypothesizing that a scaffold-free TEBV may be advantageous, we constructed a tubular structure (1 mm i.d.) from aligned human mesenchymal cell sheets (hMSC) as the wall and human endothelial progenitor cell (hEPC) coating as the lumen. The burst pressure of the scaffold-free TEBV was above 200 mmHg after three weeks of sequential culture in a rotating wall bioreactor and perfusion at 6.8 dynes/cm(2). The interwoven organization of the cell layers and extensive extracellular matrix (ECM) formation of the hMSC-based TEBV resembled that of native blood vessels. The TEBV exhibited flow-mediated vasodilation, vasoconstriction after exposure to 1 μM phenylephrine and released nitric oxide in a manner similar to that of porcine femoral vein. HL-60 cells attached to the TEBV lumen after TNF-α activation to suggest a functional endothelium. This study demonstrates the potential of a hEPC endothelialized hMSC-based TEBV for drug screening.

Polyploidy and Mitotic Cell Death are Two Distinct HIV-1 Vpr-Driven Outcomes in Renal Tubule Epithelial Cells

Given the emerging epidemic of renal disease in HIV+ patients and the fact that HIV DNA and RNA persist in the kidneys of HIV+ patients despite therapy, it is necessary to understand the role of direct HIV-1 infection of the kidney. HIV-associated kidney disease pathogenesis is attributed in large part to viral proteins. Expression of Vpr in renal tubule epithelial cells (RTECs) induces G2 arrest, apoptosis and polyploidy. The ability of a subset of cells to overcome the G2/M block and progress to polyploidy is not well understood. Polyploidy frequently associates with a bypass of cell death and disease pathogenesis. Given the ability of the kidney to serve as a unique compartment for HIV-1 infection, and the observed occurrence of polyploid cells in HIV+ renal cells, it is critical to understand the mechanisms and consequences of Vpr-induced polyploidy.

Here I determined effects of HIV-1 Vpr expression in renal cells using highly efficient transduction with VSV.G pseudotyped lentiviral vectors expressing Vpr in the HK2 human tubule epithelial cell line. Using FACS, fluorescence microscopy, and live cell imaging I show that G2 escape immediately precedes a critical junction between two distinct outcomes in Vpr+ RTECs: mitotic cell death and polyploidy. Vpr+ cells that evade aberrant mitosis and become polyploid have a substantially higher survival rate than those that undergo complete mitosis, and this survival correlates with enrichment for polyploidy in cell culture over time. Further, I identify a novel role for ATM kinase in promoting G2 arrest escape and polyploidy in this context. In summary, my work identifies ATM-dependent override of Vpr-mediated G2/M arrest as a critical determinant of cell fate Vpr+ RTECs. Further, our work highlights how a poorly understood HIV mechanism, ploidy increase, may offer insight into key processes of reservoir establishment and disease pathogenesis in HIV+ kidneys.

Spontaneous Unfolding and Refolding of FNIII Domains Assayed by Thiol Exchange

Fibronectin (FN) is a large extracellular matrix (ECM) protein that is made up of

type I (FNI), type II (FNII), & type III (FNIII) domains. It assembles into an insoluble

supra-­‐‑molecular structure: the fibrillar FN matrix. FN fibrillogenesis is a cell‐‑mediated process, which is initiated when FN binds to integrins on the cell surface. The FN matrix plays an important role in cell migration, proliferation, signaling & adhesion. Despite decades of research, the FN matrix is one of the least understood supra-­‐‑molecular protein assemblies. There have been several attempts to elucidate the exact mechanism of matrix assembly resulting in significant progress in the field but it is still unclear as to what are FN-­‐‑FN interactions, the nature of these interactions and the domains of FN that

are in contact with each other. FN matrix fibrils are elastic in nature. Two models have been proposed to explain the elasticity of the fibrils. The first model: the ‘domain unfolding’ model postulates that the unraveling of FNIII domains under tension explains fibril elasticity.

The second model relies on the conformational change of FN from compact to extended to explain fibril elasticity. FN contain 15 FNIII domains, each a 7-­‐‑strand beta sandwich. Earlier work from our lab used the technique of labeling a buried Cys to study the ‘domain unfolding’ model. They used mutant FNs containing a buried Cys in a single FNIII domain and found that 6 of the 15 FNIII domains label in matrix fibrils. Domain unfolding due to tension, matrix associated conformational changes or spontaneous folding and unfolding are all possible explanation for labeling of the buried Cys. The present study also uses the technique of labeling a buried Cys to address whether it is spontaneous folding and unfolding that labels FNIII domains in cell culture. We used thiol reactive DTNB to measure the kinetics of labeling of buried Cys in eleven FN III domains over a wide range of urea concentrations (0-­‐‑9M). The kinetics data were globally fit using Mathematica. The results are equivalent to those of H-­‐‑D exchange, and

provide a comprehensive analysis of stability and unfolding/folding kinetics of each

domain. For two of the six domains spontaneous folding and unfolding is possibly the reason for labeling in cell culture. For the rest of the four domains it is probably matrix associated conformational changes or tension induced unfolding.

A long-­‐‑standing debate in the protein-­‐‑folding field is whether unfolding rate

constants or folding rate constants correlate to the stability of a protein. FNIII domains all have the same ß sandwich structure but very different stabilities and amino acid sequences. Our study analyzed the kinetics of unfolding and folding and stabilities of eleven FNIII domains and our results show that folding rate constants for FNIII domains are relatively similar and the unfolding rates vary widely and correlate to stability. FN forms a fibrillar matrix and the FN-­‐‑FN interactions during matrix fibril formation are not known. FNI 1-­‐‑9 or the N-­‐‑ terminal region is indispensible for matrix formation and its major binding partner has been shown to be FNIII 2. Earlier work from our lab, using FRET analysis showed that the interaction of FNI 1-­‐‑9 with a destabilized FNIII 2 (missing the G strand, FNIII 2ΔG) reduces the FRET efficiency. This efficiency is restored in the presence of FUD (bacterial adhesion from S. pyogenes) that has been known to interact with FNI 1-­‐‑9 via a tandem ß zipper. In the present study we

use FRET analysis and a series of deletion mutants of FNIII 2ΔG to study the shortest fragment of FNIII 2ΔG that is required to bind FNI 1-­‐‑9. Our results presented here are qualitative and show that FNIII 2ΔC’EFG is the shortest fragment required to bind FNI 1-­‐‑9. Deletion of one more strand abolishes the interaction with FNI 1-­‐‑9.