DNA binding properties of oligodeoxynucleotides containing pyrrolidino C-nucleosides

We have incorporated pyrrolidino-C-nucleosides (pyrrolidino-pseudonucleosides) containing the base uracil and N-1-methyl uracil into oligodeoxynucleotides and compared their thermal duplex and triplex stabilities with unmodified or pseudouridine-containing oligodeoxynucleotides. We find relative destabilizations of triplex formation by ca. -13 to -1 degrees C per modification (relative to thymidine) in a strongly sequence dependent mode. Duplex formation is less destabilizing and more homogeneous with -4 to -2 degrees C per modification

Stability and kinetics of nucleic acid triplexes with chimaeric DNA/RNA third strands

A series of chimaeric DNA/RNA triplex-forming oligonucleotides (TFOs) with identical base-sequence but varying sequential composition of the sugar residues were prepared. The structural, kinetic and thermodynamic properties of triplex formation with their corresponding double-helical DNA target were investigated by spectroscopic methods. Kinetic and thermodynamic data were obtained from analysis of non-equilibrium UV-melting- and annealing curves in the range of pH 5.1 to 6.7 in a 10 mM citrate/phosphate buffer containing 0.1M NaCl and 1 mM EDTA. It was found that already single substitutions of ribo- for deoxyribonucleotides in the TFOs greatly affect stability and kinetics of triplex formation in a strongly sequence dependent manner. Within the sequence context investigated, triplex stability was found to increase when deoxyribonucleotides were present at the 5'-side and ribonucleotides in the center of the TFO. Especially the substitution of thymidines for uridines in the TFO was found to accelerate both, the association and dissociation process, in a strongly position-dependent way. Differential structural information on triplexes and TFO single-strands was obtained from CD-spectroscopy and gel mobility experiments. Only minor changes were observed in the CD spectra of the triplexes at all pH values investigated, and the electrophoretic mobility was nearly identical in all cases, indicating a high degree of structural similarity. In contrast, the single-stranded TFOs showed high structural variability as determined in the same way. The results are discussed in the context of the design of TFOs for therapeutic or biochemical applications.

Conformational diversity versus nucleic acid triplex stability, a combinatorial study

The stability of a triple helix formed between a DNA duplex and an incoming oligonucleotide strand strongly depends on the solvent conditions and on intrinsic chemical and conformational factors. Attempts to increase triple helix stability in the past included chemical modification of the backbone, sugar ring, and bases in the third strand. However, the predictive power of such modifications is still rather poor. We therefore developed a method that allows for rapid screening of conformationally diverse third strand oligonucleotides for triplex stability in the parallel pairing motif to a given DNA double helix sequence. Combinatorial libraries of oligonucleotides of the requisite (fixed) base composition and length that vary in their sugar unit (ribose or deoxyribose) at each position were generated. After affinity chromatography against their corresponding immobilized DNA target duplex, utilizing a temperature gradient as the selection criterion, the oligonucleotides forming the most stable triple helices were selected and characterized by physicochemical methods. Thus, a series of oligonucleotides were identified that allowed us to define basic rules for triple helix stability in this conformationally diverse system. It was found that ribocytidines in the third strand increase triplex stability relative to deoxyribocytidines independently of the neighboring bases and position along the strand. However, remarkable sequence-dependent differences in stability were found for (deoxy)thymidines and uridines

AsaGEI2b: a new variant of a genomic island identified in the Aeromonas salmonicida subsp. salmonicida JF3224 strain isolated from a wild fish in Switzerland.

Aeromonas salmonicida subsp. salmonicida is the causal agent of furunculosis in salmonids. We recently identified a group of genomic islands (AsaGEI) in this bacterium. AsaGEI2a, one of these genomic islands, has almost exclusively been identified in isolates from North America. To date, Aeromonas salmonicida subsp. salmonicida JF3224, a strain isolated from a wild brown trout (Salmo trutta) caught in Switzerland, was the only European isolate that appeared to bear AsaGEI2a. We analyzed the genome of JF3224 and showed that the genomic island in JF3224 is a new variant of AsaGEI, which we have called AsaGEI2b. While AsaGEI2b shares the same integrase gene and insertion site as AsaGEI2a, it is very different in terms of many other features. Additional genomic investigations combined with PCR genotyping revealed that JF3224 is sensitive to growth at 25°C, leading to insertion sequence-dependent rearrangement of the locus on the pAsa5 plasmid that encodes a type three secretion system, which is essential for the virulence of the bacterium. The analysis of the JF3224 genome confirmed that AsaGEIs are accurate indicators of the geographic origins of A. salmonicida subsp. salmonicida isolates and is another example of the susceptibility of the pAsa5 plasmid to DNA rearrangements.

Estimates of in situ Larval Development Time for the Lobster, Homarus Americanus

Larval development time is a critical factor in assessing the potential for larval transport, mortality. and subsequently, the connectivity of marine populations through larval exchange. Most estimates of larval duration are based on laboratory studies and may not reflect development times in nature. For larvae of the American lobster (Homarus americanus), temperature-dependent development times have been established in previous laboratory studies. Here, we used the timing of seasonal abundance curves for newly hatched larvae (stage 1) and the final plankonic instar (postlarva), coupled with a model of temperature-dependent development to assess development time in the field. We were unable to reproduce the timing of the seasonal abundance curves using laboratory development rates in our model. Our results suggest that larval development in situ may be twice as fast as reported laboratory rates. This will result in reduced estimates of larval transport potential, and increased estimates of instantaneous mortality rate and production.

Mg/Ca and stable oxygen measurement of the past 27 Ma of ODP Hole 120-747A

We explore the applicability of paired Mg/Ca and 18O/16O measurements on benthic foraminifera from Southern Ocean site 747 to paleoceanographic reconstructions on pre-Pleistocene timescales. We focus on the late Oligocene through Pleistocene (27-0 Ma) history of paleotemperatures and the evolution of the d18O values of seawater (d18Osw) at a temporal resolution of ~100-200 kyr. Absolute paleotemperature estimates depend on assumptions of how Mg/Ca ratios of seawater have changed over the past 27 Myr, but relative changes that occur on geologically brief timescales are robust. Results indicate that at the Oligocene to Miocene boundary (23.8 Ma), temperatures lag the increase in global ice-volume deduced from benthic foraminiferal d18O values, but the smaller-scale Miocene glaciations are accompanied by ocean cooling of -1°C. During the mid-Miocene phase of Antarctic ice sheet growth (~15-13 Ma), water temperatures cool by ~3°C. Unlike the benthic foraminiferal d18O values, which remain relatively constant thereafter, temperatures vary (by 3°C) and reach maxima at ~12 and ~8.5 Ma. The onset of significant Northern Hemisphere glaciation during the late Pliocene is synchronous with an ~4°C cooling at site 747. A comparison of our d18Osw curve to the Haq et al. (1987, doi:10.1126/science.235.4793.1156 ) sea level curve yields excellent agreement between sequence boundaries and times of increasing seawater 18O/16O ratios. At ~12-11 Ma in particular, when benthic foraminiferal d18O values do not support a further increase in ice volume, the d18Osw curve comes to a maximum that corresponds to a major mid-Miocene sea level regression. The agreement between the character of our Mg/Ca-based d18Osw curve and sequence stratigraphy demonstrates that benthic foramaniferal Mg/Ca ratios can be used to trace the d18Osw on pre-Pleistocene timescales despite a number of uncertainties related to poorly constrained temperature calibrations and paleoseawater Mg/Ca ratios. The Mg/Ca record also highlights that deep ocean temperatures can vary independently and unexpectedly from ice volume changes, which can lead to misinterpretations of the d18O record.

Architectural specificity in chromatin structure at the TATA box in vivo: Nucleosome displacement upon β-phaseolin gene activation

Extensive studies of the β-phaseolin (phas) gene in transgenic tobacco have shown that it is highly active during seed embryogenesis but is completely silent in leaf and other vegetative tissues. In vivo footprinting revealed that the lack of even basal transcriptional activity in vegetative tissues is associated with the presence of a nucleosome that is rotationally positioned with base pair precision over three phased TATA boxes present in the phas promoter. Positioning is sequence-dependent because an identical rotational setting is obtained upon nucleosome reconstitution in vitro. A comparison of DNase I and dimethyl sulfate footprints in vivo and in vitro strongly suggests that this repressive chromatin architecture is remodeled concomitant with gene activation in the developing seed. This leads to the disruption of histone-mediated DNA wrapping and the assembly of the TATA boxes into a transcriptionally competent nucleoprotein complex.

Folding protein models with a simple hydrophobic energy function: The fundamental importance of monomer inside/outside segregation

The present study explores a “hydrophobic” energy function for folding simulations of the protein lattice model. The contribution of each monomer to conformational energy is the product of its “hydrophobicity” and the number of contacts it makes, i.e., E(h⃗, c⃗) = −Σi=1N cihi = −(h⃗.c⃗) is the negative scalar product between two vectors in N-dimensional cartesian space: h⃗ = (h1, … , hN), which represents monomer hydrophobicities and is sequence-dependent; and c⃗ = (c1, … , cN), which represents the number of contacts made by each monomer and is conformation-dependent. A simple theoretical analysis shows that restrictions are imposed concomitantly on both sequences and native structures if the stability criterion for protein-like behavior is to be satisfied. Given a conformation with vector c⃗, the best sequence is a vector h⃗ on the direction upon which the projection of c⃗ − c̄⃗ is maximal, where c̄⃗ is the diagonal vector with components equal to c̄, the average number of contacts per monomer in the unfolded state. Best native conformations are suggested to be not maximally compact, as assumed in many studies, but the ones with largest variance of contacts among its monomers, i.e., with monomers tending to occupy completely buried or completely exposed positions. This inside/outside segregation is reflected on an apolar/polar distribution on the corresponding sequence. Monte Carlo simulations in two dimensions corroborate this general scheme. Sequences targeted to conformations with large contact variances folded cooperatively with thermodynamics of a two-state transition. Sequences targeted to maximally compact conformations, which have lower contact variance, were either found to have degenerate ground state or to fold with much lower cooperativity.

Antitumor activity and pharmacokinetics of a mixed-backbone antisense oligonucleotide targeted to the RIα subunit of protein kinase A after oral administration

Overexpression of the RIα subunit of cAMP-dependent protein kinase (PKA) has been demonstrated in various human cancers. PKA has been suggested as a potential target for cancer therapy. The goal of the present study was to evaluate an anti-PKA antisense oligonucleotide (mixed-backbone oligonucleotide) as a therapeutic approach to human cancer treatment. The identified oligonucleotide inhibited the growth of cell lines of human colon cancer (LS174T, DLD-1), leukemia (HL-60), breast cancer (MCF-7, MDA-MB-468), and lung cancer (A549) in a time-, concentration-, and sequence-dependent manner. In a dose-dependent manner, the oligonucleotide displayed in vivo antitumor activity in severe combined immunodeficient and nude mice bearing xenografts of human cancers of the colon (LS174T), breast (MDA-MB-468), and lung (A549). The routes of drug administration were intraperitoneal and oral. Synergistic effects were found when the antisense oligonucleotide was used in combination with the cancer chemotherapeutic agent cisplatin. The pharmacokinetics of the oligonucleotide after oral administration of 35S-labeled oligonucleotide into tumor-bearing mice indicated an accumulation and retention of the oligonucleotide in tumor tissue. This study further provides a basis for clinical studies of the antisense oligonucleotide targeted to the RIα subunit of PKA (GEM 231) as a cancer therapeutic agent used alone or in combination with conventional chemotherapy.

Staphylococcal protein A: Unfolding pathways, unfolded states, and differences between the B and E domains

We present multiple native and denaturation simulations of the B and E domains of the three-helix bundle protein A, totaling 60 ns. The C-terminal helix (H3) consistently denatures later than either of the other two helices and contains residual helical structure in the denatured state. These results are consistent with experiments suggesting that the isolated H3 fragment is more stable than H1 and H2 and that H3 forms early in folding. Interestingly, the denatured state of the B domain is much more compact than that of the E domain. This sequence-dependent effect on the dimensions of the denatured state and the lack of correlation with structure suggest that the radius of gyration can be a misleading reaction coordinate for unfolding/folding. Various unfolding and refolding events are observed in the denaturation simulations. In some cases, the transitions are facilitated through interactions with other portions of the protein—contact-assisted helix formation. In the native simulations, the E domain is very stable: after 6 ns, the Cα root-mean-square deviation from the starting structure is less than 1.4 Å. In contrast, the native state of the B domain deviates more and its inter-helical angles fluctuate. In apparent contrast, we note that the B domain is thermodynamically more stable than the E domain. The simulations suggest that the increased stability of the B domain may be due to heightened mobility, and therefore entropy, in the native state and decreased mobility/entropy in the more compact denatured state.

p53-induced DNA bending and twisting: p53 tetramer binds on the outer side of a DNA loop and increases DNA twisting

DNA binding activity of p53 is crucial for its tumor suppressor function. Our recent studies have shown that four molecules of the DNA binding domain of human p53 (p53DBD) bind the response elements with high cooperativity and bend the DNA. By using A-tract phasing experiments, we find significant differences between the bending and twisting of DNA by p53DBD and by full-length human wild-type (wt) p53. Our data show that four subunits of p53DBD bend the DNA by 32–36°, whereas wt p53 bends it by 51–57°. The directionality of bending is consistent with major groove bends at the two pentamer junctions in the consensus DNA response element. More sophisticated phasing analyses also demonstrate that p53DBD and wt p53 overtwist the DNA response element by ≈35° and ≈70°, respectively. These results are in accord with molecular modeling studies of the tetrameric complex. Within the constraints imposed by the protein subunits, the DNA can assume a range of conformations resulting from correlated changes in bend and twist angles such that the p53–DNA tetrameric complex is stabilized by DNA overtwisting and bending toward the major groove at the CATG tetramers. This bending is consistent with the inherent sequence-dependent anisotropy of the duplex. Overall, the four p53 moieties are placed laterally in a staggered array on the external side of the DNA loop and have numerous interprotein interactions that increase the stability and cooperativity of binding. The novel architecture of the p53 tetrameric complex has important functional implications including possible p53 interactions with chromatin.

Dynamical principles in biological processes: A model of charge migration in proteins and DNA

The generalized master equations (GMEs) that contain multiple time scales have been derived quantum mechanically. The GME method has then been applied to a model of charge migration in proteins that invokes the hole hopping between local amino acid sites driven by the torsional motions of the floppy backbones. This model is then applied to analyze the experimental results for sequence-dependent long-range hole transport in DNA reported by Meggers et al. [Meggers, E., Michel-Beyerle, M. E., & Giese, B. (1998) J. Am. Chem. Soc. 120, 12950–12955]. The model has also been applied to analyze the experimental results of femtosecond dynamics of DNA-mediated electron transfer reported by Zewail and co-workers [Wan, C., Fiebig, T., Kelley, S. O., Treadway, C. R., Barton, J. K. & Zewail, A. H. (1999) Proc. Natl. Acad. Sci. USA 96, 6014–6019]. The initial events in the dynamics of protein folding have begun to attract attention. The GME obtained in this paper will be applicable to this problem.

Solution structure and dynamics of GCN4 cognate DNA: NMR investigations

A 12 bp long GCN4-binding, self-complementary duplex DNA d(CATGACGTCATG)2 has been investigated by NMR spectroscopy to study the structure and dynamics of the molecule in aqueous solution. The NMR structure of the DNA obtained using simulated annealing and iterative relaxation matrix calculations compares quite closely with the X-ray structure of ATF/CREB DNA in complex with GCN4 protein (DNA-binding domain). The DNA is also seen to be curved in the free state and this has a significant bearing on recognition by the protein. The dynamic characteristics of the molecule have been studied by 13C relaxation measurements at natural abundance. A correlation has been observed between sequence-dependent dynamics and recognition by GCN4 protein.

Mapping the intrinsic curvature and flexibility along the DNA chain

The energy of DNA deformation plays a crucial and active role in its packaging and its function in the cell. Considerable effort has gone into developing methodologies capable of evaluating the local sequence-directed curvature and flexibility of a DNA chain. These studies thus far have focused on DNA constructs expressly tailored either with anomalous flexibility or curvature tracts. Here we demonstrate that these two structural properties can be mapped also along the chain of a “natural” DNA with any sequence on the basis of its scanning force microscope (SFM) images. To know the orientation of the sequence of the investigated DNA molecules in their SFM images, we prepared a palindromic dimer of the long DNA molecule under study. The palindromic symmetry also acted as an internal gauge of the statistical significance of the analysis carried out on the SFM images of the dimer molecules. It was found that although the curvature modulus is not efficient in separating static and dynamic contributions to the curvature of the population of molecules, the curvature taken with its direction (its sign in two dimensions) permits the direct separation of the intrinsic curvature from the flexibility contributions. The sequence-dependent flexibility seems to vary monotonically with the chain's intrinsic curvature; the chain rigidity was found to modulate as its local thermodynamic stability and does not correlate with the dinucleotide chain rigidities evaluation made from x-ray data by other authors.

Identification of the yeast cytidine deaminase CDD1 as an orphan C→U RNA editase

Yeast co-expressing rat APOBEC-1 and a fragment of human apolipoprotein B (apoB) mRNA assembled functional editosomes and deaminated C6666 to U in a mooring sequence-dependent fashion. The occurrence of APOBEC-1-complementing proteins suggested a naturally occurring mRNA editing mechanism in yeast. Previously, a hidden Markov model identified seven yeast genes encoding proteins possessing putative zinc-dependent deaminase motifs. Here, only CDD1, a cytidine deaminase, is shown to have the capacity to carry out C→U editing on a reporter mRNA. This is only the second report of a cytidine deaminase that can use mRNA as a substrate. CDD1-dependent editing was growth phase regulated and demonstrated mooring sequence-dependent editing activity. Candidate yeast mRNA substrates were identified based on their homology with the mooring sequence-containing tripartite motif at the editing site of apoB mRNA and their ability to be edited by ectopically expressed APOBEC-1. Naturally occurring yeast mRNAs edited to a significant extent by CDD1 were, however, not detected. We propose that CDD1 be designated an orphan C→U editase until its native RNA substrate, if any, can be identified and that it be added to the CDAR (cytidine deaminase acting on RNA) family of editing enzymes.