Inferring phylogenies from DNA sequences of unequal base compositions.

A new method for computing evolutionary distances between DNA sequences is proposed. Contrasting with classical methods, the underlying model does not assume that sequence base compositions (A, C, G, and T contents) are at equilibrium, thus allowing unequal base compositions among compared sequences. This makes the method more efficient than the usual ones in recovering phylogenetic trees from sequence data when base composition is heterogeneous within the data set, as we show by using both simulated and empirical data. When applied to small-subunit ribosomal RNA sequences from several prokaryotic or eukaryotic organisms, this method provides evidence for an early divergence of the microsporidian Vairimorpha necatrix in the eukaryotic lineage.

A single-stranded DNA binding protein binds the origin of replication of the duplex kinetoplast DNA.

Replication of the kinetoplast DNA (kDNA) minicircle of trypanosomatids initiates at a conserved 12-nt sequence, 5'-GGGGTTGGTGTA-3', termed the universal minicircle sequence (UMS). A sequence-specific single-stranded DNA-binding protein from Crithidia fasciculata binds the heavy strand of the 12-mer UMS. Whereas this UMS-binding protein (UMSBP) does not bind a duplex UMS dodecamer, it binds the double-stranded kDNA minicircle as well as a duplex minicircle fragment containing the origin-associated UMS. Binding of the minicircle origin region by the single-stranded DNA binding protein suggested the local unwinding of the DNA double helix at this site. Modification of thymine residues at this site by KMnO4 revealed that the UMS resides within an unwound or otherwise sharply distorted DNA at the minicircle origin region. Computer analysis predicts the sequence-directed curving of the minicircle origin region. Electrophoresis of a minicircle fragment containing the origin region in polyacrylamide gels revealed a significantly lower electrophoretic mobility than expected from its length. The fragment anomalous electrophoretic mobility is displayed only in its native conformation and is dependent on temperature and gel porosity, indicating the local curving of the DNA double helix. We suggest that binding of UMSBP at the minicircle origin of replication is possible through local unwinding of the DNA double helix at the UMS site. It is hypothesized here that this local melting is initiated through the untwisting of unstacked dinucleotide sequences at the bent origin site.

Thyroid-stimulating hormone (TSH)-directed induction of the CREM gene in the thyroid gland participates in the long-term desensitization of the TSH receptor.

Thyroid gland function is regulated by the hypothalamic-pituitary axis via the secretion of TSH, according to environmental, developmental, and circadian stimuli. TSH modulates both the secretion of thyroid hormone and gland trophism through interaction with a specific guanine nucleotide-binding protein-coupled receptor (TSH receptor; TSH-R), which elicits the activation of the cAMP-dependent signaling pathway. After TSH stimulation, the levels of TSH-R RNA are known to decrease dramatically within a few hours. This phenomenon ultimately leads to homologous long-term desensitization of the TSH-R. Here we show that TSH drives the induction of the inducible cAMP early repressor (ICER) isoform of the cAMP response element (CRE) modulator gene both in rat thyroid gland and in the differentiated thyroid cell line FRTL-5. The kinetics of ICER protein induction mirrors the down-regulation of TSH-R mRNA. ICER binds to a CRE-like sequence in the TSH-R promoter and represses its expression. Thus, ICER induction by TSH in the thyroid gland represents a paradigm of the molecular mechanism by which pituitary hormones elicit homologous long-term desensitization.

Protein crosslinking studies suggest that Rhizobium meliloti C4-dicarboxylic acid transport protein D, a sigma 54-dependent transcriptional activator, interacts with sigma 54 and the beta subunit of RNA polymerase.

Rhizobium meliloti C4-dicarboxylic acid transport protein D (DCTD) activates transcription by a form of RNA polymerase holoenzyme that has sigma 54 as its sigma factor (referred to as E sigma 54). DCTD catalyzes the ATP-dependent isomerization of closed complexes between E sigma 54 and the dctA promoter to transcriptionally productive open complexes. Transcriptional activation probably involves specific protein-protein interactions between DCTD and E sigma 54. Interactions between sigma 54-dependent activators and E sigma 54 are transient, and there has been no report of a biochemical assay for contact between E sigma 54 and any activator to date. Heterobifunctional crosslinking reagents were used to examine protein-protein interactions between the various subunits of E sigma 54 and DCTD. DCTD was crosslinked to Salmonella typhimurium sigma 54 with the crosslinking reagents succinimidyl 4-(N-maleimidomethyl)cyclohexane-1-carboxylate and N-hydroxysulfosuccinimidyl-4-azidobenzoate. Cys-307 of sigma 54 was identified by site-directed mutagenesis as the residue that was crosslinked to DCTD. DCTD was also crosslinked to the beta subunit of Escherichia coli core RNA polymerase with succinimidyl 4-(N-maleimidomethyl)cyclohexane-1-carboxylate, but not with N-hydroxysulfosuccinimidyl-4-azidobenzoate. These data suggest that interactions of DCTD with sigma 54 and the beta subunit may be important for transcriptional activation and offer evidence for interactions between a sigma 54-dependent activator and sigma 54, as well as the beta subunit of RNA polymerase.

Transcription termination at intrinsic terminators: the role of the RNA hairpin.

Intrinsic termination of transcription in Escherichia coli involves the formation of an RNA hairpin in the nascent RNA. This hairpin plays a central role in the release of the transcript and polymerase at intrinsic termination sites on the DNA template. We have created variants of the lambda tR2 terminator hairpin and examined the relationship between the structure and stability of this hairpin and the template positions and efficiencies of termination. The results were used to test the simple nucleic acid destabilization model of Yager and von Hippel and showed that this model must be modified to provide a distinct role for the rU-rich sequence in the nascent RNA, since a perfect palindromic sequence that is sufficiently long to form an RNA hairpin that could destabilize the entire putative 12-bp RNA-DNA hybrid does not trigger termination at the expected positions. Rather, our results show that both a stable terminator hairpin and the run of 6-8 rU residues that immediately follows are required for effective intrinsic termination and that termination occurs at specific and invariant template positions relative to these two components. Possible structural or kinetic modifications of the simple model are proposed in the light of these findings and of recent results implicating "inchworming" and possible conformational heterogeneity of transcription complexes in intrinsic termination. Thus, these findings argue that the structure and dimensions of the hairpin are important determinants of the termination-elongation decision and suggest that a complete mechanism is likely to involve specific interactions of the polymerase, the RNA terminator hairpin, and, perhaps, the dT-rich template sequence that codes for the run of rU residues at the 3' end of the nascent transcript.

Phosphorylation-dependent conformational changes in OmpR, an osmoregulatory DNA-binding protein of Escherichia coli.

Osmoregulated porin gene expression in Escherichia coli is controlled by the two-component regulatory system EnvZ and OmpR. EnvZ, the osmosensor, is an inner membrane protein and a histidine kinase. EnvZ phosphorylates OmpR, a cytoplasmic DNA-binding protein, on an aspartyl residue. Phospho-OmpR binds to the promoters of the porin genes to regulate the expression of ompF and ompC. We describe the use of limited proteolysis by trypsin and ion spray mass spectrometry to characterize phospho-OmpR and the conformational changes that occur upon phosphorylation. Our results are consistent with a two-domain structure for OmpR, an N-terminal phosphorylation domain joined to a C-terminal DNA-binding domain by a flexible linker region. In the presence of acetyl phosphate, OmpR is phosphorylated at only one site. Phosphorylation induces a conformational change that is transmitted to the C-terminal domain via the central linker. Previous genetic analysis identified a region in the C-terminal domain that is required for transcriptional activation. Our results indicate that this region is within a surface-exposed loop. We propose that this loop contacts the alpha subunit of RNA polymerase to activate transcription. Mass spectrometry also reveals an unusual dephosphorylated form of OmpR, the potential significance of which is discussed.

Toward the therapeutic editing of mutated RNA sequences.

If RNA editing could be rationally directed to mutated RNA sequences, genetic diseases caused by certain base substitutions could be treated. Here we use a synthetic complementary RNA oligonucleotide to direct the correction of a premature stop codon mutation in dystrophin RNA. The complementary RNA oligonucleotide was hybridized to a premature stop codon and the hybrid was treated with nuclear extracts containing the cellular enzyme double-stranded RNA adenosine deaminase. When the treated RNAs were translated in vitro, a dramatic increase in expression of a downstream luciferase coding region was observed. The cDNA sequence data are consistent with deamination of the adenosine in the UAG stop codon to inosine by double-stranded RNA adenosine deaminase. Injection of oligonucleotide-mRNA hybrids into Xenopus embryos also resulted in an increase in luciferase expression. These experiments demonstrate the principle of therapeutic RNA editing.

A yeast gene required for DNA replication encodes a protein with homology to DNA helicases.

A yeast gene has been identified by screening for DNA replication mutants using a permeabilized cell replication assay. The mutant is temperature sensitive for growth and shows a cell cycle phenotype typical of DNA replication mutants. RNA synthesis is normal in the mutant but DNA synthesis ceases upon shift to the nonpermissive temperature. The DNA2 gene was cloned by complementation of the dna2ts gene phenotype. The gene is essential for viability. The gene encodes a 172-kDa protein with characteristic DNA helicase motifs. A hemagglutinin epitope-Dna2 fusion protein was prepared and purified by conventional and immunoaffinity chromatography. The purified protein is a DNA-dependent ATPase and has 3' to 5' DNA helicase activity specific for forked substrates. A nuclease activity that endonucleolytically cleaves DNA molecules having a single-stranded 5' tail adjacent to a duplex region copurifies through all steps with the fusion protein.

Interaction of herpes simplex virus 1 origin-binding protein with DNA polymerase alpha.

The herpes simplex virus 1 (HSV-1) genome encodes seven polypeptides that are required for its replication. These include a heterodimeric DNA polymerase, a single-strand-DNA-binding protein, a heterotrimeric helicase/primase, and a protein (UL9 protein) that binds specifically to an HSV-1 origin of replication (oris). We demonstrate here that UL9 protein interacts specifically with the 180-kDa catalytic subunit of the cellular DNA polymerase alpha-primase. This interaction can be detected by immunoprecipitation with antibodies directed against either of these proteins, by gel mobility shift of an oris-UL9 protein complex, and by stimulation of DNA polymerase activity by the UL9 protein. These findings suggest that enzymes required for cellular DNA replication also participate in HSV-1 DNA replication.

Transcription regulation by inflexibility of promoter DNA in a looped complex.

The gal operon of Escherichia coli is negatively regulated by repressor binding to bipartite operators separated by 11 helical turns of DNA. Synergistic binding of repressor to separate sites on DNA results in looping, with the intervening DNA as a topologically closed domain containing the two promoters. A closed DNA loop of 11 helical turns, which is in-flexible to torsional changes, disables the promoters either by resisting DNA unwinding needed for open complex formation or by impeding the processive DNA contacts by an RNA polymerase in flux during transcription initiation. Interaction between two proteins bound to different sites on DNA modulating the activity of the intervening segment toward other proteins by allostery may be a common mechanism of regulation in DNA-multiprotein complexes.

Identification of the human prostatic carcinoma oncogene PTI-1 by rapid expression cloning and differential RNA display.

Elucidating the relevant genomic changes mediating development and evolution of prostate cancer is paramount for effective diagnosis and therapy. A putative dominant-acting nude mouse prostatic carcinoma tumor-inducing gene, PTI-1, has been cloned that is expressed in patient-derived human prostatic carcinomas but not in benign prostatic hypertrophy or normal prostate tissue. PTI-1 was detected by cotransfecting human prostate carcinoma DNA into CREF-Trans 6 cells, inducing tumors in nude mice, and isolating genes displaying increased expression in tumor-derived cells by using differential RNA display (DD). Screening a human prostatic carcinoma (LNCaP) cDNA library with a 214-bp DNA fragment found by DD permitted the cloning of a full-length 2.0-kb PTI-1 cDNA. Sequence analysis indicates that PTI-1 is a gene containing a 630-bp 5' sequence and a 3' sequence homologous to a truncated and mutated form of human elongation factor 1 alpha. In vitro translation demonstrates that the PTI-1 cDNA encodes a predominant approximately 46-kDa protein. Probing Northern blots with a DNA fragment corresponding to the 5' region of PTI-1 identifies multiple PTI-1 transcripts in RNAs from human carcinoma cell lines derived from the prostate, lung, breast, and colon. In contrast, PTI-1 RNA is not detected in human melanoma, neuroblastoma, osteosarcoma, normal cerebellum, or glioblastoma multiforme cell lines. By using a pair of primers recognizing a 280-bp region within the 630-bp 5' PTI-1 sequence, reverse transcription-PCR detects PTI-1 expression in patient-derived prostate carcinomas but not in normal prostate or benign hypertrophic prostate tissue. In contrast, reverse transcription-PCR detects prostate-specific antigen expression in all of the prostate tissues. These results indicate that PTI-1 may be a member of a class of oncogenes that could affect protein translation and contribute to carcinoma development in human prostate and other tissues. The approaches used, rapid expression cloning with the CREF-Trans 6 system and the DD strategy, should prove widely applicable for identifying and cloning additional human oncogenes.

Binding studies of cationic thymidyl deoxyribonucleic guanidine to RNA homopolynucleotides.

Deoxyribonucleic guanidine is a potential antisense agent that is generated via the replacement of the negative phosphodiester linkages of DNA [--O--(PO2-)--O--] with positively-charged guanidinium (g) linkages [--NH--C(==NH2+)--NH--]. A pentameric thymidyl deoxyribonucleic guanidine molecule [d(Tg)4T-azido] has been shown to base pair specifically to poly(rA) with an unprecedented affinity. Both double and triple strands consisting of one and two equivalents of d(Tg)4T-azido paired with one equivalent of poly(rA) are indicated by thermal denaturation experiments. At an ionic strength of 0.22, the five bases of d(Tg)4T-azido are estimated to dissociate from a double helix with poly(rA) at > 100 degrees C! The effect of ionic strength on thermal denaturation is very pronounced, with stability greatest at low ionic strengths. The method of continuous variation indicates that there is an equilibrium complex with a molar ratio of d(Tg) to r(Ap) or d(Ap) of 2:1. Based on this evidence, models of the structures of d(Tg)9T-azido bound to r(Ap)9A are proposed.

DNA polymerase beta bypasses in vitro a single d(GpG)-cisplatin adduct placed on codon 13 of the HRAS gene.

We have examined the capacity of calf thymus DNA polymerases alpha, beta, delta, and epsilon to perform in vitro translesion synthesis on a substrate containing a single d(GpG)-cisplatin adduct placed on codon 13 of the human HRAS gene. We found that DNA synthesis catalyzed by DNA polymerases alpha, delta, and epsilon was blocked at the base preceding the lesion. Addition of proliferating cell nuclear antigen to DNA polymerase delta and replication protein A to DNA polymerase alpha did not restore their capacity to elongate past the adduct. On the other hand, DNA polymerase beta efficiently bypassed the cisplatin adduct. Furthermore, we observed that DNA polymerase beta was the only polymerase capable of primer extension of a 3'-OH located opposite the base preceding the lesion. Likewise, DNA polymerase beta was able to elongate the arrested replication products of the other three DNA polymerases, thus showing its capacity to successfully compete with polymerases alpha, delta, and epsilon in the stalled replication complex. Our data suggest (i) a possible mechanism enabling DNA polymerase beta to bypass a d(GpG)-cisplatin adduct in vitro and (ii) a role for this enzyme in processing DNA damage in vivo.

A class of single-stranded telomeric DNA-binding proteins required for Rap1p localization in yeast nuclei.

We have identified a class of proteins that bind single-stranded telomeric DNA and are required for the nuclear organization of telomeres and/or telomere-associated proteins. Rlf6p was identified by its sequence similarity to Gbp1p, a single-stranded telomeric DNA-binding protein from Chlamydomonas reinhardtii. Rlf6p and Gbp1p bind yeast single-stranded G-strand telomeric DNA. Both proteins include at least two RNA recognition motifs, which are found in many proteins that interact with single-stranded nucleic acids. Disruption of RLF6 alters the distribution of repressor/activator protein 1 (Rap1p), a telomere-associated protein. In wild-type yeast cells, Rap1p localizes to a small number of perinuclear spots, while in rlf6 cells Rap1p appears diffuse and nuclear. Interestingly, telomere position effect and telomere length control, which require RAP1, are unaffected by rlf6 mutations, demonstrating that Rap1p localization can be uncoupled from other Rap1p-dependent telomere functions. In addition, expression of Chlamydomonas GBP1 restores perinuclear, punctate Rap1p localization in rlf6 mutant cells. The functional complementation of a fungal gene by an algal gene suggests that Rlf6p and Gbp1p are members of a conserved class of single-stranded telomeric DNA-binding proteins that influence nuclear organization. Furthermore, it demonstrates that, despite their unusual codon bias, C. reinhardtii genes can be efficiently translated in Saccharomyces cerevisiae cells.

Rapid RNA polymerase genetics: one-day, no-column preparation of reconstituted recombinant Escherichia coli RNA polymerase.

We present a simple, rapid procedure for reconstitution of Escherichia coli RNA polymerase holoenzyme (RNAP) from individual recombinant alpha, beta, beta', and sigma 70 subunits. Hexahistidine-tagged recombinant alpha subunit purified by batch-mode metal-ion-affinity chromatography is incubated with crude recombinant beta, beta', and sigma 70 subunits from inclusion bodies, and the resulting reconstituted recombinant RNAP is purified by batch-mode metal-ion-affinity chromatography. RNAP prepared by this procedure is indistinguishable from RNAP prepared by conventional methods with respect to subunit stoichiometry, alpha-DNA interaction, catabolite gene activator protein (CAP)-independent transcription, and CAP-dependent transcription. Experiments with alpha (1-235), an alpha subunit C-terminal deletion mutant, establish that the procedure is suitable for biochemical screening of subunit lethal mutants.