Fos and Jun do not bend the AP-1 recognition site.

We have used a solution-based DNA cyclization assay and a gel-phasing method to show that contrary to previous reports [Kerppola, T. K. & Curran, T. (1991) Cell 66, 317-326], basic region leucine zipper proteins Fos and Jun do not significantly bend their AP-1 recognition site. We have constructed two sets of DNA constructs that contain the 7-bp 5'-TGACTCA-3' AP-1 binding site, from either the yeast or the human collagenase gene, which is well separated from and phased by 3-4 helical turns against an A tract-directed bend. The cyclization probabilities of DNAs with altered phasings are not significantly affected by Fos-Jun binding. Similarly, Fos-Jun and Jun-Jun bound to differently phased DNA constructs show insignificant variations in gel mobilities. Both these methods independently indicate that Fos and Jun bend their AP-1 target site by <5 degrees, an observation that has important implications in understanding their mechanism of transcriptional regulation.

Simultaneous assessment of loss of heterozygosity at multiple microsatellite loci using semi-automated fluorescence-based detection: subregional mapping of chromosome 4 in cervical carcinoma.

Detection of loss of heterozygosity (LOH) by comparison of normal and tumor genotypes using PCR-based microsatellite loci provides considerable advantages over traditional Southern blotting-based approaches. However, current methodologies are limited by several factors, including the numbers of loci that can be evaluated for LOH in a single experiment, the discrimination of true alleles versus "stutter bands," and the use of radionucleotides in detecting PCR products. Here we describe methods for high throughput simultaneous assessment of LOH at multiple loci in human tumors; these methods rely on the detection of amplified microsatellite loci by fluorescence-based DNA sequencing technology. Data generated by this approach are processed by several computer software programs that enable the automated linear quantitation and calculation of allelic ratios, allowing rapid ascertainment of LOH. As a test of this approach, genotypes at a series of loci on chromosome 4 were determined for 58 carcinomas of the uterine cervix. The results underscore the efficacy, sensitivity, and remarkable reproducibility of this approach to LOH detection and provide subchromosomal localization of two regions of chromosome 4 commonly altered in cervical tumors.

DNA sequencing by delayed extraction-matrix-assisted laser desorption/ionization time of flight mass spectrometry.

Matrix-assisted laser desorption/ionization (MALDI) time of flight mass spectrometry was used to detect and order DNA fragments generated by Sanger dideoxy cycle sequencing. This was accomplished by improving the sensitivity and resolution of the MALDI method using a delayed ion extraction technique (DE-MALDI). The cycle sequencing chemistry was optimized to produce as much as 100 fmol of each specific dideoxy terminated fragment, generated from extension of a 13-base primer annealed on 40- and 50-base templates. Analysis of the resultant sequencing mixture by DE-MALDI identified the appropriate termination products. The technique provides a new non-gel-based method to sequence DNA which may ultimately have considerable speed advantages over traditional methodologies.

A unified theory of carcinogenesis based on order–disorder transitions in DNA structure as studied in the human ovary and breast

Fourier transform-infrared/statistics models demonstrate that the malignant transformation of morphologically normal human ovarian and breast tissues involves the creation of a high degree of structural modification (disorder) in DNA, before restoration of order in distant metastases. Order–disorder transitions were revealed by methods including principal components analysis of infrared spectra in which DNA samples were represented by points in two-dimensional space. Differences between the geometric sizes of clusters of points and between their locations revealed the magnitude of the order–disorder transitions. Infrared spectra provided evidence for the types of structural changes involved. Normal ovarian DNAs formed a tight cluster comparable to that of normal human blood leukocytes. The DNAs of ovarian primary carcinomas, including those that had given rise to metastases, had a high degree of disorder, whereas the DNAs of distant metastases from ovarian carcinomas were relatively ordered. However, the spectra of the metastases were more diverse than those of normal ovarian DNAs in regions assigned to base vibrations, implying increased genetic changes. DNAs of normal female breasts were substantially disordered (e.g., compared with the human blood leukocytes) as were those of the primary carcinomas, whether or not they had metastasized. The DNAs of distant breast cancer metastases were relatively ordered. These findings evoke a unified theory of carcinogenesis in which the creation of disorder in the DNA structure is an obligatory process followed by the selection of ordered, mutated DNA forms that ultimately give rise to metastases.

The solution structure of the Zα domain of the human RNA editing enzyme ADAR1 reveals a prepositioned binding surface for Z-DNA

Double-stranded RNA deaminase I (ADAR1) contains the Z-DNA binding domain Zα. Here we report the solution structure of free Zα and map the interaction surface with Z-DNA, confirming roles previously assigned to residues by mutagenesis. Comparison with the crystal structure of the (Zα)2/Z-DNA complex shows that most Z-DNA contacting residues in free Zα are prepositioned to bind Z-DNA, thus minimizing the entropic cost of binding. Comparison with homologous (α+β)helix–turn–helix/B-DNA complexes suggests that binding of Zα to B-DNA is disfavored by steric hindrance, but does not eliminate the possibility that related domains may bind to both B- and Z-DNA.

Nucleosomes, linker DNA, and linker histone form a unique structural motif that directs the higher-order folding and compaction of chromatin

The compaction level of arrays of nucleosomes may be understood in terms of the balance between the self-repulsion of DNA (principally linker DNA) and countering factors including the ionic strength and composition of the medium, the highly basic N termini of the core histones, and linker histones. However, the structural principles that come into play during the transition from a loose chain of nucleosomes to a compact 30-nm chromatin fiber have been difficult to establish, and the arrangement of nucleosomes and linker DNA in condensed chromatin fibers has never been fully resolved. Based on images of the solution conformation of native chromatin and fully defined chromatin arrays obtained by electron cryomicroscopy, we report a linker histone-dependent architectural motif beyond the level of the nucleosome core particle that takes the form of a stem-like organization of the entering and exiting linker DNA segments. DNA completes ≈1.7 turns on the histone octamer in the presence and absence of linker histone. When linker histone is present, the two linker DNA segments become juxtaposed ≈8 nm from the nucleosome center and remain apposed for 3–5 nm before diverging. We propose that this stem motif directs the arrangement of nucleosomes and linker DNA within the chromatin fiber, establishing a unique three-dimensional zigzag folding pattern that is conserved during compaction. Such an arrangement with peripherally arranged nucleosomes and internal linker DNA segments is fully consistent with observations in intact nuclei and also allows dramatic changes in compaction level to occur without a concomitant change in topology.

Attenuated T2 relaxation by mutual cancellation of dipole–dipole coupling and chemical shift anisotropy indicates an avenue to NMR structures of very large biological macromolecules in solution

Fast transverse relaxation of 1H, 15N, and 13C by dipole-dipole coupling (DD) and chemical shift anisotropy (CSA) modulated by rotational molecular motions has a dominant impact on the size limit for biomacromolecular structures that can be studied by NMR spectroscopy in solution. Transverse relaxation-optimized spectroscopy (TROSY) is an approach for suppression of transverse relaxation in multidimensional NMR experiments, which is based on constructive use of interference between DD coupling and CSA. For example, a TROSY-type two-dimensional 1H,15N-correlation experiment with a uniformly 15N-labeled protein in a DNA complex of molecular mass 17 kDa at a 1H frequency of 750 MHz showed that 15N relaxation during 15N chemical shift evolution and 1HN relaxation during signal acquisition both are significantly reduced by mutual compensation of the DD and CSA interactions. The reduction of the linewidths when compared with a conventional two-dimensional 1H,15N-correlation experiment was 60% and 40%, respectively, and the residual linewidths were 5 Hz for 15N and 15 Hz for 1HN at 4°C. Because the ratio of the DD and CSA relaxation rates is nearly independent of the molecular size, a similar percentagewise reduction of the overall transverse relaxation rates is expected for larger proteins. For a 15N-labeled protein of 150 kDa at 750 MHz and 20°C one predicts residual linewidths of 10 Hz for 15N and 45 Hz for 1HN, and for the corresponding uniformly 15N,2H-labeled protein the residual linewidths are predicted to be smaller than 5 Hz and 15 Hz, respectively. The TROSY principle should benefit a variety of multidimensional solution NMR experiments, especially with future use of yet somewhat higher polarizing magnetic fields than are presently available, and thus largely eliminate one of the key factors that limit work with larger molecules.

Knowledge-based analysis of microarray gene expression data by using support vector machines

We introduce a method of functionally classifying genes by using gene expression data from DNA microarray hybridization experiments. The method is based on the theory of support vector machines (SVMs). SVMs are considered a supervised computer learning method because they exploit prior knowledge of gene function to identify unknown genes of similar function from expression data. SVMs avoid several problems associated with unsupervised clustering methods, such as hierarchical clustering and self-organizing maps. SVMs have many mathematical features that make them attractive for gene expression analysis, including their flexibility in choosing a similarity function, sparseness of solution when dealing with large data sets, the ability to handle large feature spaces, and the ability to identify outliers. We test several SVMs that use different similarity metrics, as well as some other supervised learning methods, and find that the SVMs best identify sets of genes with a common function using expression data. Finally, we use SVMs to predict functional roles for uncharacterized yeast ORFs based on their expression data.

Solution structure and dynamics of GCN4 cognate DNA: NMR investigations

A 12 bp long GCN4-binding, self-complementary duplex DNA d(CATGACGTCATG)2 has been investigated by NMR spectroscopy to study the structure and dynamics of the molecule in aqueous solution. The NMR structure of the DNA obtained using simulated annealing and iterative relaxation matrix calculations compares quite closely with the X-ray structure of ATF/CREB DNA in complex with GCN4 protein (DNA-binding domain). The DNA is also seen to be curved in the free state and this has a significant bearing on recognition by the protein. The dynamic characteristics of the molecule have been studied by 13C relaxation measurements at natural abundance. A correlation has been observed between sequence-dependent dynamics and recognition by GCN4 protein.

Fluorescent quenching-based quantitative detection of specific DNA/RNA using a BODIPY® FL-labeled probe or primer

We have developed a simple method for the quantitative detection of specific DNA or RNA molecules based on the finding that BODIPY® FL fluorescence was quenched by its interaction with a uniquely positioned guanine. This approach makes use of an oligonucleotide probe or primer containing a BODIPY® FL-modified cytosine at its 5′-end. When such a probe was hybridized with a target DNA, its fluorescence was quenched by the guanine in the target, complementary to the modified cytosine, and the quench rate was proportional to the amount of target DNA. This widely applicable technique will be used directly with larger samples or in conjunction with the polymerase chain reaction to quantify small DNA samples.

Simultaneous A8344G heteroplasmy and mitochondrial DNA copy number quantification in Myoclonus Epilepsy and Ragged-Red Fibers (MERRF) syndrome by a multiplex Molecular Beacon based real-time fluorescence PCR

The association of a particular mitochondrial DNA (mtDNA) mutation with different clinical phenotypes is a well-known feature of mitochondrial diseases. A simple genotype–phenotype correlation has not been found between mutation load and disease expression. Tissue and intercellular mosaicism as well as mtDNA copy number are thought to be responsible for the different clinical phenotypes. As disease expression of mitochondrial tRNA mutations is mostly in postmitotic tissues, studies to elucidate disease mechanisms need to be performed on patient material. Heteroplasmy quantitation and copy number estimation using small patient biopsy samples has not been reported before, mainly due to technical restrictions. In order to resolve this problem, we have developed a robust assay that utilizes Molecular Beacons to accurately quantify heteroplasmy levels and determine mtDNA copy number in small samples carrying the A8344G tRNALys mutation. It provides the methodological basis to investigate the role of heteroplasmy and mtDNA copy number in determining the clinical phenotypes.

Non-Arrhenius kinetics for the loop closure of a DNA hairpin

Intramolecular chain diffusion is an elementary process in the conformational fluctuations of the DNA hairpin-loop. We have studied the temperature and viscosity dependence of a model DNA hairpin-loop by FRET (fluorescence resonance energy transfer) fluctuation spectroscopy (FRETfs). Apparent thermodynamic parameters were obtained by analyzing the correlation amplitude through a two-state model and are consistent with steady-state fluorescence measurements. The kinetics of closing the loop show non-Arrhenius behavior, in agreement with theoretical prediction and other experimental measurements on peptide folding. The fluctuation rates show a fractional power dependence (β = 0.83) on the solution viscosity. A much slower intrachain diffusion coefficient in comparison to that of polypeptides was derived based on the first passage time theory of SSS [Szabo, A., Schulten, K. & Schulten, Z. (1980) J. Chem. Phys. 72, 4350–4357], suggesting that intrachain interactions, especially stacking interaction in the loop, might increase the roughness of the free energy surface of the DNA hairpin-loop.

Solution structure of the 2-amino-1- methyl-6-phenylimidazo[4,5-b]pyridine C8-deoxyguanosine adduct in duplex DNA

The carcinogenic heterocyclic amine (HA) 2-amino-1-methyl-6-phenylimidazo[4,5-b]pyridine (PhIP) is formed during the cooking of various meats. To enable structure/activity studies aimed at understanding how DNA damaged by a member of the HA class of compounds can ultimately lead to cancer, we have determined the first solution structure of an 11-mer duplex containing the C8-dG adduct formed by reaction with N-acetoxy-PhIP. A slow conformational exchange is observed in which the PhIP ligand either intercalates into the DNA helix by denaturing and displacing the modified base pair (main form) or is located outside the helix in a minimally perturbed B-DNA duplex (minor form). In the main base-displaced intercalation structure, the minor groove is widened, and the major groove is compressed at the lesion site because of the location of the bulky PhIP-N-methyl and phenyl ring in the minor groove; this distortion causes significant bending of the helix. The PhIP phenyl ring interacts with the phosphodiester-sugar ring backbone of the complementary strand and its fast rotation with respect to the intercalated imidazopyridine ring causes substantial distortions at this site, such as unwinding and bulging-out of the strand. The glycosidic torsion angle of the [PhIP]dG residue is syn, and the displaced guanine base is directed toward the 3′ end of the modified strand. This study contributes, to our knowledge, the first structural information on the biologically relevant HA class to a growing body of knowledge about how conformational similarities and differences for a variety of types of lesions can influence protein interactions and ultimately biological outcome.

Multiple DNA and protein sequence alignment based on segment-to-segment comparison.

In this paper, a new way to think about, and to construct, pairwise as well as multiple alignments of DNA and protein sequences is proposed. Rather than forcing alignments to either align single residues or to introduce gaps by defining an alignment as a path running right from the source up to the sink in the associated dot-matrix diagram, we propose to consider alignments as consistent equivalence relations defined on the set of all positions occurring in all sequences under consideration. We also propose constructing alignments from whole segments exhibiting highly significant overall similarity rather than by aligning individual residues. Consequently, we present an alignment algorithm that (i) is based on segment-to-segment comparison instead of the commonly used residue-to-residue comparison and which (ii) avoids the well-known difficulties concerning the choice of appropriate gap penalties: gaps are not treated explicity, but remain as those parts of the sequences that do not belong to any of the aligned segments. Finally, we discuss the application of our algorithm to two test examples and compare it with commonly used alignment methods. As a first example, we aligned a set of 11 DNA sequences coding for functional helix-loop-helix proteins. Though the sequences show only low overall similarity, our program correctly aligned all of the 11 functional sites, which was a unique result among the methods tested. As a by-product, the reading frames of the sequences were identified. Next, we aligned a set of ribonuclease H proteins and compared our results with alignments produced by other programs as reported by McClure et al. [McClure, M. A., Vasi, T. K. & Fitch, W. M. (1994) Mol. Biol. Evol. 11, 571-592]. Our program was one of the best scoring programs. However, in contrast to other methods, our protein alignments are independent of user-defined parameters.

Differences in water release for the binding of EcoRI to specific and nonspecific DNA sequences.

The free energy difference between complexes of the restriction nuclease EcoRI with nonspecific DNA and with the enzyme's recognition sequence is linearly dependent on the water chemical potential of the solution, set using several very different solutes, ranging from glycine and glycerol to triethylene glycol and sucrose. This osmotic dependence indicates that the nonspecific complex sequesters some 110 waters more than the specific complex with the recognition sequence. The insensitivity of the difference in number of waters released to the solute identity further indicates that this water is sequestered in a space that is sterically inaccessible to solutes, most likely at the protein-DNA interface of the nonspecific complex. Calculations based on the structure of the specific complex suggest that the apposing DNA and protein surfaces in the nonspecific complex retain approximately a full hydration layer of water.