Analysis of DNA sequence transformations on grids

Study of the evolution of species or organisms is essential for various biological applications. Evolution is typically studied at the molecular level by analyzing the mutations of DNA sequences of organisms. Techniques have been developed for building phylogenetic or evolutionary trees for a set of sequences. Though phylogenetic trees capture the overall evolutionary relationships among the sequences, they do not reveal fine-level details of the evolution. In this work, we attempt to resolve various fine-level sequence transformation details associated with a phylogenetic tree using cellular automata. In particular, our work tries to determine the cellular automata rules for neighbor-dependent mutations of segments of DNA sequences. We also determine the number of time steps needed for evolution of a progeny from an ancestor and the unknown segments of the intermediate sequences in the phylogenetic tree. Due to the existence of vast number of cellular automata rules, we have developed a grid system that performs parallel guided explorations of the rules on grid resources. We demonstrate our techniques by conducting experiments on a grid comprising machines in three countries and obtaining potentially useful statistics regarding evolutions in three HIV sequences. In particular, our work is able to verify the phenomenon of neighbor-dependent mutations and find that certain combinations of neighbor-dependent mutations, defined by a cellular automata rule, occur with greater than 90% probability. We also find the average number of time steps for mutations for some branches of phylogenetic tree over a large number of possible transformations with standard deviations less than 2.

Role of DNA sequence based structural features of promoters in transcription initiation and gene expression

D Regulatory information for transcription initiation is present in a stretch of genomic DNA, called the promoter region that is located upstream of the transcription start site (TSS) of the gene. The promoter region interacts with different transcription factors and RNA polymerase to initiate transcription and contains short stretches of transcription factor binding sites (TFBSs), as well as structurally unique elements. Recent experimental and computational analyses of promoter sequences show that they often have non-B-DNA structural motifs, as well as some conserved structural properties, such as stability, bendability, nucleosome positioning preference and curvature, across a class of organisms. Here, we briefly describe these structural features, the differences observed in various organisms and their possible role in regulation of gene expression.

Sex Determination - A Hypothesis Based On Noncoding Dna

Certain recent models of sex determination in mammals, Drosophila melanogaster, Caenorhabditis elegans, and snakes are examined in the light of the hypothesis that the relevant genetic regulatory mechanisms are similar and interrelated. The proposed key element in each of these instances is a noncoding DNA sequence, which serves as a high-affinity binding site for a repressor-like molecule regulating the activity of a major "sex-determining" gene. On this basis it is argued that, in several eukaryotes, (i) certain DNA sequences that are sex-determining are noncoding, in the sense that they are not the structural genes of a sex-determining protein; (ii) in some species these noncoding sequences are present in one sex and absent in the other, while in others their copy number or accessibility to regulatory molecules is significantly unequal between the two sexes; and (iii) this inequality determines whether the embryo develops into a male or a female.

Chemistry and biology of DNA methyltransferases

Recognition of a specific DNA sequence by a protein is probably the best example of macromolecular interactions leading to various events. It is a prerequisite to understanding the basis of protein-DNA interactions to obtain a better insight into fundamental processes such as transcription, replication, repair, and recombination. DNA methyltransferases with varying sequence specificities provide an excellent model system for understanding the molecular mechanism of specific DNA recognition. Sequence comparison of cloned genes, along with mutational analyses and recent crystallographic studies, have clearly defined the functions of various conserved motifs. These enzymes access their target base in an elegant manner by flipping it out of the DNA double helix. The drastic protein-induced DNA distortion, first reported for HhaI DNA methyltransferase, appears to be a common mechanism employed by various proteins that need to act on bases. A remarkable feature of the catalytic mechanism of DNA (cytosine-5) methyltransferases is the ability of these enzymes to induce deamination of the target cytosine in the absence of S-adenosyl-L-methionine or its analogs. The enzyme-catalyzed deamination reaction is postulated to be the major cause of mutational hotspots at CpG islands responsible for various human genetic disorders. Methylation of adenine residues in Escherichia coli is known to regulate various processes such as transcription, replication, repair, recombination, transposition, and phage packaging.

DNA duplex with the potential to change handedness after every half a turn

Polymorphic forms of the DNA duplex with long stretches of structural monotony are known. Several alternating purine-pyrimidine sequences have been shown to adopt left-handed Z-conformation. We report a DNA sequence d(CGCGCGATCGAT)n exhibiting alternating right-handed B and left-handed Z helical conformation after every half a turn. Further, this unusual conformation with change in handedness after every six base pairs was induced at physiological superhelical density.

DNA Compaction by a Dendrimer

At physiological pH, a PAMAM dendrimer is positively charged and can effectively bind negatively charged DNA. Currently, there has been great interest in understanding this complexation reaction both for fundamental (as a model for complex biological reactions) as well as for practical (as a gene delivery material and probe for sensing DNA sequence) reasons. Here, we have studied the complexation between double-stranded DNA (dsDNA) and various generations of PAMAM dendrimers (G3-05) through atomistic molecular dynamics simulations in the presence of water and ions. We report the compaction of DNA on a nanosecond time scale. This is remarkable, given the fact that such a short DNA duplex with a length close to 13 nm is otherwise thought to be a rigid rod. Using several nanoseconds long MD simulations, we have observed various binding modes of dsDNA and dendrimers for various generations of PAMAM dendrimers at varying charge ratios, and it confirms some of the binding modes proposed earlier. The binding is driven by the electrostatic interaction, and the larger the dendrimer charge, the stronger the binding affinity. As DNA wraps/binds to the dendrimer, counterions originally condensed onto DNA (Na+) and the dendrimer (Cl-) get released. We calculate the entropy of counterions and show that there is gain in entropy due to counterion release during the complexation. MD simulations demonstrate that, when the charge ratio is greater than 1 (as in the case of the G5 dendrimer), the optimal wrapping of DNA is observed. Calculated binding energies of the complexation follow the trend G5 > 04 > 03, in accordance with the experimental data. For a lower-generation dendrimer, such as G3, and, to some extent, for G4 also, we see considerable deformation in the dendrimer structure due to their flexible nature. We have also calculated the various helicoidal parameters of DNA to study the effect of dendrimer binding on the structure of DNA. The B form of the DNA is well preserved in the complex, as is evident from various helical parameters, justifying the use of the PAMAM dendrimer as a suitable delivery vehicle.

An ensemble of B-DNA dinucleotide geometries lead to characteristic nucleosomal DNA structure and provide plasticity required for gene expression

Background: A nucleosome is the fundamental repeating unit of the eukaryotic chromosome. It has been shown that the positioning of a majority of nucleosomes is primarily controlled by factors other than the intrinsic preference of the DNA sequence. One of the key questions in this context is the role, if any, that can be played by the variability of nucleosomal DNA structure. Results: In this study, we have addressed this question by analysing the variability at the dinucleotide and trinucleotide as well as longer length scales in a dataset of nucleosome X-ray crystal structures. We observe that the nucleosome structure displays remarkable local level structural versatility within the B-DNA family. The nucleosomal DNA also incorporates a large number of kinks. Conclusions: Based on our results, we propose that the local and global level versatility of B-DNA structure may be a significant factor modulating the formation of nucleosomes in the vicinity of high-plasticity genes, and in varying the probability of binding by regulatory proteins. Hence, these factors should be incorporated in the prediction algorithms and there may not be a unique `template' for predicting putative nucleosome sequences. In addition, the multimodal distribution of dinucleotide parameters for some steps and the presence of a large number of kinks in the nucleosomal DNA structure indicate that the linear elastic model, used by several algorithms to predict the energetic cost of nucleosome formation, may lead to incorrect results.

Nonspecific Interaction between DNA and Protein allows for Cooperativity: A Case Study with Mycobacterium DNA Binding Protein

Different DNA-binding proteins have different interaction modes with DNA. Sequence-specific DNA protein interaction has been mostly associated with regulatory processes inside a cell, and as such extensive studies have been made. Adequate data is also available on nonspecific DNA protein interaction, as an intermediate to protein's search for its cognate partner. Multidomain nonspecific DNA protein interaction involving physical sequestering of DNA has often been implicated to regulate gene expression indirectly. However, data available on this type of interaction is limited. One such interaction is the binding of DNA with mycobacterium DNA binding proteins. We have used the Langmuir-Blodgett technique to evaluate for the first time the kinetics and thermodynamics of Mycobacterium smegmatis Dps 1 binding to DNA. By immobilizing one of the interacting partners, we have shown that, when a kinetic bottleneck is applied, the binding mechanism showed cooperative binding (n = 2.72) at lower temperatures, but the degree of cooperativity gradually reduces (n = 1.38) as the temperature was increased We have also compared the kinetics and thermodynamics of sequence-specific and nonspecific DNA protein interactions under the same set of conditions.

Stabilization and Structural Alteration of the G-Quadruplex DNA Made from the Human Telomeric Repeat Mediated by Troger's Base Based Novel Benzimidazole Derivatives

Ligand-induced stabilization of the G-quadruplex DNA structure derived from the single-stranded 3'-overhang of the telomeric DNA is an attractive strategy for the inhibition of the telomerase activity. The agents that can induce/stabilize a DNA sequence into a G-quadruplex structure are therefore potential anticancer drugs. Herein we present the first report of the interactions of two novel bisbenzimidazoles (TBBz1 and TBBz2) based on Troger's base skeleton with the G-quadruplex DNA (G4DNA). These Troger's base molecules stabilize the G4DNA derived from a human telomeric sequence. Evidence of their strong interaction with the G4DNA has been obtained from CD spectroscopy, thermal denaturation, and UV-vis titration studies. These ligands also possess significantly higher affinity toward the G4DNA over the duplex DNA. The above results obtained are in excellent agreement with the biological activity, measured in vitro using a modified TRAP assay. Furthermore, the ligands are selectively more cytotoxic toward the cancerous cells than the corresponding noncancerous cells. Computational studies suggested that the adaptive scaffold might allow these ligands to occupy not only the G-quartet planes but also the grooves of the G4DNA.

Detection of G-Quadruplex DNA Using Primer Extension as a Tool

DNA sequence and structure play a key role in imparting fragility to different regions of the genome. Recent studies have shown that non-B DNA structures play a key role in causing genomic instability, apart from their physiological roles at telomeres and promoters. Structures such as G-quadruplexes, cruciforms, and triplexes have been implicated in making DNA susceptible to breakage, resulting in genomic rearrangements. Hence, techniques that aid in the easy identification of such non-B DNA motifs will prove to be very useful in determining factors responsible for genomic instability. In this study, we provide evidence for the use of primer extension as a sensitive and specific tool to detect such altered DNA structures. We have used the G-quadruplex motif, recently characterized at the BCL2 major breakpoint region as a proof of principle to demonstrate the advantages of the technique. Our results show that pause sites corresponding to the non-B DNA are specific, since they are absent when the G-quadruplex motif is mutated and their positions change in tandem with that of the primers. The efficiency of primer extension pause sites varied according to the concentration of monovalant cations tested, which support G-quadruplex formation. Overall, our results demonstrate that primer extension is a strong in vitro tool to detect non-B DNA structures such as G-quadruplex on a plasmid DNA, which can be further adapted to identify non-B DNA structures, even at the genomic level.

Is Human-X Chromosome Inactivation A Sex-Determining Device

The evolutionary function of X chromosome inactivation is thought to be dosage compensation. However, there is, at present, little evidence to suggest that most X chromosome-linked genes require such compensation. Another view--that X chromosome inactivation may be related to sex determination--is examined here. Consider a hypothetical DNA sequence regulating a major structural gene concerned with the determination of maleness. If this regulatory sequence occurs in both X and Y chromosomes and if its copy number in the Y chromosome is significantly greater than in the X chromosome, then the male-determining properties of the Y chromosome could be attributed to this higher copy number. On the other hand, if the Y chromosome has the same copy number of this sequence as the X chromosome, it is difficult to see how determination of two sexes would occur under such circumstances because XX and XY genomes would then be indistinguishable in this regard. Such a situation seems to occur in the human species with respect to the banded krait minor satellite, a repetitious DNA sequence associated with sex determination. This apparent difficulty may be resolved if X chromosome inactivation renders regulatory as well as structural genes nonfunctional and thereby brings about a significant reduction in the effective copy number of X chromosome-linked DNA sequences concerned with sex determination. It is suggested that X chromosome inactivation brings about, in this manner, a critical inequality between XX and XY embryos and that sex determination in humans is a consequence of this inequality. An analogous situation appears to exist in certain insects in which inactivation of a haploid set of chromosomes (and presumably, therefore, a 50% reduction in the effective copy number of most genes) is associated with maleness. If this line of reasoning is correct, it would suggest that sex determination may be the primary function of X chromosome inactivation.

Construction of a cDNA clone for a nuclear-coded subunit of cytochrome c oxidase from rat liver

A cDNA library for 6S–9S poly(A)-containing RNA from rat liver was constructed in Image . Initial screening of the clones was carried out using single stranded 32P-labeled cDNA prepared against poly(A)-containing RNA isolated from immunoadsorbed polyribosomes enriched for the nuclear-coded subunit messenger RNAs of cytochrome c oxidase. One of the clones, pCO89, was found to hybridize with the messenger RNA for subunit VIC. The DNA sequence of the insert in pCO89 was carried out and it has got extensive homology with the C-terminal 33 amino acids of subunit VIC from beef heart cytochrome c oxidase. In addition, the insert contained 146 bp, corresponding to a portion of the 3′-non-coding region. Northern blot analysis of rat liver RNA with the nick-translated insert of pCO89 revealed that the messenger RNA for subunit VI would contain around 510 bases.

Polypurine/polypyrimidine sequences as cis-acting transcriptional regulators

Genome sequence information has generated increasing evidence for the claim that repetitive DNA sequences present within and around genes could play a important role in the regulation of gene expression. Polypurine/polypyrimidine sequences [poly(Pu/Py)] have been observed in the vicinity of promoters and within the transcribed regions of many genes. To understand whether such sequences influence the level of gene expression, we constructed several prokaryotic and eukaryotic expression vectors incorporating poly(Pu/Py) repeats both within and upstream of a reporter gene, lacZ (encoding β-galactosidase), and studied its expression in vivo. We find that, in contrast to the situation in Escherichia coli, the presence of poly(Pu/Py) sequences within the gene does not significantly inhibit gene expression in mammalian cells. On the other hand, the presence of such sequences upstream of lacZ leads to a several-fold reduction of gene expression in mammalian cells. Similar down-regulation was observed when a structural cassette containing poly(Pu/Py) sequences upstream of lacZ was integrated into yeast chromosome V. Sequence analysis of the nine totally sequenced yeast chromosomes shows that a large number of such sequences occur upstream of ORFs. On the basis of our experimental results and DNA sequence analysis, we propose that these sequences can function as cis-acting transcriptional regulators.

Cloning and sequencing of winged bean (Psophocarpus tetragonolobus) basic agglutinin (WBA I): presence of second glycosylation site and its implications in quaternary structure

We report cloning of the DNA encoding winged bean basic agglutinin (WBA I). Using oligonucleotide primers corresponding to N- and C-termini of the mature lectin, the complete coding sequence for WBA I could be amplified from genomic DNA. DNA sequence determination by the chain termination method revealed the absence of any intervening sequences in the gene. The DNA deduced amino acid sequence of WBA I displayed some differences with its primary structure established previously by chemical means. Comparison of the sequence of WBA I with that of other legume lectins highlighted several interesting features, including the existence of the largest specificity determining loop which might account for its oligosaccharide-binding specificity and the presence of an additional N-glycosylation site. These data also throw some light on the relationship between the primary structure of the protein and its probable mode of dimerization.

Simple repetitive sequences in the genome: structure and functional significance

The current explosion of DNA sequence information has generated increasing evidence for the claim that noncoding repetitive DNA sequences present within and around different genes could play an important role in genetic control processes, although the precise role and mechanism by which these sequences function are poorly understood. Several of the simple repetitive sequences which occur in a large number of loci throughout the human and other eukaryotic genomes satisfy the sequence criteria for forming non-B DNA structures in vitro. We have summarized some of the features of three different types of simple repeats that highlight the importance of repetitive DNA in the control of gene expression and chromatin organization. (i) (TG/CA)n repeats are widespread and conserved in many loci. These sequences are associated with nucleosomes of varying linker length and may play a role in chromatin organization. These Z-potential sequences can help absorb superhelical stress during transcription and aid in recombination. (ii) Human telomeric repeat (TTAGGG)n adopts a novel quadruplex structure and exhibits unusual chromatin organization. This unusual structural motif could explain chromosome pairing and stability. (iii) Intragenic amplification of (CTG)n/(CAG)n trinucleotide repeat, which is now known to be associated with several genetic disorders, could down-regulate gene expression in vivo. The overall implications of these findings vis-à-vis repetitive sequences in the genome are summarized.