Molecular dynamics studies of trinucleotide repeat DNA involved in neurodegenerative disorders.

Expansion of trinucleotide repeat DNA of the classes CAG�·CTG, CGG�·CCG and GAA�·TTC are found to be associated with several neurodegenerative disorders. Different mechanisms have been attributed to the expansion of triplets, mainly involving the formation of alternate secondary structures by such repeats. This paper reports the molecular dynamics simulation of triplet repeat DNA sequences to study the basic structural features of DNA that are responsible for the formation of structures such as hairpins and slip-strand DNA leading to expansion. All the triplet repeat sequences studied were found to be more flexible compared to the control sequence unassociated with disease. Moreover, flexibility was found to be in the order CAG�·CTG > CGG�·CCG = GAA�·TTC, the highly flexible CAG�·CTG repeat being the most common cause of neurodegenerative disorders. In another simulation, a single G�·C to T�·A mutation at the 9th position of the CAG�·CTG repeat exhibited a reduction in bending compared to the pure 15-mer CAGâ�¢CTG repeat. EPM1 dodecamer repeat associated with the pathogenesis of progressive myoclonus epilepsy was also simulated and showed flexible nature suggesting a similar expansion mechanism.

The unstable (CTG)$\sb{\rm n}$ trinucleotide repeat associated with myotonic dystrophy: What it does and where it came from

Myotonic dystrophy (DM), an autosomal dominant disorder mapping to human chromosome 19q13.3, is the most common neuromuscular disease in human adults.^ Following the identification of the mutation underlying the DM phenotype, an unstable (CTG)$\sb{n}$ trinucleotide repeat in the 3$\prime$ untranslated region (UTR) of a gene encoding a ser/thr protein kinase named DM protein kinase (DMPK), the study was targeted at two questions: (1) the identification of the disease-causing mechanism(s) of the unstable repeat, and at a more basic level, (2) the identification of the origin and the mechanism(s) involved in repeat instability. The first goal was to identify the pathophysiological mechanisms of the (CTG)$\sb{n}$ repeat.^ The normal repeat is transcribed but not translated; therefore, initial studies centered on the effect on RNA transcript levels. The vast majority of DM affecteds are heterozygous for the mutant expansion, so that the normal allele interferes with the analysis of the mutant allele. A quantitative allele-specific RT-PCR procedure was developed and applied to a spectrum of patient tissue samples and cell lines. Equal levels of unprocessed pre-mRNA were determined for the wild type (+) and disease (DM) alleles in skeletal muscle and cell lines of heterozygous DM patients, indicating that any nucleosome binding has no effect at the level of transcriptional initiation and transcription of the mutant DMPK locus. In contrast, processed mRNA levels from the DM allele were reduced relative to the + allele as the size of the expansion increased. The unstable repeat, therefore, impairs post-transcriptional processing of DM allele transcripts. This phenomenon has profound effects on overall DMPK locus steady-state transcript levels in cells missing a wild type allele and does not appear to be mediated by imprinting, decreased mRNA stability, generation of aberrant splice forms, or absence of polyadenylation of the mutant allele.^ In Caucasian DM subjects, the unstable repeat is in complete linkage disequlibrium with a single haplotype composed of nine alleles within and flanking DMPK over a physical distance of 30 kb. A detailed haplotype analysis of the DM region was conducted on a Nigerian (Yoruba) DM family, the only indigenous sub-Saharan DM case reported to date. Each affected member of this family had an expanded (CTG)$\sb{n}$ repeat in one of their DMPK alleles. However, unlike all other DM populations studied thus far, disassociation of the (CTG)$\sb{n}$ repeat expansion from other alleles of the putative predisposing haplotype was found. Thus, the expanded (CTG)$\sb{n}$ repeat in this family was the result of an independent mutational event. Consequently, the origin of DM is unlikely the result of a single mutational event, and the hypothesis that a single ancestral haplotype predisposes to repeat expansion is not compelling. (Abstract shortened by UMI.) ^

The in vivo and in vitro formation of sticky DNA and the GAA.TTC trinucleotide repeat associated with Friedreich's ataxia

Friedreich's ataxia is caused by the expansion of the GAA•TTC trinucleotide repeat sequence located in intron 1 of the frataxin gene. The long GAA•TTC repeats are known to form several non-B DNA structures including hairpins, triplexes, parallel DNA and sticky DNA. Therefore it is believed that alternative DNA structures play a role in the loss of mRNA transcript and functional frataxin protein in FRDA patients. We wanted to further elucidate the characteristics for formation and stability of sticky DNA by evaluating the structure in a plasmid based system in vitro and in vivo in Escherichia coli. The negative supercoil density of plasmids harboring different lengths of GAA•TTC repeats, as well as either one or two repeat tracts were studied in E. coli to determine if plasmids containing two long tracts (≥60 repeats) in a direct repeat orientation would have a different topological effect in vivo compared to plasmids that harbored only one GAA•TTC tract or two tracts of < 60 repeats. The experiments revealed that, in fact, sticky DNA forming plasmids had a lower average negative supercoil density (-σ) compared to all other control plasmids used that had the potential to form other non-B DNA structures such as triplexes or Z-DNA. Also, the requirements for in vitro dissociation and reconstitution of the DNA•DNA associated region of sticky DNA were evaluated. Results conclude that the two repeat tracts associate in the presence of negative supercoiling and MgCl 2 or MnCl2 in a time and concentration-dependent manner. Interaction of the repeat sequences was not observed in the absence of negative supercoiling and/or MgCl2 or in the presence of other monovalent or divalent cations, indicating that supercoiling and quite specific cations are needed for the association of sticky DNA. These are the first experiments studying a more specific role of supercoiling and cation influence on this DNA conformation. To support our model of the topological effects of sticky DNA in plasmids, changes in sticky DNA band migration was measured with reference to the linear DNA after treatment with increasing concentrations of ethidium bromide (EtBr). The presence of independent negative supercoil domains was confirmed by this method and found to be segregated by the DNA-DNA associated region. Sequence-specific polyamide molecules were used to test the effect of binding of the ligands to the GAA•TTC repeats on the inhibition of sticky DNA. The destabilization of the sticky DNA conformation in vitro through this binding of the polyamides demonstrated the first conceptual therapeutic approach for the treatment of FRDA at the DNA molecular level. ^ Thus, examining the properties of sticky DNA formed by these long repeat tracts is important in the elucidation of the possible role of sticky DNA in Friedreich's ataxia. ^

Orientation-dependent and sequence-specific expansions of CTG/CAG trinucleotide repeats in Saccharomyces cerevisiae

A quantitative and selective genetic assay was developed to monitor expansions of trinucleotide repeats (TNRs) in yeast. A promoter containing 25 repeats allows expression of a URA3 reporter gene and yields sensitivity to the drug 5-fluoroorotic acid. Expansion of the TNR to 30 or more repeats turns off URA3 and provides drug resistance. When integrated at either of two chromosomal loci, expansion rates were 1 × 10−5 to 4 × 10−5 per generation if CTG repeats were replicated on the lagging daughter strand. PCR analysis indicated that 5–28 additional repeats were present in 95% of the expanded alleles. No significant changes in CTG expansion rates occurred in strains deficient in the mismatch repair gene MSH2 or the recombination gene RAD52. The frequent nature of CTG expansions suggests that the threshold number for this repeat is below 25 in this system. In contrast, expansions of the complementary repeat CAG occurred at 500- to 1,000-fold lower rates, similar to a randomized (C,A,G) control sequence. When the reporter plasmid was inverted within the chromosome, switching the leading and lagging strands of replication, frequent expansions were observed only when CTG repeats resided on the lagging daughter strand. Among the rare CAG expansions, the largest gain in tract size was 38 repeats. The control repeats CTA and TAG showed no detectable rate of expansions. The orientation-dependence and sequence-specificity data support the model that expansions of CTG and CAG tracts result from aberrant DNA replication via hairpin-containing Okazaki fragments.

Expansion of a CUG trinucleotide repeat in the 3′ untranslated region of myotonic dystrophy protein kinase transcripts results in nuclear retention of transcripts

Expansion of a CTG trinucleotide repeat in the 3′ untranslated region (UTR) of DMPK, the gene encoding myotonic dystrophy protein kinase, induces the dominantly inherited neuromuscular disorder myotonic dystrophy (DM). Transcripts containing the expanded trinucleotide are abundant in differentiated cultured myoblasts, and they are spliced and polyadenylylated normally. However, mutant transcripts never reach the cytoplasm in these nonmitotic cells; instead, they form stable clusters that are tightly linked to the nuclear matrix, which can prevent effective biochemical purification of these transcripts. In DM patients, reduced DMPK protein levels, consequent to nuclear retention of mutant transcripts, are probably a cause of disease development. Formation of nuclear foci is a novel mechanism for preventing transcript export and effecting a loss of gene function.

Evolution of the Friedreich’s ataxia trinucleotide repeat expansion: Founder effect and premutations

Friedreich’s ataxia, the most frequent inherited ataxia, is caused, in the vast majority of cases, by large GAA repeat expansions in the first intron of the frataxin gene. The normal sequence corresponds to a moderately polymorphic trinucleotide repeat with bimodal size distribution. Small normal alleles have approximately eight to nine repeats whereas a more heterogeneous mode of large normal alleles ranges from 16 to 34 GAA. The latter class accounts for ≈17% of normal alleles. To identify the origin of the expansion mutation, we analyzed linkage disequilibrium between expansion mutations or normal alleles and a haplotype of five polymorphic markers within or close to the frataxin gene; 51% of the expansions were associated with a single haplotype, and the other expansions were associated with haplotypes that could be related to the major one by mutation at a polymorphic marker or by ancient recombination. Of interest, the major haplotype associated with expansion is also the major haplotype associated with the larger alleles in the normal size range and was almost never found associated with the smaller normal alleles. The results indicate that most if not all large normal alleles derive from a single founder chromosome and that they represent a reservoir for larger expansion events, possibly through “premutation” intermediates. Indeed, we found two such alleles (42 and 60 GAA) that underwent cataclysmic expansion to pathological range in a single generation. This stepwise evolution to large trinucleotide expansions already was suggested for myotonic dystrophy and fragile X syndrome and may relate to a common mutational mechanism, despite sequence motif differences.

Trinucleotide repeats and long homopeptides in genes and proteins associated with nervous system disease and development.

Several human neurological disorders are associated with proteins containing abnormally long runs of glutamine residues. Strikingly, most of these proteins contain two or more additional long runs of amino acids other than glutamine. We screened the current human, mouse, Drosophila, yeast, and Escherichia coli protein sequence data bases and identified all proteins containing multiple long homopeptides. This search found multiple long homopeptides in about 12% of Drosophila proteins but in only about 1.7% of human, mouse, and yeast proteins and none among E. coli proteins. Most of these sequences show other unusual sequence features, including multiple charge clusters and excessive counts of homopeptides of length > or = two amino acid residues. Intriguingly, a large majority of the identified Drosophila proteins are essential developmental proteins and, in particular, most play a role in central nervous system development. Almost half of the human and mouse proteins identified are homeotic homologs. The role of long homopeptides in fine-tuning protein conformation for multiple functional activities is discussed. The relative contributions of strand slippage and of dynamic mutation are also addressed. Several new experiments are proposed.

Oxidative DNA Damage Modulates Trinucleotide Repeat Instability Via DNA Base Excision Repair

Trinucleotide repeat (TNR) expansion is the cause of more than 40 types of human neurodegenerative diseases such as Huntington’s disease. Recent studies have linked TNR expansion with oxidative DNA damage and base excision repair (BER). In this research, we provided the first evidence that oxidative DNA damage can induce CAG repeat deletion/contraction via BER. We found that BER of an oxidized DNA base lesion, 8-oxoguanine in a CAG repeat tract, resulted in the formation of a CTG hairpin at the template strand. DNA polymerase β (pol b) then skipped over the hairpin creating a 5’-flap that was cleaved by flap endonuclease 1 (FEN1) leading to CAG repeat deletion. To further investigate whether BER may help to shorten an expanded TNR tract, we examined BER in a CAG repeat hairpin loop. We found that 8-oxoguanine DNA glycosylase removed the oxidized base located in the loop region of the hairpin leaving an abasic site. Apurinic/apyrimidinic (AP) endonuclease 1 then incised the 5’-end of the abasic site leaving a nick in the loop. This further converted the hairpin into an intermediate with a 3’-flap and a 5’-flap. As a 5’-3’ endonuclease, FEN1 cleaved the 5’-flap, whereas a 3’-5’ endonuclease, Mus81/Eme1, removed the 3’-flap. The coordination between FEN1 and Mus81/Eme1 ultimately resulted in removal of a CAG repeat hairpin attenuating or preventing TNR expansion. To further explore if pol β bypass of an oxidized base lesion, 5’,8-cyclodeoxyadenosine, may affect TNR instability, we examined pol β DNA synthesis in bypassing this base lesion and found that the lesion preferentially induced TNR deletion during BER and Okazaki fragment maturation. The repeat deletion was mediated by the formation of a loop in the template strand induced specifically by the damage. Pol β then skipped over the loop structure creating a 5’-flap that was efficiently removed by FEN1 leading to repeat deletion. Our study demonstrates that pol β-mediated BER plays an important role in mediating TNR deletion and removing a TNR hairpin to prevent TNR expansion. Our research provides a molecular basis for further developing BER as a target for prevention and treatment of neurodegenerative diseases caused by TNR expansion.

Immunological studies on nucleic acids and their components.8.Antibodies specific to a deoxyribotrinucleotide sequence

Bovine serum albumin conjugates of two trinucleotides, dpTpTpA and dTpTpAp, were prepared by linking the trinucleotides through their end phosphates by the ‘carbodiimide method’. Antibodies were raised in rabbits by injecting the trinucleotide-bovine serum albumin conjugates. Analysis by double diffusion in agar gel, quantitative precipitin reaction and its inhibition by haptens showed clearly the presence of antibodies specific to the whole trinucleotide molecule. The titre of antibodies obtained by the trinucleotide-rabbit serum albumin conjugates with their respective antisera was approximately the same, indicating that linking the trinucleotide through either 5′ or 3′ phosphate does not have an appreciable effect on the titre of antibodies. The results also demonstrate that the nucleotide(s) away from the carrier protein is more immunodominant than the one linked directly to the protein.

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.

On self-interacting oligoribonucleotides. I. Absorption and optical rotatory dispersion of 2′-5′- and 3′-5′-oligoguanylic acids

A series of 2′-5′-oligoguanylic acids are prepared by reacting G(cyclic)p with takadiastase T1 ribonuclease and separating the products chromatographically. The 3′-5′-oligoguanylic acids are obtained by separating the products of alkaline degradation of 3′-5′-poly(G). The optical rotatory dispersion and hypochromism of both 2′-5′- and 3′-5′-oligoguanylic acids are studied at two different pH. The optical rotatory dispersion spectrum of 2′-5′-GpG is significantly different from that of 3′-5′-GpG. The magnitude of rotation of the long-wavelength peak of 2′-5′-GpG is larger than that of 3′-5′-GpG. This finding contradicts the explanation that the extra stability and more intense circular dichroism band of other 3′-5′-dinucleoside monophosphates is due to H-bond formation between 2′-OH and either the base or the phosphate oxygen. The end phosphate group has a marked effect on the spectrum of GpG between 230 and 250 mμ. In addition the optical rotatory dispersion spectra of 2′-5′ exhibit strong pH, temperature, and solvent dependence between 230 and 250 mμ. ΔH and AS for order ⇌ disorder transition is estimated to be 9.7 kcal/mole and 35.2 eu, respectively. The optical rotatory dispersion spectra of guanine-rich oligoribonucleotides, GpGpC, GpGpU, GpGpGpC, and GpGpGpU are compared to the calculated optical rotatory dispersion from the semiempirical expression of Cantor and Tinoco, using measured optical rotatory dispersion of dimers. Contrary to previous studies, agreement is found not at all satisfactory. However, optical rotatory dispersion of 3′-5′-GpGpGpC and GpGpGpU can be estimated from the semiempirical expression, if a next-nearest interaction parameter is introduced empirically. Such interaction parameter can be calculated from the measured properties of trinucleotide sequences like GpGpG, GpGpC, and GpGpU, assuming that only the nearest-neighbor interaction is important. The optical rotatory dispersion of single-stranded poly(G) is also predicted. The importance of syn-anti equilibrium and next-nearest-neighbor interaction in oligoguanylic acids is suggested as a probable explanation.

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.

Structural study of CUG-repeating small RNAs complexed with silencing suppressor P19

The study is focused on structural aspects of interaction between silencing suppressor p19 and CUG-repeating small RNAs. The work involves crystal structure determination of a protein-unbound RNA form and RNA fragments of various lengths (19, 20, 21 nucleotides) complexed with p19-suppressor. Results prove the ability of silencing suppressor p19 to bind CUG-repeating small RNAs, as well as reveal features of U•U mismatches flanked by Watson-Crick C•G base pairs in p19-bound and p19-unbound states. In addition, structural data reveal a p19 specific site for anchoring extra nucleotides in small RNAs. In general, the study extends our knowledge about the mechanism of small RNA recognition by silencing suppressor p19.

Studies on 4S and 5S RNA of HeLa cells

SECTION I

Section I is concerned with a partial sequence analysis conducted on 5S RNA from HeLa cells. Analysis of the oligonucleotide pattern after pancreatic ribonuclease digestion of a highly-purified preparation of 5S RNA gave results which were in general agreement with those published for KB cells, both with respect to the identity and the frequency of the partial sequences. However, the presence of a trinucleotide not found in the KB 5S pattern, together with the reproducibly much lower than expected molar yield of the larger oligonucleotides strongly suggested the occurrence of alternate sequences at various sites in the 5S molecules of human cells. The presence of ppGp and pppGp at the 5'-terminus of HeLa 5S RNA was clearly demonstrated. The implications of this finding with regard to the origin of 5S RNA are discussed.

SECTION II

In Section II the proportion of the HeLa cell genome complementary to tRNA was investigated by using RNA- DNA hybridization. The value for saturation of the HeLa DNA by tRNA was found to be 1.1 x 10^-5, which corresponds to about 4900 sites for tRNA per HeLa cell in an exponentially growing culture. Analysis of the nucleotide composition of the hybridized tRNA revealed significant differences from the nucleotide composition of the input tRNA, with the purine to pyrimidine ratio indicating, however, that these differences were not produced by excessive RNase attack of the hybrid. The size of the hybridized tRNA was only moderately smaller than that of the input RNA; the average S value in formaldehyde was 2.7 (corresponding to a length of about 65 nucleotides), suggesting that a relatively small portion near the ends of the hybridized 4S chains had been removed by RNase.

SECTION III

The proportion of the HeLa cell genome complementary to 5S RNA was investigated by using RNA-DNA hybridization. The value for saturation of the HeLa DNA by 5S RNA was found to be 2.3 x 10^-5, which corresponds to about 7,000 sites for 5S RNA per HeLa cell in an exponentially growing culture. Analysis of the nucleotide composition of the hybridized 5S RNA revealed no significant difference from the nucleotide composition of the input RNA. At the RNA to DNA input ratio of 1:1000, the average S value in formaldehyde of the hybridized 5S RNA corresponded to a polynucleotide chain about two-thirds the size of the input RNA.