Biochemical Characterization of Human Exonuclease 1 Protein-protein Interactions

Human Exonuclease 1 (Exo1) plays important roles in numerous DNA metabolic/repair pathways including DNA mismatch repair, DNA double strand break repair, Okazaki fragment maturation etc. The nuclease activity of Exo1 is tightly regulated in vivo. The regulation of Exo1 in different pathways is achieved by interactions with different protein partners. The focus of this dissertation will be on characterization of Exo1 interactions with traditional protein partners and providing experimental evidences for new Exo1 interactions.

Molecular cloning, biochemical assays, collaborative nuclear magnetic resonance and X-ray crystallography have been employed to study Exo1 interactions with protein partners. This work contains: (i) the experimental evidence for new Exo1 interactions, and (ii) the detailed characterization of Exo1 interactions with PCNA, MLH1 and MutSα/β.

Taken together, the research progress presented in this dissertation further advances our understanding of traditional Exo1 interaction network and probably provides new insights to new functions and new regulations of Exo1.

Development and Application of Covalent-Labeling Strategies for the Large-Scale Thermodynamic Analysis of Protein Folding and Ligand Binding

Thermodynamic stability measurements on proteins and protein-ligand complexes can offer insights not only into the fundamental properties of protein folding reactions and protein functions, but also into the development of protein-directed therapeutic agents to combat disease. Conventional calorimetric or spectroscopic approaches for measuring protein stability typically require large amounts of purified protein. This requirement has precluded their use in proteomic applications. Stability of Proteins from Rates of Oxidation (SPROX) is a recently developed mass spectrometry-based approach for proteome-wide thermodynamic stability analysis. Since the proteomic coverage of SPROX is fundamentally limited by the detection of methionine-containing peptides, the use of tryptophan-containing peptides was investigated in this dissertation. A new SPROX-like protocol was developed that measured protein folding free energies using the denaturant dependence of the rate at which globally protected tryptophan and methionine residues are modified with dimethyl (2-hydroxyl-5-nitrobenzyl) sulfonium bromide and hydrogen peroxide, respectively. This so-called Hybrid protocol was applied to proteins in yeast and MCF-7 cell lysates and achieved a ~50% increase in proteomic coverage compared to probing only methionine-containing peptides. Subsequently, the Hybrid protocol was successfully utilized to identify and quantify both known and novel protein-ligand interactions in cell lysates. The ligands under study included the well-known Hsp90 inhibitor geldanamycin and the less well-understood omeprazole sulfide that inhibits liver-stage malaria. In addition to protein-small molecule interactions, protein-protein interactions involving Puf6 were investigated using the SPROX technique in comparative thermodynamic analyses performed on wild-type and Puf6-deletion yeast strains. A total of 39 proteins were detected as Puf6 targets and 36 of these targets were previously unknown to interact with Puf6. Finally, to facilitate the SPROX/Hybrid data analysis process and minimize human errors, a Bayesian algorithm was developed for transition midpoint assignment. In summary, the work in this dissertation expanded the scope of SPROX and evaluated the use of SPROX/Hybrid protocols for characterizing protein-ligand interactions in complex biological mixtures.

Sequence and Structural Determinants of Specificity Differences between Paralogous Transcription Factors

Transcription factors (TFs) control the temporal and spatial expression of target genes by interacting with DNA in a sequence-specific manner. Recent advances in high throughput experiments that measure TF-DNA interactions in vitro and in vivo have facilitated the identification of DNA binding sites for thousands of TFs. However, it remains unclear how each individual TF achieves its specificity, especially in the case of paralogous TFs that recognize distinct target genomic sites despite sharing very similar DNA binding motifs. In my work, I used a combination of high throughput in vitro protein-DNA binding assays and machine-learning algorithms to characterize and model the binding specificity of 11 paralogous TFs from 4 distinct structural families. My work proves that even very closely related paralogous TFs, with indistinguishable DNA binding motifs, oftentimes exhibit differential binding specificity for their genomic target sites, especially for sites with moderate binding affinity. Importantly, the differences I identify in vitro and through computational modeling help explain, at least in part, the differential in vivo genomic targeting by paralogous TFs. Future work will focus on in vivo factors that might also be important for specificity differences between paralogous TFs, such as DNA methylation, interactions with protein cofactors, or the chromatin environment. In this larger context, my work emphasizes the importance of intrinsic DNA binding specificity in targeting of paralogous TFs to the genome.

Identification of the endophilins (SH3p4/p8/p13) as novel binding partners for the beta1-adrenergic receptor.

Several G-protein coupled receptors, such as the beta1-adrenergic receptor (beta1-AR), contain polyproline motifs within their intracellular domains. Such motifs in other proteins are known to mediate protein-protein interactions such as with Src homology (SH)3 domains. Accordingly, we used the proline-rich third intracellular loop of the beta1-AR either as a glutathione S-transferase fusion protein in biochemical "pull-down" assays or as bait in the yeast two-hybrid system to search for interacting proteins. Both approaches identified SH3p4/p8/p13 (also referred to as endophilin 1/2/3), a SH3 domain-containing protein family, as binding partners for the beta1-AR. In vitro and in human embryonic kidney (HEK) 293 cells, SH3p4 specifically binds to the third intracellular loop of the beta1-AR but not to that of the beta2-AR. Moreover, this interaction is mediated by the C-terminal SH3 domain of SH3p4. Functionally, overexpression of SH3p4 promotes agonist-induced internalization and modestly decreases the Gs coupling efficacy of beta1-ARs in HEK293 cells while having no effect on beta2-ARs. Thus, our studies demonstrate a role of the SH3p4/p8/p13 protein family in beta1-AR signaling and suggest that interaction between proline-rich motifs and SH3-containing proteins may represent a previously underappreciated aspect of G-protein coupled receptor signaling.

Development and Characterization of Monovalent and Bivalent RNA Aptamers Targeting the Common Pathway of Coagulation

Anticoagulant agents are commonly used drugs to reduce blood coagulation in acute and chronic clinical settings. Many of these drugs target the common pathway of coagulation because it is critical for thrombin generation and disruption of this portion of the pathway has profound effects on the hemostatic process. Currently available drugs for these indications struggle with balancing desired activity with immunogenicity and poor reversibility or irreversibility in the event of hemorrhage. While improvements are being made with the current drugs, new drugs with better therapeutic indices are needed for surgical intervention and chronic indications to prevent thrombosis from occurring.

A class of therapeutics known as aptamers may be able to meet the need for safer anticoagulant agents. Aptamer are short single-stranded RNA oligonucleotides that adopt specific secondary and tertiary structures based upon their sequence. They can be generated to both enzymes and cofactors because they derive their inhibitory activity by blocking protein-protein interactions, rather than active site inhibition. They inhibit their target proteins with a high level of specificity and bind with high affinity to their target. Additionally, they can be reversed using two different antidote approaches, specific oligonucleotide antidotes, or with cationic, “universal” antidotes. The reversal of their activity is both rapid and durable.

The ability of aptamers to be generated to cofactors has been conclusively proven by generating an aptamer targeting the common pathway coagulation cofactor, Factor V (FV). We developed two aptamers with anticoagulant ability that bind to both FV and FVa, the active cofactor. Both aptamers were truncated to smaller functional sizes and had specific point mutant aptamers developed for use as controls. The anticoagulant activity of both aptamer-mutant pairs was characterized using plasma-based clotting assays and whole blood assays. The mechanism of action resulting in anticoagulant activity was assessed for one aptamer. The aptamer was found to block FVa docking to membrane surfaces, a mechanism not previously observed in any of our other anticoagulant aptamers.

To explore development of aptamers as anticoagulant agents targeting the common pathway for surgical interventions, we fused two anticoagulant aptamers targeting Factor X and prothrombin into a single molecule. The bivalent aptamer was truncated to a minimal size while maintaining robust anticoagulant activity. Characterization of the bivalent aptamer in plasma-based clotting assays indicated we had generated a very robust anticoagulant therapeutic. Furthermore, we were able to simultaneously reverse the activity of both aptamers with a single oligonucleotide antidote. This rapid and complete reversal of anticoagulant activity is not available in the antithrombotic agents currently used in surgery.

Occurrence and Function of Hoogsteen Base Pairs in Nucleic Acids

Nucleic acids (DNA and RNA) play essential roles in the central dogma of biology for the storage and transfer of genetic information. The unique chemical and conformational structures of nucleic acids – the double helix composed of complementary Watson-Crick base pairs, provide the structural basis to carry out their biological functions. DNA double helix can dynamically accommodate Watson-Crick and Hoogsteen base-pairing, in which the purine base is flipped by ~180° degrees to adopt syn rather than anti conformation as in Watson-Crick base pairs. There is growing evidence that Hoogsteen base pairs play important roles in DNA replication, recognition, damage or mispair accommodation and repair. Here, we constructed a database for existing Hoogsteen base pairs in DNA duplexes by a structure-based survey from the Protein Data Bank, and structural analyses based on the resulted Hoogsteen structures revealed that Hoogsteen base pairs occur in a wide variety of biological contexts and can induce DNA kinking towards the major groove. As there were documented difficulties in modeling Hoogsteen or Watson-Crick by crystallography, we collaborated with the Richardsons’ lab and identified potential Hoogsteen base pairs that were mis-modeled as Watson-Crick base pairs which suggested that Hoogsteen can be more prevalent than it was thought to be. We developed solution NMR method combined with the site-specific isotope labeling to characterize the formation of, or conformational exchange with Hoogsteen base pairs in large DNA-protein complexes under solution conditions, in the absence of the crystal packing force. We showed that there are enhanced chemical exchange, potentially between Watson-Crick and Hoogsteen, at a sharp kink site in the complex formed by DNA and the Integration Host Factor protein. In stark contrast to B-form DNA, we found that Hoogsteen base pairs are strongly disfavored in A-form RNA duplex. Chemical modifications N1-methyl adenosine and N1-methyl guanosine that block Watson-Crick base-pairing, can be absorbed as Hoogsteen base pairs in DNA, but rather potently destabilized A-form RNA and caused helix melting. The intrinsic instability of Hoogsteen base pairs in A-form RNA endows the N1-methylation as a functioning post-transcriptional modification that was known to facilitate RNA folding, translation and potentially play roles in the epitranscriptome. On the other hand, the dynamic property of DNA that can accommodate Hoogsteen base pairs could be critical to maintaining the genome stability.

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.

Binding of MetJ repressor to specific and nonspecific DNA and effect of S-adenosylmethionine on these interactions.

We have used analytical ultracentrifugation to characterize the binding of the methionine repressor protein, MetJ, to synthetic oligonucleotides containing zero to five specific recognition sites, called metboxes. For all lengths of DNA studied, MetJ binds more tightly to repeats of the consensus sequence than to naturally occurring metboxes, which exhibit a variable number of deviations from the consensus. Strong cooperative binding occurs only in the presence of two or more tandem metboxes, which facilitate protein-protein contacts between adjacent MetJ dimers, but weak affinity is detected even with DNA containing zero or one metbox. The affinity of MetJ for all of the DNA sequences studied is enhanced by the addition of SAM, the known cofactor for MetJ in the cell. This effect extends to oligos containing zero or one metbox, both of which bind two MetJ dimers. In the presence of a large excess concentration of metbox DNA, the effect of cooperativity is to favor populations of DNA oligos bound by two or more MetJ dimers rather than a stochastic redistribution of the repressor onto all available metboxes. These results illustrate the dynamic range of binding affinity and repressor assembly that MetJ can exhibit with DNA and the effect of the corepressor SAM on binding to both specific and nonspecific DNA.

Receptor and G betagamma isoform-specific interactions with G protein-coupled receptor kinases.

The G protein-coupled receptor (GPCR) kinases (GRKs) phosphorylate and desensitize agonist-occupied GPCRs. GRK2-mediated receptor phosphorylation is preceded by the agonist-dependent membrane association of this enzyme. Previous in vitro studies with purified proteins have suggested that this translocation may be mediated by the recruitment of GRK2 to the plasma membrane by its interaction with the free betagamma subunits of heterotrimeric G proteins (G betagamma). Here we demonstrate that this mechanism operates in intact cells and that specificity is imparted by the selective interaction of discrete pools of G betagamma with receptors and GRKs. Treatment of Cos-7 cells transiently overexpressing GRK2 with a beta-receptor agonist promotes a 3-fold increase in plasma membrane-associated GRK2. This translocation of GRK2 is inhibited by the carboxyl terminus of GRK2, a known G betagamma sequestrant. Furthermore, in cells overexpressing both GRK2 and G beta1 gamma2, activation of lysophosphatidic acid receptors leads to the rapid and transient formation of a GRK/G betagamma complex. That G betagamma specificity exists at the level of the GPCR and the GRK is indicated by the observation that a GRK2/G betagamma complex is formed after agonist occupancy of the lysophosphatidic acid and beta-adrenergic but not thrombin receptors. In contrast to GRK2, GRK3 forms a G betagamma complex after stimulation of all three GPCRs. This G betagamma binding specificity of the GRKs is also reflected at the level of the purified proteins. Thus the GRK2 carboxyl terminus binds G beta1 and G beta2 but not G beta3, while the GRK3 fusion protein binds all three G beta isoforms. This study provides a direct demonstration of a role for G betagamma in mediating the agonist-stimulated translocation of GRK2 and GRK3 in an intact cellular system and demonstrates isoform specificity in the interaction of these components.

Interactions of chromatin context, binding site sequence content, and sequence evolution in stress-induced p53 occupancy and transactivation.

Cellular stresses activate the tumor suppressor p53 protein leading to selective binding to DNA response elements (REs) and gene transactivation from a large pool of potential p53 REs (p53REs). To elucidate how p53RE sequences and local chromatin context interact to affect p53 binding and gene transactivation, we mapped genome-wide binding localizations of p53 and H3K4me3 in untreated and doxorubicin (DXR)-treated human lymphoblastoid cells. We examined the relationships among p53 occupancy, gene expression, H3K4me3, chromatin accessibility (DNase 1 hypersensitivity, DHS), ENCODE chromatin states, p53RE sequence, and evolutionary conservation. We observed that the inducible expression of p53-regulated genes was associated with the steady-state chromatin status of the cell. Most highly inducible p53-regulated genes were suppressed at baseline and marked by repressive histone modifications or displayed CTCF binding. Comparison of p53RE sequences residing in different chromatin contexts demonstrated that weaker p53REs resided in open promoters, while stronger p53REs were located within enhancers and repressed chromatin. p53 occupancy was strongly correlated with similarity of the target DNA sequences to the p53RE consensus, but surprisingly, inversely correlated with pre-existing nucleosome accessibility (DHS) and evolutionary conservation at the p53RE. Occupancy by p53 of REs that overlapped transposable element (TE) repeats was significantly higher (p<10-7) and correlated with stronger p53RE sequences (p<10-110) relative to nonTE-associated p53REs, particularly for MLT1H, LTR10B, and Mer61 TEs. However, binding at these elements was generally not associated with transactivation of adjacent genes. Occupied p53REs located in L2-like TEs were unique in displaying highly negative PhyloP scores (predicted fast-evolving) and being associated with altered H3K4me3 and DHS levels. These results underscore the systematic interaction between chromatin status and p53RE context in the induced transactivation response. This p53 regulated response appears to have been tuned via evolutionary processes that may have led to repression and/or utilization of p53REs originating from primate-specific transposon elements.

Potential To Streamline Heterologous DNA Prime and NYVAC/Protein Boost HIV Vaccine Regimens in Rhesus Macaques by Employing Improved Antigens.

UNLABELLED: In a follow-up to the modest efficacy observed in the RV144 trial, researchers in the HIV vaccine field seek to substantiate and extend the results by evaluating other poxvirus vectors and combinations with DNA and protein vaccines. Earlier clinical trials (EuroVacc trials 01 to 03) evaluated the immunogenicity of HIV-1 clade C GagPolNef and gp120 antigens delivered via the poxviral vector NYVAC. These showed that a vaccination regimen including DNA-C priming prior to a NYVAC-C boost considerably enhanced vaccine-elicited immune responses compared to those with NYVAC-C alone. Moreover, responses were improved by using three as opposed to two DNA-C primes. In the present study, we assessed in nonhuman primates whether such vaccination regimens can be streamlined further by using fewer and accelerated immunizations and employing a novel generation of improved DNA-C and NYVAC-C vaccine candidates designed for higher expression levels and more balanced immune responses. Three different DNA-C prime/NYVAC-C+ protein boost vaccination regimens were tested in rhesus macaques. All regimens elicited vigorous and well-balanced CD8(+)and CD4(+)T cell responses that were broad and polyfunctional. Very high IgG binding titers, substantial antibody-dependent cellular cytotoxicity (ADCC), and modest antibody-dependent cell-mediated virus inhibition (ADCVI), but very low neutralization activity, were measured after the final immunizations. Overall, immune responses elicited in all three groups were very similar and of greater magnitude, breadth, and quality than those of earlier EuroVacc vaccines. In conclusion, these findings indicate that vaccination schemes can be simplified by using improved antigens and regimens. This may offer a more practical and affordable means to elicit potentially protective immune responses upon vaccination, especially in resource-constrained settings. IMPORTANCE: Within the EuroVacc clinical trials, we previously assessed the immunogenicity of HIV clade C antigens delivered in a DNA prime/NYVAC boost regimen. The trials showed that the DNA prime crucially improved the responses, and three DNA primes with a NYVAC boost appeared to be optimal. Nevertheless, T cell responses were primarily directed toward Env, and humoral responses were modest. The aim of this study was to assess improved antigens for the capacity to elicit more potent and balanced responses in rhesus macaques, even with various simpler immunization regimens. Our results showed that the novel antigens in fact elicited larger numbers of T cells with a polyfunctional profile and a good Env-GagPolNef balance, as well as high-titer and Fc-functional antibody responses. Finally, comparison of the different schedules indicates that a simpler regimen of only two DNA primes and one NYVAC boost in combination with protein may be very efficient, thus showing that the novel antigens allow for easier immunization protocols.