A Computer Graphics Study of Sequence-Directed Bending in DNA

Abstract

We present theoretical results to account for the unusual physical properties of a 423 bp DNA restriction fragment isolated from the kinetoplast of the trypanosomatid Leishmania tarentolae. This fragment has an anomalously low electrophoretic mobility in Polyacrylamide gels and a rotational relaxation time smaller than that of normally-behaved control fragments of the same molecular weight. Our earlier work (Proc. Natl. Acad. Sci. USA 79, 7664, 1982) has attributed these anomalies to the highly periodic distribution of the dinucleotide ApA in the DNA sequence. As originally proposed by Trifonov and Sussman (Proc. Natl. Acad. Sci. USA 77, 3816,1980) local features of the DNA structure such as a small bend at ApA, if repeated with the periodicity of the helix, will cause systematic bending of the molecule.

Computer graphics representations of DNA chain trajectories are presented for different structural models. It is shown that the structural model of Calladine (J. Mol. Biol. 161, 343, 1982) which is based on crystallographic data, is unsuccessful in predicting the systematic bending of DNA in solution.     

Maize chloroplast DNA encodes a protein sequence homologous to the bacterial ribosome assembly protein S4.

A cloned restriction fragment of maize chloroplast DNA (Bam H1 fragment 5) is shown to contain an open reading frame which encodes a basic protein of 201 amino acid residues with 40-50% sequence homology to E. coli ribosomal protein S4. Based on the experimentally determined sequence homology between the highly conserved bacterial ribosomal protein L12 and its chloroplast homologue (Bartsch M., Kimura, M. and Subramanian, A.R. (1982) Proc. Natl. Acad. Sci. USA 79, 6871), we conclude that this reading frame represents the maize chloroplast S4 gene. The nucleotide sequence of a 1100 base pair DNA segment containing the putative gene is presented.     

Bending and curvature calculations in B-DNA.

A simple program, BEND, has been written to calculate the magnitude of local bending and macroscopic curvature at each point along an arbitrary B-DNA sequence, using any desired bending model that specifies values of twist, roll and tilt as a function of sequence. The program has been used to evaluate six different DNA bending models in three categories. Two are bent non-A-tract models: (a) A new model based on the nucleosome positioning data of Satchwell et al 1986 (J. Mol. Biol. 191, 659-675), (b) The model of Calladine et al 1988 (J. Mol. Biol. 201, 127-137). Three are bent A-tract models: (c) The wedge model of Bolshoy et al 1991 (Proc. Natl. Acad. Sci. USA 88, 2312-2316), (d) The model of Cacchione et al 1989 (Biochem. 28, 8706-8713), (e) A reversed version of model (b). The last is a junction model: (f) The model of Koo & Crothers 1988 (Proc. Natl. Acad. Sci. USA 85, 1763-1767). Although they have widely different assumptions and values for twist, roll and tilt, all six models correctly predict experimental A-tract curvature as measured by gel retardation and cyclization kinetics, but only the new nucleosome positioning model is successful in predicting curvature in regions containing phased GGGCCC sequences. This model--showing local bending at mixed sequence DNA, strong bends at the sequence GGC, and straight, rigid A-tracts--is the only model consistent with both solution data from gel retardation and cyclization kinetics and structural data from x-ray crystallography.     

Decreased rate of fatty acid synthesis in brown adipose tissue of hypothalamic obese rats

Intranuclear coinjection of the late SV40 KpnI/BclI DNA fragment and the early promotor/enhancer HpaII/BglI DNA segment into permissive monkey and non-permissive mouse cells allows late SV40 gene expression without T-antigen synthesis and DNA replication. These conditions were chosen to analyse the effect of DNA methylation on V-antigen synthesis detached from the process of DNA replication. We found that in vitro methylation of a single cytosine nucleotide proximal to the major late mRNA cap site by the HpaII methylase does not block capsid protein synthesis. This result is in contrast to reported data obtained in Xenopus laevis oocyte injection experiments [(1982) Proc. Natl. Acad. Sci. USA 79, 5142-5146].     

Chromosomal locations of members of a family of novel endogenous human retroviral genomes.

Human cellular DNA contains two distinguishable families of retroviral related sequences. One family shares extensive nucleotide sequence homology with infectious mammalian type C retroviral genomes (T. I. Bonner, C. O'Connell, and M. Cohen, Proc. Natl. Acad. Sci. USA 79:4709-4713, 1982; M. A. Martin, T. Bryan, S. Rasheed, and A. S. Khan, Proc. Natl. Acad. Sci. USA 78:4892-4896, 1981). The other family contains major regions of homology with the pol genes of infectious type A and B and avian type C and D retroviral genomes (R. Callahan, W. Drohan, S. Tronick, and J. Schlom, Proc. Natl. Acad. Sci. USA 79:5503-5507, 1982; I. M. Chiu, R. Callahan, S. R. Tronick, J. Schlom, and S. A. Aaronson, Science 223:364-370, 1984). Analysis of the human recombinant clone HLM-2 has shown that the pol gene in the latter family is located within an endogenous proviral genome (R. Callahan, I. M. Chiu, J. F. H. Wong, S. R. Tronick, B. A. Roe, S. A. Aaronson, and J. Schlom, Science 228:1208-1211, 1985). We show that the proviral genome in HLM-2 and the related recombinant clone HLM-25 are located, respectively, on human chromosomes 1 and 5. Other related proviral genomes are located on chromosomes 7, 8, 11, 14, and 17.     

Nonrandom DNA sequencing of exonuclease III-deleted complementary DNA

The nonrandom DNA sequence analysis procedure of Poncz et al. [Proc. Natl. Acad. Sci. USA 79, 4298-4302 (1982)] was extensively modified to permit the determination of complementary DNA (cDNA) sequences containing G-C homopolymer regions. The recombinant cDNA plasmid was cleaved at a unique restriction enzyme site close to the cDNA and treated with Exonuclease III under controlled conditions to generate a set of overlapping fragments having deletions 50-1500 bases in length at the free 3' termini. After removal of single-stranded DNA regions by Bal31 and DNA polymerase I large fragment, the unique restriction enzyme site was recreated by blunt end ligation of synthetic oligonucleotides to the deleted DNA fragments and restriction enzyme digestion. The cDNA fragment was excised from the cloning vector using a second different restriction enzyme having a unique site that flanks the cDNA fragment and subsequently force-cloned into either M13 mp10 or mp11. This method should also be particularly useful for the sequencing of other types of DNA molecules with lengths 1500 bp or smaller.     

Sites of termination of in vitro DNA synthesis on cis-diamminedichloroplatinum(II) treated single-stranded DNA: a comparison between E. coli DNA polymerase I and eucaryotic DNA polymerases alpha.

We have compared the capacity of the large fragment of E. coli DNA polymerase I and highly purified DNA polymerases alpha from either Drosophila melanogaster embryos or calf thymus to replicate single-stranded M13 mp10 DNA treated with the antitumoral drug cis-diamminedichloroplatinum(II) (cis-DDP). We report that: a) although prokaryotic and eukaryotic enzymes have different structural complexity and dissimilar in vivo functions, their synthesis was blocked in vitro at similar sites on cis-DDP treated DNA; b) this inhibition occurred not only at d(G)n sequences, as previously reported for E. coli DNA polymerase I, (Pinto & Lippard (1985) Proc. Natl. Acad. Sci. USA, 82, 4616-4619) but also at other sequences which may represent putative cis-DDP-DNA adducts.     

Consistencies of individual DNA base-amino acid interactions in structures and sequences.

Amino acid-amino acid interaction energies have been derived from crystal structure data for a number of years. Here is reported the first derivation of normalized relative interaction from binding data for each of the four bases interacting with a specific amino acid, utilizing data from combinatorial multiplex DNA binding of zinc finger domains [Desjarlais, J. R. and Berg, J. M. (1994) Proc. Natl. Acad. Sci. USA, 91, 11099-11103]. The five strongest interactions are observed for lysine-guanine, lysine-thymine, arginine-guanine, aspartic acid-cytosine and asparagine-adenine. These rankings for interactions with the four bases appear to be related to base-amino acid partial charges. Also, similar normalized relative interaction energies are derived by using DNA binding data for Cro and lambda repressors and the R2R3 c-Myb protein domain [Takeda, Y., Sarai, A. and Rivera, V. M. (1989) Proc. Natl. Acad. Sci. USA, 86, 439-443; Sarai, A. and Takeda, Y. (1989) Proc. Natl. Acad. Sci. USA, 86, 6513-6517; Ogata, K. et al. (1995) submitted]. These energies correlate well with the combinatorial multiplex energies, and the strongest cases are similar between the two sets. They also correlate well with similar relative interaction energies derived directly from frequencies of bases in the bacteriophage lambda operator sequences. These results suggest that such potentials are general and that extensive combinatorial binding studies can be used to derive potential energies for DNA-protein interactions.     

Characterization of a mouse DNA clone containing an immunoglobulin variable region gene.

A 4.8 kilobase mouse embryo DNA fragment was inserted into a phage lambda genome and was subsequently characterized by electron microscopy, restriction enzyme mapping and partial DNA sequencing. This fragment contains a 400 base sequence which is homologous to that of an immunoglobulin light lambda chain mRNA which spans 3.3 to 3.7 kilobases from one end of the fragment. Restriction enzyme mapping as well as partial nucleotide sequencing of the 3' terminal of the homology region confirm the previous conclusion [Tonegawa, Brack, Hozumi and Schuller, Proc. Natl. Acad. Sci. USA. 74, 3518-3522 (1977)] that the cloned DNA fragment contains a Vlambda gene sequence which is separate from any Clambda sequence.     

Translocation of a lysosomal enzyme across the microsomal membrane requires signal recognition particle

Signal recognition particle (SRP), an 11S ribonucleoprotein (Walter and Blobel (1982) Nature 299, 691-698), is required for translocation of secretory proteins across microsomal membranes (Walter and Blobel (1980) Proc. Natl. Acad. Sci. USA 77, 7112-7116) and for asymmetric integration into microsomal membranes of a transmembrane protein (Anderson et al., (1982) J. Cell Biol. 93, 501-506). We demonstrate here that SRP is also required for translocation of the lysosomal protease cathepsin D across microsomal membranes.     

Computed spatial homology between the L12 protein of chloroplast ribosome and 1.7 A structure of Escherichia coli L12 domain

A computer-graphic model of the tertiary structure of a functional domain in an organelle ribosomal protein was generated using the amino acid sequence of chloroplast ribosomal protein L12 from spinach (Bartsch, Kimura and Subramanian, Proc. Natl. Acad. Sci. USA 79, 6871-6875, 1982) and 1.7 A resolution coordinates of the E. coli L12 C-terminal fragment crystal (Leijonmarck, Eriksson and Liljas, Nature 286, 824-826, 1980). A comparison between the model and the experimentally derived structure shows that although 40% of the primary structure of this part of the two proteins has undergone amino acid replacements, the gross spatial structure of the domain is maintained and the character of the surfaces of possible functional importance are not significantly altered.     

DNA sequence motifs which direct adeno-associated virus site-specific integration in a model system

The DNA sequence motifs which direct adeno-associated virus type 2 site-specific integration are being investigated using a shuttle vector, propagated as a stable episome in cultured cell lines, as the target for integration. Previously, we reported that the minimum episomal targeting elements comprise a 16-bp binding motif (Rep binding site [RBS]) for a viral regulatory protein (Rep) separated by a short DNA spacer from a sequence (terminal resolution site [TRS]) that can serve as a substrate for Rep-mediated nicking activity (R. M. Linden, P. Ward, C. Giraud, E. Winocour, and K. I. Berns, Proc. Natl. Acad. Sci. USA 93:11288-11294, 1996; R. M. Linden, E. Winocour, and K. I. Berns, Proc. Natl. Acad. Sci. USA 93:7966-7972, 1996). We now report that episomal integration depends upon both the sequence and the position of the spacer DNA separating the RBS and TRS motifs. The spacer thus constitutes a third element required for site-specific episomal integration.     

Horizontal gene transfer in microbial genome evolution

Horizontal gene transfer is the collective name for processes that permit the exchange of DNA among organisms of different species. Only recently has it been recognized as a significant contribution to inter-organismal gene exchange. Traditionally, it was thought that microorganisms evolved clonally, passing genes from mother to daughter cells with little or no exchange of DNA among diverse species. Studies of microbial genomes, however, have shown that genomes contain genes that are closely related to a number of different prokaryotes, sometimes to phylogenetically very distantly related ones. (Doolittle et al., 1990, J. Mol. Evol. 31, 383-388; Karlin et al., 1997, J. Bacteriol. 179, 3899-3913; Karlin et al., 1998, Annu. Rev. Genet. 32, 185-225; Lawrence and Ochman, 1998, Proc. Natl. Acad. Sci. USA 95, 9413-9417; Rivera et al., 1998, Proc. Natl. Acad. Sci. USA 95, 6239-6244; Campbell, 2000, Theor. Popul. Biol. 57 71-77; Doolittle, 2000, Sci. Am. 282, 90-95; Ochman and Jones, 2000, Embo. J. 19, 6637-6643; Boucher et al. 2001, Curr. Opin., Microbiol. 4, 285-289; Wang et al., 2001, Mol. Biol. Evol. 18, 792-800). Whereas prokaryotic and eukaryotic evolution was once reconstructed from a single 16S ribosomal RNA (rRNA) gene, the analysis of complete genomes is beginning to yield a different picture of microbial evolution, one that is wrought with the lateral movement of genes across vast phylogenetic distances. (Lane et al., 1988, Methods Enzymol. 167, 138-144; Lake and Rivera, 1996, Proc. Natl. Acad. Sci. USA 91, 2880-2881; Lake et al., 1999, Science 283, 2027-2028).     

DNA folding: structural and mechanical properties of the two-angle model for chromatin

We present a theoretical analysis of the structural and mechanical properties of the 30-nm chromatin fiber. Our study is based on the two-angle model introduced by Woodcock et al. (Woodcock, C. L., S. A. Grigoryev, R. A. Horowitz, and N. Whitaker. 1993. Proc. Natl. Acad. Sci. USA. 90:9021-9025) that describes the chromatin fiber geometry in terms of the entry-exit angle of the nucleosomal DNA and the rotational setting of the neighboring nucleosomes with respect to each other. We analytically explore the different structures that arise from this building principle, and demonstrate that the geometry with the highest density is close to the one found in native chromatin fibers under physiological conditions. On the basis of this model we calculate mechanical properties of the fiber under stretching. We obtain expressions for the stress-strain characteristics that show good agreement with the results of recent stretching experiments (Cui, Y., and C. Bustamante. 2000. Proc. Natl. Acad. Sci. USA. 97:127-132) and computer simulations (Katritch, V., C. Bustamante, and W. K. Olson. 2000. J. Mol. Biol. 295:29-40), and which provide simple physical insights into correlations between the structural and elastic properties of chromatin.     

Herpesvirus papio contains a plasmid origin of replication that acts in cis interspecies with an Epstein-Barr virus trans-acting function.

Herpesvirus papio (HVP) and Epstein-Barr virus (EBV) are closely related biologically and biochemically; lymphoblastoid cells infected with either virus contain episomal viral DNA. The putative origin of replication for EBV plasmids (oriP) has been assigned to a 1,790-base-pair fragment (cis) in the short unique region of the genome which requires a viral function supplied in trans from elsewhere in the genome (J. Yates, N. Warren, D. Reisman, and B. Sugden, Proc. Natl. Acad. Sci. USA 81:3806-3810, 1984). We report here the identification of the putative origin of replication (cis) in HVP; we assigned it to the HVP EcoRI K fragment. The results indicate that the HVP replication process requires both a cis and a trans-acting function, analogous to that found in EBV.     

The yeast DNA repair proteins RAD1 and RAD7 share similar putative functional domains

Sequence information on eukaryotic DNA repair proteins provided so far only few clues concerning possible functional domains. Since the DNA repair process involves a strict sequential complex formation of several proteins [1988) FASEB J. 2, 2696-2701], we searched for special protein-protein interacting domains, which consist of tandemly repeated leucine rich motifs (LRM). Search algorithms, capable of detecting even largely divergent repeats by assessing their significance due to the tandem repetitivity, revealed that the yeast DNA repair proteins RAD1 and RAD7 contain 9 and 12 tandem LRM repeats, respectively. These results represent the first clues concerning specific domains in these proteins and assign them to the LRM superfamily, which includes such members as yeast adenylate cyclase, cell surface protein receptors and ribonuclease/angiogenin inhibitor, all exerting their function by specific protein-protein interactions involving LRM domains [( 1988) EMBO J. 7, 4151-4156; (1990) Proc. Natl. Acad. Sci. USA 87, 8711-8715; (1989) Science 245, 494-499; (1990) Mol. Cell. Biol. 10, 6436-6444; (1989) Proc. Natl. Acad. Sci. USA 86, 6773-6777].     

Immunochemical detection of the alpha-subunit of the G-protein, GZ, in membranes and cytosols of mammalian cells

An antiserum (AS 98) was raised against a synthetic peptide deduced from published cDNA sequences of the alpha-subunit of the putative G-protein, GZ (Fong et al. Proc. Natl. Acad. Sci. USA 85, 3066-3070, 1988; Matsuoka et al. Proc. Natl. Acad. Sci. USA 85, 5384-5388, 1988). In membrane and cytosol preparations of many but not all tested mammalian tissues, AS 98 predominantly recognized two proteins of 40 and 43 kDa Mr. Whereas high levels of a 40 kDa GZ alpha-subunit were found in rat liver membranes and in brain cytosol, AS 98 failed to detect the alpha-subunit of GZ in brain membranes.     

Identification of the gene for DNA helicase II of Escherichia coli

Using a modification of the solid-phase radioimmune assay of Broome and Gilbert [Proc. Natl Acad. Sci. USA, 75, 2746 (1978)] to screen the plaques of lambda recombinant phages for the presence of an elevated level of helicase-II-specific antigen, we have identified the gene for helicase II in a library of Escherichia coli DNA. The DNA selected was subcloned from lambda into plasmid vectors; restriction analysis located the DNA region encoding helicase II in a PvuII fragment identical in size (2900 base pairs) and restriction pattern to that which contains the uvrD gene. Plasmids carrying this DNA fragment complemented the increased sensitivity to ultraviolet irradiation and the mutator phenotype of uvrD mutants. Furthermore, uvrD502 mutant cells were found to liberate no helicase II activity upon extraction. Following transformation with the cloned DNA, active helicase II was recovered from the mutant cells. These results support the view that helicase II is encoded by uvrD.     

The 5'-terminal sequence of TMV RNA. Question on the polymorphism found in vulgare strain

The complete nucleotide sequence of TMV RNA (common strain) reported in [Proc. Natl. Acad. Sci. USA (1982) 79, 5818] its 5'-end to be represented by two variants which differed in length. We have tested that result and sequenced the 5'-terminal regions of two strains of TMV RNA (common strain OM and tomato strain L) using cloned cDNA copies. The results showed that the 5'-terminal region of the TMV genome is not polymorphic and that one of the two variants cited above represents a tomato strain but not the common strain.