Rescue from apoptosis in early (CD34 selected) versus late (CD34 unselected) human hematopoietic cells by VLA-4 and VCAM-1 dependent adhesion to bone marrow stromal cells

The interaction of hematopoietic precursor cell with bone marrow stromal cells is assumed to be important to the survival of hematopoietic precursor cells during hematopoietic cell long-term culture. Early hematopoietic stem cells are preferentially found within the stromal adherent cell fraction in primary long-term bone marrow cultures. The purpose of this dissertation was to understand the molecular mechanisms that govern these interactions for the regulation of survival and proliferation of early versus late hematopoietic cells.^ Monoclonal antibodies to the VLA-4 recognize the alpha4 beta1 integrin receptor on human hematopoietic cells. This monoclonal antibody blocks the adhesion between early hematopoietic progenitor cells (CD34 selected cells) and stromal cells when added to cultures of these cells. Addition of the VLA-4 monoclonal antibody to cultures of stromal cells and CD34 selected cells was shown to induce apoptosis of CD34 selected cells in these CD34 selected cell/stromal cell cocultures, as measured by the terminal deoxynucleotidyl transferase (TdT)-mediated dUTP-biotin nick end-labeling method. In contrast to these experiments with early hematopoietic progenitor cells (CD34+), the level of adhesion between more differentiated cells (unfractionated hematopoietic cells) and stromal cells was not significantly altered by addition of the anti-VLA-4 monoclonal antibody. Similarly, the level of apoptosis of unfractionated hematopoietic cells was not significantly increased by the addition of anti-VLA-4 monoclonal antibody to cultures of the latter cells with stromal cells. The binding of the unfractionated cells is less than that of the CD34 selected. Since there is no difference between the alpha4 beta1 integrin expression level of the early and late myeloid cells, there may be a difference in the functional state of the integrin between the early and late myeloid cells. We also show that CD34+ selected precursor cells proliferate at a higher rate when these cells are plated on recombinant VCAM-1 molecules. These data indicate that the alpha4beta1 integrin receptor (VLA-4) plays a central role in the apoptosis rescue function which results from the anchorage-dependent growth of the CD34 selected early hematopoietic cells on stromal cells. The data suggest that these apoptosis rescue pathways have less significance as the cells mature and become anchorage-independent in their growth. These data should assist in the design of systems for the ex vivo proliferation and transduction of early hematopoietic cells for genetic therapy. ^

STUDIES ON ENERGY METABOLISM IN THE HIPPOCAMPUS AFTER ISCHEMIA

The gerbil model of ischemia was used to determine the effect of carotid occlusion on energy metabolites in cellular layers of discrete regions of the hippocampus and dentate gyrus. Levels of glucose, glycogen, ATP and phosphocreatine (PCr) were unchanged after 1 minute of ischemia. However, 3 minutes of ischemia produced a dramatic decrease in net levels of all metabolites. No additional decrease was observed after 15 minutes of ischemia. Re-establishment of the blood flow for 5 minutes after a 15 minute ischemic episode returned all metabolites to pre-ischemia levels. Concentrations of glucose and glycogen were elevated in sham-operated animals as a function of the pentobarbital anesthetic employed. In other studies, elevated GABA levels (produced by inhibiting GABA-transaminase with (gamma)-vinyl-GABA (GVG)) were found to decrease the rate of utilization of the high-energy phosphate metabolites ATP and PCr in the mouse cortex. In addition, glucose and glycogen levels were increased. Thus, tonic inhibition by GABA produced decreased cellular activity. Additional experiments demonstrated the attenuation of ischemia-induced metabolite depletion in cellular layers of regions of the hippocampus, dentate gyrus and cortex after GVG administration. Under ether, 1 minute of bilateral carotid occlusion produced a dramatic decrease in metabolite levels. After GVG treatment, the decrease was blocked completely for glucose, glycogen and ATP, and partially for PCr. Therefore, GABA-transaminase inhibition produced increased levels of GABA which subsequently decreased cellular activity. The protection against ischemia may have been due to (a)decreased metabolic rate; the available energy stores were utilized at a slower rate, and (b)increased levels of energy substrates; additional supplies available to maintain viability. These data suggest that the functional state of neural tissue can determine the response to metabolic stress. ^

BIOCHEMICAL CHARACTERIZATION OF BINDING PARTNERS OF TWO HSP70 CO-CHAPERONES IN SACCHAROMYCES CEREVISIAE

Cells are exposed to a variety of environmental and physiological changes including temperature, pH and nutrient availability. These changes cause stress to cells, which results in protein misfolding and altered cellular protein homeostasis. How proteins fold into their three-dimensional functional structure is a fundamental biological process with important relevance to human health. Misfolded and aggregated proteins are linked to multiple neurodegenerative diseases, cardiovascular disease and cystic fibrosis. To combat proteotoxic stress, cells deploy an array of molecular chaperones that assist in the repair or removal of misfolded proteins. Hsp70, an evolutionarily conserved molecular chaperone, promotes protein folding and helps maintain them in a functional state. Requisite co-chaperones, including nucleotide exchange factors (NEFs) strictly regulate and serve to recruit Hsp70 to distinct cellular processes or locations. In yeast and human cells, three structurally non-related cytosolic NEFs are present: Sse1 (Hsp110), Fes1 (HspBP1) and Snl1 (Bag-1). Snl1 is unique among the cytosolic NEFs as it is localized at the ER membrane with its Hsp70 binding (BAG) domain exposed to the cytosol. I discovered that Snl1 distinctly interacts with assembled ribosomes and several lines of evidence indicate that this interaction is both independent of and concurrent with binding to Hsp70 and is not dependent on membrane localization. The ribosome-binding site is identified as a short lysine-rich motif within the amino terminus of the Snl1 BAG domain distinct from the Hsp70 interaction region. In addition, I demonstrate ribosome association with the Snl1 homolog in the pathogenic fungus, Candida albicans and localize this putative NEF to a perinuclear/ER membrane, suggesting functional conservation in fungal BAG domain-containing proteins. As a first step in determining specific domain architecture in fungal BAG proteins, I present the preliminary steps of protein purification and analysis of the minimal Hsp70 binding region in in both S.cerevisiae and C. albicans Snl1. Contrary to previous in vitro evidence which showed the Fes1 NEF to interact with both cytosolic Hsp70s, Ssa and Ssb, Fes1 is shown to interact specifically with Ssa when expressed under normal cellular conditions in S. cerevisiae. This is the first reported evidence of Hsp70 binding selectivity for a cytosolic NEF, and suggests a possible mechanism to achieve specificity in Hsp70-dependent functions. Taken together, the work presented in this dissertation highlights the striking divergence among Hsp70 co-chaperones in selecting binding partners, which may correlate with their specific roles in the cell.

AMPA RECEPTOR ALLOSTERISM: MEASUREMENT OF THE CONFORMATIONAL CHANGES IN THE LIGAND BINDING DOMAIN OF A FUNCTIONAL RECEPTOR

Ionotropic glutamate receptors are important excitatory neurotransmitter receptors in the mammalian central nervous system that have been implicated in a number of neuropathologies such as epilepsy, ischemia, and amyotrophic lateral sclerosis. Glutamate binding to an extracellular ligand binding domain initiates a series of structural changes that leads to the formation of a cation selective transmembrane channel, which consequently closes due to desensitization of the receptor. The crystal structures of the AMPA subtype of the glutamate receptor have been particularly useful in providing initial insight into the conformational changes in the ligand binding domain; however, these structures are limited by crystallographic constraint. To gain a clear picture of how agonist binding is coupled to channel activation and desensitization, it is essential to study changes in the ligand binding domain in a dynamic, physiological state. In this dissertation, a technique called Luminescence Resonance Energy Transfer was used to determine the conformational changes associated with activation and desensitization in a functional AMPA receptor (ÄN*-AMPA) that contains the ligand binding domain and transmembrane segments; ÄN*-AMPA has been modified such that fluorophores can be introduced at specific sites to serve as a readout of cleft closure or to establish intersubunit distances. Previous structural studies of cleft closure of the isolated ligand binding domain in conjunction with functional studies of the full receptor suggest that extent of cleft closure correlates with extent of activation. Here, LRET has been used to show that a similar relationship between cleft closure and activation is observed in the “full length” receptor showing that the isolated ligand binding domain is a good model of the domain in the full length receptor for changes within a subunit. Similar LRET investigations were used to study intersubunit distances specifically to probe conformational changes between subunits within a dimer in the tetrameric receptor. These studies show that the dimer interface is coupled in the open state, and decoupled in the desensitized state, similar to the isolated ligand binding domain crystal structure studies. However, we show that the apo state dimer interface is not pre-formed as in the crystal structure, hence suggesting a mechanism for functional transitions within the receptor based on LRET distances obtained.

Extended kalman filter for estimation of parameters in nonlinear state-space models of biochemical networks.

It is system dynamics that determines the function of cells, tissues and organisms. To develop mathematical models and estimate their parameters are an essential issue for studying dynamic behaviors of biological systems which include metabolic networks, genetic regulatory networks and signal transduction pathways, under perturbation of external stimuli. In general, biological dynamic systems are partially observed. Therefore, a natural way to model dynamic biological systems is to employ nonlinear state-space equations. Although statistical methods for parameter estimation of linear models in biological dynamic systems have been developed intensively in the recent years, the estimation of both states and parameters of nonlinear dynamic systems remains a challenging task. In this report, we apply extended Kalman Filter (EKF) to the estimation of both states and parameters of nonlinear state-space models. To evaluate the performance of the EKF for parameter estimation, we apply the EKF to a simulation dataset and two real datasets: JAK-STAT signal transduction pathway and Ras/Raf/MEK/ERK signaling transduction pathways datasets. The preliminary results show that EKF can accurately estimate the parameters and predict states in nonlinear state-space equations for modeling dynamic biochemical networks.

Functional requirements of histone deacetylases in Tup1-Ssn6 repression

The corepressor complex Tup1-Ssn6 regulates many classes of genes in yeast including cell type specific, glucose repressible, and DNA damage inducible. Tup1 and Ssn6 are recruited to target promoters through their interactions with specific DNA binding proteins such as α2, Mig1, and Crt1. Most promoters that are repressed by this corepressor complex exhibit a high degree of nucleosomal organization. This chromatin domain occludes transcription factor access to the promoter element resulting in gene repression. Previous work indicated that Tup1 interacts with underacetylated isoforms of H3 and H4, and that mutation of these histones synergistically compromises repression. These studies predict that Tup1-hypoacetyalted histone interaction is important to the repression mechanism, and in vivo hyperacetylation might compromise the corepressors ability to repress target genes. ^ One way to alter histone acetylation levels in vivo is to alter the balance between histone acetyltransferases and histone deacetylases. To date five histone deacetylases (HDACs) have been identified in yeast Rpd3, Hos1, Hos2, Hos3 and Hda1. Deletion of single or double HDAC genes had little to no effect on Tup1-Ssn6 repression, but simultaneous deletion of three specific activities Rpd3, Hos1, and Hos2 abolished repression in vivo. Promoter regions of Tup1-Ssn6 target genes in these triple deacetylase mutant cells are dramatically hyperacetylated in both H3 and H4. Examination of bulk histone acetylation levels showed that this specific HDAC triple mutant combination (rpd3 hos1 hos2) caused a dramatic and concomitant hyperacetylation of both H3 and H4. The loss of repression in the rpd3 hos1 hos2 cells, but not in other mutants, is consistent with previous observations, which indicate that histones provide redundant functions in the repression mechanism and that high levels of acetylation are required to prevent Tup1 binding. Investigation into a potential direct interaction between the Tup1-Ssn6 corepressor complex and one or more HDAC activities showed that both Rpd3 and Hos2 interact with the corepressor complex in vivo. These findings indicate that Tup1-Ssn6 repression involves the recruitment of histone deacetylase activities to target promoters, where they locally deacetylate histone residues promoting Tup1-histone tail interaction to initiate and/or maintain the repressed state. ^

Functional data analysis approaches for genotype-phenotype association studies from next-generation sequencing

Next-generation DNA sequencing platforms can effectively detect the entire spectrum of genomic variation and is emerging to be a major tool for systematic exploration of the universe of variants and interactions in the entire genome. However, the data produced by next-generation sequencing technologies will suffer from three basic problems: sequence errors, assembly errors, and missing data. Current statistical methods for genetic analysis are well suited for detecting the association of common variants, but are less suitable to rare variants. This raises great challenge for sequence-based genetic studies of complex diseases.^ This research dissertation utilized genome continuum model as a general principle, and stochastic calculus and functional data analysis as tools for developing novel and powerful statistical methods for next generation of association studies of both qualitative and quantitative traits in the context of sequencing data, which finally lead to shifting the paradigm of association analysis from the current locus-by-locus analysis to collectively analyzing genome regions.^ In this project, the functional principal component (FPC) methods coupled with high-dimensional data reduction techniques will be used to develop novel and powerful methods for testing the associations of the entire spectrum of genetic variation within a segment of genome or a gene regardless of whether the variants are common or rare.^ The classical quantitative genetics suffer from high type I error rates and low power for rare variants. To overcome these limitations for resequencing data, this project used functional linear models with scalar response to develop statistics for identifying quantitative trait loci (QTLs) for both common and rare variants. To illustrate their applications, the functional linear models were applied to five quantitative traits in Framingham heart studies. ^ This project proposed a novel concept of gene-gene co-association in which a gene or a genomic region is taken as a unit of association analysis and used stochastic calculus to develop a unified framework for testing the association of multiple genes or genomic regions for both common and rare alleles. The proposed methods were applied to gene-gene co-association analysis of psoriasis in two independent GWAS datasets which led to discovery of networks significantly associated with psoriasis.^

THE EFFECT OF PERSON-ENVIRONMENT FIT ON THE PSYCHOLOGICAL WELL-BEING OF PSYCHIATRIC AIDES WORKING IN A STATE HOSPITAL (COPING RESOURCES, SOCIAL SUPPORT, OCCUPATIONAL STRESS, LOCUS-OF-CONTROL)

This study addressed two purposes: (1) to determine the effect of person-environment fit on the psychological well-being of psychiatric aides and (2) to determine what role the coping resources of social support and control have on the above relationship. Two hundred and ten psychiatric aides working in a state hospital in Texas responded to a questionnaire pertaining to these issues.^ Person-environment fit, as a measure of occupational stress, was assessed through a modified version of the Work Environment Scale (WES). The WES subscales used in this study were: involvement, autonomy, job pressure, job clarity, and physical comfort. Psychological well-being was measured with the General Well-Being Schedule which was developed by the National Center for Health Statistics. Co-worker and supervisor support were measured through the WES and finally, control was assessed through Rotter's Locus of Control Scale.^ The results of this study were as follows: (1) all person-environment (p-e) dimensions appeared to have linear relationships with psychological well-being; (2) the p-e fit - well-being relationship did not appear to be confounded by demographic factors; (3) all p-e fit dimensions were significantly related to well-being except for autonomy; (4) p-e fit was more strongly related to well-being than the environmental measure alone; (5) supervisor support and non-work related support were found to have additive effects on the relationship between p-e fit and well-being, however no interaction or buffering effects were observed; (6) locus of control was found to have additive effects in the prediction of well-being and showed interactive effects with work pressure, involvement and physical comfort; and (7) the testing of the overall study model which included many of the components mentioned above yielded an R('2) = .27.^ Implications of these findings are discussed, future research suggested and applications proposed. ^

NONLINEAR TIME SERIES/SYSTEM ANALYSIS (STATE-SPACE, DYNAMICS, INTERVENTION, RESILIENCE, KALMAN)

A discussion of nonlinear dynamics, demonstrated by the familiar automobile, is followed by the development of a systematic method of analysis of a possibly nonlinear time series using difference equations in the general state-space format. This format allows recursive state-dependent parameter estimation after each observation thereby revealing the dynamics inherent in the system in combination with random external perturbations.^ The one-step ahead prediction errors at each time period, transformed to have constant variance, and the estimated parametric sequences provide the information to (1) formally test whether time series observations y(,t) are some linear function of random errors (ELEM)(,s), for some t and s, or whether the series would more appropriately be described by a nonlinear model such as bilinear, exponential, threshold, etc., (2) formally test whether a statistically significant change has occurred in structure/level either historically or as it occurs, (3) forecast nonlinear system with a new and innovative (but very old numerical) technique utilizing rational functions to extrapolate individual parameters as smooth functions of time which are then combined to obtain the forecast of y and (4) suggest a measure of resilience, i.e. how much perturbation a structure/level can tolerate, whether internal or external to the system, and remain statistically unchanged. Although similar to one-step control, this provides a less rigid way to think about changes affecting social systems.^ Applications consisting of the analysis of some familiar and some simulated series demonstrate the procedure. Empirical results suggest that this state-space or modified augmented Kalman filter may provide interesting ways to identify particular kinds of nonlinearities as they occur in structural change via the state trajectory.^ A computational flow-chart detailing computations and software input and output is provided in the body of the text. IBM Advanced BASIC program listings to accomplish most of the analysis are provided in the appendix. ^

Functional Characterization of AtxA, the Bacillus anthracis Virulence Regulator

Coordinated expression of virulence genes in Bacillus anthracis occurs via a multi-faceted signal transduction pathway that is dependent upon the AtxA protein. Intricate control of atxA gene transcription and AtxA protein function have become apparent from studies of AtxA-induced synthesis of the anthrax toxin proteins and the poly-D-glutamic acid capsule, two factors with important roles in B. anthracis pathogenesis. The amino-terminal region of the AtxA protein contains winged-helix (WH) and helix-turn-helix (HTH) motifs, structural features associated with DNA-binding. Using filter binding assays, I determined that AtxA interacted non-specifically at a low nanomolar affinity with a target promoter (Plef) and AtxA-independent promoters. AtxA also contains motifs associated with phosphoenolpyruvate: sugar phosphotransferase system (PTS) regulation. These PTS-regulated domains, PRD1 and PRD2, are within the central amino acid sequence. Specific histidines in the PRDs serve as sites of phosphorylation (H199 and H379). Phosphorylation of H199 increases AtxA activity; whereas, H379 phosphorylation decreases AtxA function. For my dissertation, I hypothesized that AtxA binds target promoters to activate transcription and that DNA-binding activity is regulated via structural changes within the PRDs and a carboxy-terminal EIIB-like motif that are induced by phosphorylation and ligand binding. I determined that AtxA has one large protease-inaccessible domain containing the PRDs and the carboxy-terminal end of the protein. These results suggest that AtxA has a domain that is distinct from the putative DNA-binding region of the protein. My data indicate that AtxA activity is associated with AtxA multimerization. Oligomeric AtxA was detected when co-affinity purification, non-denaturing gel electrophoresis, and bis(maleimido)hexane (BMH) cross-linking techniques were employed. I exploited the specificity of BMH for cysteine residues to show that AtxA was cross-linked at C402, implicating the carboxy-terminal EIIB-like region in protein-protein interactions. In addition, higher amounts of the cross-linked dimeric form of AtxA were observed when cells were cultured in conditions that promote toxin gene expression. Based on the results, I propose that AtxA multimerization requires the EIIB-like motif and multimerization of AtxA positively impacts function. I investigated the role of the PTS in the function of AtxA and the impact of phosphomimetic residues on AtxA multimerization. B. anthracis Enzyme I (EI) and HPr did not facilitate phosphorylation of AtxA in vitro. Moreover, markerless deletion of ptsHI in B. anthracis did not perturb AtxA function. Taken together, these results suggest that proteins other than the PTS phosphorylate AtxA. Point mutations mimicking phosphohistidine (H to D) and non-phosphorylated histidine (H to A) were tested for an impact on AtxA activity and multimerization. AtxA H199D, AtxA H199A, and AtxA H379A displayed multimerization phenotypes similar to that of the native protein, whereas AtxA H379D was not susceptible to BMH cross-linking or co-affinity purification with AtxA-His. These data suggest that phosphorylation of H379 may decrease AtxA activity by preventing AtxA multimerization. Overall, my data support the following model of AtxA function. AtxA binds to target gene promoters in an oligomeric state. AtxA activity is increased in response to the host-related signal bicarbonate/CO2 because this signal enhances AtxA multimerization. In contrast, AtxA activity is decreased by phosphorylation at H379 because multimerization is inhibited. Future studies will address the interplay between bicarbonate/CO2 signaling and phosphorylation on AtxA function.

CONFORMATIONAL DYNAMICS OF K-RAS AND H-RAS PROTEINS: IS THERE FUNCTIONAL SPECIFICITY AT THE CATALYTIC DOMAIN?

Ras proteins serve as crucial signaling modulators in cell proliferation through their ability to hydrolyze GTP and exist in a GTP “on” state and GTP “off” state. There are three different human Ras isoforms: H-ras, N-ras and K-ras (4A and 4B). Although their sequence identity is very high at the catalytic domain, these isoforms differ in their ability to activate different effectors and hence different signaling pathways. Much of the previous work on this topic has attributed this difference to the hyper variable region of Ras proteins, which contains most of the sequence variance among the isoforms and encodes specificity for differential distribution in the membrane. However, we hypothesize that sequence variation on lobe II of Ras catalytic domain alters dynamics and leads to differential preference for different effectors or modulators. In this work, we used all atom molecular dynamics to analyze the dynamics in the catalytic domain of H-ras and K-ras. We have also analyzed the dynamics of a transforming mutant of H-ras and K-ras and further studied the dynamics of an effectorselective mutant of H-ras. Collectively we have determined that wild type K-ras is more dynamic than H-ras and that the structure of the effector binding loop more closely resembles that of the T35S Raf-selective mutant, possibly giving us a new view and insight into the v mode of effector specificity. Furthermore we have determined that specific mutations at the same location perturb the conformational equilibrium differently in H-ras and K-ras and that an enhanced oncogenic potential may arise from different structural perturbations for each point mutation of a specific isoform.

Functional analysis of GNA -1, a Galphai superfamily member found in the filamentous fungus Neurospora crassa

Heterotrimeric GTP-binding proteins, G proteins, are integral components of eukaryotic signaling systems linking extracellular signals to intracellular responses. Through coupling to seven-transmembrane helix receptors, G proteins convey primary signaling events into multi-leveled cascades of intracellular activity by regulating downstream enzymes, collectively called effectors. The effector enzymes regulated by G proteins include adenylyl cyclase, cAMP phosphodiesterase, phospolipase C-β, mitogen-activated protein kinases, and ion channels. ^ Neurospora crassa is a multicellular, filamentous fungus that is capable of both asexual and sexual reproduction by elaboration of specialized, developmentally controlled structures that give rise to either asexual or sexual spores, respectively. N. crassa possesses at least three heterotrimeric Gα proteins (GNA-1–3) and one Gβ subunit (GNB-1). GNA-1 was the first microbial protein that could be classified in the Gαi superfamily based on its amino acid identity and demonstration that it is a substrate for ADP-ribosylation by pertussis toxin. ^ Experiments were designed to identify the signal transduction pathways and the effector enzymes regulated by GNA-1. Targeted gene-replacement of gna-1 revealed that GNA-1 controls multiple developmental pathways including both asexual and sexual reproduction, maintenance of growth, and resistance to osmotic stress. The Gαi and Gαz members of the Gαi superfamily negatively regulate adenylyl cyclase activity in mammalian cells; therefore, adenylyl cyclase and cAMP levels were measured in Δgna-1 strains and also in strains that were deleted for both gna-1 and gna-2, a second Gα in N. crassa shown to have overlapping functions with GNA-1. Direct measurements of adenylyl cyclase activity revealed that GNA-1, but not GNA-2, was responsible for GTP-stimulated adenylyl cyclase activity in N. crassa. Furthermore, anti-GNA-1 IgG could specifically inhibit GTP-stimulated adenylyl cyclase activity in wild-type strain extracts. These studies also provided evidence that N. crassa possesses feedback mechanisms that control steady-state cAMP levels through indirect regulation of cAMP-phosphodiesterase activity; mutations in gna-1 and gna-2 were additive in their effect on lowering cAMP-phosphodiesterase activity under growth conditions where steady-state cAMP levels were normal but GTP-stimulated adenylyl cyclase activity was reduced 90% in comparison to control strains. ^ Genetic and biochemical epistasis experiments utilizing a Δ gna-1 cr-1 mutant suggest that GNA-1 is essential for female fertility in a cAMP-independent pathway. Furthermore, deletion of gna-1 in a cr-1 background exacerbated many of the defects already observed in the cr-1 strain including more severe growth restriction and developmental defects. However, deletion of gna-1 had no effect on the increased thermotolerance of cr-1, which has been attributed to loss of cAMP. cr-1 possesses GNA-1 protein, and crude membrane fractions from this strain reconstituted GTP-stimulated adenylyl cyclase activity in Δgna-1 membrane fractions. These studies provide direct evidence for the involvement of Gα proteins in the regulation of adenylyl cyclase activity in eukaryotic microbes. ^