Sequence-Dependent Biophysical Modeling of DNA Amplification

A theoretical framework for prediction of the dynamic evolution of chemical species in DNA amplification reactions, for any specified sequence and operating conditions, is reported. Using the polymerase chain reaction (PCR) as an example, we developed a sequence- and temperature-dependent kinetic model for DNA amplification using first-principles biophysical modeling of DNA hybridization and polymerization. We compare this kinetic model with prior PCR models and discuss the features of our model that are essential for quantitative prediction of DNA amplification efficiency for arbitrary sequences and operating conditions. Using this model, the kinetics of PCR is analyzed. The ability of the model to distinguish between the dynamic evolution of distinct DNA sequences in DNA amplification reactions is demonstrated. The kinetic model is solved for a typical PCR temperature protocol to motivate the need for optimization of the dynamic operating conditions of DNA amplification reactions. It is shown that amplification efficiency is affected by dynamic processes that are not accurately represented in the simplified models of DNA amplification that form the basis of conventional temperature cycling protocols. Based on this analysis, a modified temperature protocol that improves PCR efficiency is suggested. Use of this sequence-dependent kinetic model in a control theoretic framework to determine the optimal dynamic operating conditions of DNA amplification reactions, for any specified amplification objective, is discussed.     

Sequence-Dependent Biophysical Modeling of DNA Amplification

A theoretical framework for prediction of the dynamic evolution of chemical species in DNA amplification reactions, for any specified sequence and operating conditions, is reported. Using the polymerase chain reaction (PCR) as an example, we developed a sequence- and temperature-dependent kinetic model for DNA amplification using first-principles biophysical modeling of DNA hybridization and polymerization. We compare this kinetic model with prior PCR models and discuss the features of our model that are essential for quantitative prediction of DNA amplification efficiency for arbitrary sequences and operating conditions. Using this model, the kinetics of PCR is analyzed. The ability of the model to distinguish between the dynamic evolution of distinct DNA sequences in DNA amplification reactions is demonstrated. The kinetic model is solved for a typical PCR temperature protocol to motivate the need for optimization of the dynamic operating conditions of DNA amplification reactions. It is shown that amplification efficiency is affected by dynamic processes that are not accurately represented in the simplified models of DNA amplification that form the basis of conventional temperature cycling protocols. Based on this analysis, a modified temperature protocol that improves PCR efficiency is suggested. Use of this sequence-dependent kinetic model in a control theoretic framework to determine the optimal dynamic operating conditions of DNA amplification reactions, for any specified amplification objective, is discussed.     

Sequence-Dependent Stability of Intramolecular DNA Triple Helices

Abstract

A set of 21 oligodeoxynucleotides were designed to fold into intramolecular triple helices of the pyrimidine motif under appropriate conditions. UV melting experiments on the triplexes which only differ in the number and distribution of third strand cytosines reveal the influence of sequence and pH on triplex stability and can be summarized as follows: (1) increasing the cytosine content in the third strand results in a higher thermal stability of the triplex at acidic pH but lowers the triplex to duplex melting temperature at neutral pH; (2) cytosines at terminal positions destabilize the triple helical structure as compared to non-terminal positions; (3) contiguous cytosines lead to a pH dependent destabilization of the triplex, the destabilizing effect being more pronounced at higher pH. Analysis of these effects in terms of the various interactions within a triple helical complex indicate that the sequence-dependent stabilities are largely determined by the extent of protonation for individual third strand cytosines.     

A Model of Sequence-Dependent Protein Diffusion Along DNA

We introduce a probabilistic model for protein sliding motion along DNA during the search of a target sequence. The model accounts for possible effects due to sequence-dependent interaction between the nonspecific DNA and the protein. Hydrogen bonds formed at the target site are used as the main sequence-dependent interaction between protein and DNA. The resulting dynamical properties and the possibility of an experimental verification are discussed in details. We show that, while at large times the process reaches a linear diffusion regime, it initially displays a sub-diffusive behavior. The sub-diffusive regime can last sufficiently long to be of biological interest.     

Nucleosomes Shape DNA Polymorphism and Divergence

An estimated 80% of genomic DNA in eukaryotes is packaged as nucleosomes, which, together with the remaining interstitial linker regions, generate higher order chromatin structures [1]. Nucleosome sequences isolated from diverse organisms exhibit ∼10 bp periodic variations in AA, TT and GC dinucleotide frequencies. These sequence elements generate intrinsically curved DNA and help establish the histone-DNA interface. We investigated an important unanswered question concerning the interplay between chromatin organization and genome evolution: do the DNA sequence preferences inherent to the highly conserved histone core exert detectable natural selection on genomic divergence and polymorphism? To address this hypothesis, we isolated nucleosomal DNA sequences from Drosophila melanogaster embryos and examined the underlying genomic variation within and between species. We found that divergence along the D. melanogaster lineage is periodic across nucleosome regions with base changes following preferred nucleotides, providing new evidence for systematic evolutionary forces in the generation and maintenance of nucleosome-associated dinucleotide periodicities. Further, Single Nucleotide Polymorphism (SNP) frequency spectra show striking periodicities across nucleosomal regions, paralleling divergence patterns. Preferred alleles occur at higher frequencies in natural populations, consistent with a central role for natural selection. These patterns are stronger for nucleosomes in introns than in intergenic regions, suggesting selection is stronger in transcribed regions where nucleosomes undergo more displacement, remodeling and functional modification. In addition, we observe a large-scale (∼180 bp) periodic enrichment of AA/TT dinucleotides associated with nucleosome occupancy, while GC dinucleotide frequency peaks in linker regions. Divergence and polymorphism data also support a role for natural selection in the generation and maintenance of these super-nucleosomal patterns. Our results demonstrate that nucleosome-associated sequence periodicities are under selective pressure, implying that structural interactions between nucleosomes and DNA sequence shape sequence evolution, particularly in introns.     

Sequence-Dependent Elasticity and Electrostatics of Single-Stranded DNA: Signatures of Base-Stacking

Base-stacking is a key factor in the energetics that determines nucleic acid structure. We measure the tensile response of single-stranded DNA as a function of sequence and monovalent salt concentration to examine the effects of base-stacking on the mechanical and thermodynamic properties of single-stranded DNA. By comparing the elastic response of highly stacked poly(dA) and that of a polypyrimidine sequence with minimal stacking, we find that base-stacking in poly(dA) significantly enhances the polymer’s rigidity. The unstacking transition of poly(dA) at high force reveals that the intrinsic electrostatic tension on the molecule varies significantly more weakly on salt concentration than mean-field predictions. Further, we provide a model-independent estimate of the free energy difference between stacked poly(dA) and unstacked polypyrimidine, finding it to be ∼−0.25 k_BT/base and nearly constant over three orders of magnitude in salt concentration.     

Sequence-Dependent Structural Variations of DNA Revealed by Chemical Ligation

Abstract

The comparative study of nick-sealing efficiency under the action of BrCN or water-soluble carbodiimide was carried out for 14 dinucleotide combinations in double helix. The difference between the lowest (17%, GpG) and the highest (94%, CpT) coupling yields was found to be more than five fold, both condensing reagents showing a similar sequence-specific trend. A strong correlation observed between coupling yields at different dinucleotide combinations and ³¹P NMR parameters supports the idea that variations in chemical ligation efficiency arise from sequence-specific modulations of the helix geometry and confirms close similarity of the intrinsic fine structure of intact and nicked DNAs.     

The Duplexing of the Genetic Code and Sequence-Dependent DNA Geometry

It is well known that sequences of bases in DNA are translated into sequences of amino acids in cells via the genetic code. More recently, it has been discovered that the sequence of DNA bases also influences the geometry and deformability of the DNA. These two correspondences represent a naturally arising example of duplexed codes, providing two different ways of interpreting the same DNA sequence. This paper will set up the notation and basic results necessary to mathematically investigate the relationship between these two natural DNA codes. It then undertakes two very different such investigations: one graphical approach based only on expected values and another analytic approach incorporating the deformability of the DNA molecule and approximating the mutual information of the two codes. Special emphasis is paid to whether there is evidence that pressure to maximize the duplexing efficiency influenced the evolution of the genetic code. Disappointingly, the results fail to support the hypothesis that the genetic code was influenced in this way. In fact, applying both methods to samples of realistic alternative genetic codes shows that the duplexing of the genetic code found in nature is just slightly less efficient than average. The implications of this negative result are considered in the final section of the paper.     

Sequence-Dependent Elasticity and Electrostatics of Single-Stranded DNA: Signatures of Base-Stacking

Base-stacking is a key factor in the energetics that determines nucleic acid structure. We measure the tensile response of single-stranded DNA as a function of sequence and monovalent salt concentration to examine the effects of base-stacking on the mechanical and thermodynamic properties of single-stranded DNA. By comparing the elastic response of highly stacked poly(dA) and that of a polypyrimidine sequence with minimal stacking, we find that base-stacking in poly(dA) significantly enhances the polymer’s rigidity. The unstacking transition of poly(dA) at high force reveals that the intrinsic electrostatic tension on the molecule varies significantly more weakly on salt concentration than mean-field predictions. Further, we provide a model-independent estimate of the free energy difference between stacked poly(dA) and unstacked polypyrimidine, finding it to be ∼−0.25 k_BT/base and nearly constant over three orders of magnitude in salt concentration.     

Sequence-Dependent DNA Separation by Anion-Exchange High-Performance Liquid Chromatography

High-performance liquid chromatography (HPLC) system with a new nonporous anion-exchange resin, DNA–NPR, made it possible to rapidly separate DNA fragments up to 20 kbp with high resolution. In order to further characterize this chromatographic DNA separation system, we prepared a mixture of double-stranded DNAs of constant length carrying a fully degenerated 50-bp region and analyzed their chromatographic behavior on the DNA–NPR column. The results indicated that the separation of DNA fragments on the anion-exchange HPLC was governed not only by size, but also by nucleotide sequence: even DNA fragments with the same size and the same base content could be separated on this column. Taking advantage of this characteristic feature of the anion-exchange HPLC, we could readily fractionate human cDNAs with practically acceptable recovery and high resolution. Furthermore, the combination of HPLC and gel electrophoresis realized separation of a mixture of DNA fragments in a two-dimensional pattern.     

Nucleosomes protect DNA from DNA methylation in vivo and in vitro

Positioned nucleosomes limit the access of proteins to DNA. However, the impact of nucleosomes on DNA methylation in vitro and in vivo is poorly understood. Here, we performed a detailed analysis of nucleosome binding and nucleosomal DNA methylation by the de novo methyltransferases. We show that compared to linker DNA, nucleosomal DNA is largely devoid of CpG methylation. ATP-dependent chromatin remodelling frees nucleosomal CpG dinucleotides and renders the remodelled nucleosome a 2-fold better substrate for Dnmt3a methyltransferase compared to free DNA. These results reflect the situation in vivo, as quantification of nucleosomal DNA methylation levels in HeLa cells shows a 2-fold decrease of nucleosomal DNA methylation levels compared to linker DNA. Our findings suggest that nucleosomal positions are stably maintained in vivo and nucleosomal occupancy is a major determinant of global DNA methylation patterns in vivo.     

Placeholder Nucleosomes Underlie Germline-to-Embryo DNA Methylation Reprogramming

1. Download high-res image (222KB)
2. Download full-size image

     

Initial stages of DNA Base Excision Repair in Nucleosomes

Molecular Biology - In mammalian cells, base excision repair (BER) is the main pathway responsible for the correction of a variety of chemically modified DNA bases. DNA packaging in chromatin...     

Herpes Simplex Virus 1 DNA Is in Unstable Nucleosomes throughout the Lytic Infection Cycle,and the Instability of the Nucleosomes Is Independent of DNA Replication

Herpes simplex virus 1 (HSV-1) DNA is chromatinized during latency and consequently regularly digested by micrococcal nuclease (MCN) to nucleosome-size fragments. In contrast, MCN digests HSV-1 DNA in lytically infected cells to mostly heterogeneous sizes. Yet HSV-1 DNA coimmunoprecipitates with histones during lytic infections. We have shown that at 5 h postinfection, most nuclear HSV-1 DNA is in particularly unstable nucleoprotein complexes and consequently is more accessible to MCN than DNA in cellular chromatin. HSV-1 DNA was quantitatively recovered at this time in complexes with the biophysical properties of mono- to polynucleosomes following a modified MCN digestion developed to detect potential unstable intermediates. We proposed that most HSV-1 DNA is in unstable nucleosome-like complexes during lytic infections. Physiologically, nucleosome assembly typically associates with DNA replication, although DNA replication transiently disrupts nucleosomes. It therefore remained unclear whether the instability of the HSV-1 nucleoprotein complexes was related to the ongoing viral DNA replication. Here we tested whether HSV-1 DNA is in unstable nucleosome-like complexes before, during, or after the peak of viral DNA replication or when HSV-1 DNA replication is inhibited. HSV-1 DNA was quantitatively recovered in complexes fractionating as mono- to polynucleosomes from nuclei harvested at 2, 5, 7, or 9 h after infection, even if viral DNA replication was inhibited. Therefore, most HSV-1 DNA is in unstable nucleosome-like complexes throughout the lytic replication cycle, and the instability of these complexes is surprisingly independent of HSV-1 DNA replication. The specific accessibility of nuclear HSV-1 DNA, however, varied at different times after infection.     

Nucleosomes positioned by ORC facilitate the initiation of DNA replication

The packaging of eukaryotic DNA into nucleosomes is a critical regulator of nuclear events. To address the interplay between chromatin and replication initiation, we have assessed the determinants and function of the nucleosomal configuration of S. cerevisiae replication origins. Using in vitro and in vivo assays, we demonstrate that the yeast initiator, the origin recognition complex (ORC), is required to maintain the nucleosomal configuration adjacent to origins. Disruption of the ORC-directed nucleosomal arrangement at an origin interferes with initiation of replication, but does not alter the association of ORC with the origin. Instead, the nucleosomes positioned by ORC are important for prereplicative complex formation. These findings suggest that origin-proximal nucleosomes facilitate replication initiation, and that local chromatin structure affects origin function.     

Selective Anchoring of DNA Methyltransferases 3A and 3B to Nucleosomes Containing Methylated DNA

Proper DNA methylation patterns are essential for mammalian development and differentiation. DNA methyltransferases (DNMTs) primarily establish and maintain global DNA methylation patterns; however, the molecular mechanisms for the generation and inheritance of methylation patterns are still poorly understood. We used sucrose density gradients of nucleosomes prepared by partial and maximum micrococcal nuclease digestion, coupled with Western blot analysis to probe for the interactions between DNMTs and native nucleosomes. This method allows for analysis of the in vivo interactions between the chromatin modification enzymes and their actual nucleosomal substrates in the native state. We show that little free DNA methyltransferase 3A and 3B (DNMT3A/3B) exist in the nucleus and that almost all of the cellular contents of DNMT3A/3B, but not DNMT1, are strongly anchored to a subset of nucleosomes. This binding of DNMT3A/3B does not require the presence of other well-known chromatin-modifying enzymes or proteins, such as proliferating cell nuclear antigen, heterochromatin protein 1, methyl-CpG binding protein 2, Enhancer of Zeste homolog 2, histone deacetylase 1, and UHRF1, but it does require an intact nucleosomal structure. We also show that nucleosomes containing methylated SINE and LINE elements and CpG islands are the main sites of DNMT3A/3B binding. These data suggest that inheritance of DNA methylation requires cues from the chromatin component in addition to hemimethylation.Proper DNA methylation patterns are essential for mammalian development and differentiation. More than three decades ago, de novo cytosine DNA methylation and its maintenance were proposed to exist in eukaryotic cells (29, 54); however, the molecular mechanisms for the generation and inheritance of methylation patterns are still poorly understood. DNA methyltransferases (DNMTs) DNMT1, DNMT3A, and DNMT3B primarily establish and maintain global DNA methylation patterns (39, 48). DNMT1 preferentially methylates hemimethylated DNA in vitro (7) and is tethered to replication foci during S phase (38). In contrast, DNMT3A and DNMT3B (DNMT3A/3B) have no preference for hemimethylated DNA (49) and are required for de novo methylation of genomic DNA (48). It has been thought that DNMT1 acts mainly as a “maintenance methyltransferase” during DNA synthesis and that DNMT3A and DNMT3B act as “de novo” enzymes. However, more recent studies indicate that DNMT1 may also be required for de novo methylation of genomic DNA (17, 30) and that DNMT3A/3B are also required for maintenance functions (11, 40, 55). Furthermore, the different DNMTs cooperate in maintaining the methylation of some regions of the genome, particularly repetitive elements (40, 53).Recruitment of individual DNMT enzymes to different regions of chromatin in vivo, particularly to gene regulatory regions, may require interaction with auxiliary factors (28, 36). DNMT1, which is diffusely localized throughout nuclei in non-S-phase cells (38), is targeted to replication foci by interacting with proliferating cell nuclear antigen (PCNA) (15) and also physically interacts with UHRF1 (ubiquitinlike, containing PHD and RING finger domains 1) that binds to hemimethylated DNA (3, 4, 8, 27, 62). DNMT3 enzymes are usually found localized to heterochromatin regions in most transient-expression assays (5, 12). As genomic DNA in chromatin is packaged into nucleosomes which might limit the accessibility of target sites to the enzymes, the interaction of DNMTs with nucleosomes in a chromatin context is important for the regulation of genomic methylation.Genetic and biochemical studies have provided many insights into the distinct and cooperative functions of the DNMT enzymes; however, few of these studies have addressed how they interact with chromatin in vivo. Recombinant DNMT1 and DNMT3 enzymes can methylate the CpG sites on nucleosomes assembled in vitro (26, 50, 56, 65). Recently DNMT3L has been found to connect DNMT3A2 to nucleosomes in embryonic stem cells (52). However, DNMT3L is expressed only during gametogenesis and embryonic stages (1, 9), suggesting that other mechanisms might be necessary for directing the enzyme to specific chromatin regions in somatic cells.In the present study, we investigated how different DNMT enzymes interact with chromatin at the nucleosomal level in somatic cell lines. Micrococcal nuclease (MNase) treatment of nuclei in a low-ionic-strength buffer digests nucleosomal linker DNA regions, thereby minimizing the disruption of protein complexes on the nucleosomes. We prepared nucleosomes from partial or maximum MNase-digested nuclei and resolved them on sucrose density gradients to analyze their interactions with chromatin proteins. The results indicate that while DNMT1 interacts primarily with linker DNA, DNMT3A/3B enzymes interact strongly with nucleosomes containing methylated repetitive elements and also containing methylated CpG islands (CGIs) and may not require additional proteins for this strong binding. These data are particularly intriguing in that they provide insights into the mechanisms of the interaction of DNMTs with chromatin and maintenance of DNA methylation in somatic cells.     

Bending and flexibility of kinetoplast DNA

We have evaluated the extent of bending at an anomalous locus in DNA restriction fragments from the kinetoplast body of Leishmania tarentolae using transient electric dichroism to measure the rate of rotational diffusion of DNA fragments in solution. We compare the rate of rotational diffusion of two fragments identical in sequence except for circular permutation, which places the bend near the center in one case and near one end of the molecule in the other. Hydrodynamic theory was used to conclude that the observed 20% difference in rotational relaxation times is a consequence of an overall average bending angle of 84 +/- 6 degrees between the end segments of the fragment that contains the bending locus near its center. If it is assumed that bending results from structural dislocations at the junctions between oligo(dA).oligo(dT) tracts and adjacent segments of B DNA, a bend angle of 9 +/- 0.5 degrees at each junction is required to explain the observations. The extent of bending is little affected by ionic conditions and is weakly dependent on temperature. Comparison of one of the anomalous fragments with an electrophoretically normal control fragment leads to the conclusion that they differ measurably in apparent stiffness, consistent with a significantly increased persistence length or contour length in the kinetoplast fragments.     

Characterization of Dnmt1 Binding and DNA Methylation on Nucleosomes and Nucleosomal Arrays

The packaging of DNA into nucleosomes and the organisation into higher order structures of chromatin limits the access of sequence specific DNA binding factors to DNA. In cells, DNA methylation is preferentially occuring in the linker region of nucleosomes, suggesting a structural impact of chromatin on DNA methylation. These observations raise the question whether DNA methyltransferases are capable to recognize the nucleosomal substrates and to modify the packaged DNA. Here, we performed a detailed analysis of nucleosome binding and nucleosomal DNA methylation by the maintenance DNA methyltransferase Dnmt1. Our binding studies show that Dnmt1 has a DNA length sensing activity, binding cooperatively to DNA, and requiring a minimal DNA length of 20 bp. Dnmt1 needs linker DNA to bind to nucleosomes and most efficiently recognizes nucleosomes with symmetric DNA linkers. Footprinting experiments reveal that Dnmt1 binds to both DNA linkers exiting the nucleosome core. The binding pattern correlates with the efficient methylation of DNA linkers. However, the enzyme lacks the ability to methylate nucleosomal CpG sites on mononucleosomes and nucleosomal arrays, unless chromatin remodeling enzymes create a dynamic chromatin state. In addition, our results show that Dnmt1 functionally interacts with specific chromatin remodeling enzymes to enable complete methylation of hemi-methylated DNA in chromatin.