Electronic couplings and on-site energies for hole transfer in DNA: systematic quantum mechanical/molecular dynamic study

The electron hole transfer (HT) properties of DNA are substantially affected by thermal fluctuations of the π stack structure. Depending on the mutual position of neighboring nucleobases, electronic coupling V may change by several orders of magnitude. In the present paper, we report the results of systematic QM/molecular dynamic (MD) calculations of the electronic couplings and on-site energies for the hole transfer. Based on 15 ns MD trajectories for several DNA oligomers, we calculate the average coupling squares 〈 V2 〉 and the energies of basepair triplets X G+ Y and X A+ Y, where X, Y=G, A, T, and C. For each of the 32 systems, 15 000 conformations separated by 1 ps are considered. The three-state generalized Mulliken-Hush method is used to derive electronic couplings for HT between neighboring basepairs. The adiabatic energies and dipole moment matrix elements are computed within the INDO/S method. We compare the rms values of V with the couplings estimated for the idealized B -DNA structure and show that in several important cases the couplings calculated for the idealized B -DNA structure are considerably underestimated. The rms values for intrastrand couplings G-G, A-A, G-A, and A-G are found to be similar, ∼0.07 eV, while the interstrand couplings are quite different. The energies of hole states G+ and A+ in the stack depend on the nature of the neighboring pairs. The X G+ Y are by 0.5 eV more stable than X A+ Y. The thermal fluctuations of the DNA structure facilitate the HT process from guanine to adenine. The tabulated couplings and on-site energies can be used as reference parameters in theoretical and computational studies of HT processes in DNA

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.

Assessment of microbial community structure changes by amplified ribosomal DNA restriction analysis (ARDRA)

Amplified ribosomal DNA restriction analysis (ARDRA) is a simple method based on restriction endonuclease digestion of the amplified bacterial 16S rDNA. In this study we have evaluated the suitability of this method to detect differences in activated sludge bacterial communities fed on domestic or industrial wastewater, and subject to different operational conditions. The ability of ARDRA to detect these differences has been tested in modified Ludzack-Ettinger (MLE) configurations. Samples from three activated sludge wastewater treatment plants (WWTPs) with the MLE configuration were collected for both oxic and anoxic reactors, and ARDRA patterns using double enzyme digestions AluI+MspI were obtained. A matrix of Dice similarity coefficients was calculated and used to compare these restriction patterns. Differences in the community structure due to influent characteristics and temperature could be observed, but not between the oxic and anoxic reactors of each of the three MLE configurations. Other possible applications of ARDRA for detecting and monitoring changes in activated sludge systems are also discussed

Elastic properties and secondary structure formation of single-stranded DNA at monovalent and divalent salt conditions

Single-stranded DNA (ssDNA) plays a major role in several biological processes. It is therefore of fundamental interest to understand how the elastic response and the formation of secondary structures are modulated by the interplay between base pairing and electrostatic interactions. Here we measure force-extension curves (FECs) of ssDNA molecules in optical tweezers set up over two orders of magnitude of monovalent and divalent salt conditions, and obtain its elastic parameters by fitting the FECs to semiflexible models of polymers. For both monovalent and divalent salts, we find that the electrostatic contribution to the persistence length is proportional to the Debye screening length, varying as the inverse of the square root of cation concentration. The intrinsic persistence length is equal to 0.7 nm for both types of salts, and the effectivity of divalent cations in screening electrostatic interactions appears to be 100-fold as compared with monovalent salt, in line with what has been recently reported for single-stranded RNA. Finally, we propose an analysis of the FECs using a model that accounts for the effective thickness of the filament at low salt condition and a simple phenomenological description that quantifies the formation of non-specific secondary structure at low forces.

Elastic properties and secondary structure formation of single-stranded DNA at monovalent and divalent salt conditions

Single-stranded DNA (ssDNA) plays a major role in several biological processes. It is therefore of fundamental interest to understand how the elastic response and the formation of secondary structures are modulated by the interplay between base pairing and electrostatic interactions. Here we measure force-extension curves (FECs) of ssDNA molecules in optical tweezers set up over two orders of magnitude of monovalent and divalent salt conditions, and obtain its elastic parameters by fitting the FECs to semiflexible models of polymers. For both monovalent and divalent salts, we find that the electrostatic contribution to the persistence length is proportional to the Debye screening length, varying as the inverse of the square root of cation concentration. The intrinsic persistence length is equal to 0.7 nm for both types of salts, and the effectivity of divalent cations in screening electrostatic interactions appears to be 100-fold as compared with monovalent salt, in line with what has been recently reported for single-stranded RNA. Finally, we propose an analysis of the FECs using a model that accounts for the effective thickness of the filament at low salt condition and a simple phenomenological description that quantifies the formation of non-specific secondary structure at low forces.

Geographic patterns of genetic variation in a broadly distributed marine vertebrate: new insights into loggerhead turtle stock structure from expanded mitochondrial DNA sequences

Previous genetic studies have demonstrated that natal homing shapes the stock structure of marine turtle nesting populations. However, widespread sharing of common haplotypes based on short segments of the mitochondrial control region often limits resolution of the demographic connectivity of populations. Recent studies employing longer control region sequences to resolve haplotype sharing have focused on regional assessments of genetic structure and phylogeography. Here we synthesize available control region sequences for loggerhead turtles from the Mediterranean Sea, Atlantic, and western Indian Ocean basins. These data represent six of the nine globally significant regional management units (RMUs) for the species and include novel sequence data from Brazil, Cape Verde, South Africa and Oman. Genetic tests of differentiation among 42 rookeries represented by short sequences (380 bp haplotypes from 3,486 samples) and 40 rookeries represented by long sequences (~800 bp haplotypes from 3,434 samples) supported the distinction of the six RMUs analyzed as well as recognition of at least 18 demographically independent management units (MUs) with respect to female natal homing. A total of 59 haplotypes were resolved. These haplotypes belonged to two highly divergent global lineages, with haplogroup I represented primarily by CC-A1, CC-A4, and CC-A11 variants and haplogroup II represented by CC-A2 and derived variants. Geographic distribution patterns of haplogroup II haplotypes and the nested position of CC-A11.6 from Oman among the Atlantic haplotypes invoke recent colonization of the Indian Ocean from the Atlantic for both global lineages. The haplotypes we confirmed for western Indian Ocean RMUs allow reinterpretation of previous mixed stock analysis and further suggest that contemporary migratory connectivity between the Indian and Atlantic Oceans occurs on a broader scale than previously hypothesized. This study represents a valuable model for conducting comprehensive international cooperative data management and research in marine ecology.

Solid-phase synthesis and NMR studies of modified oligonucleotides forming triplex helix and of oligonucleopeptides mimicking the topoisomerase I-DNA covalent complex

Report for the scientific sojourn carried out at the Institut de Biologia Molecular de Barcelona of the CSIC –state agency – from april until september 2007. Topoisomerase I is an essential nuclear enzyme that modulates the topological status of DNA, facilitating DNA helix unwinding during replication and transcription. We have prepared the oligonucleotide-peptide conjugate Ac-NLeu-Asn-Tyr(p-3’TTCAGAAGC5’)-LeuC-CONH-(CH2)6-OH as model compound for NMR studies of the Topoisomerase I- DNA complex. Special attention was made on the synthetic aspects for the preparation of this challenging compound especially solid supports and protecting groups. The desired peptide was obtained although we did not achieve the amount of the conjugate needed for NMR studies. Most probably the low yield is due to the intrinsic sensitive to hydrolysis of the phosphate bond between oligonucleotide and tyrosine. We have started the synthesis and the structural characterization of oligonucleotides carrying intercalating compounds. At the present state we have obtained model duplex and quadruplex sequences modified with acridine and NMR studies are underway. In addition to this project we have successfully resolved the structure of a fusion peptide derived from hepatitis C virus envelope synthesized by the group of Dr. Haro and we have synthesized and started the characterization of a modified G-quadruplex.

An assessment of gene prediction accuracy in large DNA sequences

One of the first useful products from the human genome will be a set of predicted genes. Besides its intrinsic scientific interest, the accuracy and completeness of this data set is of considerable importance for human health and medicine. Though progress has been made on computational gene identification in terms of both methods and accuracy evaluation measures, most of the sequence sets in which the programs are tested are short genomic sequences, and there is concern that these accuracy measures may not extrapolate well to larger, more challenging data sets. Given the absence of experimentally verified large genomic data sets, we constructed a semiartificial test set comprising a number of short single-gene genomic sequences with randomly generated intergenic regions. This test set, which should still present an easier problem than real human genomic sequence, mimics the approximately 200kb long BACs being sequenced. In our experiments with these longer genomic sequences, the accuracy of GENSCAN, one of the most accurate ab initio gene prediction programs, dropped significantly, although its sensitivity remained high. Conversely, the accuracy of similarity-based programs, such as GENEWISE, PROCRUSTES, and BLASTX was not affected significantly by the presence of random intergenic sequence, but depended on the strength of the similarity to the protein homolog. As expected, the accuracy dropped if the models were built using more distant homologs, and we were able to quantitatively estimate this decline. However, the specificities of these techniques are still rather good even when the similarity is weak, which is a desirable characteristic for driving expensive follow-up experiments. Our experiments suggest that though gene prediction will improve with every new protein that is discovered and through improvements in the current set of tools, we still have a long way to go before we can decipher the precise exonic structure of every gene in the human genome using purely computational methodology.

Counter-ion effect on surfactant-DNA gel particles as controlled DNA delivery systems

The ability to entrap drugs within vehicles and subsequently release them has led to new treatments for a number of diseases. Based on an associative phase separation and interfacial diffusion approach, we developed a way to prepare DNA gel particles without adding any kind of cross-linker or organic solvent. Among the various agents studied, cationic surfactants offered particularly efficient control for encapsulation and DNA release from these DNA gel particles. The driving force for this strong association is the electrostatic interaction between the two components, as induced by the entropic increase due to the release of the respective counter-ions. However, little is known about the influence of the respective counter-ions on this surfactant-DNA interaction. Here we examined the effect of different counter-ions on the formation and properties of the DNA gel particles by mixing DNA (either single- (ssDNA) or double-stranded (dsDNA)) with the single chain surfactant dodecyltrimethylammonium (DTA). In particular, we used as counter-ions of this surfactant the hydrogen sulfate and trifluoromethane sulfonate anions and the two halides, chloride and bromide. Effects on the morphology of the particles obtained, the encapsulation of DNA and its release, as well as the haemocompatibility of these particles, are presented, using the counter-ion structure and the DNA conformation as controlling parameters. Analysis of the data indicates that the degree of counter-ion dissociation from the surfactant micelles and the polar/hydrophobic character of the counter-ion are important parameters in the final properties of the particles. The stronger interaction with amphiphiles for ssDNA than for dsDNA suggests the important role of hydrophobic interactions in DNA.

Current directions in DNA gel particles

A general understanding of interactions between DNA andoppositely charged compounds forms the basis for developing novelDNA-based materials, including gel particles. The association strength,which is altered by varying the chemical structure of the cationiccosolute, determines the spatial homogeneity of the gelation process,creating DNA reservoir devices and DNA matrix devices that can bedesigned to release either single- (ssDNA) or double-stranded(dsDNA) DNA. This paper reviews the preparation of DNA gelparticles using surfactants, proteins and polysaccharides. Particlemorphology, swelling/dissolution behaviour, degree of DNAentrapment and DNA release responses as a function of the nature ofthe cationic agent used are discussed. Current directions in thehaemocompatible and cytotoxic characterization of these DNA gelparticles have been also included.

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.

Allelic diversity and population structure in Vibrio cholerae O139 Bengal based on nucleotide sequence analysis

Comparative analysis of gene fragments of six housekeeping loci, distributed around the two chromosomes of Vibrio cholerae, has been carried out for a collection of 29 V. cholerae O139 Bengal strains isolated from India during the first epidemic period (1992 to 1993). A toxigenic O1 ElTor strain from the seventh pandemic and an environmental non-O1/non-O139 strain were also included in this study. All loci studied were polymorphic, with a small number of polymorphic sites in the sequenced fragments. The genetic diversity determined for our O139 population is concordant with a previous multilocus enzyme electrophoresis study in which we analyzed the same V. cholerae O139 strains. In both studies we have found a higher genetic diversity than reported previously in other molecular studies. The results of the present work showed that O139 strains clustered in several lineages of the dendrogram generated from the matrix of allelic mismatches between the different genotypes, a finding which does not support the hypothesis previously reported that the O139 serogroup is a unique clone. The statistical analysis performed in the V. cholerae O139 isolates suggested a clonal population structure. Moreover, the application of the Sawyer's test and split decomposition to detect intragenic recombination in the sequenced gene fragments did not indicate the existence of recombination in our O139 population.