Sequence-dependent DNA Research Articles

Transcriptional Regulation of the Clostridium cellulolyticum cip-cel Operon: a Complex Mechanism Involving a Catabolite-Responsive Element

The cip-cel cluster of genes plays an important role in the catabolism of the substrate cellulose by Clostridium cellulolyticum. It encodes several key components of the cellulosomes, including the scaffolding protein CipC and the major cellulase Cel48F. All the genes of this cluster display linked transcription, focusing attention on the promoter upstream from the first gene, cipC. We analyzed the regulation of the cipC promoter using a transcriptional fusion approach. A single promoter is located between nucleotides -671 and -643 with respect to the ATG start codon, and the large mRNA leader sequence is processed at position -194. A catabolite-responsive element (CRE) 414 nucleotides downstream from the transcriptional start site has been shown to be involved in regulating this operon by a carbon catabolite repression mechanism. This CRE is thought to bind a CcpA-like regulator complexed with a P-Ser-Crh-like protein. Sequences surrounding the promoter sequence may also be involved in direct (sequence-dependent DNA curvature) or indirect (unknown regulator binding) regulation.

Indirect readout in drug-DNA recognition: role of sequence-dependent DNA conformation.

DNA-binding drugs have numerous applications in the engineered gene regulation. However, the drug-DNA recognition mechanism is poorly understood. Drugs can recognize specific DNA sequences not only through direct contacts but also indirectly through sequence-dependent conformation, in a similar manner to the indirect readout mechanism in protein-DNA recognition. We used a knowledge-based technique that takes advantage of known DNA structures to evaluate the conformational energies. We built a dataset of non-redundant free B-DNA crystal structures to calculate the distributions of adjacent base-step and base-pair conformations, and estimated the effective harmonic potentials of mean force (PMF). These PMFs were used to calculate the conformational energy of drug-DNA complexes, and the Z-score as a measure of the binding specificity. Comparing the Z-scores for drug-DNA complexes with those for free DNA structures with the same sequence, we observed that in several cases the Z-scores became more negative upon drug binding. Furthermore, the specificity is position-dependent within the drug-bound region of DNA. These results suggest that DNA conformation plays an important role in the drug-DNA recognition. The presented method provides a tool for the analysis of drug-DNA recognition and can facilitate the development of drugs for targeting a specific DNA sequence.

Chemical biology that controls DNA structure and function

All genetic information, which is necessary for life, is written in 3 billion base pairs of DNA. Sequence-dependent local DNA conformations are also thought to play an important role in gene expression. Many diseases, including cancer, and hereditary and viral diseases, can now be understood at the DNA sequence and structure level. Therefore, direct control of the expression level of a specific gene would provide a promising approach for knowledge-based therapy. In this presentation I focus on our long-term efforts in controlling of structure and function of DNA by chemical biological approach.

Functional Specificity of a Hox Protein Mediated by the Recognition of Minor Groove Structure

The recognition of specific DNA-binding sites by transcription factors is a critical yet poorly understood step in the control of gene expression. Members of the Hox family of transcription factors bind DNA by making nearly identical major groove contacts via the recognition helices of their homeodomains. In vivo specificity, however, often depends on extended and unstructured regions that link Hox homeodomains to a DNA-bound cofactor, Extradenticle (Exd). Using a combination of structure determination, computational analysis, and in vitro and in vivo assays, we show that Hox proteins recognize specific Hox-Exd binding sites via residues located in these extended regions that insert into the minor groove but only when presented with the correct DNA sequence. Our results suggest that these residues, which are conserved in a paralog-specific manner, confer specificity by recognizing a sequence-dependent DNA structure instead of directly reading a specific DNA sequence.

Sequence-dependent DNA deformability studied using molecular dynamics simulations

Proteins recognize specific DNA sequences not only through direct contact between amino acids and bases, but also indirectly based on the sequence-dependent conformation and deformability of the DNA (indirect readout). We used molecular dynamics simulations to analyze the sequence-dependent DNA conformations of all 136 possible tetrameric sequences sandwiched between CGCG sequences. The deformability of dimeric steps obtained by the simulations is consistent with that by the crystal structures. The simulation results further showed that the conformation and deformability of the tetramers can highly depend on the flanking base pairs. The conformations of xATx tetramers show the most rigidity and are not affected by the flanking base pairs and the xYRx show by contrast the greatest flexibility and change their conformations depending on the base pairs at both ends, suggesting tetramers with the same central dimer can show different deformabilities. These results suggest that analysis of dimeric steps alone may overlook some conformational features of DNA and provide insight into the mechanism of indirect readout during protein–DNA recognition. Moreover, the sequence dependence of DNA conformation and deformability may be used to estimate the contribution of indirect readout to the specificity of protein–DNA recognition as well as nucleosome positioning and large-scale behavior of nucleic acids.

ABCB1 genetic polymorphism influences the pharmacology of the new pyrrolobenzodiazepine derivative SJG-136

ATP-binding cassette transporter P-glycoprotein (ABCB1) is responsible for the multidrug resistance (MDR1) phenotype observed in cancer cells. SJG-136, a new pyrrolobenzodiazepine dimer, is a sequence-dependent DNA crosslinking agent and substrate of ABCB1. We previously showed that colon cancer cell lines expressing high levels of ABCB1 showed a lower sensitivity to SJG-136. Here, we show that in 3T3 isogenic fibroblasts, ABCB1 genetic polymorphism differentially affects ABCB1 gene expression and transport function. However, this genotype-phenotype relationship was not observed in immortalized lymphocytes, which expressed 10- to 1000-fold less ABCB1 than colon cancer cell lines. Consistent with this, the cytotoxicity of SJG-136 in 3T3 fibroblasts was affected by ABCB1 genetic polymorphism but not in immortalized lymphocytes. ABCB1 genetic polymorphism is therefore likely to affect drug sensitivity in tissues expressing high levels of the transporter and in which significant variability is observed.

On Bifurcations of Equilibria of Intrinsically Curved, Electrically Charged, Rod-like Structures that Model DNA Molecules in Solution

DNA molecules in the familiar Watson–Crick double helical B form can be treated as though they have rod-like structures obtained by stacking dominoes one on top of another with each rotated by approximately one-tenth of a full turn with respect to its immediate predecessor in the stack. These “dominoes” are called base pairs. A recently developed theory of sequence-dependent DNA elasticity (Coleman, Olson, & Swigon, J. Chem. Phys. 118:7127–7140, 2003) takes into account the observation that the step from one base pair to the next can be one of several distinct types, each having its own mechanical properties that depend on the nucleotide composition of the step. In the present paper, which is based on that theory, emphasis is placed on the fact that, as each base in a base pair is attached to the sugar-phosphate backbone chain of one of the two DNA strands that have come together to form the Watson–Crick structure, and each phosphate group in a backbone chain bears one electronic charge, two such charges are associated with each base pair, which implies that each base pair is subject to not only the elastic forces and moments exerted on it by its neighboring base pairs but also to long range electrostatic forces that, because they are only partially screened out by positively charged counter ions, can render the molecule’s equilibrium configurations sensitive to changes in the concentration c of salt in the medium. When these electrostatic forces are taken into account, the equations of mechanical equilibrium for a DNA molecule with N + 1 base pairs are a system of μN non-linear equations, where μ, the number of kinematical variables describing the relative displacement and orientation of adjacent pairs is in general 6; it reduces to 3 when base-pair steps are assumed to be inextensible and non-shearable. As a consequence of the long-range electrostatic interactions of base pairs, the μN × μN Jacobian matrix of the equations of equilibrium is full. An efficient numerically stable computational scheme is here presented for solving those equations and determining the mechanical stability of the calculated equilibrium configurations. That scheme is employed to compute and analyze bifurcation diagrams in which c is the bifurcation parameter and to show that, for an intrinsically curved molecule, small changes in c can have a strong effect on stable equilibrium configurations. Cases are presented in which several stable configurations occur at a single value of c.

Structural basis for sequence-dependent DNA cleavage by nonspecific endonucleases.

Nonspecific endonucleases hydrolyze DNA without sequence specificity but with sequence preference, however the structural basis for cleavage preference remains elusive. We show here that the nonspecific endonuclease ColE7 cleaves DNA with a preference for making nicks after (at 3′O-side) thymine bases but the periplasmic nuclease Vvn cleaves DNA more evenly with little sequence preference. The crystal structure of the ‘preferred complex’ of the nuclease domain of ColE7 bound to an 18 bp DNA with a thymine before the scissile phosphate had a more distorted DNA phosphate backbone than the backbones in the non-preferred complexes, so that the scissile phosphate was compositionally closer to the endonuclease active site resulting in more efficient DNA cleavage. On the other hand, in the crystal structure of Vvn in complex with a 16 bp DNA, the DNA phosphate backbone was similar and not distorted in comparison with that of a previously reported complex of Vvn with a different DNA sequence. Taken together these results suggest a general structural basis for the sequence-dependent DNA cleavage catalyzed by nonspecific endonucleases, indicating that nonspecific nucleases could induce DNA to deform to distinctive levels depending on the local sequence leading to different cleavage rates along the DNA chain.

Operation of a DNA Robot Arm Inserted into a 2D DNA Crystalline Substrate

The success of nanorobotics requires the precise placement and subsequent operation of specific nanomechanical devices at particular locations. The structural programmability of DNA makes it a particularly attractive system for nanorobotics. We have developed a cassette that enables the placement of a robust, sequence-dependent DNA robot arm within a two-dimensional (2D) crystalline DNA array. The cassette contains the device, an attachment site, and a reporter of state. We used atomic force microscopy to demonstrate that the rotary device is fully functional after insertion. Thus, a nanomechanical device can operate within a fixed frame of reference.

Indirect Recognition in Sequence-specific DNA Binding by Escherichia coli Integration Host Factor: THE ROLE OF DNA DEFORMATION ENERGY

Integration host factor (IHF) is a bacterial histone-like protein whose primary biological role is to condense the bacterial nucleoid and to constrain DNA supercoils. It does so by binding in a sequence-independent manner throughout the genome. However, unlike other structurally related bacterial histone-like proteins, IHF has evolved a sequence-dependent, high affinity DNA-binding motif. The high affinity binding sites are important for the regulation of a wide range of cellular processes. A remarkable feature of IHF is that it employs an indirect readout mechanism to bind and wrap DNA at both the nonspecific and high affinity (sequence-dependent) DNA sites. In this study we assessed the contributions of pre-formed and protein-induced DNA conformations to the energetics of IHF binding. Binding energies determined experimentally were compared with energies predicted for the IHF-induced deformation of the DNA helix (DNA deformation energy) in the IHF-DNA complex. Combinatorial sets of de novo DNA sequences were designed to systematically evaluate the influence of sequence-dependent structural characteristics of the conserved IHF recognition elements of the consensus DNA sequence. We show that IHF recognizes pre-formed conformational characteristics of the consensus DNA sequence at high affinity sites, whereas at all other sites relative affinity is determined by the deformational energy required for nearest-neighbor base pairs to adopt the DNA structure of the bound DNA-IHF complex.

RNA used to control a DNA rotary nanomachine.

The PX-JX2 device is a robust sequence-dependent DNA nanomachine. It is controlled by the addition of DNA strands to the solution. The ability to control such devices with RNA will enable the machine to respond to signals generated by transcriptional logic circuits in vitro or, possibly, in vivo. The PX-JX2 device does not respond properly to RNA strands, so we have adapted the signaling method by using a cover-strand strategy, This approach entails using DNA control strands that contain sequence regions to set the device to either state. An RNA strand is used to mask one of these two regions, enabling the machine to assume the state corresponding to the unmasked region.

AFM imaging and theoretical modeling studies of sequence-dependent nucleosome positioning

Telomeric chromatin has different features with respect to bulk chromatin, since nucleosomal repeat along the chain is unusually short. We studied the role of telomeric DNA sequences on nucleosomal spacing in a model system. Nucleosomal arrays, assembled on a 1500-bp-long human telomeric DNA and on a DNA fragment containing 8 copies of the 601 strong nucleosome positioning sequence, have been studied at the single molecule level, by atomic force microscopy imaging. Random nucleosome positioning was found in the case of human telomeric DNA. On the contrary, nucleosome positioning on 601 DNA is characterized by preferential positions of nucleosome dyad axis each 200 bp. The AFM-derived nucleosome organization is in satisfactory agreement with that predicted by theoretical modeling, based on sequence-dependent DNA curvature and flexibility. The reported results show that DNA sequence has a main role, not only in mononucleosome thermodynamic stability, but also in the organization of nucleosomal arrays.

Modeling the Lac repressor-operator assembly: The influence of DNA looping on Lac repressor conformation

Repression of transcription of the Escherichia coli Lac operon by the Lac repressor (LacR) is accompanied by the simultaneous binding of LacR to two operators and the formation of a DNA loop. A recently developed theory of sequence-dependent DNA elasticity enables one to relate the fine structure of the LacR-DNA complex to a wide range of heretofore-unconnected experimental observations. Here, that theory is used to calculate the configuration and free energy of the DNA loop as a function of its length and base-pair sequence, its linking number, and the end conditions imposed by the LacR tetramer. The tetramer can assume two types of conformations. Whereas a rigid V-shaped structure is observed in the crystal, EM images show extended forms in which two dimer subunits are flexibly joined. Upon comparing our computed loop configurations with published experimental observations of permanganate sensitivities, DNase I cutting patterns, and loop stabilities, we conclude that linear DNA segments of short-to-medium chain length (50-180 bp) give rise to loops with the extended form of LacR and that loops formed within negatively supercoiled plasmids induce the V-shaped structure.

Electronic density of states in sequence dependent DNA molecules

We report in this work a numerical study of the electronic density of states (DOS) in π-stacked arrays of DNA single-strand segments made up from the nucleotides guanine G, adenine A, cytosine C and thymine T, forming a Rudin–Shapiro (RS) as well as a Fibonacci (FB) polyGC quasiperiodic sequences. Both structures are constructed starting from a G nucleotide as seed and following their respective inflation rules. Our theoretical method uses Dyson’s equation together with a transfer-matrix treatment, within an electronic tight-binding Hamiltonian model, suitable to describe the DNA segments modelled by the quasiperiodic chains. We compared the DOS spectra found for the quasiperiodic structure to those using a sequence of natural DNA, as part of the human chromosome Ch22, with a remarkable concordance, as far as the RS structure is concerned. The electronic spectrum shows several peaks, corresponding to localized states, as well as a striking self-similar aspect.

Ability of spermine to differentiate between DNA sequences—Preferential stabilization of A-tracts

The regulatory roles fulfilled by polyamines by governance of chromatin structure are made possible by their strong association with cellular DNA, and hence by their ability to modulate DNA structure and function. Towards this end, it is crucial to understand the manifestation of sequence-dependent polyamine binding at the secondary and tertiary structural levels of DNA. This study utilizes circular dichroism (CD) and isothermal titration calorimetry (ITC) to address this relationship by using 20 bp oligonucleotides with sequences–poly(dA):poly(dT), poly(dAdT):poly(dAdT), poly(dG):poly(dC), poly(dGdC):poly(dGdC)–that yield physiologically relevant structures, and poly(dIdC):poly(dIdC). CD studies show that at physiological ionic strength (150 mM NaCl), spermine preferentially stabilizes A-tracts, and increases flexibility of the G-tract oligomer; the latter is also suggested by the larger change in entropy (Δ S) of spermine binding to G-tracts. Given the chromatin destabilizing property of these sequences, these findings suggest a role for spermine in stabilization of non-nucleosomal A-tracts, and a compensating mechanism for incorporation of G-tracts in the chromatin structure. Other implications of these findings in sequence dependent DNA packaging are discussed.

2P173 Sequence-dependent DNA conformation and flexibility characterized by hydration pattern(36. DNA to chromatin,Poster Session,Abstract,Meeting Program of EABS & BSJ 2006)

2P173 Sequence-dependent DNA conformation and flexibility characterized by hydration pattern(36. DNA to chromatin,Poster Session,Abstract,Meeting Program of EABS & BSJ 2006)

Sequence Dependence of Charge Transport through DNA Domains

Here we examine the photooxidation of two kinetically fast electron hole traps, N4-cyclopropylcytosine (CPC) and N2-cyclopropylamine-guanosine (CPG), incorporated in DNA duplexes of various sequence using different photooxidants. DNA oxidation studies are carried out either with noncovalently bound [Ru(phen)(dppz)(bpy')]3+ (dppz = dipyridophenazine) and [Rh(phi)2(bpy)]3+ (phi = phenanthrenequinone diimine) or with anthraquinone tethered to DNA. Because the cyclopropylamine-substituted bases decompose rapidly upon oxidation, their efficiency of decomposition provides a measure of relative hole localization. Consistent with a higher oxidation potential for CPC versus CPG in DNA, CPC decomposes with photooxidation by [Rh(phi)2(bpy)]3+, while CPG undergoes ring-opening both with photoexcited [Rh(phi)2(bpy)]3+ and with [Ru(phen)(dppz)(bpy')]3+. Anthraquinone-modified DNA assemblies of identical base composition but different base sequence are also probed. Single and double base substitutions within adenine tracts modulate CPC decomposition. In fact, the entire sequence within the DNA assembly is seen to govern CPC oxidation, not simply the bases intervening between CPC and the tethered photooxidant. These data are reconciled in the context of a mechanistic model of conformationally gated charge transport through delocalized DNA domains. Photooxidations of anthraquinone-modified DNA assemblies containing both CPC and CPG, but with varied distances separating the modified bases, point to a domain size of at least three bases. Our model for DNA charge transport is distinguished from polaron models. In our model, delocalized domains within the base pair stack form transiently based upon sequence-dependent DNA structure and dynamics. Given these results, DNA charge transport is indeed remarkably sensitive to DNA sequence and structure.

Role of inter and intramolecular interactions in protein–DNA recognition

Protein–DNA recognition plays an essential role in the regulation of gene expression. Regulatory proteins are known to recognize specific DNA sequences directly through atomic contacts between protein and DNA, and/or indirectly through the conformational properties of the DNA. In this work, we have analyzed the specificity of intermolecular interactions by statistical analysis of base–amino acid interactions within protein–DNA complexes as well as the computer simulations of base–amino acid interactions. The specificity of the intramolecular interactions was studied by statistical analysis of the sequence-dependent DNA conformational parameters and the elastic properties of DNA. Systematic comparison of these specificities in a large number of protein–DNA complexes revealed that both intermolecular and intramolecular interactions contribute to the specificity of protein–DNA recognition, and their relative contributions vary depending upon the protein–DNA complex. We demonstrated that combination of the intermolecular and intramolecular energies leads to enhanced specificity and the combined energy could explain experimental data on binding affinity changes caused by base mutations. These results provided new insight into the relationship between specificity and structure in the process of protein–DNA recognition, which would lead to prediction of specific protein–DNA binding sites.

Reversed DNA Strand Cleavage Specificity in Initiation of Cre–LoxP Recombination Induced by the His289Ala Active-site Substitution

During the first steps of site-specific recombination, Cre protein cleaves and religates a specific homologous pair of LoxP strands to form a Holliday junction (HJ) intermediate. The HJ is resolved into recombination products through exchange of the second homologous strand pair. CreH289A, containing a His to Ala substitution in the conserved R-H-R catalytic motif, has a 150-fold reduced recombination rate and accumulates HJs. However, to produce these HJs, CreH289A exchanges the opposite set of strands compared to wild-type Cre (CreWT). To investigate how CreH289A and CreWT impose strand exchange order, we characterized their reactivities and strand cleavage preferences toward LoxP duplex and HJ substrates containing 8bp spacer substitutions. Remarkably, CreH289A had different and often opposite strand exchange preferences compared to CreWT with nearly all substrates. CreH289N was much less perturbed, implying that overall recombination rate and strand exchange depend more on His289 hydrogen bonding capability than on its acid/base properties. LoxP substitutions immediately 5' (S1 nucleotide) or 3' (S1' nucleotide) of the scissile phosphate had large effects on substrate utilization and strand exchange order. S1' substitutions, designed to alter base-unstacking events concomitant with Cre-induced LoxP bending, caused HJ accumulation and dramatically inverted the cleavage preferences. That pre-formed HJs were resolved via either strand in vitro suggests that inhibition of the "conformational switch" isomerization required to trigger the second strand exchange accounts for the observed HJ accumulation. Rather than reflecting CreWT behavior, CreH289A accumulates HJs of opposite polarity through a combination of its unique cleavage specificity and an HJ isomerization defect. The overall implication is that cleavage specificity is mediated by sequence-dependent DNA deformations that influence the scissile phosphate positioning and reactivity. A role of His289 may be to selectively stabilize the "activated" phosphate conformation in order to promote cleavage.

Brownian dynamics simulations of sequence-dependent duplex denaturation in dynamically superhelical DNA

The topological state of DNA in vivo is dynamically regulated by a number of processes that involve interactions with bound proteins. In one such process, the tracking of RNA polymerase along the double helix during transcription, restriction of rotational motion of the polymerase and associated structures, generates waves of overtwist downstream and undertwist upstream from the site of transcription. The resulting superhelical stress is often sufficient to drive double-stranded DNA into a denatured state at locations such as promoters and origins of replication, where sequence-specific duplex opening is a prerequisite for biological function. In this way, transcription and other events that actively supercoil the DNA provide a mechanism for dynamically coupling genetic activity with regulatory and other cellular processes. Although computer modeling has provided insight into the equilibrium dynamics of DNA supercoiling, to date no model has appeared for simulating sequence-dependent DNA strand separation under the nonequilibrium conditions imposed by the dynamic introduction of torsional stress. Here, we introduce such a model and present results from an initial set of computer simulations in which the sequences of dynamically superhelical, 147 base pair DNA circles were systematically altered in order to probe the accuracy with which the model can predict location, extent, and time of stress-induced duplex denaturation. The results agree both with well-tested statistical mechanical calculations and with available experimental information. Additionally, we find that sites susceptible to denaturation show a propensity for localizing to supercoil apices, suggesting that base sequence determines locations of strand separation not only through the energetics of interstrand interactions, but also by influencing the geometry of supercoiling.