Computational Inference of Genome-Wide Protein-DNA Interactions Using High-Throughput Genomic Data

Transcriptional regulation has been studied intensively in recent decades. One important aspect of this regulation is the interaction between regulatory proteins, such as transcription factors (TF) and nucleosomes, and the genome. Different high-throughput techniques have been invented to map these interactions genome-wide, including ChIP-based methods (ChIP-chip, ChIP-seq, etc.), nuclease digestion methods (DNase-seq, MNase-seq, etc.), and others. However, a single experimental technique often only provides partial and noisy information about the whole picture of protein-DNA interactions. Therefore, the overarching goal of this dissertation is to provide computational developments for jointly modeling different experimental datasets to achieve a holistic inference on the protein-DNA interaction landscape.

We first present a computational framework that can incorporate the protein binding information in MNase-seq data into a thermodynamic model of protein-DNA interaction. We use a correlation-based objective function to model the MNase-seq data and a Markov chain Monte Carlo method to maximize the function. Our results show that the inferred protein-DNA interaction landscape is concordant with the MNase-seq data and provides a mechanistic explanation for the experimentally collected MNase-seq fragments. Our framework is flexible and can easily incorporate other data sources. To demonstrate this flexibility, we use prior distributions to integrate experimentally measured protein concentrations.

We also study the ability of DNase-seq data to position nucleosomes. Traditionally, DNase-seq has only been widely used to identify DNase hypersensitive sites, which tend to be open chromatin regulatory regions devoid of nucleosomes. We reveal for the first time that DNase-seq datasets also contain substantial information about nucleosome translational positioning, and that existing DNase-seq data can be used to infer nucleosome positions with high accuracy. We develop a Bayes-factor-based nucleosome scoring method to position nucleosomes using DNase-seq data. Our approach utilizes several effective strategies to extract nucleosome positioning signals from the noisy DNase-seq data, including jointly modeling data points across the nucleosome body and explicitly modeling the quadratic and oscillatory DNase I digestion pattern on nucleosomes. We show that our DNase-seq-based nucleosome map is highly consistent with previous high-resolution maps. We also show that the oscillatory DNase I digestion pattern is useful in revealing the nucleosome rotational context around TF binding sites.

Finally, we present a state-space model (SSM) for jointly modeling different kinds of genomic data to provide an accurate view of the protein-DNA interaction landscape. We also provide an efficient expectation-maximization algorithm to learn model parameters from data. We first show in simulation studies that the SSM can effectively recover underlying true protein binding configurations. We then apply the SSM to model real genomic data (both DNase-seq and MNase-seq data). Through incrementally increasing the types of genomic data in the SSM, we show that different data types can contribute complementary information for the inference of protein binding landscape and that the most accurate inference comes from modeling all available datasets.

This dissertation provides a foundation for future research by taking a step toward the genome-wide inference of protein-DNA interaction landscape through data integration.

Ghrelin suppresses glucose-stimulated insulin secretion and deteriorates glucose tolerance in healthy humans.

OBJECTIVE: The orexigenic gut hormone ghrelin and its receptor are present in pancreatic islets. Although ghrelin reduces insulin secretion in rodents, its effect on insulin secretion in humans has not been established. The goal of this study was to test the hypothesis that circulating ghrelin suppresses glucose-stimulated insulin secretion in healthy subjects. RESEARCH DESIGN AND METHODS: Ghrelin (0.3, 0.9 and 1.5 nmol/kg/h) or saline was infused for more than 65 min in 12 healthy patients (8 male/4 female) on 4 separate occasions in a counterbalanced fashion. An intravenous glucose tolerance test was performed during steady state plasma ghrelin levels. The acute insulin response to intravenous glucose (AIRg) was calculated from plasma insulin concentrations between 2 and 10 min after the glucose bolus. Intravenous glucose tolerance was measured as the glucose disappearance constant (Kg) from 10 to 30 min. RESULTS: The three ghrelin infusions raised plasma total ghrelin concentrations to 4-, 15-, and 23-fold above the fasting level, respectively. Ghrelin infusion did not alter fasting plasma insulin or glucose, but compared with saline, the 0.3, 0.9, and 1.5 nmol/kg/h doses decreased AIRg (2,152 +/- 448 vs. 1,478 +/- 2,889, 1,419 +/- 275, and 1,120 +/- 174 pmol/l) and Kg (0.3 and 1.5 nmol/kg/h doses only) significantly (P < 0.05 for all). Ghrelin infusion raised plasma growth hormone and serum cortisol concentrations significantly (P < 0.001 for both), but had no effect on glucagon, epinephrine, or norepinephrine levels (P = 0.44, 0.74, and 0.48, respectively). CONCLUSIONS: This is a robust proof-of-concept study showing that exogenous ghrelin reduces glucose-stimulated insulin secretion and glucose disappearance in healthy humans. Our findings raise the possibility that endogenous ghrelin has a role in physiologic insulin secretion, and that ghrelin antagonists could improve beta-cell function.

Relationships between polybrominated diphenyl ether concentrations in house dust and serum.

Polybrominated diphenyl ethers (PBDEs) have been measured in the home environment and in humans, but studies linking environmental levels to body burdens are limited. This study examines the relationship between PBDE concentrations in house dust and serum from adults residing in these homes. We measured PBDE concentrations in house dust from 50 homes and in serum of male-female couples from 12 of the homes. Detection rates, dust-serum, and within-matrix correlations varied by PBDE congener. There was a strong correlation (r = 0.65-0.89, p < 0.05) between dust and serum concentrations of several predominant PBDE congeners (BDE 47, 99, and 100). Dust and serum levels of BDE 153 were not correlated (r < 0.01). The correlation of dust and serum levels of BDE 209 could not be evaluated due to low detection rates of BDE 209 in serum. Serum concentrations of the sum of BDE 47, 99, and 100 were also strongly correlated within couples (r = 0.85, p = 0.0005). This study provides evidence that house dust is a primary exposure pathway of PBDEs and supports the use of dust PBDE concentrations as a marker for exposure to PBDE congeners other than BDE 153.

Osmolyte-induced folding of an intrinsically disordered protein: folding mechanism in the absence of ligand.

Understanding the interconversion between thermodynamically distinguishable states present in a protein folding pathway provides not only the kinetics and energetics of protein folding but also insights into the functional roles of these states in biological systems. The protein component of the bacterial RNase P holoenzyme from Bacillus subtilis (P protein) was previously shown to be unfolded in the absence of its cognate RNA or other anionic ligands. P protein was used in this study as a model system to explore general features of intrinsically disordered protein (IDP) folding mechanisms. The use of trimethylamine N-oxide (TMAO), an osmolyte that stabilizes the unliganded folded form of the protein, enabled us to study the folding process of P protein in the absence of ligand. Transient stopped-flow kinetic traces at various final TMAO concentrations exhibited multiphasic kinetics. Equilibrium "cotitration" experiments were performed using both TMAO and urea during the titration to produce a urea-TMAO titration surface of P protein. Both kinetic and equilibrium studies show evidence of a previously undetected intermediate state in the P protein folding process. The intermediate state is significantly populated, and the folding rate constants are relatively slow compared to those of intrinsically folded proteins similar in size and topology. The experiments and analysis described serve as a useful example for mechanistic folding studies of other IDPs.

Adenovirus F protein as a delivery vehicle for botulinum B.

BACKGROUND: Immunization with recombinant carboxyl-terminal domain of the heavy chain (Hc domain) of botulinum neurotoxin (BoNT) stimulates protective immunity against native BoNT challenge. Most studies developing a botulism vaccine have focused on the whole Hc; however, since the principal protective epitopes are located within beta-trefoil domain (Hcbetatre), we hypothesize that immunization with the Hcbetatre domain is sufficient to confer protective immunity. In addition, enhancing its uptake subsequent to nasal delivery prompted development of an alternative vaccine strategy, and we hypothesize that the addition of targeting moiety adenovirus 2 fiber protein (Ad2F) may enhance such uptake during vaccination. RESULTS: The Hcbetatre serotype B immunogen was genetically fused to Ad2F (Hcbetatre/B-Ad2F), and its immunogenicity was tested in mice. In combination with the mucosal adjuvant, cholera toxin (CT), enhanced mucosal IgA and serum IgG Ab titers were induced by nasal Hcbetatre-Ad2F relative to Hcbetatre alone; however, similar Ab titers were obtained upon intramuscular immunization. These BoNT/B-specific Abs induced by nasal immunization were generally supported in large part by Th2 cells, as opposed to Hcbetatre-immunized mice that showed more mixed Th1 and Th2 cells. Using a mouse neutralization assay, sera from animals immunized with Hcbetatre and Hcbetatre-Ad2F protected mice against 2.0 LD50. CONCLUSION: These results demonstrate that Hcbetatre-based immunogens are highly immunogenic, especially when genetically fused to Ad2F, and Ad2F can be exploited as a vaccine delivery platform to the mucosa.

Interaction of Cryptococcus neoformans Rim101 and protein kinase A regulates capsule.

Cryptococcus neoformans is a prevalent human fungal pathogen that must survive within various tissues in order to establish a human infection. We have identified the C. neoformans Rim101 transcription factor, a highly conserved pH-response regulator in many fungal species. The rim101 multiply sign in circle mutant strain displays growth defects similar to other fungal species in the presence of alkaline pH, increased salt concentrations, and iron limitation. However, the rim101 multiply sign in circle strain is also characterized by a striking defect in capsule, an important virulence-associated phenotype. This capsular defect is likely due to alterations in polysaccharide attachment to the cell surface, not in polysaccharide biosynthesis. In contrast to many other C. neoformans capsule-defective strains, the rim101 multiply sign in circle mutant is hypervirulent in animal models of cryptococcosis. Whereas Rim101 activation in other fungal species occurs through the conserved Rim pathway, we demonstrate that C. neoformans Rim101 is also activated by the cAMP/PKA pathway. We report here that C. neoformans uses PKA and the Rim pathway to regulate the localization, activation, and processing of the Rim101 transcription factor. We also demonstrate specific host-relevant activating conditions for Rim101 cleavage, showing that C. neoformans has co-opted conserved signaling pathways to respond to the specific niche within the infected host. These results establish a novel mechanism for Rim101 activation and the integration of two conserved signaling cascades in response to host environmental conditions.

G protein signaling and vein graft intimal hyperplasia: reduction of intimal hyperplasia in vein grafts by a Gbetagamma inhibitor suggests a major role of G protein signaling in lesion development.

Vein grafting results in the development of intimal hyperplasia with accompanying changes in guanine nucleotide-binding (G) protein expression and function. Several serum mitogens that act through G protein-coupled receptors, such as lysophosphatidic acid, stimulate proliferative pathways that are dependent on the G protein betagamma subunit (Gbetagamma)-mediated activation of p21ras. This study examines the role of Gbetagamma signaling in intimal hyperplasia by targeting a gene encoding a specific Gbetagamma inhibitor in an experimental rabbit vein graft model. This inhibitor, the carboxyl terminus of the beta-adrenergic receptor kinase (betaARK(CT)), contains a Gbetagamma-binding domain. Vein graft intimal hyperplasia was significantly reduced by 37% (P<0.01), and physiological studies demonstrated that the normal alterations in G protein coupling phenotypically seen in this model were blocked by betaARK(CT) treatment. Thus, it appears that Gbetagamma-mediated pathways play a major role in intimal hyperplasia and that targeting inhibitors of Gbetagamma signaling offers novel intraoperative therapeutic modalities to inhibit the development of vein graft intimal hyperplasia and subsequent vein graft failure.

Impaired formation of beta-adrenergic receptor-nucleotide regulatory protein complexes in pseudohypoparathyroidism.

Decreased activity of the guanine nucleotide regulatory protein (N) of the adenylate cyclase system is present in cell membranes of some patients with pseudohypoparathyrodism (PHP-Ia) whereas others have normal activity of N (PHP-Ib). Low N activity in PHP-Ia results in a decrease in hormone (H)-stimulatable adenylate cyclase in various tissues, which might be due to decreased ability to form an agonist-specific high affinity complex composed of H, receptor (R), and N. To test this hypothesis, we compared beta-adrenergic agonist-specific binding properties in erythrocyte membranes from five patients with PHP-Ia (N = 45% of control), five patients with PHP-Ib (N = 97%), and five control subjects. Competition curves that were generated by increasing concentrations of the beta-agonist isoproterenol competing with [125I]pindolol were shallow (slope factors less than 1) and were computer fit to a two-state model with corresponding high and low affinity for the agonist. The agonist competition curves from the PHP-Ia patients were shifted significantly (P less than 0.02) to the right as a result of a significant (P less than 0.01) decrease in the percent of beta-adrenergic receptors in the high affinity state from 64 +/- 22% in PHP-Ib and 56 +/- 5% in controls to 10 +/- 8% in PHP-Ia. The agonist competition curves were computer fit to a "ternary complex" model for the two-step reaction: H + R + N in equilibrium HR + N in equilibrium HRN. The modeling was consistent with a 60% decrease in the functional concentration of N, and was in good agreement with the biochemically determined decrease in erythrocyte N protein activity. These in vitro findings in erythrocytes taken together with the recent observations that in vivo isoproterenol-stimulated adenylate cyclase activity is decreased in patients with PHP (Carlson, H. E., and A. S. Brickman, 1983, J. Clin. Endocrinol. Metab. 56:1323-1326) are consistent with the notion that N is a bifunctional protein interacting with both R and the adenylate cyclase. It may be that in patients with PHP-Ia a single molecular and genetic defect accounts for both decreased HRN formation and decreased adenylate cyclase activity, whereas in PHP-Ib the biochemical lesion(s) appear not to affect HRN complex formation.