Germ-line transcripts of the immunoglobulin lambda J-C clusters in the mouse: characterization of the initiation sites and regulatory elements.

Transcription of unrearranged immunoglobulin gene segments strongly correlates with their accessibility to the V(D)J recombination machinery. The regulatory mechanisms governing this germ-line transcription are still poorly defined. In order to identify new regulatory elements, we first carried out a detailed characterization of the transcription initiation sites for the J-C germ-line transcripts, using rapid amplification of 5' cDNA ends, assisted by a template switching mechanism at the 5'-end of the RNA. Transcripts were observed that initiated heterogeneously, starting up to 293 (lambda1), 116 bp (lambda2) and 79 bp (lambda3) upstream from the respective Jlambda gene segment. Additional RT-PCR analysis revealed the existence of germ-line transcripts of lambda and also of kappa that initiate even more upstream of these transcription initiation sites, although their frequencies were low. Promoter activity was detected in vitro 5' of Jlambda2, with the minimal promoter activity mapping to the region between positions -35 and -120. In addition, computer analysis allowed the prediction of a nuclear scaffold/matrix attachment (S/MAR) region between the two J-C gene clusters at each hemi-locus. This region between the lambda1/lambda3 clusters binds to the nuclear matrix in vitro, and J-C lambda1 germ-line transcription initiates a short distance downstream from this S/MAR element.

AT-rich islands in genomic DNA as a novel target for AT-specific DNA-reactive antitumor drugs.

Interstrand cross-links at T(A/T)4A sites in cellular DNA are associated with hypercytotoxicity of an anticancer drug, bizelesin. Here we evaluated whether these lethal effects reflect targeting critical genomic regions. An in silico analysis of human sequences showed that T(A/T)4A motifs are on average scarce and scattered. However, significantly higher local motif densities were identified in distinct minisatellite regions (200-1000 base pairs of approximately 85-100% AT), herein referred to as "AT islands." Experimentally detected bizelesin lesions agree with these in silico predictions. Actual bizelesin adducts clustered within the model AT island naked DNA, whereas motif-poor sequences were only sparsely adducted. In cancer cells, bizelesin produced high levels of lesions (approximately 4.7-7.1 lesions/kilobase pair/microM drug) in several prominent AT islands, compared with markedly lower lesion levels in several motif-poor loci and in bulk cellular DNA (approximately 0.8-1.3 and approximately 0.9 lesions/kilobase pair/microM drug, respectively). The identified AT islands exhibit sequence attributes of matrix attachment regions (MARs), domains that organize DNA loops on the nuclear matrix. The computed "MAR potential" and propensity for supercoiling-induced duplex destabilization (both predictive of strong MARs) correlate with the total number of bizelesin binding sites. Hence, MAR-like AT-rich non-coding domains can be regarded as a novel class of critical targets for anticancer drugs.

Stress-induced DNA duplex destabilization in transcriptional initiation.

Stress-induced destabilization of the DNA double helix (SIDD) is involved in several mechanisms by which transcription is regulated. This paper describes a computational method for predicting the locations and extents of destabilization as functions of DNA sequence and imposed superhelical stress. This method is used to investigate several transcriptional regulatory events. These include IHF-mediated activation of gene expression in E. coli, the bimodal control of the initiation of transcription from the human c-myc gene, and the determination of the minimal requirements for transcriptional activity in yeast. Collaborations with experimental groups have established the central role of SIDD in each of these processes.

Development of ribozymes that target stathmin, a major regulator of the mitotic spindle.

Stathmin is a major cytosolic phosphoprotein that plays an important role in the control of cellular proliferation by regulating the dynamics of the microtubules that make up the mitotic spindle. Because stathmin is expressed at high levels in all human cancers, it is an attractive molecular target for anticancer interventions. We had shown previously that antisense stathmin inhibition results in marked abrogation of the transformed phenotype of leukemic cells in vitro and in vivo. Unlike the antisense approach, ribozymes can catalytically cleave several molecules of target RNA. This may provide a more efficient strategy for downregulating genes, such as stathmin, that are expressed at very high levels in cancer cells. We designed several antistathmin hammerhead ribozymes and tested their cleavage activity against short synthetic stathmin RNA substrates. In vitro cleavage studies demonstrated site-specific cleavage of stathmin RNA that was dependent on ribozyme concentration and duration of exposure to ribozyme. The most active antistathmin ribozyme was capable of cleaving >90% stathmin RNA in a catalytic manner, cleaving multiple substrate molecules per ribozyme molecule. We also demonstrated that the designed antistathmin ribozymes are capable of selectively cleaving native stathmin RNA in a mixture of total RNA isolated from leukemic cells. These antistathmin ribozymes may provide a novel and effective form of gene therapy that may be applicable to a wide variety of human cancers.

An initiation element in the yeast CUP1 promoter is recognized by RNA polymerase II in the absence of TATA box-binding protein if the DNA is negatively supercoiled.

Purified RNA polymerase II initiated transcription from the yeast CUP1 promoter fused to a C-less cassette if the DNA was negatively supercoiled. Relaxed plasmid was not transcribed. Transcription did not require addition of any other transcription factors. TATA box-binding protein (TBP) was not detectable in the polymerase preparation and the TATA box was not required. Deletion analysis of the CUP1 promoter revealed that a 25-bp element containing the initiation region was sufficient for recognition by polymerase. Two transcription start sites were mapped, one of which is identical to one of the two major start sites observed in vivo. Our observations can be accounted for by using a theoretical analysis of the probability of DNA melting within the plasmid as a function of superhelix density: the CUP1 initiation element is intrinsically unstable to superhelical stress, permitting entry of the polymerase, which then scans the DNA to locate the start site. In support of this analysis, the CUP1 promoter was sensitive to mung bean nuclease. These observations and a previous theoretical analysis of yeast genes support the idea that promoters are stress points within the DNA superhelix. The role of transcription factors might be to mark the promoter and to regulate specific melting of promoter DNA.

Sequential determination of ligands binding to discrete components in heterogeneous mixtures by iterative panning and blocking (IPAB).

Biopanning has been used extensively in conjunction with purified components, but there are also examples in which mixtures of targets have been investigated. This study introduces a methodological innovation, termed iterative panning and blocking (IPAB), to extend the range of specific interactions that can be probed in mixtures. Here this procedure is used to probe a mixture of high molecular mass components of human cord blood with phage-peptide display libraries. The initial panning recovered phage that bore the consensus motif Gly-Pro-Arg-Pro, a known fibrinogen-binding motif. These phage bound specifically to purified fibrinogen. A series of peptides containing the Gly-Pro-Arg-Pro motif efficiently blocked the binding of phage having the same motif, presumably by binding to their common target. A second round of panning was performed against the same target mixture in the presence of this blocking peptide. Phage recovered from this second panning exhibited a motif (Ser-His-Tyr) that was subsequently shown to bind specifically to complement component C1q. A second peptide containing this motif specifically blocked the interaction of the phage with C1q. A third round of panning performed in the presence of both the fibrinogen- and the C1q- blocking peptides yielded phage with a new peptide motif (Asn-Pro-Phe) that also bound specifically to C1q, apparently at a new site. The three motifs isolated through this iterative process were distinct in that each was blocked only by its corresponding peptide. This IPAB strategy can be applied to many high diversity selection procedures that target complex mixtures.

Loss of FBP function arrests cellular proliferation and extinguishes c-myc expression.

The c-myc regulatory region includes binding sites for a large set of transcription factors. The present studies demonstrate that in the absence of FBP [far upstream element (FUSE)-binding protein], which binds to the single-stranded FUSE, the remainder of the set fails to sustain endogenous c-myc expression. A dominant-negative FBP DNA-binding domain lacking effector activity or an antisense FBP RNA, expressed via replication-defective adenovirus vectors, arrested cellular proliferation and extinguished native c-myc transcription from the P1 and P2 promoters. The dominant-negative FBP initially augmented the single-stranded character of FUSE; however, once c-myc expression was abolished, melting at FUSE could no longer be supported. In contrast, with antisense FBP RNA, the single-stranded character of FUSE decreased monotonically as the transcription of endogenous c-myc declined. Because transcription is the major source of super-coiling in vivo, we propose that by binding torsionally strained DNA, FBP measures promoter activity directly. We also show that FUSE is predicted to behave as a torsion-regulated switch poised to regulate c-myc and to confer a higher order regulation on a large repertoire of factors.

Inhibition of DNA supercoiling-dependent transcriptional activation by a distant B-DNA to Z-DNA transition.

Negative DNA superhelicity can destabilize the local B-form DNA structure and can drive transitions to other conformations at susceptible sites. In a molecule containing multiple susceptible sites, superhelicity can couple these alternatives together, causing them to compete. In principle, these superhelically driven local structural transitions can be either facilitated or inhibited by proteins that bind at or near potential transition sites. If a DNA region that is susceptible to forming a superhelically induced alternate structure is stabilized in the B-form by a DNA-binding protein, its propensity for transition will be transferred to other sites within the same domain. If one of these secondary sites is in a promoter region, this transfer could facilitate open complex formation and thereby activate gene expression. We previously proposed that a supercoiling-dependent, DNA structural transmission mechanism of this type is responsible for the integration host factor-mediated activation of transcription from the ilvPG promoter of Escherichia coli (Sheridan, S. D., Benham, C. J. & Hatfield, G. W. (1998) J. Biol. Chem. 273, 21298-21308). In this report we confirm the validity of this mechanism by demonstrating the ability of a distant Z-DNA-forming site to compete with the superhelical destabilization that is required for integration host factor-mediated transcriptional activation, and thereby delay its occurrence.

Activation of gene expression by a novel DNA structural transmission mechanism that requires supercoiling-induced DNA duplex destabilization in an upstream activating sequence.

We have previously demonstrated that integration host factor (IHF)-mediated activation of transcription from the ilvPG promoter of Escherichia coli requires a supercoiled DNA template and occurs in the absence of specific interactions between IHF and RNA polymerase. In this report, we describe a novel, supercoiling-dependent, DNA structural transmission mechanism for this activation. We provide theoretical evidence for a supercoiling-induced DNA duplex destabilized (SIDD) structure in the A + T-rich, ilvPG regulatory region between base pair positions +1 and -160. We show that the region of this SIDD sequence immediately upstream of an IHF binding site centered at base pair position -92 is, in fact, destabilized by superhelical stress and that this duplex destabilization is inhibited by IHF binding. Thus, in the presence of IHF, the negative superhelical twist normally absorbed by this DNA structure in the promoter distal half of the SIDD sequence is transferred to the downstream portion of the SIDD sequence containing the ilvPG promoter site. This IHF-mediated translocation of superhelical energy facilitates duplex destabilization in the -10 region of the downstream ilvPG promoter and activates transcription by increasing the rate of open complex formation.

Analysis of the structure of a natural alternating d(TA)n sequence in yeast chromatin.

We address here the question of the in vivo structure of a natural alternating d(TA)n sequence found at the 3' region of the Saccharomyces cerevisiae FBP1 gene. This sequence consists of 13 TA pairs interrupted by a TT dinucleotide in the middle of the tract. Previous experiments with cruciform-specific nucleases S1 and Endonuclease VII demonstrated the presence in vitro of a cruciform in this region. We also showed this region to be part of a nuclease hypersensitive site flanked by nucleosomes in yeast chromatin. Here we demonstrate, by means of S1 in vivo footprinting, that in yeast plasmids also adopts in vivo a non B-DNA structure where is not a cruciform. A theoretical analysis of this region that it contains a site susceptible to superhelical stress duplex destabilization. The locations and conditions under which alternative structures form in the wild-type sequence and in deletion mutants agree with these theoretical predictions, suggesting that some kind of denaturation is the alternative structure adopted by the sequence in vivo. This suggests that negative superhelical stress sufficient for local denaturation exists in nucleosomal DNA. We also demonstrate by micrococcal nuclease digestions that the deletion of the alternating d(TA)n sequence modifies the chromatin hypersensitive site but does not affect nucleosome positioning.

Computation of DNA structural variability--a new predictor of DNA regulatory regions.

MOTIVATION: Local separation of the two strands of the DNA duplex is an essential step in important biological activities, including the initiation of transcription and replication. This conformational transition occurs in response to imposed stresses, which are stringently regulated in vivo. RESULTS: This paper describes two computational methods to analyse this phenomenon: an approximate statistical mechanical method and a Monte Carlo sampling technique. Analysis of genomic DNA sequences shows that sites of predicted duplex destabilization are closely associated with regions regulating transcription, with the most destabilized sites coinciding with 3' gene termini. Experimental results supporting this conclusion are described. The incorporation of this technique into computational searches for regulatory regions is discussed. AVAILABILITY: Because the programs implementing these calculations are complex to use and not completely developed, they have not been released yet. However, the author will analyze any DNA sequences upon request. Sequences may be submitted electronically to benham/msvax.mssm.edu.

Duplex destabilization in superhelical DNA is predicted to occur at specific transcriptional regulatory regions.

Analytic methods that accurately calculate the extent of duplex destabilization induced in each base-pair of a DNA molecule by superhelical stresses are used to analyze several genomic DNA sequences. Sites predicted to be susceptible to stress-induced duplex destabilization (SIDD) are found to be closely associated with specific transcriptional regulatory regions. Operators within the promoters of SOS-regulated Escherichia coli genes are destabilized by superhelical stresses, whereas closely related sequences present elsewhere on that genome are not. Analysis of genomic sequences from the budding yeast Saccharomyces cerevisiae finds a distinctive tripartite pattern, in which the 3' and 5' termini of genes are destabilized, but the sequence encoding the primary transcript is not. Three rDNA genes from higher eukaryotes exhibit a similar pattern. Implications of these results regarding possible mechanisms of activity of the regions involved are discussed. A strategy is presented for designing experiments in which the susceptibility to SIDD of a local region is altered without changing its local base sequence. The occurrence of the observed SIDD patterns provides a new approach to searching uncharacterized genomic sequences for transcriptionally active regions.

Energetics of coupled twist and writhe changes in closed circular pSM1 DNA.

The extent of local denaturation in closed circular pSM1 DNA depends upon the linking difference, delta Lk, and the temperature, t. We have determined the denaturation profiles, using gel electrophoresis, over the ranges -37 < or = delta Lk < or = +16 and 25 degrees C < or = t < or = 65 degrees C. We have applied statistical mechanical methods to these data to evaluate the free energies of superhelix formation, of the twisting of single strands around each other, and of the initration of local denaturation. Because the complete nucleotide sequence is needed for this analysis, the complete pSM1 DNA sequence was determined and is reported here. The values of the free energy parameters found in this work agree closely with those previously obtained from experiments with pBR322 DNA, suggesting that there is little dependence of these values on the particular DNA sequence. We find the temperature dependence of these free energies by the appropriate statistical mechanical analysis of the temperature-dependent denaturation profiles produced by supercoiling. Calculations of the transition probability profiles indicate that the course of local denaturation in pSM1 DNA involves a complex competition among several sites of comparable susceptibility. This contrasts with the melting of pBR322 DNA, in which one principal site dominates. In both molecules the sites of predicted denaturation occur at or near regulatory regions, suggesting that duplex destabilization may be associated with their biological activities.

The free energy, enthalpy and entropy of native and of partially denatured closed circular DNA.

We have used gel electrophoresis to measure the progress of local denaturation in closed circular pBR322 DNA as a function of temperature and linking deficiency, delta Lk. Local denaturation is closely coupled to supercoiling in closed DNA, requiring statistical mechanical methods for analysis. We have applied these methods to the experimental data to evaluate the free energies for three associated molecular processes. These processes are changes in the residual linking deficiency, delta Lkr, initiation of local denaturation, and twisting of denatured strands about one another. Our results confirm the quadratic dependence of the supercoiling free energy upon delta Lk, with a free energy coefficient of 740/N kcal/mol at 37 degrees C, where N is the number of base-pairs. The free energy of initiation of denaturation is 10.2(+/- 0.7) kcal/mol. The free energy of interstrand twisting of denatured regions varies with the square of the twist density, with proportionality coefficient C tau = 1.62 (+/- 0.11) kcal/rad2 at 37 degrees C. We have also calculated the entropy and enthalpy of these three processes, using the temperature dependence of the respective free energies. We find that both the entropy and the enthalpy of supercoiling are positive and vary quadratically with delta Lk. The free energy of initiation of denaturation is independent of temperature, hence arises primarily from a change in enthalpy. The entropy and enthalpy of interstrand twisting of denatured regions are both positive, and the twisting force constant decreases with temperature. These results differ considerably from expectations based solely upon considerations of chain configuration in vacuo, indicating the importance of solvent-dependent factors in determining the structure of closed circular DNA.

Sites of predicted stress-induced DNA duplex destabilization occur preferentially at regulatory loci.

This paper describes a computational method to predict the sites on a DNA molecule where imposed superhelical stresses destabilize the duplex. Several DNA sequences are analyzed in this way, including the pBR322 and ColE1 plasmids, bacteriophage f1, and the polyoma and bovine papilloma virus genomes. Superhelical destabilization in these molecules is predicted to occur at small numbers of discrete sites, most of which are within regulatory regions. The most destabilized sites include the terminator and promoter regions of specific plasmid operons, the LexA binding sites of genes under SOS control, the intergenic control region of bacteriophage f1, and the polyadenylylation sites in eukaryotic viruses. These results demonstrate the existence of close correspondences between sites of predicted superhelical duplex destabilization and specific types of regulatory regions. The use of these correspondences to supplement string-matching techniques in the search for regulatory loci is discussed.

Disulfide bonding patterns and protein topologies.

This paper examines the topological properties of protein disulfide bonding patterns. First, a description of these patterns in terms of partially directed graphs is developed. The topologically distinct disulfide bonding patterns available to a polypeptide chain containing n disulfide bonds are enumerated, and their symmetry and reducibility properties are examined. The theoretical probabilities are calculated that a randomly chosen pattern of n bonds will have any combination of symmetry and reducibility properties, given that all patterns have equal probability of being chosen. Next, the National Biomedical Research Foundation protein sequence and Brookhaven National Laboratories protein structure (PDB) databases are examined, and the occurrences of disulfide bonding patterns in them are determined. The frequencies of symmetric and/or reducible patterns are found to exceed theoretical predictions based on equiprobable pattern selection. Kauzmann's model, in which disulfide bonds form during random encounters as the chain assumes random coil conformations, finds that bonds are more likely to form with near neighbor cysteines than with remote cysteines. The observed frequencies of occurrence of disulfide patterns are found here to be virtually uncorrelated with the predictions of this alternative random bonding model. These results strongly suggest that disulfide bond pattern formation is not the result of random factors, but instead is a directed process. Finally, the PDB structure database is examined to determine the extrinsic topologies of polypeptides containing disulfide bonds. A complete survey of all structures in the database found no instances in which two loops formed by disulfide bonds within the same polypeptide chain are topologically linked. Similarly, no instances are found in which two loops present on different polypeptide chains in a structure are catenated. Further, no examples of topologically knotted loops occur. In contrast, pseudolinking has been found to be a relatively frequent event. These results show a complete avoidance of nontrivial topological entanglements that is unlikely to be the result of chance events. A hypothesis is presented to account for some of these observations.

An algebraic representation of RNA secondary structures.

This paper develops mathematical methods for describing and analyzing RNA secondary structures. It was motivated by the need to develop rigorous yet efficient methods to treat transitions from one secondary structure to another, which we propose here may occur as motions of loops within RNAs having appropriate sequences. In this approach a molecular sequence is described as a vector of the appropriate length. The concept of symmetries between nucleic acid sequences is developed, and the 48 possible different types of symmetries are described. Each secondary structure possible for a particular nucleotide sequence determines a symmetric, signed permutation matrix. The collection of all possible secondary structures is comprised of all matrices of this type whose left multiplication with the sequence vector leaves that vector unchanged. A transition between two secondary structures is given by the product of the two corresponding structure matrices. This formalism provides an efficient method for describing nucleic acid sequences that allows questions relating to secondary structures and transitions to be addressed using the powerful methods of abstract algebra. In particular, it facilitates the determination of possible secondary structures, including those containing pseudoknots. Although this paper concentrates on RNA structure, this formalism also can be applied to DNA.

Energetics of the strand separation transition in superhelical DNA.

In this paper the values of three free energy parameters governing the superhelical strand separation transition are determined by analysis of available experimental data. These are the free energy, a, needed to initiate a run of separation, the torsional stiffness, C, associated with interstrand winding of the two single strands comprising a separated site and the coefficient, K, of the quadratic free energy associated to residual linking. The experimental data used in this analysis are the locations and relative amounts of strand separation occurring in the pBR322 DNA molecule and the measured residual linking, both evaluated over a range of negative linking differences. The analytic method used treats strand separation as a heteropolymeric, co-operative, two-state transition to a torsionally deformable alternative conformation, which takes place in a circular DNA molecule constrained by the constancy of its linking number. The values determined for these parameters under the experimental conditions (T = 310 K, pH = 7.0, monovalent cation concentration = 0.01 M) are a = 10.84(+/- 0.2) kcal/mol, C = 2.5(+/- 0.3) x 10(-13) erg/rad2 and K = 2350(+/- 80) RT/N, where N is the molecular length in base-pairs. In order to assess the accuracy of the author's theoretical methods, these free energy parameters are incorporated into the analysis of superhelical strand separation in different molecules and under other conditions than those used in their evaluation. First, the temperature dependence of transition is treated, then superhelical strand separation is analyzed in a series of DNA molecules having systematic sequence modifications, and the results of these theoretical analyses are compared with those from experiments. In all molecules, transition is predicted in the range of linking differences where it is seen experimentally. Moreover, it occurs at the specific sequence locations that the analysis predicts, and with approximately the predicted relative amounts of transition at each location. The known sensitivities of this transition to changes of temperature and to small sequence modifications are predicted in a quantitatively precise manner by the theoretical results. The demonstrated high-level precision of these theoretical methods provides a tool for the screening of DNA sequences for sites susceptible to superhelical strand separation, some of which may have regulatory or other biological significance.

Negative supercoiling increases the sensitivity of plasmid DNA to single-strand break induction by X-rays.

Negatively supercoiled topoisomers of the plasmid pIBI 30 were irradiated with 250 kV X-rays and assayed for strand scission by agarose gel electrophoresis. The survival of supercoiled molecules (Form I) decreased exponentially with increasing X-ray exposure and the dose required to reduce the fraction of DNA in Form I to 37% of its value in unirradiated controls (D37) decreased with increasing negative superhelicity. This enhanced radiation sensitivity of underwound DNA is tentatively attributed to the transient denaturation of the double helix that increases the susceptibility of individual strands to free radical attack.

The influence of tertiary structural restraints on conformational transitions in superhelical DNA.

This paper examines theoretically the effects that restraints on the tertiary structure of a superhelical DNA domain exert on the energetics of linking and the onset of conformational transitions. The most important tertiary constraint arises from the nucleosomal winding of genomic DNA in vivo. Conformational transitions are shown to occur at equilibrium at less extreme superhelicities in DNA whose tertiary structure is restrained than in unrestrained molecules where the residual linking difference alpha res (that part of the superhelical deformation which is not absorbed by transitions) may be freely partitioned between twisting and bending. In the extreme case of a rigidly held tertiary structure, this analysis predicts that the B-Z transition will occur at roughly half the superhelix density needed to drive the same transition in solution, other factors remaining fixed. This suggests that superhelical transitions may occur at more moderate superhelical deformations in vivo than in solution. The influence on transition behavior of the tertiary structural restraints imposed by gel conditions also are discussed.