Helix bending as a factor in protein/DNA recognition

DNA selection by the master transcription factor PU.1

The master transcriptional regulator PU.1/Spi-1 engages DNA sites with affinities spanning multiple orders of magnitude. To elucidate this remarkable plasticity, we have characterized 22 high-resolution co-crystallographic PU.1/DNA complexes across the addressable affinity range in myeloid gene transactivation. Over a purine-rich core (such as 5'-GGAA-3') flanked by variable sequences, affinity is negotiated by direct readout on the 5' flank via a critical glutamine (Q226) sidechain and by indirect readout on the 3' flank by sequence-dependent helical flexibility. Direct readout by Q226 dynamically specifies PU.1's characteristic preference for purines and explains the pathogenic mutation Q226E in Waldenström macroglobulinemia. The structures also reveal how disruption of Q226 mediates strand-specific inhibition by DNA methylation and the recognition of non-canonical sites, including the authentic binding sequence at the CD11b promoter. A re-synthesis of phylogenetic and structural data on the ETS family, considering the centrality of Q226 in PU.1, unifies the model of DNA selection by ETS proteins.

Structural Characterization of the Interaction of Hypoxia Inducible Factor-1 with Its Hypoxia Responsive Element at the -964G > A Variation Site of the HLA-G Promoter Region

Human Antigen Leukocyte-G (HLA-G) gene encodes an immune checkpoint molecule that has restricted tissue expression in physiological conditions; however, the gene may be induced in hypoxic conditions by the interaction with the hypoxia inducible factor-1 (HIF1). Hypoxia regulatory elements (HRE) located at the HLA-G promoter region and at exon 2 are the major HIF1 target sites. Since the G allele of the −964G > A transversion induces higher HLA-G expression when compared to the A allele in hypoxic conditions, here we analyzed HIF1-HRE complex interaction at the pair-atom level considering both −964G > A polymorphism alleles. Mouse HIF2 dimer crystal (Protein Data Bank ID: 4ZPK) was used as template to perform homology modelling of human HIF1 quaternary structure using MODELLER v9.14. Two 3D DNA structures were built from 5′GCRTG’3 HRE sequence containing the −964G/A alleles using x3DNA. Protein-DNA docking was performed using the HADDOCK v2.4 server, and non-covalent bonds were computed by DNAproDB server. Molecular dynamic simulation was carried out per 200 ns, using Gromacs v.2019. HIF1 binding in the HRE containing −964G allele results in more hydrogen bonds and van der Waals contact formation than HRE with −964A allele. Protein-DNA complex trajectory analysis revealed that HIF1-HRE-964G complex is more stable. In conclusion, HIF1 binds in a more stable and specific manner at the HRE with G allele.

CTCF binding modulates UV damage formation to promote mutation hot spots in melanoma

Somatic mutations in DNA-binding sites for CCCTC-binding factor (CTCF) are significantly elevated in many cancers. Prior analysis has suggested that elevated mutation rates at CTCF-binding sites in skin cancers are a consequence of the CTCF-cohesin complex inhibiting repair of UV damage. Here, we show that CTCF binding modulates the formation of UV damage to induce mutation hot spots. Analysis of genome-wide CPD-seq data in UV-irradiated human cells indicates that formation of UV-induced cyclobutane pyrimidine dimers (CPDs) is primarily suppressed by CTCF binding but elevated at specific locations within the CTCF motif. Locations of CPD hot spots in the CTCF-binding motif coincide with mutation hot spots in melanoma. A similar pattern of damage formation is observed at CTCF-binding sites in vitro, indicating that UV damage modulation is a direct consequence of CTCF binding. We show that CTCF interacts with binding sites containing UV damage and inhibits repair by a model repair enzyme in vitro. Structural analysis and molecular dynamic simulations reveal the molecular mechanism for how CTCF binding modulates CPD formation.

Application of portrait recognition system for emergency evacuation in mass emergencies

A portrait recognition system can play an important role in emergency evacuation in mass emergencies. This paper designed a portrait recognition system, analyzed the overall structure of the system and the method of image preprocessing, and used the Single Shot MultiBox Detector (SSD) algorithm for portrait detection. It also designed an improved algorithm combining principal component analysis (PCA) with linear discriminant analysis (LDA) for portrait recognition and tested the system by applying it in a shopping mall to collect and monitor the portrait and establish a data set. The results showed that the missing detection rate and false detection rate of the SSD algorithm were 0.78 and 2.89%, respectively, which were lower than those of the AdaBoost algorithm. Comparisons with PCA, LDA, and PCA + LDA algorithms demonstrated that the recognition rate of the improved PCA + LDA algorithm was the highest, which was 95.8%, the area under the receiver operating characteristic curve was the largest, and the recognition time was the shortest, which was 465 ms. The experimental results show that the improved PCA + LDA algorithm is reliable in portrait recognition and can be used for emergency evacuation in mass emergencies.

A Conformational Switch in the Zinc Finger Protein Kaiso Mediates Differential Readout of Specific and Methylated DNA Sequences

Recognition of the epigenetic mark 5-methylcytosine (mC) at CpG sites in DNA has emerged as a novel function of many eukaryotic transcription factors (TFs). It remains unclear why the sequence specificity of these TFs differs for CpG-methylated motifs and consensus motifs. Here, we dissect the structural and dynamic basis for this differential DNA binding specificity in the human zinc finger TF Kaiso, which exhibits high affinity for two consecutive mCpG sites in variable contexts and also for a longer, sequence-specific Kaiso binding site (KBS). By integrating structural analysis and DNA binding studies with targeted protein mutagenesis and nucleotide substitutions, we identify distinct mechanisms for readout of methylated and KBS motifs by Kaiso. We show that a key glutamate residue (E535), critical for mCpG site recognition, adopts different conformations in complexes with specific and methylated DNA. These conformational differences, together with intrinsic variations in DNA flexibility and/or solvation at TpG versus mCpG sites, contribute to the different DNA affinity and sequence specificity. With methylated DNA, multiple direct contacts between E535 and the 5' mCpG site dominate the binding affinity, allowing for tolerance of different flanking DNA sequences. With KBS, Kaiso employs E535 as part of an indirect screen of the 5' flanking sequence, relying on key tyrosine-DNA interactions to stabilize an optimal DNA conformation and select against noncognate sites. These findings demonstrate how TFs use conformational adaptation and exploit variations in DNA flexibility to achieve distinct DNA readout outcomes and target a greater variety of regulatory and epigenetic sites than previously appreciated.

Energy Landscape and Pathways for Transitions between Watson-Crick and Hoogsteen Base Pairing in DNA

The recent discovery that Hoogsteen (HG) base pairs are widespread in DNA across diverse sequences and positional contexts could have important implications for understanding DNA replication and DNA-protein recognition. While evidence is emerging that the Hoogsteen conformation could be a thermodynamically accessible conformation of the DNA duplex and provide a means to expand its functionality, relatively little is known about the molecular mechanism underlying the Watson-Crick (WC) to HG transition. In this Perspective, we describe pathways and kinetics for this transition at an atomic level of detail, using the energy landscape perspective. We show that competition between the duplex conformations results in a double funnel landscape, which explains some recent experimental observations. The interconversion pathways feature a number of intermediates, with a variable number of WC and HG base pairs. The relatively slow kinetics, with possible deviations from two-state behavior, suggest that this conformational switch is likely to be a challenging target for both simulation and experiment.

Structural basis of interstrand cross-link repair by O6-alkylguanine DNA alkyltransferase

Conformation of the alkylene lesion may play a role in interstrand cross-link repair by O6-alkylguanine DNA alkyltransferases.

DNA polymorphism in crystals: three stable conformations for the decadeoxynucleotide d(GCATGCATGC)

High-resolution structures of DNA fragments determined using X-ray crystallography or NMR have provided descriptions of a veritable alphabet of conformations. They have also shown that DNA is a flexible molecule, with some sequences capable of adopting two different structures. Here, the first example is presented of a DNA fragment that can assume three different and distinct conformations in crystals. The decanucleotide d(GCATGCATGC) was previously reported to assume a single-stranded double-fold structure. In one of the two crystal structures described here the decamer assumes both the double-fold conformation and, simultaneously, the more conventional B-type double-helical structure. In the other crystal the sequence assumes the A-type double-helical conformation. These results, taken together with CD spectra, which were recorded as the decamer was titrated against four metal ions and spermine, indicate that the molecule may exist as a mixed population of structures in solution. Small differences in the environmental conditions, such as the concentration of metal ion, may decide which of these crystallizes out. The results also support the idea that it may be possible for DNA to change its structure to suit the binding requirements of proteins or drugs.

μABC: a systematic microsecond molecular dynamics study of tetranucleotide sequence effects in B-DNA

We present the results of microsecond molecular dynamics simulations carried out by the ABC group of laboratories on a set of B-DNA oligomers containing the 136 distinct tetranucleotide base sequences. We demonstrate that the resulting trajectories have extensively sampled the conformational space accessible to B-DNA at room temperature. We confirm that base sequence effects depend strongly not only on the specific base pair step, but also on the specific base pairs that flank each step. Beyond sequence effects on average helical parameters and conformational fluctuations, we also identify tetranucleotide sequences that oscillate between several distinct conformational substates. By analyzing the conformation of the phosphodiester backbones, it is possible to understand for which sequences these substates will arise, and what impact they will have on specific helical parameters.

Symmetric curvature descriptors for label-free analysis of DNA

High-resolution microscopy techniques such as electron microscopy, scanning tunnelling microscopy and atomic force microscopy represent well-established, powerful tools for the structural characterization of adsorbed DNA molecules at the nanoscale. Notably, the analysis of DNA contours allows mapping intrinsic curvature and flexibility along the molecular backbone. This is particularly suited to address the impact of the base-pairs sequence on the local conformation of the strands and plays a pivotal role for investigations relating the inherent DNA shape and flexibility to other functional properties. Here, we introduce novel chain descriptors aimed to characterize the local intrinsic curvature and flexibility of adsorbed DNA molecules with unknown orientation. They consist of stochastic functions that couple the curvatures of two nanosized segments, symmetrically placed on the DNA contour. We show that the fine mapping of the ensemble-averaged functions along the molecular backbone generates characteristic patterns of variation that highlight all pairs of tracts with large intrinsic curvature or enhanced flexibility. We demonstrate the practical applicability of the method for DNA chains imaged by atomic force microscopy. Our approach paves the way for the label-free comparative analysis of duplexes, aimed to detect nanoscale conformational changes of physical or biological relevance in large sample numbers.

Sequence dependent free energy profiles of localized B- to A-form transition of DNA in water

DNA carries an inherent polymorphism, which surfaces under various external conditions. While B-form remains predominant under normal physiological conditions for most of the DNA sequences, low humidity and increased ion concentration cause B- to A-form transition. Certain proteins and molecules also sometimes cause local deformation of the DNA to the specific A-form. Previous experimental and computational studies focused on the overall B- to A-form transition. Here for the first time we investigated thermodynamics and mechanism of B- to A-form transition in water for various DNA sequences at a local dinucleotide base pair level. We introduced a new reaction coordinate Zp', based on the unique order parameter Zp, to drive B- to A-form transition locally and thereby calculate free energy profiles for the same for all the ten different dinucleotide steps embedded in a twelve base pair DNA. Results show that the trend of "A" and "B" philicity observed in experiment is preserved even at this local dinucleotide level, indicating its localized origin. Higher free energy cost obtained here is attributed to the cost of creating B∕A junctions along with formation of B->A transition at dimer level. We find that while water energetically stabilizes A-form for all the ten different dinucleotide steps to various extents, entropy acts against it. Therefore, we find that the stability of B-form DNA in water is entropic in origin. Mechanism of the conversion appears to be triggered by Slide; however, backbone parameters change concertedly.

Bioinformatic analysis of the protein/DNA interface

To investigate the principles driving recognition between proteins and DNA, we analyzed more than thousand crystal structures of protein/DNA complexes. We classified protein and DNA conformations by structural alphabets, protein blocks [de Brevern, Etchebest and Hazout (2000) (Bayesian probabilistic approach for predicting backbone structures in terms of protein blocks. Prots. Struct. Funct. Genet., 41:271-287)] and dinucleotide conformers [Svozil, Kalina, Omelka and Schneider (2008) (DNA conformations and their sequence preferences. Nucleic Acids Res., 36:3690-3706)], respectively. Assembling the mutually interacting protein blocks and dinucleotide conformers into 'interaction matrices' revealed their correlations and conformer preferences at the interface relative to their occurrence outside the interface. The analyzed data demonstrated important differences between complexes of various types of proteins such as transcription factors and nucleases, distinct interaction patterns for the DNA minor groove relative to the major groove and phosphate and importance of water-mediated contacts. Water molecules mediate proportionally the largest number of contacts in the minor groove and form the largest proportion of contacts in complexes of transcription factors. The generally known induction of A-DNA forms by complexation was more accurately attributed to A-like and intermediate A/B conformers rare in naked DNA molecules.

Cations form sequence selective motifs within DNA grooves via a combination of cation-pi and ion-dipole/hydrogen bond interactions

The fine conformational subtleties of DNA structure modulate many fundamental cellular processes including gene activation/repression, cellular division, and DNA repair. Most of these cellular processes rely on the conformational heterogeneity of specific DNA sequences. Factors including those structural characteristics inherent in the particular base sequence as well as those induced through interaction with solvent components combine to produce fine DNA structural variation including helical flexibility and conformation. Cation-pi interactions between solvent cations or their first hydration shell waters and the faces of DNA bases form sequence selectively and contribute to DNA structural heterogeneity. In this paper, we detect and characterize the binding patterns found in cation-pi interactions between solvent cations and DNA bases in a set of high resolution x-ray crystal structures. Specifically, we found that monovalent cations (Tl⁺) and the polarized first hydration shell waters of divalent cations (Mg²⁺, Ca²⁺) form cation-pi interactions with DNA bases stabilizing unstacked conformations. When these cation-pi interactions are combined with electrostatic interactions a pattern of specific binding motifs is formed within the grooves.

DNA phosphate crowding correlates with protein cationic side chain density and helical curvature in protein/DNA crystal structures

Sequence-specific binding of proteins to their DNA targets involves a complex spectrum of processes that often induce DNA conformational variation in the bound complex. The forces imposed by protein binding that cause the helical deformations are intimately interrelated and difficult to parse or rank in importance. To investigate the role of electrostatics in helical deformation, we quantified the relationship between protein cationic residue density (Cpc) and DNA phosphate crowding (Cpp). The correlation between Cpc and Cpp was then calculated for a subset of 58 high resolution protein-DNA crystal structures. Those structures containing strong Cpc/Cpp correlation (>±0.25) were likely to contain DNA helical curvature. Further, the correlation factor sign predicted the direction of helical curvature with positive (16 structures) and negative (seven structures) correlation containing concave (DNA curved toward protein) and convex (DNA curved away from protein) curvature, respectively. Protein-DNA complexes without significant Cpc/Cpp (36 structures) correlation (-0.25<0<0.25) tended to contain DNA without significant curvature. Interestingly, concave and convex complexes also include more arginine and lysine phosphate contacts, respectively, whereas linear complexes included essentially equivalent numbers of Lys/Arg phosphate contacts. Together, these findings suggest an important role for electrostatic interactions in protein-DNA complexes involving helical curvature.

Label-free, atomic force microscopy-based mapping of DNA intrinsic curvature for the nanoscale comparative analysis of bent duplexes

We propose a method for the characterization of the local intrinsic curvature of adsorbed DNA molecules. It relies on a novel statistical chain descriptor, namely the ensemble averaged product of curvatures for two nanosized segments, symmetrically placed on the contour of atomic force microscopy imaged chains. We demonstrate by theoretical arguments and experimental investigation of representative samples that the fine mapping of the average product along the molecular backbone generates a characteristic pattern of variation that effectively highlights all pairs of DNA tracts with large intrinsic curvature. The centrosymmetric character of the chain descriptor enables targetting strands with unknown orientation. This overcomes a remarkable limitation of the current experimental strategies that estimate curvature maps solely from the trajectories of end-labeled molecules or palindromes. As a consequence our approach paves the way for a reliable, unbiased, label-free comparative analysis of bent duplexes, aimed to detect local conformational changes of physical or biological relevance in large sample numbers. Notably, such an assay is virtually inaccessible to the automated intrinsic curvature computation algorithms proposed so far. We foresee several challenging applications, including the validation of DNA adsorption and bending models by experiments and the discrimination of specimens for genetic screening purposes.

ParaDock: a flexible non-specific DNA--rigid protein docking algorithm

Accurate prediction of protein–DNA complexes could provide an important stepping stone towards a thorough comprehension of vital intracellular processes. Few attempts were made to tackle this issue, focusing on binding patch prediction, protein function classification and distance constraints-based docking. We introduce ParaDock: a novel ab initio protein–DNA docking algorithm. ParaDock combines short DNA fragments, which have been rigidly docked to the protein based on geometric complementarity, to create bent planar DNA molecules of arbitrary sequence. Our algorithm was tested on the bound and unbound targets of a protein–DNA benchmark comprised of 47 complexes. With neither addressing protein flexibility, nor applying any refinement procedure, CAPRI acceptable solutions were obtained among the 10 top ranked hypotheses in 83% of the bound complexes, and 70% of the unbound. Without requiring prior knowledge of DNA length and sequence, and within <2 h per target on a standard 2.0 GHz single processor CPU, ParaDock offers a fast ab initio docking solution.

Genome-wide analysis of mobile genetic element insertion sites

Mobile genetic elements (MGEs) account for a significant fraction of eukaryotic genomes and are implicated in altered gene expression and disease. We present an efficient computational protocol for MGE insertion site analysis. ELAN, the suite of tools described here uses standard techniques to identify different MGEs and their distribution on the genome. One component, DNASCANNER analyses known insertion sites of MGEs for the presence of signals that are based on a combination of local physical and chemical properties. ISF (insertion site finder) is a machine-learning tool that incorporates information derived from DNASCANNER. ISF permits classification of a given DNA sequence as a potential insertion site or not, using a support vector machine. We have studied the genomes of Homo sapiens, Mus musculus, Drosophila melanogaster and Entamoeba histolytica via a protocol whereby DNASCANNER is used to identify a common set of statistically important signals flanking the insertion sites in the various genomes. These are used in ISF for insertion site prediction, and the current accuracy of the tool is over 65%. We find similar signals at gene boundaries and splice sites. Together, these data are suggestive of a common insertion mechanism that operates in a variety of eukaryotes.

Structural and thermodynamic signatures of DNA recognition by Mycobacterium tuberculosis DnaA

An essential protein, DnaA, binds to 9-bp DNA sites within the origin of replication oriC. These binding events are prerequisite to forming an enigmatic nucleoprotein scaffold that initiates replication. The number, sequences, positions, and orientations of these short DNA sites, or DnaA boxes, within the oriCs of different bacteria vary considerably. To investigate features of DnaA boxes that are important for binding Mycobacterium tuberculosis DnaA (MtDnaA), we have determined the crystal structures of the DNA binding domain (DBD) of MtDnaA bound to a cognate MtDnaA-box (at 2.0 Å resolution) and to a consensus Escherichia coli DnaA-box (at 2.3 Å). These structures, complemented by calorimetric equilibrium binding studies of MtDnaA DBD in a series of DnaA-box variants, reveal the main determinants of DNA recognition and establish the [T/C][T/A][G/A]TCCACA sequence as a high-affinity MtDnaA-box. Bioinformatic and calorimetric analyses indicate that DnaA-box sequences in mycobacterial oriCs generally differ from the optimal binding sequence. This sequence variation occurs commonly at the first 2 bp, making an in vivo mycobacterial DnaA-box effectively a 7-mer and not a 9-mer. We demonstrate that the decrease in the affinity of these MtDnaA-box variants for MtDnaA DBD relative to that of the highest-affinity box TTGTCCACA is less than 10-fold. The understanding of DnaA-box recognition by MtDnaA and E. coli DnaA enables one to map DnaA-box sequences in the genomes of M. tuberculosis and other eubacteria.

Nuclear magnetic resonance analysis of protein-DNA interactions

Recent methodological and instrumental advances in solution-state nuclear magnetic resonance have opened up the way to investigating challenging problems in structural biology such as large macromolecular complexes. This review focuses on the experimental strategies currently employed to solve structures of protein-DNA complexes and to analyse their dynamics. It highlights how these approaches can help in understanding detailed molecular mechanisms of target recognition.

Molecular dynamics of DNA: comparison of force fields and terminal nucleotide definitions

Despite DNA being a very important target for several proteins and drugs, molecular dynamics simulations with nucleic acids still encompass many challenges, such as the reliability of the chosen force field. In this paper, we carried out molecular dynamics simulations of the Dickerson-Drew dodecamer comparing GROMOS 53A6 and AMBER 03 force fields. While the AMBER force field presents specific topologies for the 5' and 3' terminal nucleotides, the GROMOS force field considers all nucleotides in the same way. To investigate the effects of the terminal nucleotide definitions, both force fields were modified to be applied in the two possible ways: with or without specific terminal nucleotide topologies. The analysis of global stability (rmsd, number of base pairs and hydrogen bonds) showed that both systems simulated with AMBER were stable, while the system simulated with the original GROMOS topologies was very unstable after 5 ns. When specific terminal topologies were included for GROMOS force field, DNA denaturation was delayed until 15 ns, but not avoided. The alpha/gamma transitions also displayed a strong dependence on the force field, but not on the terminal nucleotide definitions: AMBER simulations mainly explored configurations corresponding to the global minimum, while GROMOS simulations exhibited, very early in the simulations, an extensive sampling of local minima that may facilitate transitions to A-DNA isoform. The epsilon/zeta sampling was dependent both on the force field and on the terminal nucleotide definitions: while the AMBER simulations displayed well-defined B-I --> B-II transitions, the GROMOS force field clearly favored the B-I conformation. Also, the system simulated with the original GROMOS topologies displayed uncoupled epsilon/zeta transitions, leading to noncanonical conformations, but this was reverted when the new terminal nucleotide topologies were applied. Finally, the GROMOS force field leads to strong geometrical deformations on the DNA (overestimated groove widths and roll and strongly underestimated twist and slide), which restrict the use of GROMOS force field in long time scale DNA simulations unless a further reparametrization is made.

Energetics, conformation, and recognition of DNA duplexes modified by monodentate Ru(II) complexes containing terphenyl arenes

We studied the thermodynamic properties, conformation, and recognition of DNA duplexes site-specifically modified by monofunctional adducts of Ru(II) complexes of the type [Ru(II)(eta(6)-arene)(Cl)(en)](+), in which arene=para-, meta-, or ortho-terphenyl (complexes 1, 2, and 3, respectively) and en=1,2-diaminoethane. It has been shown (J. Med. Chem. 2008, 51, 5310) that 1 exhibits promising cytotoxic effects in human tumor cells, whereas 2 and 3 are much less cytotoxic; concomitantly with the high cytotoxicity of 1, its DNA binding mode involves combined intercalative and monofunctional (coordination) binding modes, whereas less cytotoxic compounds 2 and 3 bind to DNA only through a monofunctional coordination to DNA bases. An analysis of conformational distortions induced in DNA by adducts of 1 and 2 revealed more extensive and stronger distortion and concomitantly greater thermodynamic destabilization of DNA by the adducts of nonintercalating 2. Moreover, affinity of replication protein A to the DNA duplex containing adduct of 1 was pronouncedly lower than to the adduct of 2. On the other hand, another damaged-DNA-binding protein, xeroderma pigmentosum protein A, did not recognize the DNA adduct of 1 or 2. Importantly, the adducts of 1 induced a considerably lower level of repair synthesis than the adducts of 2, which suggests enhanced persistence of the adducts of the more potent and intercalating 1 in comparison with the adducts of the less potent and nonintercalating 2. Also interestingly, the adducts of 1 inhibited DNA polymerization more efficiently than the adducts of 2, and they could also be bypassed by DNA polymerases with greater difficulty. Results of the present work along with those previously published support the view that monodentate Ru(II) arene complexes belong to a class of anticancer agents for which structure-pharmacological relationships might be correlated with their DNA-binding modes.

Structural insight into the sequence dependence of nucleosome positioning

Nucleosome positioning displays sequence dependency and contributes to genomic regulation in a site-specific manner. We solved the structures of nucleosome core particle composed of strong positioning TTTAA elements flanking the nucleosome center. The positioning strength of the super flexible TA dinucleotide is consistent with its observed central location within minor groove inward regions, where it can contribute maximally to energetically challenging minor groove bending, kinking and compression. The marked preference for TTTAA and positioning power of the site 1.5 double helix turns from the nucleosome center relates to a unique histone protein motif at this location, which enforces a sustained, extremely narrow minor groove via a hydrophobic "sugar clamp." Our analysis sheds light on the basis of nucleosome positioning and indicates that the histone octamer has evolved not to fully minimize sequence discrimination in DNA binding.

DNA deformations near charged surfaces: electron and atomic force microscopy views

DNA is a very important cell structural element, which determines the level of expression of genes by virtue of its interaction with regulatory proteins. We use electron (EM) and atomic force microscopy (AFM) to characterize the flexibility of double-stranded DNA ( approximately 150-950 nm long) close to a charged surface. Automated procedures for the extraction of DNA contours ( approximately 10-120 nm for EM data and approximately 10-300 nm for AFM data) combined with new statistical chain descriptors indicate a uniquely two-dimensional equilibration of the molecules on the substrate surface regardless of the procedure of molecule mounting. However, in contrast to AFM, the EM mounting leads to a noticeable decrease in DNA persistence length together with decreased kurtosis. Analysis of local bending on short length scales (down to 6 nm in the EM study) shows that DNA flexibility behaves as predicted by the wormlike chain model. We therefore argue that adhesion of DNA to a charged surface may lead to additional static bending (kinking) of approximately 5 degrees per dinucleotide step without impairing the dynamic behavior of the DNA backbone. Implications of this finding are discussed.

A systematic molecular dynamics study of nearest-neighbor effects on base pair and base pair step conformations and fluctuations in B-DNA

It is well recognized that base sequence exerts a significant influence on the properties of DNA and plays a significant role in protein–DNA interactions vital for cellular processes. Understanding and predicting base sequence effects requires an extensive structural and dynamic dataset which is currently unavailable from experiment. A consortium of laboratories was consequently formed to obtain this information using molecular simulations. This article describes results providing information not only on all 10 unique base pair steps, but also on all possible nearest-neighbor effects on these steps. These results are derived from simulations of 50–100 ns on 39 different DNA oligomers in explicit solvent and using a physiological salt concentration. We demonstrate that the simulations are converged in terms of helical and backbone parameters. The results show that nearest-neighbor effects on base pair steps are very significant, implying that dinucleotide models are insufficient for predicting sequence-dependent behavior. Flanking base sequences can notably lead to base pair step parameters in dynamic equilibrium between two conformational sub-states. Although this study only provides limited data on next-nearest-neighbor effects, we suggest that such effects should be analyzed before attempting to predict the sequence-dependent behavior of DNA.

Ets-1 p51 and p42 isoforms differentially modulate Stromelysin-1 promoter according to induced DNA bend orientation

The Stromelysin-1 gene promoter contains a palindrome of two Ets-binding sites (EBS) that bind the p51 and p42 isoforms of the human Ets-1-transcription factor. A previous study established that full gene transactivation is associated with a ternary complex consisting of two p51 bound to the two EBS on the promoter. p42, only able to bind one of the two EBS, induces only very weak activity. Here, we investigate the mechanism by which the Stromelysin-1 promoter discriminates between p51 and p42. The differential stoichiometry of the two Ets-1 isoforms arises from the Stromelysin-1 EBS palindrome. The ternary complex requires the presence of two inhibitory domains flanking the DNA-binding domain and the ability to form an intramolecular autoinhibition module. Most importantly, the p51-ternary and the p42-binary complexes induce DNA curvatures with opposite orientations. These results establish that differential DNA bending, via p51 and p42 differential binding, is correlated with the Stromelysin-1 promoter activation process.

Small local variations in B-form DNA lead to a large variety of global geometries which can accommodate most DNA-binding protein motifs

BACKGROUND

An important question of biological relevance is the polymorphism of the double-helical DNA structure in its free form, and the changes that it undergoes upon protein-binding. We have analysed a database of free DNA crystal structures to assess the inherent variability of the free DNA structure and have compared it with a database of protein-bound DNA crystal structures to ascertain the protein-induced variations.

RESULTS

Most of the dinucleotide steps in free DNA display high flexibility, assuming different conformations in a sequence-dependent fashion. With the exception of the AA/TT and GA/TC steps, which are 'A-phobic', and the GG/CC step, which is 'A-philic', the dinucleotide steps show no preference for A or B forms of DNA. Protein-bound DNA adopts the B-conformation most often. However, in certain cases, protein-binding causes the DNA backbone to take up energetically unfavourable conformations. At the gross structural level, several protein-bound DNA duplexes are observed to assume a curved conformation in the absence of any large distortions, indicating that a series of normal structural parameters at the dinucleotide and trinucleotide level, similar to the ones in free B-DNA, can give rise to curvature at the overall level.

CONCLUSION

The results illustrate that the free DNA molecule, even in the crystalline state, samples a large amount of conformational space, encompassing both the A and the B-forms, in the absence of any large ligands. A-form as well as some non-A, non-B, distorted geometries are observed for a small number of dinucleotide steps in DNA structures bound to the proteins belonging to a few specific families. However, for most of the bound DNA structures, across a wide variety of protein families, the average step parameters for various dinucleotide sequences as well as backbone torsion angles are observed to be quite close to the free 'B-like' DNA oligomer values, highlighting the flexibility and biological significance of this structural form.

Local and global effects of strong DNA bending induced during molecular dynamics simulations

DNA bending plays an important role in many biological processes, but its molecular and energetic details as a function of base sequence remain to be fully understood. Using a recently developed restraint, we have studied the controlled bending of four different B-DNA oligomers using molecular dynamics simulations. Umbrella sampling with the AMBER program and the recent parmbsc0 force field yield free energy curves for bending. Bending 15-base pair oligomers by 90 degrees requires roughly 5 kcal mol(-1), while reaching 150 degrees requires of the order of 12 kcal mol(-1). Moderate bending occurs mainly through coupled base pair step rolls. Strong bending generally leads to local kinks. The kinks we observe all involve two consecutive base pair steps, with disruption of the central base pair (termed Type II kinks in earlier work). A detailed analysis of each oligomer shows that the free energy of bending only varies quadratically with the bending angle for moderate bending. Beyond this point, in agreement with recent experiments, the variation becomes linear. An harmonic analysis of each base step yields force constants that not only vary with sequence, but also with the degree of bending. Both these observations suggest that DNA is mechanically more complex than simple elastic rod models would imply.

Structure-based prediction of transcription factor binding sites using a protein-DNA docking approach

Accurate identification of transcription factor binding sites is critical to our understanding of transcriptional regulatory networks. To overcome the issue of high false-positive predictions that trouble the sequence-based prediction techniques, we have developed a structure-based prediction method that takes into consideration of interactions between the amino acids of a transcription factor and the nucleotides of its DNA binding sequence at structural level, along with an efficient protein-DNA docking algorithm. The docked structures between a protein and a DNA are evaluated using a knowledge-based energy function, in conjunction with van der Waals energy. Our docking algorithm supports quasi-flexible docking, overcoming a number of limiting issues faced by similar docking algorithms. Our rigid-body docking algorithm is tested on a dataset of 141 nonredundant transcription factor-DNA complex structures. The test results show that 63.1% of the 141 complex structures are reconstructed with accuracies better than 1.0 A RMSDs (root mean square deviation) and 79.4% of the complexes are predicted with accuracies better than 3.0 A RMSDs when using the native DNA structures. Our quasi-flexible docking algorithm, assuming that the DNA structures are not known, is tested on a separate set of 45 transcription factor-DNA complexes, of which 57.8% of the docked complex conformations achieve better than 1.0 A RMSDs while 71.1% of the complexes have RMSDs less than 3.0 A. We have also applied our method to predict the binding motifs of the ferric uptake regulator in E. coli and showed that most of the experimentally identified sites can be predicted accurately.

Crystallographic analysis of a sex-specific enhancer element: sequence-dependent DNA structure, hydration, and dynamics

The crystal structure of a sex-specific enhancer element is described at a resolution of 1.6 A. This 16-bp site, designated Dsx(A), functions in the regulation of a genetic switch between male and female patterns of gene expression in Drosophila melanogaster. Related sites are broadly conserved in metazoans, including in the human genome. This enhancer element is unusually rich in general regulatory sequences related to DNA recognition by multiple classes of eukaryotic transcription factors, including the DM motifs, homeodomain, and high mobility group box. Whereas free DNA is often crystallized as an A-form double helix, Dsx(A) was crystallized as B-DNA and thus provides a model for the prebound conformation of diverse regulatory DNA complexes. Sequence-dependent conformational properties that extend features of shorter B-DNA fragments with respect to double helical parameters, groove widths, hydration, and binding of divalent metal ions are observed. The structure also exhibits a sequence-dependent pattern of isotropic thermal B-factors, suggesting possible variation in the local flexibility of the DNA backbone. Such fluctuations are in accord with structural variability observed in prior B-DNA structures. We speculate that sites of intrinsic flexibility within a DNA control element provide hinges for its protein-directed reorganization in a transcriptional preinitiation complex.

DNA conformations and their sequence preferences

The geometry of the phosphodiester backbone was analyzed for 7739 dinucleotides from 447 selected crystal structures of naked and complexed DNA. Ten torsion angles of a near-dinucleotide unit have been studied by combining Fourier averaging and clustering. Besides the known variants of the A-, B- and Z-DNA forms, we have also identified combined A + B backbone-deformed conformers, e.g. with alpha/gamma switches, and a few conformers with a syn orientation of bases occurring e.g. in G-quadruplex structures. A plethora of A- and B-like conformers show a close relationship between the A- and B-form double helices. A comparison of the populations of the conformers occurring in naked and complexed DNA has revealed a significant broadening of the DNA conformational space in the complexes, but the conformers still remain within the limits defined by the A- and B- forms. Possible sequence preferences, important for sequence-dependent recognition, have been assessed for the main A and B conformers by means of statistical goodness-of-fit tests. The structural properties of the backbone in quadruplexes, junctions and histone-core particles are discussed in further detail.

Magnitude and direction of DNA bending induced by screw-axis orientation: influence of sequence, mismatches and abasic sites

DNA-bending flexibility is central for its many biological functions. A new bending restraining method for use in molecular mechanics calculations and molecular dynamics simulations was developed. It is based on an average screw rotation axis definition for DNA segments and allows inducing continuous and smooth bending deformations of a DNA oligonucleotide. In addition to controlling the magnitude of induced bending it is also possible to control the bending direction so that the calculation of a complete (2-dimensional) directional DNA-bending map is now possible. The method was applied to several DNA oligonucleotides including A(adenine)-tract containing sequences known to form stable bent structures and to DNA containing mismatches or an abasic site. In case of G:A and C:C mismatches a greater variety of conformations bent in various directions compared to regular B-DNA was found. For comparison, a molecular dynamics implementation of the approach was also applied to calculate the free energy change associated with bending of A-tract containing DNA, including deformations significantly beyond the optimal curvature. Good agreement with available experimental data was obtained offering an atomic level explanation for stable bending of A-tract containing DNA molecules. The DNA-bending persistence length estimated from the explicit solvent simulations is also in good agreement with experiment whereas the adiabatic mapping calculations with a GB solvent model predict a bending rigidity roughly two times larger.

DNA sequence-specific recognition by a transcriptional regulator requires indirect readout of A-tracts

The bacteriophage Ø29 transcriptional regulator p4 binds to promoters of different intrinsic activities. The p4–DNA complex contains two identical protomers that make similar interactions with the target sequence 5′-AACTTTTT-15 bp-AAAATGTT-3′. To define how the various elements in the target sequence contribute to p4's affinity, we studied p4 binding to a series of mutated binding sites. The binding specificity depends critically on base pairs of the target sequence through both direct as well as indirect readout. There is only one specific contact between a base and an amino acid residue; other contacts take place with the phosphate backbone. Alteration of direct amino acid–base contacts, or mutation of non-contacted A·T base pairs at A-tracts abolished binding. We generated three 5 ns molecular dynamics (MD) simulations to investigate the basis for the p4–DNA complex specificity. Recognition is controlled by the protein and depends on DNA dynamic properties. MD results on protein–DNA contacts and the divergence of p4 affinity to modified binding sites reveal an inherent asymmetry, which is required for p4-specific binding and may be crucial for transcription regulation.

The transcriptional repressor TtgV recognizes a complex operator as a tetramer and induces convex DNA bending

The TtgV repressor belongs to the large but infrequently investigated IclR family of transcriptional regulators. Although members of this family usually exhibit high effector specificity, TtgV possesses multidrug binding properties. The TtgV protein regulates the expression of the ttgGHI operon encoding the main solvent extrusion pump of the extremophile Pseudomonas putida DOT-T1E strain. Here we used a multidisciplinary approach to study the functional oligomeric state of TtgV during repression and derepression events, as well as the molecular basis of TtgV-DNA operator interactions. Analytical ultracentrifugation studies (AUC) show that TtgV is a tetramer in solution and that this oligomeric state does not change in the presence of effectors. We also show that the binding of effectors leads to the dissociation of TtgV as a tetramer from the DNA-TtgV complex. Previous dimethyl sulfate and DNase I footprints revealed that TtgV protected a 42 bp region. Based on AUC, electrophorectic mobility shift assays and isothermal titration calorimetry analyses we show that TtgV recognition specificity is restricted within this operator to a 34-nucleotide stretch and that TtgV may interact with intercalated inverted repeats that share no significant DNA sequence similarities within this short 34-nucleotide segment. Binding stoichiometry is one TtgV tetramer per operator, and affinity for its target DNA is around 200 nM. Circular dichroism analysis reveals that TtgV binding causes DNA distortion and atomic force microscopy imaging of TtgV-DNA operator complexes shows that TtgV induces a 57 degrees convex bend in its operator DNA. We propose that the mechanism of TtgV repression is based on the steric occlusion of the RNA polymerase binding site reinforced by DNA-bending of the ttgV-ttgG promoter region.

DNA stretching and extreme kinking in the nucleosome core

DNA stretching in chromatin may facilitate its compaction and influence site recognition by nuclear factors. In vivo, stretching has been estimated to occur at the equivalent of one to two base-pairs (bp) per nucleosome. We have determined the crystal structure of a nucleosome core particle containing 145 bp of DNA (NCP145). Compared to the structure with 147 bp, the NCP145 displays two incidences of stretching one to two double-helical turns from the particle dyad axis. The stretching illustrates clearly a mechanism for shifting DNA position by displacement of a single base-pair while maintaining nearly identical histone-DNA interactions. Increased DNA twist localized to a short section between adjacent histone-DNA binding sites advances the rotational setting, while a translational component involves DNA kinking at a flanking region that initiates elongation by unstacking bases. Furthermore, one stretched region of the NCP145 displays an extraordinary 55 degrees kink into the minor groove situated 1.5 double-helical turns from the particle dyad axis, a hot spot for gene insertion by HIV-integrase, which prefers highly distorted substrate. This suggests that nucleosome position and context within chromatin could promote extreme DNA kinking that may influence genomic processes.

Identification of insertion hot spots for non-LTR retrotransposons: computational and biochemical application to Entamoeba histolytica

The genome of the human pathogen Entamoeba histolytica contains non-long terminal repeat (LTR) retrotransposons, the EhLINEs and EhSINEs, which lack targeted insertion. We investigated the importance of local DNA structure, and sequence preference of the element-encoded endonuclease (EN) in selecting target sites for retrotransposon insertion. Pre-insertion loci were tested computationally to detect unique features based on DNA structure, thermodynamic considerations and protein interaction measures. Target sites could readily be distinguished from other genomic sites based on these criteria. The contribution of the EhLINE1-encoded EN in target site selection was investigated biochemically. The sequence-specificity of the EN was tested in vitro with a variety of mutated substrates. It was possible to assign a consensus sequence, 5′-GCATT-3′, which was efficiently nicked between A-T and T-T. The upstream G residue enhanced EN activity, possibly serving to limit retrotransposition in the A+T-rich E.histolytica genome. Mutated substrates with poor EN activity showed structural differences compared with normal substrates. Analysis of retrotransposon insertion sites from a variety of organisms showed that, in general, regions of favorable DNA structure were recognized for retrotransposition. A combination of favorable DNA structure and preferred EN nicking sequence in the vicinity of this structure may determine the genomic hotspots for retrotransposition.

A GCUA tetranucleotide loop found in the poliovirus oriL by in vivo SELEX (un)expectedly forms a YNMG-like structure: Extending the YNMG family with GYYA

The cloverleaf structure in the 5'-untranslated region of enterovirus RNA that regulates viral RNA replication contains an evolutionarily conserved YNMG tetraloop closed by a Y-G base pair. This loop is believed to interact specifically with the viral protease 3C. To further characterize the specificity of this interaction, the tetraloop and two flanking base pairs of the poliovirus RNA were randomized, and viable viral clones were obtained using in vivo SELEX. Among many different mutants with the canonical YNMG sequences to be described elsewhere, a large-plaque-forming clone contained a deviating uGCUAg sequence. The NMR structure of a small hairpin capped with uGCUAg that we present here shows that the GCUA tetraloop adopts a novel fold, which is highly similar to that of the YNMG tetraloop with common stacking properties and hydrogen-bond interactions including an unusual syn conformation of the adenosine. Thermodynamic studies show moderate stabilities of hairpins with canonical YNMG and the novel GCUA loops, which, together with the similarity of spatial structures, illustrates that the tetraloop structure itself is crucial for the RNA-protein interaction required for the viral replication. A re-evaluation of the ribosomal secondary structure database reveals a hairpin containing a GCUA loop, which covaries with YNMG and is involved in a tertiary interaction, and in the 50S ribosomal subunit from Haloarcula marismortui the structurally comparable apex of stem-loop 35a is a recognition site for protein L2. These observations show a more general occurrence and importance of the so-far unrecognized GYYA hairpin loops.

Information-driven protein-DNA docking using HADDOCK: it is a matter of flexibility

Intrinsic flexibility of DNA has hampered the development of efficient protein−DNA docking methods. In this study we extend HADDOCK (High Ambiguity Driven DOCKing) [C. Dominguez, R. Boelens and A. M. J. J. Bonvin (2003) J. Am. Chem. Soc.125, 1731–1737] to explicitly deal with DNA flexibility. HADDOCK uses non-structural experimental data to drive the docking during a rigid-body energy minimization, and semi-flexible and water refinement stages. The latter allow for flexibility of all DNA nucleotides and the residues of the protein at the predicted interface. We evaluated our approach on the monomeric repressor−DNA complexes formed by bacteriophage 434 Cro, the Escherichia coli Lac headpiece and bacteriophage P22 Arc. Starting from unbound proteins and canonical B-DNA we correctly predict the correct spatial disposition of the complexes and the specific conformation of the DNA in the published complexes. This information is subsequently used to generate a library of pre-bent and twisted DNA structures that served as input for a second docking round. The resulting top ranking solutions exhibit high similarity to the published complexes in terms of root mean square deviations, intermolecular contacts and DNA conformation. Our two-stage docking method is thus able to successfully predict protein−DNA complexes from unbound constituents using non-structural experimental data to drive the docking.

Studies of DNA dumbbells VIII. Melting analysis of DNA dumbbells with dinucleotide repeat stem sequences

Melting curves and circular dichroism spectra were measured for a number of DNA dumbbell and linear molecules containing dinucleotide repeat sequences of different lengths. To study effects of different sequences on the melting and spectroscopic properties, six DNA dumbbells whose stems contain the central sequences (AA)(10), (AC)(10), (AG)(10), (AT)(10), (GC)(10), and (GG)(10) were prepared. These represent the minimal set of 10 possible dinucleotide repeats. To study effects of dinucleotide repeat length, dumbbells with the central sequences (AG)(n), n = 5 and 20, were prepared. Control molecules, dumbbells with a random central sequence, (RN)(n), n = 5, 10, and 20, were also prepared. The central sequence of each dumbbell was flanked on both sides by the same 12 base pairs and T(4) end-loops. Melting curves were measured by optical absorbance and differential scanning calorimetry in solvents containing 25, 55, 85, and 115 mM Na(+). CD spectra were collected from 20 to 45 degrees C and [Na(+)] from 25 to 115 mM. The spectral database did not reveal any apparent temperature dependence in the pretransition region. Analysis of the melting thermodynamics evaluated as a function of Na(+) provided a means for quantitatively estimating the counterion release with melting for the different sequences. Results show a very definite sequence dependence, indicating the salt-dependent properties of duplex DNA are also sequence dependent. Linear DNA molecules containing the (AG)(n) and (RN)(n), sequences, n = 5, 10, 20, and 30, were also prepared and studied. The linear DNA molecules had the exact sequences of the dumbbell stems. That is, the central repeat sequence in each linear duplex was flanked on both sides by the same 12-bp sequence. Melting and CD studies were also performed on the linear DNA molecules. Comparison of results obtained for the same sequences in dumbbell and linear molecular environments reveals several interesting features of the interplay between sequence-dependent structural variability, sequence length, and the unconstrained (linear) or constrained (dumbbell) molecular environments.

Simulations of DNA Coiling around a Synthetic Supramolecular Cylinder That Binds in the DNA Major Groove

In this work we present the results of a molecular simulation study of the interaction between a tetracationic bis iron(II) supramolecular cylinder, [Fe2(C25H20N4)3]4+, and DNA. This supramolecular cylinder has been shown to bind in the major groove of DNA and to induce dramatic coiling of the DNA. The simulations have been designed to elucidate the interactions that lead the cylinder to target the major groove and that drive the subsequent DNA conformational changes. Three sets of multi-nanosecond simulations have been performed: one of the uncomplexed d(CCCCCTTTTTCC) d(GGAAAAAGGGGG) dodecamer; one of this DNA complexed with the cylinder molecule; and one of this DNA complexed with a neutralised version of the cylinder. Coiling of the DNA was observed in the DNA-cylinder simulations, giving insight into the molecular level nature of the supramolecular coiling observed experimentally. The cylinder charge was found not to be essential for the DNA coiling, which implies that the DNA response is moderated by the short range interactions that define the molecular shape. Cylinder charge did, however, affect the integrity of the DNA duplex, to the extent that, under some circumstances, the tetracationic cylinder induced defects in the DNA base pairing at locations adjacent to the cylinder binding site.

Protein-nucleic acid recognition: statistical analysis of atomic interactions and influence of DNA structure

We analyzed structural features of 11,038 direct atomic contacts (either electrostatic, H-bonds, hydrophobic, or other van der Waals interactions) extracted from 139 protein-DNA and 49 protein-RNA nonhomologous complexes from the Protein Data Bank (PDB). Globally, H-bonds are the most frequent interactions (approximately 50%), followed by van der Waals, hydrophobic, and electrostatic interactions. From the protein viewpoint, hydrophilic amino acids are over-represented in the interaction databases: Positively charged amino acids mainly contact nucleic acid phosphate groups but can also interact with base edges. From the nucleotide point of view, DNA and RNA behave differently: Most protein-DNA interactions involve phosphate atoms, while protein-RNA interactions involve more frequently base edge and ribose atoms. The increased participation of DNA phosphate involves H-bonds rather than salt bridges. A statistical analysis was performed to find the occurrence of amino acid-nucleotide pairs most different from chance. These pairs were analyzed individually. Finally, we studied the conformation of DNA in the interaction sites. Despite the prevalence of B-DNA in the database, our results suggest that A-DNA is favored in the interaction sites.

Probing of HIV-1 integrase/DNA interactions using novel analogs of viral DNA

The specific activity of the human immunodeficiency virus, type 1 (HIV-1), integrase on the viral long terminal repeat requires the binding of the enzyme to certain sequences located in the U3 and U5 regions at the ends of viral DNA, but the determinants of this specific DNA-protein recognition are not yet completely understood. We synthesized DNA duplexes mimicking the U5 region and containing either 2'-modified nucleosides or 1,3-propanediol insertions and studied their interactions with HIV-1 integrase, using Mn2+ or Mg2+ ions as integrase cofactors. These DNA modifications had no strong effect on integrase binding to the substrate analogs but significantly affected 3'-end processing rate. The effects of nucleoside modifications at positions 5, 6, and especially 3 strongly depended on the cationic cofactor used. These effects were much more pronounced in the presence of Mg2+ than in the presence of Mn2+. Modifications of base pairs 7-9 affected 3'-end processing equally in the presence of both ions. Adenine from the 3rd bp is thought to form at least two hydrogen bonds with integrase that are crucial for specific DNA recognition. The complementary base, thymine, is not important for integrase activity. For other positions, our results suggest that integrase recognizes a fine structure of the sugar-phosphate backbone rather than heterocyclic bases. Integrase interactions with the unprocessed strand at positions 5-8 are more important than interactions with the processed strand for specific substrate recognition. Based on our results, we suggest a model for integrase interaction with the U5 substrate.

Crystal structure of the restriction-modification system control element C.Bcll and mapping of its binding site

Protection from DNA invasion is afforded by restriction-modification systems in many bacteria. The efficiency of protection depends crucially on the relative expression levels of restriction versus methytransferase genes. This regulation is provided by a controller protein, named C protein. Studies of the Bcll system in E. coli suggest that C.Bcll functions as a negative regulator for M.Bcll expression, implying that it plays a role in defense against foreign DNA during virus infection. C.Bcll binds (Kd = 14.3 nM) to a 2-fold symmetric C box DNA sequence that overlaps with the putative -35 promoter region upstream of the bcllM and bcllC genes. The C.Bcll fold comprises five alpha helices: two helices form a helix-turn-helix motif, and the remaining three helices form the extensive dimer interface. The C.Bcll-DNA model proposed suggests that DNA bending might play an important role in gene regulation, and that Glu27 and Asp31 in C.Bcll might function critically in the regulation.

Thermodynamic properties of damaged DNA and its recognition by xeroderma pigmentosum group A protein and replication protein A

The effects of the lesions induced by single, site-specific 1,2-GG or 1,3-GTG intrastrand adducts of cis-diamminedichloroplatinum(II) formed in oligodeoxyribonucleotide duplexes on energetics of DNA were examined by means of differential scanning calorimetry. These effects were correlated with affinity of these duplexes for damaged-DNA binding-proteins XPA and RPA; this affinity was examined by gel electrophoresis. The results confirm that rigid DNA bending is the specific determinant responsible for high-affinity interactions of XPA with damaged DNA, but that an additional important factor, which affects affinity of XPA to damaged DNA, is a change of thermodynamic stability of DNA induced by the damage. In addition, the results also confirm that RPA preferentially binds to DNA distorted so that hydrogen bonds between complementary bases are interrupted. RPA also binds to non-denaturational distortions in double-helical DNA, but affinity of RPA to these distortions is insensitive to alterations of thermodynamic stability of damaged DNA.

Automatic intrinsic DNA curvature computation from AFM images

Critical information on several biological processes such as DNA-protein interactions and DNA transcription can be derived from analysis of DNA curvature. Under thermal perturbation, the curvature is composed of static and dynamic contributions, thus, can be described as the sum of intrinsic curvature and a fluctuation contribution. Without considering thermal agitations, the DNA curvature is reducible to the intrinsic component, which is a function of the DNA nucleotide sequence only. In this paper, we present an automated algorithm to determine the DNA intrinsic curvature profiles and the molecular spatial orientations in Atomic Force Microscope images. The algorithm allows to reconstruct the intrinsic curvature profile by filtering the thermal contribution. It detects fragment orientation on atomic force microscope images without labels with a percentage of correct molecular-orientation detection of 96.79% in computer-generated benchmarks, for molecules with a high curvature peak. The automated algorithm reconstructs the intrinsic curvature profile of DNA molecules with a mean square error of 3.8122 x 10(-4) rads over a profile with a central peak value of 0.196 rads, and 6.1 x 10(-3) rads over a curvature profile with two symmetric peaks of about 0.08 rads. Moreover, it correctly detects the location of the peaks in the molecules with a deviation of about 1% of molecule length.

Principal component analysis of DNA oligonucleotide structural data

The microstructure of a DNA helix is characterized by several base pair and base step parameters such as twist, rise, roll, propeller twist, etc., in addition to conformational parameters such as the backbone and the glycosidic torsion angles. Among these only a few, which are independent of all others and of each other, may be used to precisely characterize the helix. The problem however is to identify these independent parameters. We have used principal component analysis to identify a relatively small set of independent parameters, with which to characterize each DNA helix. We show that these principal components clearly discriminate between A and B DNA helical types. The calculations further suggest that the microstructure of a DNA helix is better characterized using dinucleotides.