Nuclear proteins interacting with an AP-1/CRE-like promoter sequence in the human TNF$\alpha$ gene

Monocyte developmental heterogeneity is reflected at the cellular level by differential activation competence, at the molecular level by differential regulation of gene expression. LPS activates monocytes to produce tumor necrosis factor-$\alpha$ (TNF). Events occurring at the molecular level necessary for TNF regulation have not been elucidated, but depend both on activation signals and the maturation state of the cell: Peripheral blood monocytes produce TNF upon LPS stimulation, but only within the first 72 hours of culture. Expression of c-fos is associated with monocytic differentiation and activation; the fos-associated protein, c-jun, is also expressed during monocyte activation. Increased cAMP levels are associated with down regulation of macrophage function, including LPS-induced TNF transcription. Due to these associations, we studied a region of the TNF promoter which resembles the binding sites for both AP-1(fos/jun) and CRE-binding protein (or ATF) in order to identify potential molecular markers defining activation competent populations of monocytic cells.^ Nuclear protein binding studies using extracts from THP-1 monocytic cells stimulated with LPS, which stimulates, or dexamethasone (Dex) or pentoxyfilline (PTX), which inhibit TNF production, respectively, suggest that a low mobility doublet complex may be involved in regulation through this promoter region. PTX or Dex increase binding of these complexes equivalently over untreated cells; approximately two hours after LPS induction, the upper complex is undetectable. The upper complex is composed of ATF2 (CRE-BP1); the lower is a heterodimer of jun/ATF2. LPS induces c-jun and thus may enhance formation of jun-ATF2 complexes. The simultaneous presence of both complexes may reduce the amount of TNF transcription through competitive binding, while a loss of the upper (ATF2) and/or gain of the lower (jun-ATF2) allow increased transcription. AP-1 elements generally transduce signals involving PKC; the CRE mediates a cAMP response, involving PKA. Thus, this element has the potential of receiving signals through divergent signalling pathways. Our findings also suggest that cAMP-induced inhibition of macrophage functions may occur via down regulation of activation-associated genes through competitive binding of particular cAMP-responsive nuclear protein complexes. ^

Models of DNA sequence evolution

Models of DNA sequence evolution and methods for estimating evolutionary distances are needed for studying the rate and pattern of molecular evolution and for inferring the evolutionary relationships of organisms or genes. In this dissertation, several new models and methods are developed.^ The rate variation among nucleotide sites: To obtain unbiased estimates of evolutionary distances, the rate heterogeneity among nucleotide sites of a gene should be considered. Commonly, it is assumed that the substitution rate varies among sites according to a gamma distribution (gamma model) or, more generally, an invariant+gamma model which includes some invariable sites. A maximum likelihood (ML) approach was developed for estimating the shape parameter of the gamma distribution $(\alpha)$ and/or the proportion of invariable sites $(\theta).$ Computer simulation showed that (1) under the gamma model, $\alpha$ can be well estimated from 3 or 4 sequences if the sequence length is long; and (2) the distance estimate is unbiased and robust against violations of the assumptions of the invariant+gamma model.^ However, this ML method requires a huge amount of computational time and is useful only for less than 6 sequences. Therefore, I developed a fast method for estimating $\alpha,$ which is easy to implement and requires no knowledge of tree. A computer program was developed for estimating $\alpha$ and evolutionary distances, which can handle the number of sequences as large as 30.^ Evolutionary distances under the stationary, time-reversible (SR) model: The SR model is a general model of nucleotide substitution, which assumes (i) stationary nucleotide frequencies and (ii) time-reversibility. It can be extended to SRV model which allows rate variation among sites. I developed a method for estimating the distance under the SR or SRV model, as well as the variance-covariance matrix of distances. Computer simulation showed that the SR method is better than a simpler method when the sequence length $L>1,000$ bp and is robust against deviations from time-reversibility. As expected, when the rate varies among sites, the SRV method is much better than the SR method.^ The evolutionary distances under nonstationary nucleotide frequencies: The statistical properties of the paralinear and LogDet distances under nonstationary nucleotide frequencies were studied. First, I developed formulas for correcting the estimation biases of the paralinear and LogDet distances. The performances of these formulas and the formulas for sampling variances were examined by computer simulation. Second, I developed a method for estimating the variance-covariance matrix of the paralinear distance, so that statistical tests of phylogenies can be conducted when the nucleotide frequencies are nonstationary. Third, a new method for testing the molecular clock hypothesis was developed in the nonstationary case. ^