The corrected nucleotide sequences of the TaqI restriction and modification enzymes reveal a thirteen-codon overlap.

The nucleotide sequence of the genes encoding methyltransferase TaqI (M.TaqI) and restriction endonuclease TaqI (R.TaqI) with the recognition sequence, TCGA, were analyzed in clones isolated from independent libraries. The genes, originally reported as 363 and 236 codons long [Slatko et al., Nucleic Acids Res. 15 (1987) 9781-9796] were redetermined as 421 and 263 codons long, respectively. The C terminus of the taqIM gene overlaps the N terminus of the taqIR gene by 13 codons, as observed with the isoschizomeric TthHB8I restriction-modification system [Barany et al., Gene 112 (1992) 13-20]. Removal of the overlapping codons did not interfere with in vivo M.TaqI activity. We postulate the overlap plays a role in regulating taqIR expression.     

Identification of the binding site for the extrahelical target base in N6-adenine DNA methyltransferases by photo-cross-linking with duplex oligodeoxyribonucleotides containing 5-iodouracil at the target position.

DNA methyltransferases flip their target bases out of the DNA double helix for catalysis. Base flipping of C5-cytosine DNA methyltransferases was directly observed in the protein-DNA cocrystal structures of M.HhaI and M.HaeIII. Indirect structural evidence for base flipping of N6-adenine and N4-cytosine DNA methyltransferases was obtained by modeling DNA into the three-dimensional structures of M.TaqI and M.PvuII in complex with the cofactor. In addition, biochemical evidence of base flipping was reported for different N6-adenine DNA methyltransferases. As no protein-DNA cocrystal structure for the related N6-adenine and N4-cytosine DNA methyltransferases is available, we used light-induced photochemical cross-linking to identify the binding site of the extrahelical target bases. The N6-adenine DNA methyltransferases M.TaqI and M.CviBIII, which both methylate adenine within the double-stranded 5'-TCGA-3' DNA sequence, were photo-cross-linked to duplex oligodeoxyribonucleotides containing 5-iodouracil at the target position in 50-60% and almost quantitative yield, respectively. Proteolytic fragmentation of the M. CviBIII-DNA complex followed by Edman degradation and electrospray ionization mass spectrometry indicates photo-cross-linking to tyrosine 122. In addition, the mutant methyltransferases M. TaqI/Y108A and M.TaqI/F196A were photo-cross-linked with 6-fold and 2-fold reduced efficiency, respectively, which suggests that tyrosine 108 is the primary site of modification in M.TaqI. Our results indicate a close proximity between the extrahelical target base and tyrosine 122 in M.CviBIII or tyrosine 108 in M.TaqI. As both residues belong to the conserved motif IV ((N/D/S)(P/I)P(Y/F/W)) found in all N6-adenine and N4-cytosine DNA as well as in N6-adenine RNA methyltransferases, a similar spatial relationship between the target bases and the aromatic amino acid residue within motif IV is expected for all these methyltransferases.     

Photofootprinting of drug-binding sites on DNA using diazo- and azido-9-aminoacridine derivatives

It is demonstrated that DNA photofootprinting analysis of the intercalating depsipeptide echinomycin, and the minor groove-binders distamicyn, 4',6-diamidino-2-phenylindole (DAPI) and Hoechst 33258 can be performed using 9-[6-(2-diazocyclopentadienylcarbonyloxy)hexylamino]acridine (DHA) [Nielsen et al. (1988) Nucleic Acids Res. 16, 3877-3888] or 2-methoxy-6-azido-9-aminoacridine (MAA) [Jeppesen et al. (1988) Nucleic Acids Res. 16, 5755-5770]. Both the extent of the drug-binding sites and their relative strength can be determined with either reagent. DNA has the advantage of giving virtually sequence-uniform DNA photocleavage. On the other hand, structural changes in the DNA are detected by MAA. Using the 232-base-pair EcoRI-PvuII pUC19 restriction fragment, it is found that cleavage protection by distamycin, DAPI and Hoechst 33258 all require an (A.T)4 sequence, whereas protection by echinomycin was confined to a G + C-rich 8-base-pair region.     

Characterization of two novel human retropseudogenes related to the calmodulin-encoding gene, CaMII

Two human genomic clones (lambda hg22 and lambda hg29), containing two novel calmodulin (CaM) retropseudogenes, were isolated and characterized. The two pseudogenes show high similarity with the human CaMII cDNA, hCE1 [SenGupta et al., J. Biol. Chem. 262 (1987) 16663-16670] and the CaMII-type retropseudogene, hCE2 (CaMII-psi 1) [SenGupta et al., Nucleic Acids Res. 17 (1989) 2868]. One of them, in clone lambda hg22 (CaMII-psi 2), shows all the characteristics of a processed pseudogene. In clone lambda hg29 (CaMII-psi 3), however, an Alu repetitive sequence was detected immediately upstream from the ancestral 5'-untranslated region. Downstream from the truncated 3'-untranslated region, three additional copies of Alu repetitive sequences flanking about 750 nucleotides of unknown origin were found. Such a processed retropseudogene flanked by multiple Alu repeats may be a target for further recombination events. The three retropseudogenes CaMII-psi 1, CaMII-psi 2 and CaMII-psi 3 are estimated to be about 49, 21 and 25 million years old, respectively.     

The Chinese hamster Alu-equivalent sequence: a conserved highly repetitious, interspersed deoxyribonucleic acid sequence in mammals has a structure suggestive of a transposable element.

A consensus sequence has been determined for a major interspersed deoxyribonucleic acid repeat in the genome of Chinese hamster ovary cells (CHO cells). This sequence is extensively homologous to (i) the human Alu sequence (P. L. Deininger et al., J. Mol. Biol., in press), (ii) the mouse B1 interspersed repetitious sequence (Krayev et al., Nucleic Acids Res. 8:1201-1215, 1980) (iii) an interspersed repetitious sequence from African green monkey deoxyribonucleic acid (Dhruva et al., Proc. Natl. Acad. Sci. U.S.A. 77:4514-4518, 1980) and (iv) the CHO and mouse 4.5S ribonucleic acid (this report; F. Harada and N. Kato, Nucleic Acids Res. 8:1273-1285, 1980). Because the CHO consensus sequence shows significant homology to the human Alu sequence it is termed the CHO Alu-equivalent sequence. A conserved structure surrounding CHO Alu-equivalent family members can be recognized. It is similar to that surrounding the human Alu and the mouse B1 sequences, and is represented as follows: direct repeat-CHO-Alu-A-rich sequence-direct repeat. A composite interspersed repetitious sequence has been identified. Its structure is represented as follows: direct repeat-residue 47 to 107 of CHO-Alu-non-Alu repetitious sequence-A-rich sequence-direct repeat. Because the Alu flanking sequences resemble those that flank known transposable elements, we think it likely that the Alu sequence dispersed throughout the mammalian genome by transposition.     

Changing the hydrogen-bonding potential in the DNA binding site of EcoRI by site-directed mutagenesis drastically reduces the enzymatic activity, not, however, the preference of this restriction endonuclease for cleavage within the site-GAATTC-

According to the X-ray structure analysis of an EcoRI-oligodeoxynucleotide complex [McClarin et al. (1986) Science 234, 1526], sequence specificity is mediated by 12 hydrogen bonds, 6 from each of the two identical subunits of the dimeric enzyme to the recognition site -GAATTC-: Arg200 forms two hydrogen bonds with guanine, while Glu144 and Arg145 form four hydrogen bonds to adjacent adenine residues. Changing the hydrogen-bonding potential at the recognition site without perturbing the rest of the interface should lead to the recognition of degenerate sequences [Rosenberg et al. (1987) in Protein Engineering (Oxender, D. L., & Fox, C. F., Eds.) pp 237-250, Liss, New York]. We have shown previously that replacing Glu144 by Gln and Arg145 by Lys affects the activity of the enzyme, not, however, its specificity [Wolfes et al. (1986) Nucleic Acids Res. 14, 9063]. We show now that also the mutation of Arg200 to Lys, the double mutation Glu144Arg145 to GlnLys, and the triple mutation Glu144Arg145Arg200 to GlnLysLys do not lead to a detectable degeneracy of the specificity of cleavage by EcoRI but significantly impair the catalytic activity of this enzyme. A detailed analysis of the steady-state kinetics of cleavage of pUC8 DNA and a tridecadeoxynucleotide substrate demonstrates that the reduction in activity for all DNA binding site mutants investigated so far is mainly due to a decrease in kcat, with the exception of the Arg200 to Lys mutant, which is only impaired in its KM.(ABSTRACT TRUNCATED AT 250 WORDS)     

DNA binding properties in vivo and target recognition domain sequence alignment analyses of wild-type and mutant RsrI [N6-adenine] DNA methyltransferases

A genetic selection method, the P22 challenge-phage assay, was used to characterize DNA binding in vivo by the prokaryotic beta class [N:6-adenine] DNA methyltransferase M.RSR:I. M.RSR:I mutants with altered binding affinities in vivo were isolated. Unlike the wild-type enzyme, a catalytically compromised mutant, M.RSR:I (L72P), demonstrated site-specific DNA binding in vivo. The L72P mutation is located near the highly conserved catalytic motif IV, DPPY (residues 65-68). A double mutant, M.RSR:I (L72P/D173A), showed less binding in vivo than did M.RSR:I (L72P). Thus, introduction of the D173A mutation deleteriously affected DNA binding. D173 is located in the putative target recognition domain (TRD) of the enzyme. Sequence alignment analyses of several beta class MTases revealed a TRD sequence element that contains the D173 residue. Phylogenetic analysis suggested that divergence in the amino acid sequences of these methyltransferases correlated with differences in their DNA target recognition sequences. Furthermore, MTases of other classes (alpha and gamma) having the same DNA recognition sequence as the beta class MTases share related regions of amino acid sequences in their TRDs.     

The Stanford Microarray Database

The Stanford Microarray Database (SMD) stores raw and normalized data from microarray experiments, and provides web interfaces for researchers to retrieve, analyze and visualize their data. The two immediate goals for SMD are to serve as a storage site for microarray data from ongoing research at Stanford University, and to facilitate the public dissemination of that data once published, or released by the researcher. Of paramount importance is the connection of microarray data with the biological data that pertains to the DNA deposited on the microarray (genes, clones etc.). SMD makes use of many public resources to connect expression information to the relevant biology, including SGD [Ball,C.A., Dolinski,K., Dwight,S.S., Harris,M.A., Issel-Tarver,L., Kasarskis,A., Scafe,C.R., Sherlock,G., Binkley,G., Jin,H. et al. (2000) Nucleic Acids Res., 28, 77-80], YPD and WormPD [Costanzo,M.C., Hogan,J.D., Cusick,M.E., Davis,B.P., Fancher,A.M., Hodges,P.E., Kondu,P., Lengieza,C., Lew-Smith,J.E., Lingner,C. et al. (2000) Nucleic Acids Res., 28, 73-76], Unigene [Wheeler,D.L., Chappey,C., Lash,A.E., Leipe,D.D., Madden,T.L., Schuler,G.D., Tatusova,T.A. and Rapp,B.A. (2000) Nucleic Acids Res., 28, 10-14], dbEST [Boguski,M.S., Lowe,T.M. and Tolstoshev,C.M. (1993) Nature Genet., 4, 332-333] and SWISS-PROT [Bairoch,A. and Apweiler,R. (2000) Nucleic Acids Res., 28, 45-48] and can be accessed at http://genome-www.stanford.edu/microarray.     

Proton NMR studies of [N-MeCys3,N-MeCys7]TANDEM binding to DNA oligonucleotides: sequence-specific binding at the TpA site.

[N-MeCys3,N-MeCys7]TANDEM, an undermethylated analogue of Triostin A, contains two N-methyl groups on the cysteine residues only. Footprinting results showed that [N-MeCys3,N-MeCys7]TANDEM binds strongly to DNA rich in A.T residues [Low, C. M. L., Fox, K. R., Olsen, R. K., & Waring, M. J. (1986) Nucleic Acids Res. 14, 2015-2033]. However, it was not known whether specific binding of [N-MeCys3,N-MeCys7]TANDEM requires a TpA step or an ApT step. In 1:1 saturated complexes with the octamers [d(GGATATCC)]2 and [d(GGTTAACC)]2, [N-MeCys3,N-MeCys7]TANDEM binds to each octamer as a bis-intercalator bracketing the TpA step. The octadepsipeptide ring binds in the minor groove of the DNA. Analysis of sugar coupling constants from the phase-sensitive COSY data indicates that the sugar of the thymine in the TpA binding site adopts predominantly an N-type sugar conformation, while the remaining sugars on the DNA adopt an S-type conformation, as has been observed in other Triostin A and echinomycin complexes. The drug does not bind to the octamer [d(GGAATTCC)]2 as a bis-intercalator. Only weak nonintercalative binding is observed to this DNA octamer. These results show unambiguously that [N-MeCys3,N-MeCys7]TANDEM binds sequence specifically at TpA sites in DNA. The factors underlying the sequence specificity of [N-MeCys3,N-MeCys7]TANDEM binding to DNA are discussed.     

HOMSTRAD: adding sequence information to structure-based alignments of homologous protein families

summary: We describe an extension to the Homologous Structure Alignment Database (HOMSTRAD; Mizuguchi et al., Protein Sci., 7, 2469-2471, 1998a) to include homologous sequences derived from the protein families database Pfam (Bateman et al., Nucleic Acids Res., 28, 263-266, 2000). HOMSTRAD is integrated with the server FUGUE (Shi et al., submitted, 2001) for recognition and alignment of homologues, benefitting from the combination of abundant sequence information and accurate structure-based alignments. AVAILABILITY The HOMSTRAD database is available at: http://www-cryst.bioc.cam.ac.uk/homstrad/. Query sequences can be submitted to the homology recognition/alignment server FUGUE at: http://www-cryst.bioc.cam.ac.uk/fugue/.     

Bh-DNA: variations of the poly[d(A)].poly[d(T)] structure within the framework of the fibre diffraction studies

A refinement of the recent results for poly[d(A)].poly[d(T)] (Alexeev et al., J. Biomol. Struct. Dyn. 4,989 (1987)) involving additional parameters of the base-pair structure and of the sugar-phosphate backbone expands the conformational potential of this polynucleotide of the B type to include the possibility of bifurcated hydrogen bonds of the kind recently discovered in crystalline deoxyoligonucleotide with lone d(A)n.d(T)n stretch (Nelson et al., Nature 330, 221 (1987)). Still, analysis of the available data and energy calculations do not seem to indicate that the bifurcated H-bonds are a crucial factor responsible for the anomalous structure of the d(A)n.d(T)n sequence. The unique structural properties of poly[d(A)].poly[d(T)] can hardly be explained without taking into account its interactions with the double-layer hydration spine in the minor groove. In view of the hydration mechanism stabilizing poly[d(A)].poly[d(T)] and of the polynucleotide's heteronomous prehistory (Arnott et al., Nucleic Acids Res. 11,4141 (1983)) we suggest that this B-type structure be called Bh.     

Aspartate and asparagine tRNA genes in wheat mitochondrial DNA: a cautionary note on the isolation of tRNA genes from plants.

We have identified genes encoding a "native" tRNA(Asp) (trnD-GTC) and a "chloroplast-like" tRNA(Asn) (trnN-GTT) on opposite strands and 633 bp apart within a sequenced 1640 bp RsaI restriction fragment of wheat mtDNA. The trnD gene has been found previously at a different location in wheat mtDNA (P.B.M. Joyce et al. (1988) Piant Mol. Biol. 11, 833-843); the duplicate copies of this gene are identical within the coding and immediate flanking regions (9 bp downstream and at least 68 bp upstream), after which obvious sequence similarity abruptly disappears. The trnN gene is identical to its homolog in maize ctDNA; continuation of sequence similarity beyond the coding region suggests that this gene originated as promiscuous ctDNA that is now part of the wheat mitochondrial genome. In the course of this work, we have encountered some unexpected similarities between tRNA gene regions from wheat mitochondria and other sources. Detailed analysis of these similarities leads us to suggest that trnN genes reportedly from petunia nuclear DNA (N. Bawnik et al. (1983) Nucleic Acids Res. 11, 1117-1122) and lupine mtDNA (B. Karpińska and H. Augustyniak (1988) Nucleic Acids Res. 16, 6239) are, in fact, from petunia mtDNA and lupine ctDNA, respectively, whereas a putative wheat nuclear tRNA(Ser) (trnS-TGA) gene (Z. Szwekowska-Kulińska et al. (1989) Gene 77, 163-167) is actually from wheat mtDNA. In these instances, it seems probable that the DNA samples used for cloning contained trace amounts of DNA from another sub-cellular compartment, leading to the inadvertent selection of spurious clones.     

Purification and complete primary structures of the heparin-, cell-, and DNA-binding domains of bovine plasma fibronectin

The complete amino acid sequences of the heparin-, cell- and DNA-binding domains of bovine plasma fibronectin have been determined. The fragments were generated from the 170-kDa central plasmic fragment by extensive digestion with chymotrypsin, and they contain 268, 300 and 269 amino acid residues, respectively. No half-cystines or cysteines were found in these sequences. A glucosamine-based oligosaccharide group is attached to Asn-108 in the sequence of the DNA-binding domain. Only one of the three types of internal homology found in fibronectin [Petersen et al. (1983) Proc. Natl Acad. Sci. USA 80, 137-141], namely the type III homology, occurs in these three fragments, and each of them consists of approximately three stretches of this type III homology. Part of the arrangement of peptides was derived by comparison with the partial cDNA sequence for human fibronectin recently reported [Kornblihtt et al. (1984) Nucleic Acids Res. 12, 5853-5868].     

Systematic sequencing of the Escherichia coli genome: analysis of the 2.4-4.1 min (110,917-193,643 bp) region.

The complete sequence analysis of the E. coli genome was initiated as a collaborative study in Japan. Following the initial analysis of the 0-2.4 min region (Yura, T. et al. (1992) Nucleic Acids Res. 20, 3305-3308), a contiguous sequence of 82,727 bp corresponding to the 2.4-4.1 min region (110,917-193,643 bp as counted from 0 min) was determined. The resulting sequence was found to contain at least 33 known genes and 24 putative genes predicted from protein sequence homology.     

A peptide to DNA conversion program.

A modification and extension of the computer program REVCUT (Blumenthal et al, Nucl.Acids Res. 10, 91-101 (1982) is described. The new program searches for restriction endonuclease recognition sites that are not coding DNA sequences of a protein of known aminoacid sequence using bit patterns. The modifications make the program more accurate and extend the range of the restriction endonucleases.     

Multiplex strand displacement amplification (SDA) and detection of DNA sequences from Mycobacterium tuberculosis and other mycobacteria.

Strand Displacement Amplification (SDA) is an isothermal, in vitro method of amplifying a DNA target sequence prior to detection [Walker et al (1992) Nucleic Acids Res., 20, 1691-1693]. Here we describe a multiplex form of SDA that allows two target sequences and an internal amplification control to be co-amplified by a single pair of primers after common priming sequences are spontaneously appended to the ends of target fragments. Multiplex SDA operates at a single temperature, under the same simple protocol previously developed for single-target SDA. We applied multiplex SDA to co-amplification of a target sequence (IS6110) that is specific to members of the Mycobacterium tuberculosis-complex and a target (16S ribosomal gene) that is common to most clinically relevant species of mycobacteria. Both targets are amplified 10(8)-fold during a 2 hour, single temperature incubation. The relative sensitivity of the system was evaluated across a number of clinically relevant mycobacteria and checked for crossreactivity against organisms that are closely related to mycobacteria.