Sequence preference for BI/BII conformations in DNA: MD and crystal structure data analysis

Deciphering sequence information from sugar-phosphate backbone is finely tuned through the conformational substates of DNA. BII conformation, one of the conformational substates of B-DNA, is known to play a key role in DNA-protein recognition. BI and BII are identified by the epsilon-zeta difference, which is negative in BI and positive in BII. Our analysis of MD and crystal structures shows that BII conformation is sequence specific and dinucleotides GC, CG, CA, TG, TA show high preference to take up BII conformation, while TT, TC, CT, CC dinucleotides rarely take up this conformation. Significant changes were observed in the dinucleotide parameters viz. twist, roll, and slide for the steps having BII conformation. Interestingly, the magnitude of variation in the dinucleotide parameters is seen to depend mainly on two factors, the magnitude of epsilon-zeta difference and the presence or absence of BII conformation in the second strand, across the WC base-paired dinucleotide step. Based on these two factors, the conformational substate of a dinucleotide step can be further classified as BI.BI (BI conformation in both strands), BI.BII (BI conformation in one strand and BII conformation in the other), and BII.BII (BII conformation in both strands). The occurrence of BII in both strands was found to be quite rare and thus, it can be concluded that BI.BI and BI.BII hybrid steps are more favorable than a BII.BII step. In conformity with the sequence preference seen for dinucleotides in each strand, BII.BII combination of backbone conformation was observed only for GC, CG, CA, and TG containing dinucleotide steps. We further classified BII.BII step as strong BII and weak BII depending on the magnitude of the average epsilon-zeta difference. The dinucleotide steps which belong to the category of strong BII, have large twist, high positive slide and negative roll values, while those in the weak BII group have roll, twist, and slide values similar to that of hybrid BI.BII steps. This conformational property could be contributing to the groove opening/closing and thus can modulate protein-DNA interaction.     

The covariation between TpA deficiency, CpG deficiency, and G+C content of human isochores is due to a mathematical artifact

CpG and TpA dinucleotides are underrepresented in the human genome. The CpG deficiency is due to the high mutation rate from C to T in methylated CpG's. The TpA suppression was thought to reflect a counterselection against TpA's destabilizing effect in RNA. Unexpectedly, the TpA and CpG deficiencies vary according to the G+C contents of sequences. It has been proposed that the variation in CpG suppression was correlated with a particular chromatin organization in G+C-rich isochores. Here, we present an improved model of dinucleotide evolution accounting for the overlap between successive dinucleotides. We show that an increased mutation rate from CpG to TpG or CpA induces both an apparent TpA deficiency and a correlation between CpG and TpA deficiencies and G+C content. Moreover, this model shows that the ratio of observed over expected CpG frequency underestimates the real CpG deficiency in G+C-rich sequences. The predictions of our model fit well with observed frequencies in human genomic data. This study suggests that previously published selectionist interpretations of patterns of dinucleotide frequencies should be taken with caution. Moreover, we propose new criteria to identify unmethylated CpG islands taking into account this bias in the measure of CpG depletion.     

DNA fine structure and dynamics in crystals and in solution: the impact of BI/BII backbone conformations

Sugar-phosphate backbone conformations are an important structural element for a complete understanding of specific recognition in nucleic acid-protein interactions. They can be involved both in early stages of target discrimination and in structural adaptation upon binding. In the first part of this study, we have analyzed high-resolution structures of double-stranded B-DNA either isolated or bound to proteins, and explored the impact of both the standard BI and the unusual BII phosphate backbone conformations on neighboring sugar puckers and on selected helical parameters. Correlations are found to be similar for free and bound DNA, and in both categories, the possible facing backbone conformations (BI.BI, BI.BII, and BII.BII) define well-characterized substates in the B-DNA conformational space. Notably, BII.BII steps are characterized by specific, and sequence-independent, structural effects involving reduced standard deviations for almost all conformational parameters. In the second part of this work, we analyze four 10 ns molecular dynamics simulations in explicit solvent on the DNA targets of NF-kappaB and bovine papillomavirus E2 proteins, highlighting the multiplicity of backbone dynamical behavior. These results show sequence effects on the percentages of BI and BII conformers, the preferential state of facing backbones, the occurrence of coupled transitions. The backbone states can consequently be seen as a mechanism for transmitting information from the bases to the phosphate groups and thus for modulating the overall structural properties of the target DNA.     

Genetic organization of the polymorphic equine alpha globin locus and sequence of the BII alpha 1 gene.

The equine alpha globin gene complex comprises two functional alpha genes and an alpha-like pseudogene arranged in the order 5'-alpha 2-(5kb)-alpha 1-(3kb)-psi alpha-3'. A single (embryonic) zeta-like sequence lies within a 12 kb region 5' to the alpha 2 gene. We have determined the sequence of the alpha 1 gene of the BII haplotype, one of two most common haplotypes (the other being BI) which encode alpha globins with either Tyr (BI) or Phe (BII) at codon 24 in both linked alpha genes. In BI and BII the non-allelic alpha 2 and alpha 1 genes respectively code for Gln or Lys at codon 60, thus accounting for the 4 alpha globin types seen in BI/BII heterozygotes. Genomic restriction enzyme maps of the BII alpha complex (24Phe/60Lys,Gln) and the allelic BI (24Tyr/60Lys,Gln) are identical to each other, and to those of a rarer normal haplotype, A, which encodes only alpha 24Tyr/60Gln globin, and a low expression mutant of BII which encodes only 24Phe/60Lys globin. These two latter haplotypes must therefore have a linked pair of alpha genes, as in BI and BII, but with identical coding properties, and it is suggested that this has arisen by gene conversion.     

Water and ion binding around r(UpA)12 and d(TpA)12 oligomers--comparison with RNA and DNA (CpG)12 duplexes

The structural and dynamic properties of the water and ion first coordination shell of the r(A-U) and d(A-T) base-pairs embedded within the r(UpA)12 and d(TpA)12 duplexes are described on the basis of two 2.4 ns molecular dynamics simulations performed in a neutralizing aqueous environment with 0.25 M added KCl. The results are compared to previous molecular dynamics simulations of the r(CpG)12 and d(CpG)12 structures performed under similar conditions. It can be concluded that: (i) RNA helices are more rigid than DNA helices of identical sequence, as reflected by the fact that RNA duplexes keep their initial A-form shape while DNA duplexes adopt more sequence-specific shapes. (ii) Around these base-pairs, the water molecules occupy 21 to 22 well-defined hydration sites, some of which are partially occupied by potassium ions. (iii) These hydration sites are occupied by an average of 21.9, 21.0, 20.1, and 19.8 solvent molecules (water and ions) around the r(G=C), r(A-U), d(G=C), and d(A-T) pairs, respectively. (iv) From a dynamic point of view, the stability of the hydration shell is the strongest for the r(G=C) pairs and the weakest for the d(A-T) pairs. (v) For RNA, the observed long-lived hydration patterns are essentially non-sequence dependent and involve water bridges located in the deep groove and linking OR atoms of adjacent phosphate groups. Maximum lifetimes are close to 400 ps. (vi) In contrast, for DNA, long-lived hydration patterns are sequence dependent and located in the minor groove. For d(CpG)12, water bridges linking the (G)N3 and (C)O2 with the O4' atoms of adjacent nucleotides with 400 ps maximum lifetimes are characterized while no such bridges are observed for d(TpA)12. (vii) Potassium ions are observed to bind preferentially to deep/major groove atoms at RpY steps, essentially d(GpC), r(GpC), and r(ApU), by forming ion-bridges between electronegative atoms of adjacent base-pairs. On average, about half an ion is observed per base-pair. Positive ion-binding determinants are related to the proximity of two or more electronegative atoms. Negative binding determinants are associated with the electrostatic and steric hindrance due to the proximity of electropositive amino groups and neutral methyl groups. Potassium ions form only transient contacts with phosphate groups.     

DNA structure influences sequence specific cleavage by bleomycin.

We have examined the cleavage of several synthetic DNA sequences by iron(II)-bleomycin. We find that, although bleomycin cuts mixed sequence DNAs with a preference for GC = GT > GA >> GG, it efficiently cleaves regions of (AT)n cutting exclusively at ApT, not TpA. Isolated ApT steps show very little cleavage while blocks of three or more contiguous ATs are cut as efficiently as GpT. This cleavage is specific for (AT)n, since sequences of the type (TAA)n.(TTA)n and (ATT)n.(AAT)n are hardly cut at all. No cleavage is observed at ApC or CpA within sequences of the type (AC)n.(GT)n; regions of An.Tn are also not cut. Although the cobalt-bleomycin complex (which binds to but does not cleave DNA) yields good DNase I footprints at GT and GC sites, no footprints are observed within (AT)n, suggesting that although the cleavage reaction is efficient, the binding affinity is relatively weak. We propose a model in which bleomycin cleavage is determined by local DNA structure, while strong binding requires the presence of a guanine residue.     

DNA-binding proteins in cells and membrane blebs of Neisseria gonorrhoeae.

Naturally elaborated membrane bleb fractions BI and BII of Neisseria gonorrhoeae contain both linear and circular DNAs. Because little is known about the interactions between DNA and blebs, studies were initiated to identify specific proteins that bind DNA in elaborated membrane blebs. Western immunoblots of whole-cell and bleb proteins from transformation-competent and DNA-uptake-deficient (dud) mutants were probed with single- or double-stranded gonococcal DNA, pBR322, or synthetic DNA oligomers containing intact or altered gonococcal transformation uptake sequences. The specificity and sensitivity of a nonradioactive DNA-binding protein assay was evaluated, and the assay was used to visualize DNA-protein complexes on the blots. The complexes were then characterized by molecular mass, DNA-binding specificity, and expression in bleb fractions. The assay effectively detected blotted DNA-binding proteins. At least 17 gonococcal DNA-binding proteins were identified; unique subsets occurred in BI and BII. Certain DNA-binding proteins had varied affinities for single- and double-stranded DNA, and the intact transformation uptake sequence competitively displaced the altered sequence from a BI protein at 11 kilodaltons (kDa). A dud mutant, strain FA660, lacked DNA-binding activity at the 11-kDa protein in BI. The segregation of DNA-binding proteins within BI and BII correlates with their distinct protein profiles and suggests that these vesicles may play different roles. Although the DNA-binding proteins expressed in BII may influence the nuclease-resistant export of plasmids within BII vesicles, the BI 11-kDa protein may bind transforming DNA.     

Invariant and variable base stacking geometries in B-DNA and A-DNA

We calculated the interatomic distances between all couples of non-hydrogen atoms belonging to the neighboring Watson-Crick base pairs in the available crystal structures of DNA. Their standard deviations revealed remarkably large differences in the variability of the base stacking geometries of the particular steps. In line with experimental studies in solution, (CpA)-(TpG) and (TpA).(TpA) were identified as the most variable or flexible steps in the crystal structures of B-DNA. On the other hand, base stacking geometries of the (ApT).(ApT) steps were the most invariant, which was very surprising because all three steps composed only of C and G were much more flexible. This finding suggests that conformational stability of DNA and the rigidity have different origins. Furthermore, the nucleotide sequence dependence of the flexibility was almost reversed in A-DNA because the most flexible steps in B-DNA were the least flexible in A-DNA. The most invariant steps of B-DNA were variable in A-DNA. The (ApT).(ApT) step was a notable exception to this rule because it belonged to the most rigid steps in both B-DNA and A-DNA. The present results are fully consistent with the properties that poly(dA-dT).poly(dA-dT), poly(dA).poly(dT), poly(dAdC).poly(dG-dT) and poly(dA-dG).poly(dC-dT) exhibit in solution.     

The effects of sequence context on base dynamics at TpA steps in DNA studied by NMR.

Base dynamics, heretofore observed only at TpA steps in DNA, were investigated as a function of sequence context by NMR spectroscopy. The large amplitude conformational dynamics have been previously observed in TnAn segments where n > or = 2. In order to determine whether the dynamic characteristics occur in more general sequence contexts, we examined four self-complementary DNA sequences, [d(CTTTA-NATNTAAAG)2] (where N = A, C, T, G and N = complement of N). The anomalous broadening of the TpA adenine H2 resonance which is indicative of large amplitude base motion was observed in all nine unique four nucleotide contexts. Furthermore, all the adenine H2 resonances experienced a linewidth maximum as a function of temperature, which is a characteristic of the dynamic process. Interestingly, the temperature of the linewidth maximum varied with sequence indicating that the thermodynamics of TpA base dynamics are also sequence dependent. In one example, neither a T preceding nor an A trailing the TpA step was required for base dynamics. These results show that base dynamics, heretofore observed in only a few isolated sequences, occurs at all TpA steps which are either preceded or followed by a thymine or adenine, respectively, and may be characteristic of all TpA steps in DNA notwithstanding sequence context.     

The incidence and distribution of CpG----TpG transitions in the coagulation factor IX gene. A fresh look at CpG mutational hotspots.

The mutations of 76 haemophilia B patients representing the whole population registered with the Malm? haemophilia centre (42) and referrals from the UK, were characterised. RFLP haplotype analysis of the defective genes indicated that 51 single base pair substitutions were definitely of independent origin and 27 of these were CpG----TpG or CpA transitions. This represents a 38-fold excess over other single-base changes. Most of such transitions (82%) occur at 9 CpG sites occupying critical positions (transitions at 3 sites substitute essential arginines, while at 6 sites transition to TpG creates stop codons). Sixteen of the 18 possible transitions at these 9 sites cause clear haemophilia B and should be fully ascertained in our haemophilia B population. This allowed the direct estimate of the rate of CpG transitions. This is 1.05 x 10(-7) substitutions per base per gamete per generation. The marked excess of CpG transitions in haemophilia B appears partly due to the high proportion of CpG sites at critical positions (at least 9 out of 20). We propose that this follows from the fact that male hemizygosity makes X-linked genes particularly susceptible to selective forces that tend to fix CpG sites arising at critical positions.     

Viruses infecting a warm water picoeukaryote shed light on spatial co-occurrence dynamics of marine viruses and their hosts

The marine picoeukaryote Bathycoccus prasinos has been considered a cosmopolitan alga, although recent studies indicate two ecotypes exist, Clade BI (B. prasinos) and Clade BII. Viruses that infect Bathycoccus Clade BI are known (BpVs), but not that infect BII. We isolated three dsDNA prasinoviruses from the Sargasso Sea against Clade BII isolate RCC716. The BII-Vs do not infect BI, and two (BII-V2 and BII-V3) have larger genomes (~210 kb) than BI-Viruses and BII-V1. BII-Vs share ~90% of their proteins, and between 65% to 83% of their proteins with sequenced BpVs. Phylogenomic reconstructions and PolB analyses establish close-relatedness of BII-V2 and BII-V3, yet BII-V2 has 10-fold higher infectivity and induces greater mortality on host isolate RCC716. BII-V1 is more distant, has a shorter latent period, and infects both available BII isolates, RCC716 and RCC715, while BII-V2 and BII-V3 do not exhibit productive infection of the latter in our experiments. Global metagenome analyses show Clade BI and BII algal relative abundances correlate positively with their respective viruses. The distributions delineate BI/BpVs as occupying lower temperature mesotrophic and coastal systems, whereas BII/BII-Vs occupy warmer temperature, higher salinity ecosystems. Accordingly, with molecular diagnostic support, we name Clade BII Bathycoccus calidus sp. nov. and propose that molecular diversity within this new species likely connects to the differentiated host-virus dynamics observed in our time course experiments. Overall, the tightly linked biogeography of Bathycoccus host and virus clades observed herein supports species-level host specificity, with strain-level variations in infection parameters.Subject terms: Microbial biooceanography, Phylogenetics     

The loss of dinucleotides CpG from DNA. IV. Methylation and divergence of genes and pseudogenes of small nuclear RNA

The frequency of neighboring base pairs in nucleotide sequences of over 80 genes and pseudogenes of low molecular weight RNAs U1-U8, 4.5S and 7S in different eukaryotes was determined. The probable frequency of CpG----TpG + CpA substitutions, caused as a result of the deamination of 5-methylcytosine residues in DNA, was determined. It was found that the genes of small RNAs do not reveal a single level of CpG methylation for all the species studied. In most cases CpG in the genes of U1, 4.5S and 7S RNA are methylated, whereas in the genes of U2-U6-RNA these sites must have never been subjected to methylation. Nearly all the investigated pseudogenes of different small RNAs are strongly methylated due to a considerable lack of CpG. It was established that CpG----TpG + CpA transitions may amount to as much third of all the mutations accumulated in the genes of the same RNAs in different species. Such transitions in pseudogenes may account for 40% of all the nucleotide substitutions. This disproportionately high level of mutations in CpG dinucleotides (3-5-fold higher than in other DNA dupletes) must be the direct result of the methylation of these sites. Consequently, CpG methylation causes a dramatic acceleration of the divergence rate of DNA sequences. It has been concluded that protection of most vital genes against methylation is one of the essential conditions for sustaining the high level of stability of the macromolecular structure and for the reliability of macromolecular functioning in a cell.     

CpG Dinucleotide Frequencies Reveal the Role of Host Methylation Capabilities in Parvovirus Evolution

Parvoviruses are rapidly evolving viruses that infect a wide range of hosts, including vertebrates and invertebrates. Extensive methylation of the parvovirus genome has been recently demonstrated. A global pattern of methylation of CpG dinucleotides is seen in vertebrate genomes, compared to “fractional” methylation patterns in invertebrate genomes. It remains unknown if the loss of CpG dinucleotides occurs in all viruses of a given DNA virus family that infect host species spanning across vertebrates and invertebrates. We investigated the link between the extent of CpG dinucleotide depletion among autonomous parvoviruses and the evolutionary lineage of the infected host. We demonstrate major differences in the relative abundance of CpG dinucleotides among autonomous parvoviruses which share similar genome organization and common ancestry, depending on the infected host species. Parvoviruses infecting vertebrate hosts had significantly lower relative abundance of CpG dinucleotides than parvoviruses infecting invertebrate hosts. The strong correlation of CpG dinucleotide depletion with the gain in TpG/CpA dinucleotides and the loss of TpA dinucleotides among parvoviruses suggests a major role for CpG methylation in the evolution of parvoviruses. Our data present evidence that links the relative abundance of CpG dinucleotides in parvoviruses to the methylation capabilities of the infected host. In sum, our findings support a novel perspective of host-driven evolution among autonomous parvoviruses.     

Structure and Dynamics of Double Helical DNA in Torsion Angle Hyperspace: A Molecular Mechanics Approach

Abstract

Analysis of the conformational space populated by the torsion angles and the correlation between the conformational energy and the sequence of DNA are important for fully understanding DNA structure and function. Presence of seven variable torsion angles about single covalent bonds in DNA main chain puts a big challenge for such analysis. We have carried out restrained energy minimization studies for four representative dinucleosides, namely d(ApA):d(TpT), d(CpG):d(CpG), d(GpC):d(GpC) and d(CpA):d(TpG) to determine the energy hyperspace of DNA in context to the values of the torsion angles and the structural properties of the DNA conformations populating the favorable regions of this energy hyperspace. The torsion angles were manipulated by constraining their values at the reference points and then performing energy minimization. The energy minima obtained on the potential energy contour plots mostly correspond to the conformations populated in crystal structures of DNA. Some novel favorable conformations that are not present in crystal structure data are also found. The plots also suggest few low energy routes for conformational transitions or the associated energy barrier heights. Analyses of base pairing and stacking possibility reveal structural changes accompanying these transitions as well as the flexibility of different base steps towards variations in different torsion angles.     

BI-BII transitions in B-DNA.

Molecular modelling is used to study the conformational and energetic aspects of BI-BII transitions within the backbone of a B-DNA dodecamer d(CATGACGTCATG) whose fine structure has previously been determined by molecular modelling combined with NMR spectroscopy. It is shown that while the dodecamer under investigation does not contain any BII junctions, the central CpG step can most easily undergo the transition. More generally, it is also found that the base sequence and hence the backbone geometry of a DNA segment, strongly influences both the conformational impact of the transition, the associated energy barrier and the stability of the resulting BII state.     

M.TaqI facilitates the base flipping via an unusual DNA backbone conformation

MD simulations have been carried out to understand the dynamical behavior of the DNA substrate of the Thermus aquaticus DNA methyltransferase (M.TaqI) in the methylation process at N6 of adenine. As starting structures, an x-ray structure of M.TaqI in complex with DNA and cofactor analogue (PDB code: 1G 38) and free decamer d(GTTCGATGTC)(2) were taken. The x-ray structure shows two consecutive BII substates that are not observed in the free decamer. These consecutive BII substates are also observed during our simulation. Additionally, their facing backbones adopt the same conformations. These double facing BII substates are stable during the last 9 ns of the trajectories and result in a stretched DNA structure. On the other hand, protein-DNA contacts on 5' and 3' phosphodiester groups of the partner thymine of flipped adenine have changed. The sugar and phosphate parts of thymine have moved further into the empty space left by the flipping base without the influence of protein. Furthermore, readily high populated BII substates at the GpA step of palindromic tetrad TCGA rather than CpG step are observed in the free decamer. On the contrary, the BI substate at the GpA step is observed on the flipped adenine strand. A restrained MD simulation, reproducing the BI/BII pattern in the complex, demonstrated the influence of the unusual backbone conformation on the dynamical behavior of the target base. This finding along with the increased nearby interstrand phosphate distance is supportive to the N6-methylation mechanism.     

Homology-dependent inactivation of LTR retrotransposons in genomes of Aspergillus fumigatus and A. nidulans

Repeat-induced point mutation (RIP) is the most intriguing among the known mechanisms of repeated sequences inactivation because of its ability to produce irreversible mutation of repeated DNA. Discovered for the first time in Neurospora crassa, RIP is characterized by C:G to T:A transitions in duplicated sequences. The mechanisms and distribution of RIP are still purely investigated. Mobile elements are a common target for the processes which lead to homology-dependent silencing because of their ability to propagate themselves. We have done comparative analysis of LTR retrotransposons in genomic scale from genomes of two aspergilli fungi--Aspergillus funmigatus and A. nidulans, based on several copies we reconstructed "de-RIP" retroelements. Investigations of frequencies of CpG, CpA and TpG sites, which are potential targets for mutagenesis, showed the much lower frequencies of these sites in mobile elements in comparison with structural genes. LTR retrotransposons from A. fumigatus and A. nidulans have different ratio of types of substitutions. Our analysis indicates that two investigated fungi have or had the RIP-like processes for repeated sequences inactivation, in various modes. Whereas in A. fumigatus the context for mutagenesis consists of both CpG and CpA sites, in A. nidulans inactivation seems to proceed only on CpG dinucleotides. The present investigation gives a theoretical background for planning of experimental studying of RIP inactivation in aspergilli.     

BII nucleotides in the B and C forms of natural-sequence polymeric DNA: A new model for the C form of DNA

A combination of solid-state (31)P and (13)C NMR, X-ray diffraction, and model building is used to show that the B and C forms of fibrous macromolecular DNA consist of two distinct nucleotide conformations, which correspond closely to the BI and BII nucleotide conformations known from oligonucleotide crystals. The proportion of the BII conformation is higher in the C form than in the B form. We show structural models for a 10(1) double helix involving BI nucleotides and a 9(1) double helix involving BII nucleotides. The 10(1) BI model is similar to a previous model of B-form DNA, while the 9(1) BII model is novel. The BII model has a very deep and narrow minor groove, a shallow and wide major groove, and highly inclined bases. This work shows that the B to C transition in fibers corresponds to BI to BII conformational changes of the individual nucleotides.