The structure of human tRNALys3 anticodon bound to HIV genome is stabilized by modified nucleosides and adjacent mismatch base pairs

Replication of human immunodeficiency virus (HIV) requires base pairing of the reverse transcriptase primer, human tRNA(Lys3), to the viral RNA. Although the major complementary base pairing occurs between the HIV primer binding sequence (PBS) and the tRNA's 3'-terminus, an important discriminatory, secondary contact occurs between the viral A-rich Loop I, 5'-adjacent to the PBS, and the modified, U-rich anticodon domain of tRNA(Lys3). The importance of individual and combined anticodon modifications to the tRNA/HIV-1 Loop I RNA's interaction was determined. The thermal stabilities of variously modified tRNA anticodon region sequences bound to the Loop I of viral sub(sero)types G and B were analyzed and the structure of one duplex containing two modified nucleosides was determined using NMR spectroscopy and restrained molecular dynamics. The modifications 2-thiouridine, s(2)U(34), and pseudouridine, Psi(39), appreciably stabilized the interaction of the anticodon region with the viral subtype G and B RNAs. The structure of the duplex results in two coaxially stacked A-form RNA stems separated by two mismatched base pairs, U(162)*Psi(39) and G(163)*A(38), that maintained a reasonable A-form helix diameter. The tRNA's s(2)U(34) stabilized the interaction between the A-rich HIV Loop I sequence and the U-rich anticodon, whereas the tRNA's Psi(39) stabilized the adjacent mismatched pairs.

Estimates of electronic coupling for excess electron transfer in DNA

Electronic coupling Vda is one of the key parameters that determine the rate of charge transfer through DNA. While there have been several computational studies of Vda for hole transfer, estimates of electronic couplings for excess electron transfer (ET) in DNA remain unavailable. In the paper, an efficient strategy is established for calculating the ET matrix elements between base pairs in a π stack. Two approaches are considered. First, we employ the diabatic-state (DS) method in which donor and acceptor are represented with radical anions of the canonical base pairs adenine-thymine (AT) and guanine-cytosine (GC). In this approach, similar values of Vda are obtained with the standard 6-31 G* and extended 6-31++ G* basis sets. Second, the electronic couplings are derived from lowest unoccupied molecular orbitals (LUMOs) of neutral systems by using the generalized Mulliken-Hush or fragment charge methods. Because the radical-anion states of AT and GC are well reproduced by LUMOs of the neutral base pairs calculated without diffuse functions, the estimated values of Vda are in good agreement with the couplings obtained for radical-anion states using the DS method. However, when the calculation of a neutral stack is carried out with diffuse functions, LUMOs of the system exhibit the dipole-bound character and cannot be used for estimating electronic couplings. Our calculations suggest that the ET matrix elements Vda for models containing intrastrand thymine and cytosine bases are essentially larger than the couplings in complexes with interstrand pyrimidine bases. The matrix elements for excess electron transfer are found to be considerably smaller than the corresponding values for hole transfer and to be very responsive to structural changes in a DNA stack

Little evidence for association between the TGFBR1*6A variant and colorectal cancer: a family-based association study on non-syndromic family members from Australia and Spain.

Genome-wide linkage studies have identified the 9q22 chromosomal region as linked with colorectal cancer (CRC) predisposition. A candidate gene in this region is transforming growth factor beta receptor 1 (TGFBR1). Investigation of TGFBR1 has focused on the common genetic variant rs11466445, a short exonic deletion of nine base pairs which results in truncation of a stretch of nine alanine residues to six alanine residues in the gene product. While the six alanine (*6A) allele has been reported to be associated with increased risk of CRC in some population based study groups this association remains the subject of robust debate. To date, reports have been limited to population-based case-control association studies, or case-control studies of CRC families selecting one affected individual per family. No study has yet taken advantage of all the genetic information provided by multiplex CRC families. Methods: We have tested for an association between rs11466445 and risk of CRC using several family-based statistical tests in a new study group comprising members of non-syndromic high risk CRC families sourced from three familial cancer centres, two in Australia and one in Spain. Results: We report a finding of a nominally significant result using the pedigree-based association test approach (PBAT; p = 0.028), while other family-based tests were non-significant, but with a p-value < 0.10 in each instance. These other tests included the Generalised Disequilibrium Test (GDT; p = 0.085), parent of origin GDT Generalised Disequilibrium Test (GDT-PO; p = 0.081) and empirical Family-Based Association Test (FBAT; p = 0.096, additive model). Related-person case-control testing using the 'More Powerful' Quasi-Likelihood Score Test did not provide any evidence for association (M-QL5; p = 0.41). Conclusions: After conservatively taking into account considerations for multiple hypothesis testing, we find little evidence for an association between the TGFBR1*6A allele and CRC risk in these families. The weak support for an increase in risk in CRC predisposed families is in agreement with recent meta-analyses of case-control studies, which estimate only a modest increase in sporadic CRC risk among 6*A allele carriers.

Charge transfer in DNA: hole charge is confined to a single base pair due to solvation effects

We include solvation effects in tight-binding Hamiltonians for hole states in DNA. The corresponding linear-response parameters are derived from accurate estimates of solvation energy calculated for several hole charge distributions in DNA stacks. Two models are considered: (A) the correction to a diagonal Hamiltonian matrix element depends only on the charge localized on the corresponding site and (B) in addition to this term, the reaction field due to adjacent base pairs is accounted for. We show that both schemes give very similar results. The effects of the polar medium on the hole distribution in DNA are studied. We conclude that the effects of polar surroundings essentially suppress charge delocalization in DNA, and hole states in (GC)n sequences are localized on individual guanines

Conformations of poly{G}-poly{C} π stacks with high hole mobility

Charge transfer properties of DNA depend strongly on the π stack conformation. In the present paper, we identify conformations of homogeneous poly-{G}-poly-{C} stacks that should exhibit high charge mobility. Two different computational approaches were applied. First, we calculated the electronic coupling squared, V2, between adjacent base pairs for all 1 ps snapshots extracted from 15 ns molecular dynamics trajectory of the duplex G15. The average value of the coupling squared 〈 V2 〉 is found to be 0.0065 eV2. Then we analyze the base-pair and step parameters of the configurations in which V2 is at least an order of magnitude larger than 〈 V2 〉. To obtain more consistent data, ∼65 000 configurations of the (G:C)2 stack were built using systematic screening of the step parameters shift, slide, and twist. We show that undertwisted structures (twist<20°) are of special interest, because the π stack conformations with strong electronic couplings are found for a wide range of slide and shift. Although effective hole transfer can also occur in configurations with twist=30° and 35°, large mutual displacements of neighboring base pairs are required for that. Overtwisted conformation (twist38°) seems to be of limited interest in the context of effective hole transfer. The results may be helpful in the search for DNA based elements for nanoelectronics

Estimation of protein coding density in a corpus of DNA sequence data

A number of experimental methods have been reported for estimating the number of genes in a genome, or the closely related coding density of a genome, defined as the fraction of base pairs in codons. Recently, DNA sequence data representative of the genome as a whole have become available for several organisms, making the problem of estimating coding density amenable to sequence analytic methods. Estimates of coding density for a single genome vary widely, so that methods with characterized error bounds have become increasingly desirable. We present a method to estimate the protein coding density in a corpus of DNA sequence data, in which a ‘coding statistic’ is calculated for a large number of windows of the sequence under study, and the distribution of the statistic is decomposed into two normal distributions, assumed to be the distributions of the coding statistic in the coding and noncoding fractions of the sequence windows. The accuracy of the method is evaluated using known data and application is made to the yeast chromosome III sequence and to C.elegans cosmid sequences. It can also be applied to fragmentary data, for example a collection of short sequences determined in the course of STS mapping.

A burst of segmental duplications in the genome of the African great ape ancestor

It is generally accepted that the extent of phenotypic change between human and great apes is dissonant with the rate of molecular change. Between these two groups, proteins are virtually identical, cytogenetically there are few rearrangements that distinguish ape-human chromosomes, and rates of single-base-pair change and retrotransposon activity have slowed particularly within hominid lineages when compared to rodents or monkeys. Studies of gene family evolution indicate that gene loss and gain are enriched within the primate lineage. Here, we perform a systematic analysis of duplication content of four primate genomes (macaque, orang-utan, chimpanzee and human) in an effort to understand the pattern and rates of genomic duplication during hominid evolution. We find that the ancestral branch leading to human and African great apes shows the most significant increase in duplication activity both in terms of base pairs and in terms of events. This duplication acceleration within the ancestral species is significant when compared to lineage-specific rate estimates even after accounting for copy-number polymorphism and homoplasy. We discover striking examples of recurrent and independent gene-containing duplications within the gorilla and chimpanzee that are absent in the human lineage. Our results suggest that the evolutionary properties of copy-number mutation differ significantly from other forms of genetic mutation and, in contrast to the hominid slowdown of single-base-pair mutations, there has been a genomic burst of duplication activity at this period during human evolution.

Dugesia sicula (Platyhelminthes, Tricladida): the colonizing success of an asexual Planarian

Background Dugesia sicula is the only species of its genus not presenting an endemic or restricted distribution within the Mediterranean area. It mostly comprises fissiparous populations (asexual reproduction by body division and regeneration), most likely sexually sterile, and characterized by an extremely low genetic diversity interpreted as the consequence of a recent anthropic expansion. However, its fissiparous reproduction can result in an apparent lack of diversity within the species, since genetic variation within individuals can be as large as between them because most individuals within a population are clones. We have estimated haplotype and nucleotide diversity of cytochrome oxidase I within and among individuals along the species distribution of a broad sample of D. sicula, including asexual and the two only sexual populations known today; and predicted its potential distribution based on climatic variables. Our aim was to determine the centre of colonisation origin, whether the populations are recent, and whether the species is expanding. Results The species presents 3 most frequent haplotypes, differing in a maximum of 11 base pairs. As expected from their fissiparous mode of reproduction, in half of all the analysed localities many individuals have multiple heteroplasmic haplotypes. The distribution of haplotypes is not geographically structured; however, the distribution of haplotypes and heteroplasmic populations shows higher diversity in the central Mediterranean region. The potential distribution predicted by climatic variables based modelling shows a preference for coastal areas and fits well with the observed data. Conclusions The distribution and frequency of the most frequent haplotypes and the presence of heteroplasmic individuals allow us to gain an understanding of the recent history of the species, together with previous knowledge on its phylogenetic relationships and age: The species most probably originated in Africa and dispersed through the central Mediterranean. After one or multiple populations became triploid and fissiparous, the species colonized the Mediterranean basin, likely both by its own means and helped by human activities. Its present distribution practically fulfils its potential distribution as modelled with climatic variables. Its prevalence in coastal regions with higher water temperatures predicts a likely future expansion to northern and more interior areas following the increase in temperatures due to climate change.

A minimal i-motif stabilized by minor groove G:T:G:T tetrads

The repetitive DNA sequences found at telomeres and centromeres play a crucial role in the structure and function of eukaryotic chromosomes. This role may be related to the tendency observed in many repetitive DNAs to adopt non-canonical structures. Although there is an increasing recognition of the importance of DNA quadruplexes in chromosome biology, the co-existence of different quadruplex-forming elements in the same DNA structure is still a matter of debate. Here we report the structural study of the oligonucleotide d(TCGTTTCGT) and its cyclic analog d. Both sequences form dimeric quadruplex structures consisting of a minimal i-motif capped, at both ends, by a slipped minor groove-aligned G:T:G:T tetrad. These mini i-motifs, which do not exhibit the characteristic CD spectra of other i-motif structures, can be observed at neutral pH, although they are more stable under acidic conditions. This finding is particularly relevant since these oligonucleotide sequences do not contain contiguous cytosines. Importantly, these structures resemble the loop moiety adopted by an 11-nucleotide fragment of the conserved centromeric protein B (CENP-B) box motif, which is the binding site for the CENP-B.

A minimal i-motif stabilized by minor groove G:T:G:T tetrads

The repetitive DNA sequences found at telomeres and centromeres play a crucial role in the structure and function of eukaryotic chromosomes. This role may be related to the tendency observed in many repetitive DNAs to adopt non-canonical structures. Although there is an increasing recognition of the importance of DNA quadruplexes in chromosome biology, the co-existence of different quadruplex-forming elements in the same DNA structure is still a matter of debate. Here we report the structural study of the oligonucleotide d(TCGTTTCGT) and its cyclic analog d. Both sequences form dimeric quadruplex structures consisting of a minimal i-motif capped, at both ends, by a slipped minor groove-aligned G:T:G:T tetrad. These mini i-motifs, which do not exhibit the characteristic CD spectra of other i-motif structures, can be observed at neutral pH, although they are more stable under acidic conditions. This finding is particularly relevant since these oligonucleotide sequences do not contain contiguous cytosines. Importantly, these structures resemble the loop moiety adopted by an 11-nucleotide fragment of the conserved centromeric protein B (CENP-B) box motif, which is the binding site for the CENP-B.

Molecular dating of caprines using ancient DNA sequences of Myotragus balearicus, an extinct endemic Balearic mammal

Background: Myotragus balearicus was an endemic bovid from the Balearic Islands (Western Mediterranean) that became extinct around 6,000-4,000 years ago. The Myotragus evolutionary lineage became isolated in the islands most probably at the end of the Messinian crisis, when the desiccation of the Mediterranean ended, in a geological date established at 5.35 Mya. Thus, the sequences of Myotragus could be very valuable for calibrating the mammalian mitochondrial DNA clock and, in particular, the tree of the Caprinae subfamily, to which Myotragus belongs. Results: We have retrieved the complete mitochondrial cytochrome b gene (1,143 base pairs), plus fragments of the mitochondrial 12S gene and the nuclear 28S rDNA multi-copy gene from a well preserved Myotragus subfossil bone. The best resolved phylogenetic trees, obtained with the cytochrome b gene, placed Myotragus in a position basal to the Ovis group. Using the calibration provided by the isolation of Balearic Islands, we calculated that the initial radiation of caprines can be dated at 6.2 ± 0.4 Mya. In addition, alpine and southern chamois, considered until recently the same species, split around 1.6 ± 0.3 Mya, indicating that the two chamois species have been separated much longer than previously thought. Conclusion: Since there are almost no extant endemic mammals in Mediterranean islands, the sequence of the extinct Balearic endemic Myotragus has been crucial for allowing us to use the Messinian crisis calibration point for dating the caprines phylogenetic tree.

A novel source for miR-21 expression through the alternative polyadenylation of VMP1 gene transcripts

miR-21 is the most commonly over-expressed microRNA (miRNA) in cancer and a proven oncogene. Hsa-miR-21 is located on chromosome 17q23.2, immediately downstream of the vacuole membrane protein-1 (VMP1) gene, also known as TMEM49. VMP1 transcripts initiate ∼130 kb upstream of miR-21, are spliced, and polyadenylated only a few hundred base pairs upstream of the miR-21 hairpin. On the other hand, primary miR-21 transcripts (pri-miR-21) originate within the last introns of VMP1, but bypass VMP1 polyadenylation signals to include the miR-21 hairpin. Here, we report that VMP1 transcripts can also bypass these polyadenylation signals to include miR-21, thus providing a novel and independently regulated source of miR-21, termed VMP1–miR-21. Northern blotting, gene-specific RT-PCR, RNA pull-down and DNA branching assays support that VMP1–miR-21 is expressed at significant levels in a number of cancer cell lines and that it is processed by the Microprocessor complex to produce mature miR-21. VMP1 and pri-miR-21 are induced by common stimuli, such as phorbol-12-myristate-13-acetate (PMA) and androgens, but show differential responses to some stimuli such as epigenetic modifying agents. Collectively, these results indicate that miR-21 is a unique miRNA capable of being regulated by alternative polyadenylation and two independent gene promoters.

Complexes of Pd(II) and Pt(II) with 9-aminoacridine: reactions with DNA and study of their antiproliferative activity

Four new metal complexes {M = Pd(II) or Pt(II)} containing the ligand 9-aminoacridine (9AA) were prepared. The compounds were characterized by FT-IR and 1H, 13C, and 195Pt NMR spectroscopies. Crystal structure of the palladium complex of formulae [Pd(9AA)(μ-Cl)]2 · 2DMF was determined by X-ray diffraction. Two 9-acridine molecules in the imine form bind symmetrically to the metal ions in a bidentate fashion through the imine nitrogen atom and the C(1) atom of the aminoacridine closing a new five-membered ring. By reaction with phosphine or pyridine, the Cl bridges broke and compounds with general formulae [Pd(9AA)Cl(L)] (where L = PPh3 or py) were formed. A mononuclear complex of platinum of formulae [Pt(9AA)Cl(DMSO)] was also obtained by direct reaction of 9-aminoacridine and the complex [PtCl2(DMSO)2]. The capacity of the compounds to modify the secondary and tertiary structures of DNA was evaluated by means of circular dichroism and electrophoretic mobility. Both palladium and platinum compounds proved active in the modification of both the secondary and tertiary DNA structures. AFM images showed noticeable modifications of the morphology of the plasmid pBR322 DNA by the compounds probably due to the intercalation of the complexes between base pairs of the DNA molecule. Finally, the palladium complex was tested for antiproliferative activity against three different human tumor cell lines. The results suggest that the palladium complex of formula [Pd(9AA)(μ-Cl)]2 has significant antiproliferative activity, although it is less active than cisplatin.

Complexes of Pd(II) and Pt(II) with 9-aminoacridine: reactions with DNA and study of their antiproliferative activity

Four new metal complexes {M = Pd(II) or Pt(II)} containing the ligand 9-aminoacridine (9AA) were prepared. The compounds were characterized by FT-IR and 1H, 13C, and 195Pt NMR spectroscopies. Crystal structure of the palladium complex of formulae [Pd(9AA)(μ-Cl)]2 · 2DMF was determined by X-ray diffraction. Two 9-acridine molecules in the imine form bind symmetrically to the metal ions in a bidentate fashion through the imine nitrogen atom and the C(1) atom of the aminoacridine closing a new five-membered ring. By reaction with phosphine or pyridine, the Cl bridges broke and compounds with general formulae [Pd(9AA)Cl(L)] (where L = PPh3 or py) were formed. A mononuclear complex of platinum of formulae [Pt(9AA)Cl(DMSO)] was also obtained by direct reaction of 9-aminoacridine and the complex [PtCl2(DMSO)2]. The capacity of the compounds to modify the secondary and tertiary structures of DNA was evaluated by means of circular dichroism and electrophoretic mobility. Both palladium and platinum compounds proved active in the modification of both the secondary and tertiary DNA structures. AFM images showed noticeable modifications of the morphology of the plasmid pBR322 DNA by the compounds probably due to the intercalation of the complexes between base pairs of the DNA molecule. Finally, the palladium complex was tested for antiproliferative activity against three different human tumor cell lines. The results suggest that the palladium complex of formula [Pd(9AA)(μ-Cl)]2 has significant antiproliferative activity, although it is less active than cisplatin.

Complexes of Pd(II) and Pt(II) with 9-aminoacridine: reactions with DNA and study of their antiproliferative activity

Four new metal complexes {M = Pd(II) or Pt(II)} containing the ligand 9-aminoacridine (9AA) were prepared. The compounds were characterized by FT-IR and 1H, 13C, and 195Pt NMR spectroscopies. Crystal structure of the palladium complex of formulae [Pd(9AA)(μ-Cl)]2 · 2DMF was determined by X-ray diffraction. Two 9-acridine molecules in the imine form bind symmetrically to the metal ions in a bidentate fashion through the imine nitrogen atom and the C(1) atom of the aminoacridine closing a new five-membered ring. By reaction with phosphine or pyridine, the Cl bridges broke and compounds with general formulae [Pd(9AA)Cl(L)] (where L = PPh3 or py) were formed. A mononuclear complex of platinum of formulae [Pt(9AA)Cl(DMSO)] was also obtained by direct reaction of 9-aminoacridine and the complex [PtCl2(DMSO)2]. The capacity of the compounds to modify the secondary and tertiary structures of DNA was evaluated by means of circular dichroism and electrophoretic mobility. Both palladium and platinum compounds proved active in the modification of both the secondary and tertiary DNA structures. AFM images showed noticeable modifications of the morphology of the plasmid pBR322 DNA by the compounds probably due to the intercalation of the complexes between base pairs of the DNA molecule. Finally, the palladium complex was tested for antiproliferative activity against three different human tumor cell lines. The results suggest that the palladium complex of formula [Pd(9AA)(μ-Cl)]2 has significant antiproliferative activity, although it is less active than cisplatin.