GANN: Genetic algorithm neural networks for the detection of conserved combinations of features in DNA

Background: The multitude of motif detection algorithms developed to date have largely focused on the detection of patterns in primary sequence. Since sequence-dependent DNA structure and flexibility may also play a role in protein-DNA interactions, the simultaneous exploration of sequence-and structure-based hypotheses about the composition of binding sites and the ordering of features in a regulatory region should be considered as well. The consideration of structural features requires the development of new detection tools that can deal with data types other than primary sequence. Results: GANN ( available at http://bioinformatics.org.au/gann) is a machine learning tool for the detection of conserved features in DNA. The software suite contains programs to extract different regions of genomic DNA from flat files and convert these sequences to indices that reflect sequence and structural composition or the presence of specific protein binding sites. The machine learning component allows the classification of different types of sequences based on subsamples of these indices, and can identify the best combinations of indices and machine learning architecture for sequence discrimination. Another key feature of GANN is the replicated splitting of data into training and test sets, and the implementation of negative controls. In validation experiments, GANN successfully merged important sequence and structural features to yield good predictive models for synthetic and real regulatory regions. Conclusion: GANN is a flexible tool that can search through large sets of sequence and structural feature combinations to identify those that best characterize a set of sequences.

Teaching principles of genetics through SNP analysis

An understanding of inheritance requires comprehension of genetic processes at all levels, from molecules to populations. Frequently genetics courses are separated into molecular and organismal genetics and students may fail to see the relationships between them. This is particularly true with human genetics, because of the difficulties in designing experimental approaches which are consistent with ethical restrictions, student abilities and background knowledge, and available time and materials. During 2005 we used analysis of single nucleotide polymorphisms (SNPs) in two genetic regions to enhance student learning and provide a practical experience in human genetics. Students scanned databases to discover SNPs in a gene of interest, used software to design PCR primers and a restriction enzyme based assay for the alleles, and carried out an analysis of the SNP on anonymous individual and family DNAs. The project occupied eight to ten hours per week for one semester, with some time spent in the laboratory and some spent in database searching, reading and writing the report. In completing their projects, students acquired a knowledge of Mendel’s first law (through looking at inheritance patterns), Mendel’s second law and the exceptions (the concepts of linkage and linkage disequilibrium), DNA structure (primer design and restriction enzyme analysis) and function (SNPs in coding and non-coding regions), population genetics and the statistical analysis of allele frequencies, genomics, bioinformatics and the ethical issues associated with the use of human samples. They also developed skills in presentation of results by publication and conference participation. Deficiencies in their understanding (for example of inheritance patterns, gene structure, statistical approaches and report writing) were detected and guidance given during the project. SNP analysis was found to be a powerful approach to enhance and integrate student understanding of genetic concepts.

Homomeric ring assemblies of eukaryotic Sm proteins have affinity for both RNA and DNA - Crystal structure of an oligomeric complex of yeast SmF

Sm and Sm-like proteins are key components of small ribonucleoproteins involved in many RNA and DNA processing pathways. In eukaryotes, these complexes contain seven unique Sm or Sm-like (Lsm) proteins assembled as hetero-heptameric rings, whereas in Archaea and bacteria six or seven-membered rings are made from only a single polypeptide chain. Here we show that single Sm and Lsm proteins from yeast also have the capacity to assemble into homo-oligomeric rings. Formation of homo-oligomers by the spliceosomal small nuclear ribonucleoprotein components SmE and SmF preclude hetero-interactions vital to formation of functional small nuclear RNP complexes in vivo. To better understand these unusual complexes, we have determined the crystal structure of the homomeric assembly of the spliceosomal protein SmF. Like its archaeal/bacterial homologs, the SmF complex forms a homomeric ring but in an entirely novel arrangement whereby two heptameric rings form a co-axially stacked dimer via interactions mediated by the variable loops of the individual SmF protein chains. Furthermore, we demonstrate that the homomeric assemblies of yeast Sm and Lsm proteins are capable of binding not only to oligo(U) RNA but, in the case of SmF, also to oligo(dT) single-stranded DNA.

Flexibility revealed by the 1.85 angstrom crystal structure of the beta sliding-clamp subunit of Escherichia coli DNA polymerase III

The beta subunit of the Escherichia coli replicative DNA polymerase III holoenzyme is the sliding clamp that interacts with the alpha (polymerase) subunit to maintain the high processivity of the enzyme. The beta protein is a ring-shaped dimer of 40.6 kDa subunits whose structure has previously been determined at a resolution of 2.5 Angstrom [Kong et al. (1992), Cell, 69, 425-437]. Here, the construction of a new plasmid that directs overproduction of beta to very high levels and a simple procedure for large-scale purification of the protein are described. Crystals grown under slightly modified conditions diffracted to beyond 1.9 Angstrom at 100 K at a synchrotron source. The structure of the beta dimer solved at 1.85 Angstrom resolution shows some differences from that reported previously. In particular, it was possible at this resolution to identify residues that differed in position between the two subunits in the unit cell; side chains of these and some other residues were found to occupy alternate conformations. This suggests that these residues are likely to be relatively mobile in solution. Some implications of this flexibility for the function of beta are discussed.

Fast Structure-Based Assignment of 15N HSQC Spectra of Selectively 15N-Labeled Paramagnetic Proteins

A novel strategy for fast NMR resonance assignment of N-15 HSQC spectra of proteins is presented. It requires the structure coordinates of the protein, a paramagnetic center, and one or more residue-selectively N-15-labeled samples. Comparison of sensitive undecoupled N-15 HSQC spectra recorded of paramagnetic and diamagnetic samples yields data for every cross-peak on pseudocontact shift, paramagnetic relaxation enhancement, cross-correlation between Curie-spin and dipole-dipole relaxation, and residual dipolar coupling. Comparison of these four different paramagnetic quantities with predictions from the three-dimensional structure simultaneously yields the resonance assignment and the anisotropy of the susceptibility tensor of the paramagnetic center. The method is demonstrated with the 30 kDa complex between the N-terminal domain of the epsilon subunit and the theta subunit of Escherichia Coll DNA polymerase III. The program PLATYPUS was developed to perform the assignment, provide a measure of reliability of the assignment, and determine the susceptibility tensor anisotropy.

FRA10B structure reveals common elements in repeat expansion and chromosomal fragile site genesis

A common mechanism for chromosomal fragile site genesis is not yet apparent. Folate-sensitive fragile sites are expanded p(CCG)n repeats that arise from longer normal alleles. Distamycin A or bromodeoxyuridine-inducible fragile site FRA16B is an expanded AT-rich similar to 33 bp repeat; however, the relationship between normal and fragile site alleles is not known. Here, we report that bromodeoxyuridine-inducible, distamycin A-insensitive fragile site FRA10B is composed of expanded similar to 42 bp repeats. Differences in repeat motif length or composition between different FRA10B families indicate multiple independent expansion events. Some FRA10B alleles comprise a mixture of different expanded repeat motifs. FRA10B fragile site and long normal alleles share flanking polymorphisms. Somatic and intergenerational FRA10B repeat instability analogous to that found in expanded trinucleotide repeats supports dynamic mutation as a common mechanism for repeat expansion.

Dehydromonocrotaline generates sequence-selective N-7 guanine alkylation and heat and alkali stable multiple fragment DNA crosslinks

Monocrotaline is a pyrrolizidine alkaloid known to cause toxicity in humans and animals. Its mechanism of biological action is still unclear although DNA crosslinking has been suggested to a play a role in its activity. In this study we found that an active metabolite of monocrotaline, dehydromonocrotaline (DHM), alkylates guanines at the N7 position of DNA with a preference for 5'-GG and 5'-GA sequences; In addition, it generates piperidine- and heat-resistant multiple DNA crosslinks, as confirmed by electrophoresis and electron microscopy. On the basis of these findings, we propose that DHM undergoes rapid polymerization to a structure which is able to crosslink several fragments of DNA.

The solution structure of the N-terminal zinc finger of GATA-1 reveals a specific binding face for the transcriptional co-factor FOG

Zinc fingers (ZnFs) are generally regarded as DNA-binding motifs. However, a number of recent reports have implicated particular ZnFs in the mediation of protein-protein interactions. The N-terminal ZnF of GATA-1 (NF) is one such finger, having been shown to interact with a number of other proteins, including the recently discovered transcriptional co-factor FOG. Here we solve the three-dimensional structure of the NF in solution using multidimensional H-1/N-15 NMR spectroscopy, and we use H-1/N-15 spin relation measurements to investigate its backbone dynamics. The structure consists of two distorted beta-hairpins and a single alpha-helix, and is similar to that of the C-terminal ZnF of chicken GATA-1. Comparisons of the NF structure with those of other C-4-type zinc binding motifs, including hormone receptor and LIM domains, also reveal substantial structural homology. Finally, we use the structure to map the spatial locations of NF residues shown by mutagenesis to be essential for FOG binding, and demonstrate that these residues all lie on a single face of the NE Notably, this face is well removed from the putative DNA-binding face of the NE an observation which is suggestive of simultaneous roles for the NF; that is, stabilisation of GATA-1 DNA complexes and recruitment of FOG to GATA-1-controlled promoter regions.

Crystal and molecular structure of 2-hydroxy-1-naphthaldehyde isonicotinoyl hydrazone (NIH) and its iron(III) complex: an iron chelator with anti-tumour activity

Previous studies have demonstrated that 2-hydroxy-1-naphthaldehyde isonicotinoyl hydrazone (NIH) and several other aroylhydrazone chelators possess anti-neoplastic activity due to their ability to bind intracellular iron. In this study we have examined the structure and properties of NIH and its Fe-III complex in order to obtain further insight into its anti-tumour activity. Two tridentate NIH ligands deprotonate upon coordination to Fe-III in a meridional fashion to form a distorted octahedral, high-spin complex. Solution electrochemistry of [Fe(NIH-H)(2)](+) shows that the trivalent oxidation state is dominant over a wide potential range and that the Fe-II analogue is not a stable form of this complex. The fact that [Fe(NIH-H)(2)](+) cannot-cycle between the Fe-II and Fe-III states suggests that the production of toxic free- radical species, e.g. OH. or O2(.-),is not part of this ligand's cytotoxic action. This suggestion is supported by cell culture experiments demonstrating that the addition of Fe-III to NIH prevents its anti-proliferative effect. The chemistry of this chelator and its Fe-III complex are discussed in the context of understanding its anti-tumour activity.

A preliminary phylogeny of the scale insects (Hemiptera : Sternorrhyncha : Coccoidea) based on nuclear small-subunit ribosomal DNA

Scale insects (Hemiptera: Sternorrhyncha: Coccoidea) are a speciose and morphologically specialized group of plant-feeding bugs in which evolutionary relationships and thus higher classification are controversial. Sequences derived from nuclear small-subunit ribosomal DNA were used to generate a preliminary molecular phylogeny for the Coccoidea based on 39 species representing 14 putative families. Monophyly of the archaeococcoids (comprising Ortheziidae, Margarodidae sensu lato, and Phenacoleachia) was equivocal, whereas monophyly of the neococcoids was supported. Putoidae, represented by Puto yuccae, was found to be outside the remainder of the neococcoid clade. These data are consistent with a single origin (in the ancestor of the neococcoid clade) of a chromosome system involving paternal genome elimination in males. Pseudococcidae (mealybugs) appear to be sister to the rest of the neococcoids and there are indications that Coccidae (soft scales) and Kerriidae (lac scales) are sister taxa. The Eriococcidae (felt scales) was not recovered as a monophyletic group and the eriococcid genus Eriococcus sensu lato was polyphyletic. (C) 2002 Elsevier Science (USA). All rights reserved.

A tale of two terminators: Crystal structures sharpen the debate on DNA replication fork arrest mechanisms

The structure of the Tus-Ter DNA replication fork arrest complex of Escherichia coli reveals a novel architecture for the bound Tus protein and a new type of DNA-binding motif, The structure of the complex may explain how Tus can block movement of a replication fork approaching from one direction and not the other.

Reorganization of terminator DNA upon binding replication terminator protein: Implications for the functional replication fork arrest complex

Termination of DNA replication in Bacillus subtilis involves the polar arrest of replication forks by a specific complex formed between the replication terminator protein (RTP) and DNA terminator sites. While determination of the crystal structure of RTP has facilitated our understanding of how a single RTP dimer interacts with terminator DNA, additional information is required in order to understand the assembly of a functional fork arrest complex, which requires an interaction between two RTP dimers and the terminator site. In this study, we show that the conformation of the major B. subtilis DNA terminator, Terl, becomes considerably distorted upon binding RTP. Binding of the first dimer of RTP to the B site of Terl causes the DNA to become slightly unwound and bent by similar to 40 degrees. Binding of a second dimer of RTP to the A site causes the bend angle to increase to similar to 60 degrees. We have used this new data to construct two plausible models that might explain how the ternary terminator complex can block DNA replication in a polar manner, in the first model, polarity of action is a consequence of the two RTP-DNA half-sites having different conformations. These different conformations result from different RTP-DNA contacts at each half-site (due to the intrinsic asymmetry at the terminator DNA), as well as interactions (direct or indirect) between the RTP dimers on the DNA. In the second model, polar fork arrest activity is a consequence of the different affinities of RTP for the A and B sites of the terminator DNA, modulated significantly by direct or indirect interactions between the RTP dimers.

Genetic population structure of the catadromous perciform: Macquaria novemaculeata (Percichthyidae)

Genetic population structure in the catadromous Australian bass Macquaria novemaculeata was investigated using samples from four locations spanning 600 km along the eastern Australian coastline. Both allozymes and mtDNA control region sequences were examined. Population subdivision estimates based on allozymes revealed low levels of population structuring (G(st)=0.043, P<0.05). However, mtDNA indicated moderate levels of geographic population structure (G(st)=0.146, P<0.01). Phylogenetic analysis of mtDNA control region sequences (mean sequence divergence 1.9%) indicated little phylogeographic structuring. Results suggested that genotypic variation within each river population, while bring affected primarily by genetic drift, was also prevented from more significant divergence by homogenizing levels of gene flow-synonymous with a one-dimensional stepping-stone model of population structure. The catadromous life history of Macquaria novemaculeata was considered to br influential on the pattern of population structure displayed. Results were compared to the few population genetic studies involving catadromous fishes, indicating that catadromy alone is unlikely to be a good predictor of population structure. A more comprehensive suite of biological characteristics than simple life-history traits must be considered fully to allow reliable predictive models of population structure to be formulated. (C) 1997 The Fisheries Society of the British Isles.

Complete DNA sequence and gene organization of the mitochondrial genome of the liverfluke, Fasciola hepatica L. (Platyhelminthes; Trematoda)

The complete nucleotide sequence of the mitochondrial (mt) DNA molecule of the liverfluke, Fasciola hepatica (phylum Platyhelminthes, class Trematoda, family Fasciolidae), was determined, It comprises 14462 bp, contains 12 protein-encoding, 2 ribosomal and 22 transfer RNA genes, and is the second complete flatworm (and the first trematode) mitochondrial sequence to be described in detail. All of the genes are transcribed from the same strand. Of the genes typically found in mitochondrial genomes of eumetazoans, only atp8 is absent. The nad4L and nad4 genes overlap by 40 nt. Most intergenic sequences are very short. Two larger non-coding regions are present. The longer one (817 nt) is located between trnG and cox3 and consists of 8 identical tandem repeats of 85 nt, rich in G and C, followed by 1 imperfect repeat. The shorter non-coding region (187 nt) exhibits no special features and is separated from the longer region by trnG. The gene arrangement resembles that of some other trematodes including the eastern Asian Schistosoma species (and cyclophyllidean cestode species) but it is strikingly different from that of the African schistosomes, represented by Schistosoma mansoni. The genetic code is as inferred previously for flatworms. Transfer RNA genes range in length from 58 to 70 nt, their products producing characteristic 'clover leaf' structures, except for tRNA(S-VON) and tRNA(S-AGN) lacking the DHU arm.