Interaction and protection mechanism between Li@C(60) and nucleic acid bases (NABs): performance of PM6-DH2 on noncovalent interaction of NABs-Li@C60

To discuss the protection mechanism of DNA from radiation as well as assess the performance of PM6-DH2 on noncovalent interactions, the interaction of four nucleic acid bases (NABs) such as adenine (A), cytosine (C), guanine (G), and thymine (T), with Li@C(60) was extensively investigated with the-state-of-art theoretical methods describing noncovalent systems, like M06-2x, PBE-D, and PM6-DH2 methods. In the gas phase, the binding strength of NABs to Li@C(60) from M06-2x decreases in the sequence, G>C>A>T. As dispersion was explicitly included, PBE-D relatively enhances the binding of A and T and corrects the sequence to, G>A>C～T. PM6-DH2 predicted similar binding energies to those from PBE-D within 0.5 kcal/mol and the same binding sequence, suggesting that the PM6-DH2 method is promising for nano-scale systems. In the aqueous solution, binding of NABs-Li@C(60) is considerably decreased, and the M06-2X and PM6-D methods yield a different sequence from the gas phase, G>A>T>C. The encapsulation of Li atom results in a lower IP for Li@C(60) than those of NABs, and the dominant localization of single-occupied molecular orbital on Li@C(60) moiety of the complexes NABs-Li@C(60) further indicates that an electron would be ejected from Li@C(60) upon radiation and Li@C(60) is therefore able to protect DNA bases from radiation. In addition, it was revealed that Li prefers coordination with the hexagonal ring at Li@C(60) , which clarifies the existing controversy in this respect. Finally, Yang's reduced density gradient approach clearly shows that the weak and strong noncovalent interaction regions in the complexes, NABs-Li@C(60) and (NABs-Li@C(60) )(+).     

Theoretical Evidence for the Stronger Ability of Thymine to Disperse SWCNT than Cytosine and Adenine: self-stacking of DNA bases vs their cross-stacking with SWCNT

Self-stacking of four DNA bases, adenine (A), cytosine (C), guanine (G) and thymine (T), and their cross-stacking with (5,5) as well as (10,0) single walled carbon nanotubes (SWCNTs) were extensively investigated with a novel hybrid DFT method, MPWB1K/cc-pVDZ. The binding energies were further corrected with MP2/6-311++G(d,p) method in both gas phase and aqueous solution, where the solvent effects were included with conductor-like polarized continuum model (CPCM) model and UAHF radii. The strongest self-stacking of G and A takes displaced anti-parallel configuration, but un-displaced or "eclipsed" anti-parallel configuration is the most stable for C and T. In gas phase the self-stacking of nucleobases decreases in the sequence G>A>C>T, while because of quite different solvent effects their self-stacking in aqueous solution exhibits a distinct sequence A>G>T>C. For a given base, cross-stacking is stronger than self-stacking in both gas phase and aqueous solution. Binding energy for cross-stacking in gas phase varies as G>A>T>C for both (10,0) and (5,5) SWCNTs, and the binding of four nucleobases to (10,0) is slightly stronger than to (5,5) SWCNT by a range of 0.1-0.5 kcal/mol. The cross-stacking in aqueous solution varies differently from that gas phase: A>G>T>C for (10,0) SWCNT and G>A>T>C for (5,5) SWCNT. It is suggested that the ability of nucleobases to disperse SWCNT depends on relative strength [Formula: see text] of self-stacking and cross-stacking with SWCNT in aqueous solution. Of the four investigated nucleobases thymine (T) exhibits the highest [Formula: see text] which can well explain the experimental finding that T more efficiently functionalizes SWCNT than C and A.     

Electrophoretic analysis of sequence variability in three mitochondrial DNA regions for ascaridoid parasites of human and animal health significance

Sequence variability in three mitochondrial DNA (mtDNA) regions, namely cytochrome c oxidase subunit 1 (cox1), NADH dehydrogenase subunits 1 and 4 (nad1 and nad4), among and within Toxocara canis, T. cati, T. malaysiensis, T. vitulorum and Toxascaris leonina from different geographical origins was examined by a mutation-scanning approach. A portion of the cox1 gene (pcox1), a portion of the nad1 and nad4 genes (pnad1 and pnad4) were amplified separately from individual ascaridoid nematodes by polymerase chain reaction and the amplicons analyzed by single-strand conformation polymorphism (SSCP). Representative samples displaying sequence variation in SSCP profiles were subjected to sequencing in order to define genetic markers for their specific identification and differentiation. While the intra-specific sequence variations within each of the five ascaridoid species were 0.2-3.7% for pcox1, 0-2.8% for pnad1 and 0-2.3% for pnad4, the inter-specific sequence differences were significantly higher, being 7.9-12.9% for pcox1, 10.7-21.1% for pnad1 and 12.9-21.7% for pnad4, respectively. Phylogenetic analyses based on the combined sequences of pcox1, pnad1 and pnad4 revealed that the recently described species T. malaysiensis was more closely related to T. cati than to T. canis. These findings provided mtDNA evidence for the validity of T. malaysiensis and also demonstrated clearly the usefulness and attributes of the mutation-scanning sequencing approach for studying the population genetic structures of these and other nematodes of socio-economic importance.     

Design of allele-specific primers and detection of the human ABO genotyping to avoid the pseudopositive problem

PCR experiments using DNA primers forming mismatch pairing with template lambda DNA at the 3' end were carried out in order to develop allele-specific primers capable of detecting SNP in genomes without generating pseudopositive amplification products, and thus avoiding the so-called pseudopositive problem. Detectable amounts of PCR products were obtained when primers forming a single or two mismatch pairings at the 3' end were used. In particular, 3' terminal A/C or T/C (primer/template) mismatches tended to allow PCR amplification to proceed, resulting in pseudopositive results in many cases. While less PCR product was observed for primers forming three terminal mismatch pairings, target DNA sequences were efficiently amplified by primers forming two mismatch pairings next to the terminal G/C base pairing. These results indicate that selecting a primer having a 3' terminal nucleotide that recognizes the SNP nucleotide and the next two nucleotides that form mismatch pairings with the template sequence can be used as an allele-specific primer that eliminates the pseudopositive problem. Trials with the human ABO genes demonstrated that this primer design is also useful for detecting a single base pair difference in gene sequences with a signal-to-noise ratio of at least 45.     

High-speed and high-accuracy DNA sequencing by capillary gel electrophoresis in a simple, low cost instrument Two-color peak-height encoded sequencing at 40°C

A low-cost DNA sequencer was constructed based on a single helium-neon laser. The two-color peak-height encoded sequencing protocol, based on the use of T7 DNA polymerase in a manganese buffer, was used to generate samples. Two termination reactions were performed. In the first, a TAMRA (applied Biosystems)-labeled primer was extended in the presence of ddATP and ddCTP. The amounts of dideoxynucleotides were adjusted to produce a 3:1 peak height ratio. Similarly, a ROX (Applied Biosystems)-labeled primer was extended in the presence of ddGTP and ddTTP; the amounts of dideoxynucleotides was adjusted to produce a 3:1 peak height ratio. The pooled fragments were separated on a 4% T LongRanger gel operated at 39°C. Over 500 bases of sequence were generated in 50 min.     

Analysis and interpretation of mixture DNA using AS‐PCR of mtDNA

Semi‐nested PCR with allele‐specific (AS) primers and sequencing of mitochondrial DNA (mtDNA) were performed to analyze and interpret DNA mixtures, especially when biological materials were degraded or contained a limited amount of DNA. SNP‐STR markers were available to identify the minor DNA component using AS‐PCR; moreover, SNPs in mtDNA could be used when the degraded or limited amounts of DNA mixtures were not successful with SNP‐STR markers. Five pairs of allele‐specific primers were designed based on three SNPs (G15043A, T16362C, and T16519C). The sequence of mtDNA control region of minor components was obtained using AS‐PCR and sequencing. Sequences of the amplification fragments were aligned and compared with the sequences of known suspects or databases. When this assay was used with the T16362C and T16519C SNPs, we found it to be highly sensitive for detecting small amounts of DNA (～30 pg) and analyzing DNA mixtures of two contributors, even at an approximately 1‰ ratio of minor and major components. An exception was tests based on the SNP G15043A, which required approximately 300 pg of a 1% DNA mixture. In simulated three contributor DNA mixtures (at rate of 1:1:1), control region fragments from each contributor were detected and interpreted. AS‐PCR combined with semi‐nested PCR was successfully used to identify the mtDNA control region of each contributor, providing biological evidence for excluding suspects in forensic cases, especially when biological materials were degraded or had a limited amount of DNA.     

Targeted next-generation sequencing by specific capture of multiple genomic loci using low-volume microfluidic DNA arrays

We report a flexible method for selective capture of sequence fragments from complex, eukaryotic genome libraries for next-generation sequencing based on hybridization to DNA microarrays. Using microfluidic array architecture and integrated hardware, the process is amenable to complete automation and does not introduce amplification steps into the standard library preparation workflow, thereby avoiding bias of sequence distribution and fragment lengths. We captured a discontiguous human genomic target region of 185 kb using a tiling design with 50mer probes. Analysis by high-throughput sequencing using an Illumina/Solexa 1G Genome Analyzer revealed 2150-fold enrichment with mean per base coverage between 4.6 and 107.5-fold for the individual target regions. This method represents a flexible and cost-effective approach for large-scale resequencing of complex genomes. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users. Stephan Bau and Nadine Schracke contributed equally to this work.     

Multiple-primer DNA sequencing method.

A multiple-primer DNA sequencing approach suitable for genotyping, detection and identification of microorganisms and viruses has been developed. In this new method two or more sequencing primers, combined in a pool, are added to a DNA sample of interest. The oligonucleotide that hybridizes to the DNA sample will function as a primer during the subsequent DNA sequencing procedure. This strategy is suited for selective detection and genotyping of relevant microorganisms and samples harboring different DNA targets such as multiple variant/infected samples as well as unspecific amplification products. This method is used here in a model system for detection and typing of high-risk oncogenic human papilloma viruses (HPVs) in samples containing multiple infections/variants or unspecific amplification products. Type-specific sequencing primers were designed for four of the most oncogenic (high-risk) HPV types (HPV-16, HPV-18, HPV-33, and HPV-45). The primers were combined and added to a sample containing a mixture of one high-risk (16, 18, 33, or 45) and one or two low-risk types. The DNA samples were sequenced by the Pyrosequencing technology and the Sanger dideoxy sequencing method. Correct genotyping was achieved in all tested combinations. This multiple-sequencing primer approach also improved the sequence data quality for samples containing unspecific amplification products. The new strategy is highly suitable for diagnostic typing of relevant species/genotypes of microorganisms.     

FeII Metallohelices Stabilize DNA G-Quadruplexes and Downregulate the Expression of G-Quadruplex-Regulated Oncogenes

DNA G-quadruplexes (G4s) have been identified within the promoter regions of many proto-oncogenes. Thus, G4s represent attractive targets for cancer therapy, and the design and development of new drugs as G4 binders is a very active field of medicinal chemistry. Here, molecular biophysics and biology methods were employed to investigate the interaction of chiral metallohelices with a series of four DNA G4s (hTelo, c-myc, c-kit1, c-kit2) that are formed by the human telomeric sequence (hTelo) and in the promoter regions of c-MYC and c-KIT proto-oncogenes. We show that the investigated water-compatible, optically pure metallohelices, which are made by self-assembly of simple nonpeptidic organic components around Fe^II ions and exhibit bioactivity emulating the natural systems, bind with high affinity to G4 DNA and much lower affinity to duplex DNA. Notably, both enantiomers of a metallohelix containing a m-xylenyl bridge ( 5 b ) were found to effectively inhibit primer elongation catalyzed by Taq DNA polymerase by stabilizing G4 structures formed in the template strands containing c-myc and c-kit2 G4-forming sequences. Moreover, both enantiomers of 5 b downregulated the expression of c-MYC and c-KIT oncogenes in human embryonic kidney cells at mRNA and protein levels. As metallohelices also bind alternative nucleic acid structures, they hold promise as potential multitargeted drugs.     

Low-cost, high-sensitivity laser-induced fluorescence detection for DNA sequencing by capillary gel electrophoresis.

A low cost, 0.75-mW helium neon laser, operating in the green region at 534.5 nm, is used to excite fluorescence from tetramethylrhodamine isothiocyanate-labelled DNA fragments that have been separated by capillary gel electrophoresis. The detection limit (3 sigma) for the dye is 500 ymol [1 yoctomole (1 ymol) = 10(-24) mol] or 300 analyte molecules in capillary zone electrophoresis; the detection limit for labeled primer separated by capillary gel electrophoresis is 2 zmol [1 zeptomole (1 zmol) = 10(-21) mol]. The Richardson-Tabor peak-height encoded sequencing technique is used to prepare DNA sequencing samples. In 6% T, 5% C acrylamide, 7 M urea gels, sequencing rates of 300 bases/hour are produced at an electric field strength of 200 V/cm; unfortunately, the data are plagued by compressions. These compressions are eliminated with addition of 20% formamide to the sequencing gel; the gel runs slowly and sequencing data are generated at a rate of about 70 bases/hour.     

Identification of epigenetic DNA modifications with a protein nanopore

Two DNA bases, 5-methylcytosine (5mC) and 5-hydroxymethylcytosine (hmC), marks of epigenetic modification, are recognized in immobilized DNA strands and distinguished from G, A, T and C by nanopore current recording. Therefore, if further aspects of nanopore sequencing can be addressed, the approach will provide a means to locate epigenetic modifications in unamplified genomic DNA.     

On the Binding Strength Sequence for Nucleic Acid Bases and C(60) with Density Functional and Dispersion-corrected Density Functional Theories: Whether C(60) could protect nucleic acid bases from radiation-induced damage?

The major objective of this paper is to address a controversial binding sequence between nucleic acid bases (NABs) and C(60) by investigating adsorptions of NABs and their cations on C(60) fullerene with a variety of density functional theories including two novel hybrid meta-GGA functionals, M05-2x and M06-2x, as well as a dispersion-corrected density functional, PBE-D. The M05-2x/6-311++G** provides the same binding sequence as previously reported, guanine(G) > cytosine(C) > adenine (A) > thymine (T); however, M06-2x switches the binding strengths of A and C, and PBE-D eventually results in the following sequence, G>A>T>C, which is the same as the widely accepted hierarchy for the stacking of NABs on other carbon nanomaterials such as single-walled carbon nanotube and graphite. The results indicate that the questionable relative binding strength is due to insufficient electron correlation treatment with the M05-2x or even the M06-2x method. The binding energy of G@C(60) obtained with the M06-2x/6-311++G(d,p) and the PBE-D/cc-pVDZ is -7.10 and -8.07 kcal/mol, respectively, and the latter is only slightly weaker than that predicted by the MP2/6-31G(d,p) (-8.10kca/mol). Thus, the PDE-D performs better than the M06-2x for the observed NAB@C(60) π-stacked complexes. To discuss whether C(60) could prevent NABs from radiation-induced damage, ionization potentials of NABs and C(60), and frontier molecular orbitals of the complexes NABs@C(60) and (NABs@C(60))(+) are also extensively investigated. These results revealed that when an electron escapes from the complexes, a hole was preferentially created in C(60) for T and C complexes, while for G and A the hole delocalizes over the entire complex, rather than a localization on the C(60) moiety. The interesting finding might open a new strategy for protecting DNA from radiation-induced damage and offer a new idea for designing C(60)-based antiradiation drugs.     

The electrophoretic mobility of DNA fragments differing by a single 3′‐terminal nucleotide in an automated capillary DNA sequencer

The electrophoretic mobility of DNA fragments that differ by a single 3′‐terminal nucleotide was assessed by capillary electrophoresis. This was accomplished using dideoxy sequencing with a 5′‐fluorescently labelled primer to generate DNA fragments with 3′‐hydrogen ends. The resulting DNA fragments were electrophoresed on the ABI 3730 automated capillary sequencer, and the data were analysed with the GeneMapper software to determine the electrophoretic mobility differences on addition of a 3′‐terminal nucleotide. It was found that the 3′‐terminal nucleotide gave rise to different electrophoretic mobility profiles depending on the identity of the terminal nucleotide. The apparent electrophoretic mobility was (faster) –C > ?A > ?T > ?G (slower). The C‐terminated fragments were the fastest and the G‐terminated fragments the slowest, relative to other nucleotides. It was proposed that the terminal nucleotide effect was due to changes in partial net charges on the nucleotides that resulted in alterations in the electrophoretic mobility of the DNA fragments in the automated capillary DNA sequencer. Other alternative explanations are also discussed. Copyright © 2012 John Wiley & Sons, Ltd.     

CHARACTERISATION OF TWO MONOCLONAL ANTIBODIES SPECIFIC FOR DIMERISED AND NON-DIMERISED ADJACENT THYMIDINES IN SINGLE STRANDED DNA*

Abstract— Hybrid cell lines (hybridomas) have been isolated from fusions between P3-NS1-1-Ag4-1 mouse myeloma cells and spleen cells from BALB/c mice hyperimmunised with UV-irradiated single-stranded DNA (UVssDNA) and UV-irradiated polydeoxythymidylic acid (UVpolydT). Monoclonal antibodies from two different hybridomas are characterised in the present report by competitive inhibition with different synthetic polynucleotides and oligonucleotides. The first antibody, designated αUVssDNA-1, recognises thymidine dimers in a polynucleotide or an oligonucleotide sequence at least four nucleotides long but not isolated thymidine dimers, suggesting that it recognises the conformational change associated with thymidine dimers. The second antibody, designated αssDNA-2, recognises unirradiated or UV-irradiated tracts of thymidine, but will not crossreact with tracts of other nucleotides (A, G, C, A.T, G.C, C.U, U). Inhibition of binding of αUVssDNA-1 to [³H]-UVssDNA by calf thymus UVssDNA is dependent on UV exposure and wavelength as expected from the antigenic determinant.     

Analysis of the variations in IL-28RA gene and their association with allergic rhinitis

IL-28RA is one of the important candidate genes for complex trait of genetic diseases, but there is no published information of the genetic variation in this gene. We scanned the seven exons and their boundary introns sequence of IL-28RA including the promoter regions to analyze genetic variation sites, and identified eighteen single nucleotide polymorphisms (SNPs) and two variation sites. We chose seven SNPs (g.-1193 A>C, g.-30 C>T, g.17654 C>T, g.27798 A>G, g.31265 C>T, g.31911 C>T and g.32349 G>A) of them for large sample size genotyping, and assessed the association of genotype and allele frequencies of these SNPs between allergic rhinitis patients and non-allergic rhinitis controls. We also compared the genotype frequencies between Korean controls and Han Chinese control or Korean Chinese control. We investigated the frequencies of haplotype constructed by these SNPs between allergic rhinitis patients and non-allergic rhinitis controls. Our results suggested that the g.32349 G>A polymorphism of IL-28RA might be associated with susceptibility to allergic rhinitis (P=0.032), but seems to have no relationship with serum total IgE levels. The haplotype frequencies by these SNPs also show significant association between controls and allergic rhinitis patients.     

Identification of methylated regions with peak search based on Poisson model from massively parallel methylated DNA immunoprecipitation-sequencing data

DNA methylation is one of the most important epigenetic modification types, which plays a critical role in gene expression. High efficient surveying of whole genome DNA methylation has been aims of many researchers for long. Recently, the rapidly developed massively parallel DNA‐sequencing technologies open the floodgates to vast volumes of sequence data, enabling a paradigm shift in profiling the whole genome methylation. Here, we describe a strategy, combining methylated DNA immunoprecipitation sequencing with peak search to identify methylated regions on a whole‐genome scale. Massively parallel methylated DNA immunoprecipitation sequencing combined with methylation DNA immunoprecipitation was adopted to obtain methylated DNA sequence data from human leukemia cell line K562, and the methylated regions were identified by peak search based on Poison model. From our result, 140 958 non‐overlapping methylated regions have been identified in the whole genome. Also, the credibility of result has been proved by its strong correlation with bisulfite‐sequencing data (Pearson R²=0.92). It suggests that this method provides a reliable and high‐throughput strategy for whole genome methylation identification.     

Sequence-selective single-molecule alkylation with a pyrrole-imidazole polyamide visualized in a DNA nanoscaffold

We demonstrate a novel strategy for visualizing sequence-selective alkylation of target double-stranded DNA (dsDNA) using a synthetic pyrrole-imidazole (PI) polyamide in a designed DNA origami scaffold. Doubly functionalized PI polyamide was designed by introduction of an alkylating agent 1-(chloromethyl)-5-hydroxy-1,2-dihydro-3H-benz[e]indole (seco-CBI) and biotin for sequence-selective alkylation at the target sequence and subsequent streptavidin labeling, respectively. Selective alkylation of the target site in the substrate DNA was observed by analysis using sequencing gel electrophoresis. For the single-molecule observation of the alkylation by functionalized PI polyamide using atomic force microscopy (AFM), the target position in the dsDNA (～200 base pairs) was alkylated and then visualized by labeling with streptavidin. Newly designed DNA origami scaffold named "five-well DNA frame" carrying five different dsDNA sequences in its cavities was used for the detailed analysis of the sequence-selectivity and alkylation. The 64-mer dsDNAs were introduced to five individual wells, in which target sequence AGTXCCA/TGGYACT (XY = AT, TA, GC, CG) was employed as fully matched (X = G) and one-base mismatched (X = A, T, C) sequences. The fully matched sequence was alkylated with 88% selectivity over other mismatched sequences. In addition, the PI polyamide failed to attach to the target sequence lacking the alkylation site after washing and streptavidin treatment. Therefore, the PI polyamide discriminated the one mismatched nucleotide at the single-molecule level, and alkylation anchored the PI polyamide to the target dsDNA.     

Mutation scanning for sequence variation in three mitochondrial DNA regions for members of the Contracaecum osculatum (Nematoda: Ascaridoidea) complex

Anisakid nematodes of seals from different geographical origins, previously identified by multilocus enzyme electrophoresis as Contracaecum osculatum A (CoA), C. osculatum B (CoB), C. osculatum C (CoC), C. osculatum D (CoD), C. osculatum E (CoE) and C. osculatum baicalensis (Cob), were characterised genetically using a mutation scanning approach, in order to define genetic markers for their specific identification and differentiation. Three mitochondrial DNA (mtDNA) regions, namely cytochrome c oxidase subunit I (COI), and the small and large subunits of rRNA (ssrRNA and IsrRNA, respectively) were amplified separately from individual nematodes by polymerase chain reaction (PCR), analysed by single-strand conformation polymorphism (SSCP), and samples displaying sequence variability were subjected to sequencing. Forty-six haplotypes were defined for 62-66 individuals (representing the six members of C. osculatum). All taxa except CoD and CoE could be identified, or delineated from one another, by nucleotide differences in the COI, ssrRNA and/or IsrRNA sequences. For all three mtDNA regions, 4 (10.5%), 7 (18.4%), 15 (39.5%) and 11 (28.9%) of 38 nucleotide positions were considered diagnostic (fixed) and could thus unequivocally delineate CoA, CoB, CoC and Cob. The lack of an unequivocal nucleotide difference in any of the three mtDNA sequences between CoD and CoE was in accordance with previous ribosomal DNA sequence data but inconsistent with multilocus enzyme electrophoretic data. Using all fixed nucleotide positions, CoA, CoD/E and CoB were genetically more similar to Cob than each was to CoC, similar to previous findings. In spite of not being able to distinguish among all six taxa of C. osculatum, the present study demonstrated clearly the usefulness and attributes of the mutation scanning approach for investigating population genetic structures of species of parasitic nematodes.     

Tabu search algorithm for DNA sequencing by hybridization with isothermic libraries

In this paper, a problem of isothermic DNA sequencing by hybridization (SBH) is considered. In isothermic SBH a new type of oligonucleotide libraries is used. The library consists of oligonucleotides of different lengths depending on an oligonucleotide content. It is assumed that every oligonucleotide in such a library has an equal melting temperature. Each nucleotide adds its increment to the oligonucleotide temperature and it is assumed that A and T add 2 degrees C and C and G add 4 degrees C. The hybridization experiment using isothermic libraries should provide data with a lower number of errors due to an expected similarity of melting temperatures. From the computational point of view the problem of isothermic DNA sequencing with errors is hard, similarly like its classical counterpart. Hence, there is a need for developing heuristic algorithms that construct good suboptimal solutions. The aim of the paper is to propose a heuristic algorithm based on tabu search approach. The algorithm solves the problem with both positive and negative errors. Results of an extensive computational experiment are presented, which prove the high quality of the proposed method.