Structure of a B-DNA dodecamer. II. Influence of base sequence on helix structure

X-ray Structure Characterization of the Selective Recognition of AT Base Pair Sequences

The rational design of small molecules that target specific DNA sequences is a promising strategy to modulate gene expression. This report focuses on a diamidinobenzimidazole compound, whose selective binding to the minor groove of AT DNA sequences holds broad significance in the molecular recognition of AT-rich human promoter sequences. The objective of this study is to provide a more detailed and systematized understanding, at an atomic level, of the molecular recognition mechanism of different AT-specific sequences by a rationally designed minor groove binder. The specialized method of X-ray crystallography was utilized to investigate how the sequence-dependent recognition properties in general, A-tract, and alternating AT sequences affect the binding of diamidinobenzimidazole in the DNA minor groove. While general and A-tract AT sequences give a narrower minor groove, the alternating AT sequences intrinsically have a wider minor groove which typically constricts upon binding. A strong and direct hydrogen bond between the N-H of the benzimidazole and an H-bond acceptor atom in the minor groove is essential for DNA recognition in all sequences described. In addition, the diamidine compound specifically utilizes an interfacial water molecule for its DNA binding. DNA complexes of AATT and AAAAAA recognition sites show that the diamidine compound can bind in two possible orientations with a preference for water-assisted hydrogen bonding at either cationic end. The complex structures of AAATTT, ATAT, ATATAT, and AAAA are bound in a singular orientation. Analysis of the helical parameters shows a minor groove expansion of about 1 Å across all the nonalternating DNA complexes. The results from this systematic approach will convey a greater understanding of the specific recognition of a diverse array of AT-rich sequences by small molecules and more insight into the design of small molecules with enhanced specificity to AT and mixed DNA sequences.

Candidate SNP Markers Significantly Altering the Affinity of TATA-Binding Protein for the Promoters of Human Hub Genes for Atherogenesis, Atherosclerosis and Atheroprotection

Atherosclerosis is a systemic disease in which focal lesions in arteries promote the build-up of lipoproteins and cholesterol they are transporting. The development of atheroma (atherogenesis) narrows blood vessels, reduces the blood supply and leads to cardiovascular diseases. According to the World Health Organization (WHO), cardiovascular diseases are the leading cause of death, which has been especially boosted since the COVID-19 pandemic. There is a variety of contributors to atherosclerosis, including lifestyle factors and genetic predisposition. Antioxidant diets and recreational exercises act as atheroprotectors and can retard atherogenesis. The search for molecular markers of atherogenesis and atheroprotection for predictive, preventive and personalized medicine appears to be the most promising direction for the study of atherosclerosis. In this work, we have analyzed 1068 human genes associated with atherogenesis, atherosclerosis and atheroprotection. The hub genes regulating these processes have been found to be the most ancient. In silico analysis of all 5112 SNPs in their promoters has revealed 330 candidate SNP markers, which statistically significantly change the affinity of the TATA-binding protein (TBP) for these promoters. These molecular markers have made us confident that natural selection acts against underexpression of the hub genes for atherogenesis, atherosclerosis and atheroprotection. At the same time, upregulation of the one for atheroprotection promotes human health.

Reactive Molecular Dynamics Simulation of the Structural Damages of the B-DNA Induced by the Oxidation/Nitration of Guanine

Reactive molecular dynamics simulations (RMD) have been carried out to investigate structural alterations of the dodecamer double-strand B-DNA due to the oxidation/nitration modifications introduced to its guanine bases, including 8-oxoguanine, 8-nitroguanine, and 5-guanidino-4-nitroimidazole, considering two distribution patterns. These modifications may arise in the case of cancer treatment using oxidative/nitrosative reactive nitrogen species as anticancer agents. Results show that these mutations affect structural characteristics of the B-DNA dodecamer in the order 8-nitroguanine > 5-guanidino-4-nitroimidazole ≫ 8-oxoguanine. For instance, the base-pair per turn for these modified B-DNA are changed respectively to 9.79, 10.88 and 10.58 from 10.51 in the native defect-free B-DNA, which is compatible with the experimental value of 10.10. In addition, these mutations allow more water molecules to diffuse into the dodecamer structure and consequently increase the possibility of the penetration of reactive and nonreactive species toward constituting nucleic base-pairs. The largest variation of the B-DNA structure is observed for the mutated B-DNA with 8-nitroguanine modifications applied to its separated CG base-pairs along the dodecamer chain. The structural changes introduced by these nitro-/oxo-modified guanine bases can be considered as a critical step in the damage of the DNA structure and alterations of its function.

Polyanion order controls liquid-to-solid phase transition in peptide/nucleic acid co-assembly

The Central Dogma highlights the mutualistic functions of protein and nucleic acid biopolymers, and this synergy appears prominently in the membraneless organelles widely distributed throughout prokaryotic and eukaryotic organisms alike. Ribonucleoprotein granules (RNPs), which are complex coacervates of RNA with proteins, are a prime example of these membranelles organelles and underly multiple essential cellular functions. Inspired by the highly dynamic character of these organelles and the recent studies that ATP both inhibits and templates phase separation of the fused in sarcoma (FUS) protein implicated in several neurodegenerative diseases, we explored the RNA templated ordering of a single motif of the Aβ peptide of Alzheimer's disease. We now know that this strong cross-β propensity motif alone assembles through a liquid-like coacervate phase that can be externally templated to form distinct supramolecular assemblies. Now we provide evidence that structured phosphates, ranging from complex structures like double stranded and quadraplex DNA to simple trimetaphosphate, differentially impact the liquid to solid phase transition necessary for paracrystalline assembly. The results from this simple model illustrate the potential of ordered environmental templates in the transition to potentially irreversible pathogenic assemblies and provides insight into the ordering dynamics necessary for creating functional synthetic polymer co-assemblies.

Synthesis, characterization, interactions with the DNA duplex dodecamer d(5'-CGCGAATTCGCG-3')2 and cytotoxicity of binuclear η6-arene-Ru(II) complexes

The novel binuclear η⁶-arene-Ru(II) complexes with the general formula {[(η⁶-cym)Ru(L)]₂(μ-BL)}(PF₆)₄, and their corresponding water soluble {[(η⁶-cym)Ru(L)]₂(μ-BL)}Cl₄, where cym = p-cymene, L = 2,2'-bipyridine (bpy) and 1,10-phenanthroline (phen), BL = 4,4'-bipyridine (BL-1), 1,2-bis(4-pyridyl)ethane (BL-2) and 1,3-bis(4-pyridyl)propane (BL-3), were synthesized and characterized. The structure of {[(η⁶-cym)Ru(phen)]₂(μ-BL-1)}(PF₆)₄ was determined by X-ray single crystal methods. The interaction of {[(η⁶-cym)Ru(phen)]₂(μ-BL-i)}Cl₄ (i = 1, 2, 3; (4), (5) and (6) correspondingly) with the DNA duplex d(5'-CGCGAATTCGCG-3')₂ was studied by means of NMR techniques and fluorescence titrations. The results show that complex (4) binds with a K_b = 12.133 × 10³ M^-1 through both intercalation and groove binding, while (5) and (6) are groove binders (K_b = 2.333 × 10³ M^-1 and K_b = 3.336 × 10³ M^-1 correspondingly). Comparison with the mononuclear complex [(η⁶-cym)Ru(phen)(py)]²⁺ reveals that it binds to the d(5'-CGCGAATTCGCG-3')₂ with a K_b value two orders of magnitude lower than (4) (K_b = 0.158 × 10³ M^-1), indicating that for the binuclear complexes both ruthenium moieties participate in the binding. The complexes were found to be cytotoxic against the A2780 and A2780 res. cancer cell line with a selectivity index (SI) in the range of 3.0-5.9.

New aspects of DNA recognition by group II WRKY transcription factor revealed by structural and functional study of AtWRKY18 DNA binding domain

WRKY transcription factors (TFs) constitute one of the largest families of plant TFs. Based on the organization of domains and motifs, WRKY TFs are divided into three Groups (I-III). The WRKY subgroup IIa includes three representatives in A. thaliana, AtWRKY18, AtWRKY40, and AtWRKY60, that participate in biotic and abiotic stress responses. Here we present crystal structures of the DNA binding domain (DBD) of AtWRKY18 alone and in the complex with a DNA duplex containing the WRKY-recognition sequence, W-box. Subgroup IIa WRKY TFs are known to form homo and heterodimers. Our data suggest that the dimerization interface of the full-length AtWRKY18 involves contacts between the DBD subunits. DNA binding experiments and structural analysis point out novel aspects of DNA recognition by WRKY TFs. In particular, AtWRKY18-DBD preferentially binds an overlapping tandem of W-boxes accompanied by a quasi-W-box motif. The binding of DNA deforms the B-type double helix, which suggests that the DNA fragment must be prone to form a specific structure. This can explain why despite the short W-box consensus, WRKY TFs can precisely control gene expression. Finally, this first experimental structure of a Group II WRKY TF allowed us to compare Group I-III representatives.

DNA Structure and the Golden Ratio Revisited

B-DNA, the informational molecule for life on earth, appears to contain ratios structured around the irrational number 1.618…, often known as the “golden ratio”. This occurs in the ratio of the length:width of one turn of the helix; the ratio of the spacing of the two helices; and in the axial structure of the molecule which has ten-fold rotational symmetry. That this occurs in the information-carrying molecule for life is unexpected, and suggests the action of some process. What this process might be is unclear, but it is central to any understanding of the formation of DNA, and so life.

Intron exon boundary junctions in human genome have in-built unique structural and energetic signals

Precise identification of correct exon–intron boundaries is a prerequisite to analyze the location and structure of genes. The existing framework for genomic signals, delineating exon and introns in a genomic segment, seems insufficient, predominantly due to poor sequence consensus as well as limitations of training on available experimental data sets. We present here a novel concept for characterizing exon–intron boundaries in genomic segments on the basis of structural and energetic properties. We analyzed boundary junctions on both sides of all the exons (3 28 368) of protein coding genes from human genome (GENCODE database) using 28 structural and three energy parameters. Study of sequence conservation at these sites shows very poor consensus. It is observed that DNA adopts a unique structural and energy state at the boundary junctions. Also, signals are somewhat different for housekeeping and tissue specific genes. Clustering of 31 parameters into four derived vectors gives some additional insights into the physical mechanisms involved in this biological process. Sites of structural and energy signals correlate well to the positions playing important roles in pre-mRNA splicing.

Single-Base Lesions and Mismatches Alter the Backbone Conformational Dynamics in DNA

DNA damage has been implicated in numerous human diseases, particularly cancer, and the aging process. Single-base lesions and mismatches in DNA can be cytotoxic or mutagenic and are recognized by a DNA glycosylase during the process of base excision repair. Altered local dynamics and conformational properties in damaged DNAs have previously been suggested to assist in recognition and specificity. Herein, we use solution nuclear magnetic resonance to quantify changes in BI-BII backbone conformational dynamics due to the presence of single-base lesions in DNA, including uracil, dihydrouracil, 1,N6-ethenoadenine, and T:G mismatches. Stepwise changes to the %BII and ΔG of the BI-BII dynamic equilibrium compared to those of unmodified sequences were observed. Additionally, the equilibrium skews toward endothermicity for the phosphates nearest the lesion/mismatched base pair. Finally, the phosphates with the greatest alterations correlate with those most relevant to the repair of enzyme binding. All of these results suggest local conformational rearrangement of the DNA backbone may play a role in lesion recognition by repair enzymes.

Beyond the double helix: DNA structural diversity and the PDB

The determination of the double helical structure of DNA in 1953 remains the landmark event in the development of modern biological and biomedical science. This structure has also been the starting point for the determination of some 2000 DNA crystal structures in the subsequent 68 years. Their structural diversity has extended to the demonstration of sequence-dependent local structure in duplex DNA, to DNA bending in short and long sequences and in the DNA wound round the nucleosome, and to left-handed duplex DNAs. Beyond the double helix itself, in circumstances where DNA sequences are or can be induced to unwind from being duplex, a wide variety of topologies and forms can exist. Quadruplex structures, based on four-stranded cores of stacked G-quartets, are prevalent though not randomly distributed in the human and other genomes and can play roles in transcription, translation, and replication. Yet more complex folds can result in DNAs with extended tertiary structures and enzymatic/catalytic activity. The Protein Data Bank is the depository of all these structures, and the resource where structures can be critically examined and validated, as well as compared one with another to facilitate analysis of conformational and base morphology features. This review will briefly survey the major structural classes of DNAs and illustrate their significance, together with some examples of how the use of the Protein Data Bank by for example, data mining, has illuminated DNA structural concepts.

Thermodynamics of DNA Hybridization from Atomistic Simulations

Studying DNA hybridization equilibrium at atomistic length scales, either via molecular dynamics (MD) or through commonly used advanced sampling approaches, is notoriously difficult. In this work, we describe an order-parameter-based advanced sampling technique to calculate the free energy of hybridization, and estimate the melting temperature of DNA oligomers at atomistic resolution. The free energy landscapes are reported as a function of a native-topology-based order parameter for the Drew-Dickerson dodecamer and for a range of DNA decamer sequences of different GC content. Our estimated melting temperatures match the experimental numbers within ±15 °C. As a test of the numerical reliability of the procedures employed, it was verified that the predicted free energy surfaces and melting temperatures of the d- and l-enantiomers of the Drew-Dickerson dodecamer were indistinguishable within numerical accuracy.

Molecular mechanism and binding free energy of doxorubicin intercalation in DNA

The intercalation process of binding doxorubicin (DOX) in DNA is studied by extensive MD simulations. Many molecular factors that control the binding affinity of DOX to DNA to form a stable complex are inspected and quantified by employing continuum solvation models for estimating the binding free energy. The modified MM-PB(GB)SA methodology provides a complete energetic profile of ΔGele, ΔGvDW, ΔGpolar, ΔGnon-polar, TΔStotal, ΔGdeform, ΔGcon, and ΔGion. To identify the sequence specificity of DOX, two different DNA sequences, d(CGATCG) or DNA1 and d(CGTACG) or DNA2, with one molecule (1 : 1 complex) or two molecule (2 : 1 complex) configurations of DOX were selected in this study. Our results show that the DNA deformation energy (ΔGdeform), the energy cost from translational and rotational entropic contributions (TΔStran+rot), the total electrostatic interactions (ΔGpolar-PB/GB + ΔGele) of incorporation, the intramolecular electrostatic interactions (ΔGele) and electrostatic polar solvation interactions (ΔGpolar-PB/GB) are all unfavorable to the binding of DOX to DNA. However, they are overcome by at least five favorable interactions: the van der Waals interactions (ΔGvDW), the non-polar solvation interaction (ΔGnon-polar), the vibrational entropic contribution (TΔSvib), and the standard concentration dependent free energies of DOX (ΔGcon) and the ionic solution (ΔGion). Specifically, the van der Waals interaction appears to be the major driving force to form a stable DOX-DNA complex. We also predict that DOX has stronger binding to DNA1 than DNA2. The DNA deformation penalty and entropy cost in the 2 : 1 complex are less than those in the 1 : 1 complex, thus they indicate that the 2 : 1 complex is more stable than the 1 : 1 complex. We have calculated the total binding free energy (BFE) (ΔGt-sim) using both MM-PBSA and MM-GBSA methods, which suggests a more stable DOX-DNA complex at lower ionic concentration. The calculated BFE from the modified MM-GBSA method for DOX-DNA1 and DOX-DNA2 in the 1 : 1 complex is -9.1 and -5.1 kcal mol-1 respectively. The same quantities from the modified MM-PBSA method are -12.74 and -8.35 kcal mol-1 respectively. The value of the total BFE ΔGt-sim in the 1 : 1 complex is in reasonable agreement with the experimental value of -7.7 ± 0.3 kcal mol-1.

New Insights into the Interactions of a DNA Oligonucleotide with mPEGylated-PAMAM by Circular Dichroism and Solution NMR

Dendrimers are well-defined, highly branched, synthetic three-dimensional molecules with a large number of reactive end groups. PAMAM dendrimers form stable complexes with DNA chemistries and constitute an important class of nonviral, cationic vectors in gene delivery. The aim of this study is to examine the interactions of a 12 bp DNA oligonucletide with PAMAM-G2 and mPEG- b-PAMAM-G3 having eight surface amine groups under physiological conditions, using constant DNA concentration but varying dendrimer concentration. 1D ³¹P NMR, 2D NOESY, and CD spectroscopic methods were employed to investigate the interactions between the dendrimer and the DNA. The CD experiments carried out with a constant DNA concentration of 25 μM and dendrimer concentrations from 0 to 100 μM indicated minimal change to the chirality of the DNA for both types of dendrimers. While the PAMAM-G2 dendrimer caused aggregation of the majority of the DNA, the 2D NMR data of the DNA with an mPEG- b-PAMAM-G3 dendrimer indicated general broadening of the 1D ³¹P peaks from the DNA phosphates, a small number of ¹H chemical shift perturbations (CSPs), and reduction of specific ¹H-¹H NOE intensities. These data suggest there is minimal structural alteration of the DNA in the complex and indicate preferential binding of the dendrimer to the central AATT region of the DNA sequence. The results herein are the first such results demonstrating a soluble DNA complex with the mPEG- b-PAMAM-G3 dendrimer analyzed by multidimensional NMR.

Branched-Chain Polyamine Found in Hyperthermophiles Induces Unique Temperature-Dependent Structural Changes in Genome-Size DNA

A pentavalent branched‐chain polyamine, N⁴‐bis(aminopropyl)spermidine 3(3)(3)4, is a unique polycation found in the hyperthermophilic archaeon Thermococcus kodakarensis, which grows at temperatures between 60 and 100 °C. We studied the effects of this branched‐chain polyamine on DNA structure at different temperatures up to 80 °C. Atomic force microscopic observation revealed that 3(3)(3)4 induces a mesh‐like structure on a large DNA (166 kbp) at 24 °C. With an increase in temperature, DNA molecules tend to unwind, and multiple nano‐loops with a diameter of 10–50 nm are generated along the DNA strand at 80 °C. These results were compared to those obtained with linear‐chain polyamines, homocaldopentamine 3334 and spermidine, the former of which is a structural isomer of 3(3)(3)4. These specific effects are expected to neatly concern with its role on high‐temperature preference in hyperthermophiles.

DFT study of geometrical and vibrational features of a 3',5'-deoxydisugar-monophosphate (dDSMP) DNA model in the presence of counterions and solvent

The B3LYP/6-31++G* theoretical level was used to study the influence of various hexahydrated monovalent (Li⁺, Na⁺, K⁺) and divalent (Mg²⁺) metal counterions in interaction with the charged PO₂^- group, on the geometrical and vibrational characteristics of the DNA fragments of 3',5'-dDSMP, represented by four conformers (g⁺g⁺, g⁺t, g^-g^- and g^-t). All complexes were optimized through two solvation models [the explicit model (6H₂O) and the hybrid model (6H₂O/Continuum)]. The results obtained established that, in the hybrid model, counterions (Li⁺, Na⁺, K⁺, Mg²⁺) always remain in the bisector plane of the O1-P-O2 angle. When these counterions are explicitly hydrated, the smallest counterions (Li⁺, Na⁺) deviate from the bisector plane, while the largest counterions (K⁺ and Mg²⁺) always remain in the same plane. On the other hand, the present calculations reveal that the g⁺g⁺ conformer is the most stable in the presence of monovalent counterions, while conformers g⁺t and g^-t are the most stable in the presence of the divalent counterion Mg²⁺. Finally, the hybrid solvation model seems to be in better agreement with the available crystallographic and spectroscopic (Raman) experiments than the explicit model. Indeed, the six conformational torsions of the C4'-C3'-O3'-PO^-₂-O5'-C5'-C4' segment of all complexes of the g^-g^- conformer in 6H₂O/Continuum remain similar to the available experimental data of A- and B-DNA forms. The calculated wavenumbers of the g⁺g⁺ conformer in the presence of the monovalent counterion and of g^-t conformer in presence of the divalent counterion in the hybrid model are in good agreement with the Raman experimental data of A- and B-DNA forms. In addition, the maximum deviation between the calculated wavenumbers in the 6H₂O/Continuum for the g⁺g⁺ conformer and experimental value measured in an aqueous solution of the DMP-Na⁺ complex, is <1.07% for the PO₂^- (asymmetric and symmetric) stretching modes and <2.03% for the O5'-C5' and O3'-C3' stretching modes. Graphical abstract dDSMP-(OO)^- Mg²⁺/6W/Continuum.

Introducing a model of pairing based on base pair specific interactions between identical DNA sequences

At present, there have been suggested two types of physical mechanism that may facilitate preferential pairing between DNA molecules, with identical or similar base pair texts, without separation of base pairs. One mechanism solely relies on base pair specific patterns of helix distortion being the same on the two molecules, discussed extensively in the past. The other mechanism proposes that there are preferential interactions between base pairs of the same composition. We introduce a model, built on this second mechanism, where both thermal stretching and twisting fluctuations are included, as well as the base pair specific helix distortions. Firstly, we consider an approximation for weak pairing interactions, or short molecules. This yields a dependence of the energy on the square root of the molecular length, which could explain recent experimental data. However, analysis suggests that this approximation is no longer valid at large DNA lengths. In a second approximation, for long molecules, we define two adaptation lengths for twisting and stretching, over which the pairing interaction can limit the accumulation of helix disorder. When the pairing interaction is sufficiently strong, both adaptation lengths are finite; however, as we reduce pairing strength, the stretching adaptation length remains finite but the torsional one becomes infinite. This second state persists to arbitrarily weak values of the pairing strength; suggesting that, if the molecules are long enough, the pairing energy scales as length. To probe differences between the two pairing mechanisms, we also construct a model of similar form. However, now, pairing between identical sequences solely relies on the intrinsic helix distortion patterns. Between the two models, we see interesting qualitative differences. We discuss our findings, and suggest new work to distinguish between the two mechanisms.

Orientation-dependent FRET system reveals differences in structures and flexibilities of nicked and gapped DNA duplexes

Differences in structures and flexibilities of DNA duplexes play important roles on recognition by DNA-binding proteins. We herein describe a novel method for structural analyses of DNA duplexes by using orientation dependence of Förster resonance energy transfer (FRET). We first analyzed canonical B-form duplex and correct structural parameters were obtained. The experimental FRET efficiencies were in excellent agreement with values theoretically calculated by using determined parameters. We then investigated DNA duplexes with nick and gaps, which are key intermediates in DNA repair systems. Effects of gap size on structures and flexibilities were successfully revealed. Since our method is facile and sensitive, it could be widely used to analyze DNA structures containing damages and non-natural molecules.

Statistical mechanical model for a closed loop plectoneme with weak helix specific forces

We develop a statistical mechanical framework, based on a variational approximation, to describe closed loop plectonemes. This framework incorporates weak helix structure dependent forces into the determination of the free energy and average structure of a plectoneme. Notably, due to their chiral nature, helix structure dependent forces break the symmetry between left and right handed supercoiling. The theoretical approach, presented here, also provides a systematic way of enforcing the topological constraint of closed loop supercoiling in the variational approximation. At large plectoneme lengths, by considering correlation functions in an expansion in terms of the spatial mean twist density about its thermally averaged value, it can be argued that topological constraint may be approximated by replacing twist and writhe by their thermal averages. A Lagrange multiplier, containing the sum of average twist and writhe, can be added to the free energy to conveniently inforce this result. The average writhe can be calculated through the thermal average of the Gauss' integral in the variational approximation. Furthermore, this approach allows for a possible way to calculate finite size corrections due to the topological constraint. Using interaction energy terms from the mean-field Kornyshev-Leikin theory, for parameter values that correspond to weak helix dependent forces, we calculate the free energy, fluctuation magnitudes and mean geometric parameters for the plectoneme. We see a slight asymmetry, where interestingly, left handed supercoils have a looser structure than right handed ones, although with a lower free energy, unlike what the previous ground state calculations would suggest.

The Sequence of Nucleosomal DNA Modulates Sliding by the Chd1 Chromatin Remodeler

Chromatin remodelers are ATP-dependent enzymes that are critical for reorganizing and repositioning nucleosomes in concert with many basic cellular processes. For the chromodomain helicase DNA-binding protein 1 (Chd1) remodeler, nucleosome sliding has been shown to depend on the DNA flanking the nucleosome, transcription factor binding at the nucleosome edge, and the presence of the histone H2A/H2B dimer on the entry side. Here, we report that Chd1 is also sensitive to the sequence of DNA within the nucleosome and slides nucleosomes made with the 601 Widom positioning sequence asymmetrically. Kinetic and equilibrium experiments show that poly(dA:dT) tracts perturb remodeling reactions if within one and a half helical turns of superhelix location 2 (SHL2), where the Chd1 ATPase engages nucleosomal DNA. These sequence-dependent effects do not rely on the Chd1 DNA-binding domain and are not due to differences in nucleosome affinity. Using site-specific cross-linking, we show that internal poly(dA:dT) tracts do not block the engagement of the ATPase motor with SHL2, yet they promote multiple translational positions of DNA with respect to both Chd1 and the histone core. We speculate that Chd1 senses the sequence-dependent response of DNA as the remodeler ATPase perturbs the duplex at SHL2. These results suggest that the sequence sensitivity of histones and remodelers occur at unique segments of DNA on the nucleosome, allowing them to work together or in opposition to determine nucleosome positions throughout the genome.

Discovery of the first dual G-triplex/G-quadruplex stabilizing compound: a new opportunity in the targeting of G-rich DNA structures?

BACKGROUND

Guanine-rich DNA motifs can form non-canonical structures known as G-quadruplexes, whose role in tumorigenic processes makes them attractive drug-target candidates for cancer therapy. Recent studies revealed that the folding and unfolding pathways of G-quadruplexes proceed through a quite stable intermediate named G-triplex.

METHODS

Virtual screening was employed to identify a small set of putative G-triplex ligands. The G-triplex stabilizing properties of these compounds were analyzed by CD melting assay. DSC, non-denaturing gel electrophoresis, NMR and molecular modeling studies were performed to investigate the interaction between the selected compound 1 and G-rich DNA structures. Cytotoxic activity of 1 was evaluated by MTT cell proliferation assay.

RESULTS

The experiments led to the identification of a promising hit that was shown to bind preferentially to G-triplex and parallel-stranded G-quadruplexes over duplex and antiparallel G-quadruplexes. Molecular modeling results suggested a partial end-stacking of 1 to the external G-triad/G-tetrads as a binding mode. Biological assays showed that 1 is endowed with cytotoxic effect on human osteosarcoma cells.

CONCLUSIONS

A tandem application of virtual screening along with the experimental investigation was employed to discover a G-triplex-targeting ligand. Experiments revealed that the selected compound actually acts as a dual G-triplex/G-quadruplex stabilizer, thus stimulating further studies aimed at its optimization.

GENERAL SIGNIFICANCE

The discovery of molecules able to bind and stabilize G-triplex structures is highly appealing, but their transient state makes challenging their recognition. These findings suggest that the identification of ligands with dual G-triplex/G-quadruplex stabilizing properties may represent a new route for the design of anticancer agents targeting the G-rich DNA structures. This article is part of a Special Issue entitled "G-quadruplex" Guest Editor: Dr. Concetta Giancola and Dr. Daniela Montesarchio.

Interactions of a biocompatible water-soluble anthracenyl polymer derivative with double-stranded DNA

We have studied the interaction of a polymeric water soluble anthracenyl derivative () with salmon testes DNA. The results from UV-Vis, fluorescence, Fourier transform infrared (FT-IR) and circular dichroism spectroscopies indicate that the groove binding process regulates the interaction between and DNA. The binding constants, calculated by absorption spectroscopy at 298, 304 and 310 K, were equal to 3.2 × 10(5) M(-1), 4.7 × 10(5) M(-1), and 6.6 × 10(5) M(-1) respectively, proving a relatively high affinity of for salmon testes DNA. Results of Hoechst 33258 displacement assays strongly support the groove binding mode of to DNA. The association stoichiometry of the :DNA adduct was found to be 1 for every 5 base pairs. FT-IR spectra, recorded at different /DNA molar ratios, indicate the involvement of the phosphate groups and adenine and thymine DNA bases in the association process. Thermodynamic results suggest that hydrophobic forces regulate the binding of with DNA without excluding some extent of involvement of van der Waals forces and hydrogen bonding arising due to surface binding between the hydrophilic polymeric arms of the ligand and the functional groups positioned on the edge of the groove. The resulting composite biomaterial could constitute a valuable candidate for future biological and/or photonic applications.

Structure elucidation of the Pribnow box consensus promoter sequence by racemic DNA crystallography

It has previously been shown that the use of racemic mixtures of naturally chiral macromolecules such as protein and DNA can significantly aid the crystallogenesis process, thereby addressing one of the major bottlenecks to structure determination by X-ray crystallographic methods-that of crystal growth. Although previous studies have provided convincing evidence of the applicability of the racemic crystallization technique to DNA through the study of well-characterized DNA structures, we sought to apply this method to a historically challenging DNA sequence. For this purpose we chose a non-self-complementary DNA duplex containing the biologically-relevant Pribnow box consensus sequence 'TATAAT'. Four racemic crystal structures of this previously un-crystallizable DNA target are reported (with resolutions in the range of 1.65-2.3 Å), with further crystallographic studies and structural analysis providing insight into the racemic crystallization process as well as structural details of this highly pertinent DNA sequence.

Evidence of protein-free homology recognition in magnetic bead force-extension experiments

Earlier theoretical studies have proposed that the homology-dependent pairing of large tracts of dsDNA may be due to physical interactions between homologous regions. Such interactions could contribute to the sequence-dependent pairing of chromosome regions that may occur in the presence or the absence of double-strand breaks. Several experiments have indicated the recognition of homologous sequences in pure electrolytic solutions without proteins. Here, we report single-molecule force experiments with a designed 60 kb long dsDNA construct; one end attached to a solid surface and the other end to a magnetic bead. The 60 kb constructs contain two 10 kb long homologous tracts oriented head to head, so that their sequences match if the two tracts fold on each other. The distance between the bead and the surface is measured as a function of the force applied to the bead. At low forces, the construct molecules extend substantially less than normal, control dsDNA, indicating the existence of preferential interaction between the homologous regions. The force increase causes no abrupt but continuous unfolding of the paired homologous regions. Simple semi-phenomenological models of the unfolding mechanics are proposed, and their predictions are compared with the data.

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.

Stacking interactions in RNA and DNA: Roll-slide energy hyperspace for ten unique dinucleotide steps

Understanding dinucleotide sequence directed structures of nuleic acids and their variability from experimental observation remained ineffective due to unavailability of statistically meaningful data. We have attempted to understand this from energy scan along twist, roll, and slide degrees of freedom which are mostly dependent on dinucleotide sequence using ab initio density functional theory. We have carried out stacking energy analysis in these dinucleotide parameter phase space for all ten unique dinucleotide steps in DNA and RNA using DFT-D by ωB97X-D/6-31G(2d,2p), which appears to satisfactorily explain conformational preferences for AU/AU step in our recent study. We show that values of roll, slide, and twist of most of the dinucleotide sequences in crystal structures fall in the low energy region. The minimum energy regions with large twist values are associated with the roll and slide values of B-DNA, whereas, smaller twist values correspond to higher stability to RNA and A-DNA like conformations. Incorporation of solvent effect by CPCM method could explain the preference shown by some sequences to occur in B-DNA or A-DNA conformations. Conformational preference of BII sub-state in B-DNA is preferentially displayed mainly by pyrimidine-purine steps and partly by purine-purine steps. The purine-pyrimidine steps show largest effect of 5-methyl group of thymine in stacking energy and the introduction of solvent reduces this effect significantly. These predicted structures and variabilities can explain the effect of sequence on DNA and RNA functionality.

Transitions of Double-Stranded DNA Between the A- and B-Forms

The structure of double-stranded DNA (dsDNA) is sensitive to solvent conditions. In solution, B-DNA is the favored conformation under physiological conditions, while A-DNA is the form found under low water activity. The A-form is induced locally in some protein-DNA complexes, and repeated transitions between the B- and A-forms have been proposed to generate the forces used to drive dsDNA into viral capsids during genome packaging. Here, we report analyses on previous molecular dynamics (MD) simulations on B-DNA, along with new MD simulations on the transition from A-DNA to B-DNA in solution. We introduce the A-B Index (ABI), a new metric along the A-B continuum, to quantify our results. When A-DNA is placed in an equilibrated solution at physiological ionic strength, there is no energy barrier to the transition to the B-form, which begins within about 1 ns. The transition is essentially complete within 5 ns, although occasionally a stretch of a few base pairs will remain A-like for up to ∼10 ns. A comparison of four sequences with a range of predicted A-phobicities shows that more A-phobic sequences make the transition more rapidly than less A-phobic sequences. Simulations on dsDNA with a region of roughly one turn locked in the A-form allow us to characterize the A/B junction, which has an average bend angle of 20-30°. Fluctuations in this angle occur with characteristic times of about 10 ns.

A historical account of Hoogsteen base-pairs in duplex DNA

In 1957, a unique pattern of hydrogen bonding between N3 and O4 on uracil and N7 and N6 on adenine was proposed to explain how poly(rU) strands can associate with poly(rA)-poly(rU) duplexes to form triplexes. Two years later, Karst Hoogsteen visualized such a noncanonical A-T base-pair through X-ray analysis of co-crystals containing 9-methyladenine and 1-methylthymine. Subsequent X-ray analyses of guanine and cytosine derivatives yielded the expected Watson-Crick base-pairing, but those of adenine and thymine (or uridine) did not yield Watson-Crick base-pairs, instead favoring "Hoogsteen" base-pairing. More than two decades ensued without experimental "proof" for A-T Watson-Crick base-pairs, while Hoogsteen base-pairs continued to surface in AT-rich sequences, closing base-pairs of apical loops, in structures of DNA bound to antibiotics and proteins, damaged and chemically modified DNA, and in polymerases that replicate DNA via Hoogsteen pairing. Recently, NMR studies have shown that base-pairs in duplex DNA exist as a dynamic equilibrium between Watson-Crick and Hoogsteen forms. There is now little doubt that Hoogsteen base-pairs exist in significant abundance in genomic DNA, where they can expand the structural and functional versatility of duplex DNA beyond that which can be achieved based only on Watson-Crick base-pairing. Here, we provide a historical account of the discovery and characterization of Hoogsteen base-pairs, hoping that this will inform future studies exploring the occurrence and functional importance of these alternative base-pairs.

Sequence-dependent structural changes in a self-assembling DNA oligonucleotide

DNA has proved to be a remarkable molecule for the construction of sophisticated two-dimensional and three-dimensional architectures because of its programmability and structural predictability provided by complementary Watson-Crick base pairing. DNA oligonucleotides can, however, exhibit a great deal of local structural diversity. DNA conformation is strongly linked to both environmental conditions and the nucleobase identities inherent in the oligonucleotide sequence, but the exact relationship between sequence and local structure is not completely understood. This study examines how a single-nucleotide addition to a class of self-assembling DNA 13-mers leads to a significantly different overall structure under identical crystallization conditions. The DNA 13-mers self-assemble in the presence of Mg(2+) through a combination of Watson-Crick and noncanonical base-pairing interactions. The crystal structures described here show that all of the predicted Watson-Crick base pairs are present, with the major difference being a significant rearrangement of noncanonical base pairs. This includes the formation of a sheared A-G base pair, a junction of strands formed from base-triple interactions, and tertiary interactions that generate structural features similar to tandem sheared G-A base pairs. The adoption of this alternate noncanonical structure is dependent in part on the sequence in the Watson-Crick duplex region. These results provide important new insights into the sequence-structure relationship of short DNA oligonucleotides and demonstrate a unique interplay between Watson-Crick and noncanonical base pairs that is responsible for crystallization fate.

DNA Structural Correlation in Short and Long Ranges

Recent single-molecule measurements have revealed the DNA allostery in protein/DNA binding. MD simulations showed that this allosteric effect is associated with the deformation properties of DNA. In this study, we used MD simulations to further investigate the mechanism of DNA structural correlation, its dependence on DNA sequence, and the chemical modification of the bases. Besides a random sequence, poly d(AT) and poly d(GC) are also used as simpler model systems, which show the different bending and twisting flexibilities. The base-stacking interactions and the methyl group on the 5-carbon site of thymine causes local structures and flexibility to be very different for the two model systems, which further lead to obviously different tendencies of the conformational deformations, including the long-range allosteric effects.

Which way up? Recognition of homologous DNA segments in parallel and antiparallel alignments

Homologous gene shuffling between DNA molecules promotes genetic diversity and is an important pathway for DNA repair. For this to occur, homologous genes need to find and recognize each other. However, despite its central role in homologous recombination, the mechanism of homology recognition has remained an unsolved puzzle of molecular biology. While specific proteins are known to play a role at later stages of recombination, an initial coarse grained recognition step has, however, been proposed. This relies on the sequence dependence of the DNA structural parameters, such as twist and rise, mediated by intermolecular interactions, in particular, electrostatic ones. In this proposed mechanism, sequences that have the same base pair text, or are homologous, have lower interaction energy than those sequences with uncorrelated base pair texts. The difference between the two energies is termed the "recognition energy." Here, we probe how the recognition energy changes when one DNA fragment slides past another, and consider, for the first time, homologous sequences in antiparallel alignment. This dependence on sliding is termed the "recognition well." We find there is a recognition well for anti-parallel, homologous DNA tracts, but only a very shallow one, so that their interaction will differ little from the interaction between two nonhomologous tracts. This fact may be utilized in single molecule experiments specially targeted to test the theory. As well as this, we test previous theoretical approximations in calculating the recognition well for parallel molecules against MC simulations and consider more rigorously the optimization of the orientations of the fragments about their long axes upon calculating these recognition energies. The more rigorous treatment affects the recognition energy a little, when the molecules are considered rigid. When torsional flexibility of the DNA molecules is introduced, we find excellent agreement between the analytical approximation and simulations.

Achiral, acyclic nucleic acids: synthesis and biophysical studies of a possible prebiotic polymer

The search for prebiotic, nucleic acid precursors is, at its best, a speculative undertaking.

Towards biochemically relevant QM computations on nucleic acids: controlled electronic structure geometry optimization of nucleic acid structural motifs using penalty restraint functions

Recent developments in dispersion-corrected density functional theory methods allow for the first time the description of large fragments of nucleic acids (hundreds of atoms) with an accuracy clearly surpassing the accuracy of common biomolecular force fields. Such calculations can significantly improve the description of the potential energy surface of nucleic acid molecules, which may be useful for studies of molecular interactions and conformational preferences of nucleic acids, as well as verification and parameterization of other methods. The first of such studies, however, demonstrated that successful applications of accurate QM calculations to larger nucleic acid building blocks are hampered by difficulties in obtaining geometries that are biochemically relevant and are not biased by non-native structural features. We present an approach that can greatly facilitate large-scale QM studies on nucleic acids, namely electronic structure geometry optimization of nucleic acid fragments utilizing a penalty function to restrain key internal coordinates with a specific focus on the torsional backbone angles. This work explores the viability of these restraint optimizations for DFT-D3, PM6-D3H and HF-3c optimizations on a set of examples (a UpA dinucleotide, a DNA G-quadruplex and a B-DNA fragment). Evaluation of different penalty function strengths reveals only a minor system-dependency and reasonable restraint values range from 0.01 to 0.05 Eh rad(-2) for the backbone torsions. Restraints are crucial to perform the QM calculations on biochemically relevant conformations in implicit solvation and gas phase geometry optimizations. The reasons for using restrained instead of constrained or unconstrained optimizations are explained and an open-source external optimizer is provided.

Type II restriction endonucleases--a historical perspective and more

This article continues the series of Surveys and Summaries on restriction endonucleases (REases) begun this year in Nucleic Acids Research. Here we discuss 'Type II' REases, the kind used for DNA analysis and cloning. We focus on their biochemistry: what they are, what they do, and how they do it. Type II REases are produced by prokaryotes to combat bacteriophages. With extreme accuracy, each recognizes a particular sequence in double-stranded DNA and cleaves at a fixed position within or nearby. The discoveries of these enzymes in the 1970s, and of the uses to which they could be put, have since impacted every corner of the life sciences. They became the enabling tools of molecular biology, genetics and biotechnology, and made analysis at the most fundamental levels routine. Hundreds of different REases have been discovered and are available commercially. Their genes have been cloned, sequenced and overexpressed. Most have been characterized to some extent, but few have been studied in depth. Here, we describe the original discoveries in this field, and the properties of the first Type II REases investigated. We discuss the mechanisms of sequence recognition and catalysis, and the varied oligomeric modes in which Type II REases act. We describe the surprising heterogeneity revealed by comparisons of their sequences and structures.

O6-Carboxymethylguanine in DNA forms a sequence context-dependent wobble base-pair structure with thymine

N-Nitrosation of glycine and its derivatives generates potent alkylating agents that can lead to the formation ofO6-carboxymethylguanine (O6-CMG) in DNA.O6-CMG has been identified in DNA derived from human colon tissue and its occurrence has been linked to diets high in red and processed meats, implying an association with the induction of colorectal cancer. By analogy toO6-methylguanine,O6-CMG is expected to be mutagenic, inducing G-to-A mutations that may be the molecular basis of increased cancer risk. Previously, the crystal structure of the DNA dodecamer d(CGCG[O6-CMG]ATTCGCG) has been reported, in whichO6-CMG forms a Watson–Crick-type pair with thymine similar to the canonical A:T pair. In order to further investigate the versatility ofO6-CMG in base-pair formation, the structure of the DNA dodecamer d(CGC[O6-CMG]AATTTGCG) containingO6-CMG at a different position has been determined by X-ray crystallography using four crystal forms obtained under conditions containing different solvent ions (Sr2+, Ba2+, Mg2+, K+or Na+) with and without Hoechst 33258. The most striking finding is that the pairing modes ofO6-CMG with T are quite different from those previously reported. In the present dodecamer, the T bases are displaced (wobbled) into the major groove to form a hydrogen bond between the thymine N3N—H and the carboxyl group ofO6-CMG. In addition, a water molecule is bridged through two hydrogen bonds between the thymine O2atom and the 2-amino group ofO6-CMG to stabilize the pairing. These interaction modes commonly occur in the four crystal forms, regardless of the differences in crystallization conditions. The previous and the present results show thatO6-CMG can form a base pair with T in two alternative modes: the Watson–Crick type and a high-wobble type, the nature of which may depend on the DNA-sequence context.

Coupling between hydrogen atoms transfer and stacking interaction in adenine-thymine/guanine-cytosine complexes: a theoretical study

Four different complexes of two base pairs, an adenine-thymine and a guanine-cytosine one, have been studied in order to understand the modifications induced by the staking interaction between the two base pairs on the hydrogen atoms transfers between the bases in either base pair. The inclusion of these two kinds of interactions allows us to clarify if some properties, as the mechanism of hydrogen transfer, is exclusively a local effect of a base pair or can be modified by a more long-range interaction between the base pairs. The results on these four complexes are compared with those of the monomeric systems, the A-T and G-C base pair, and with those of the A-T and G-C dimers. The specificity of each complex and of each hydrogen bond has been analyzed.

Structures of DNA duplexes containing O6-carboxymethylguanine, a lesion associated with gastrointestinal cancer, reveal a mechanism for inducing pyrimidine transition mutations

N-nitrosation of glycine and its derivatives generates potent alkylating agents that can lead to the formation of O(6)-carboxymethylguanine (O(6)-CMG) in DNA. O(6)-CMG has been identified in DNA derived from human colon tissue, and its occurrence has been linked to diets high in red and processed meats. By analogy to O(6)-methylguanine, O(6)-CMG is expected to be highly mutagenic, inducing G to A mutations during DNA replication that can increase the risk of gastrointestinal and other cancers. Two crystal structures of DNA dodecamers d(CGCG[O(6)-CMG]ATTCGCG) and d(CGC[O(6)-CMG]AATTCGCG) in complex with Hoechst33258 reveal that each can form a self-complementary duplex to retain the B-form conformation. Electron density maps clearly show that O(6)-CMG forms a Watson-Crick-type pair with thymine similar to the canonical A:T pair, and it forms a reversed wobble pair with cytosine. In situ structural modeling suggests that a DNA polymerase can accept the Watson-Crick-type pair of O(6)-CMG with thymine, but might also accept the reversed wobble pair of O(6)-CMG with cytosine. Thus, O(6)-CMG would permit the mis-incorporation of dTTP during DNA replication. Alternatively, the triphosphate that would be formed by carboxymethylation of the nucleotide triphosphate pool d[O(6)-CMG]TP might compete with dATP incorporation opposite thymine in a DNA template.

Theoretical investigation of the coupling between hydrogen-atom transfer and stacking interaction in adenine-thymine dimers

Three different dimers of the adenine-thymine (A-T) base pair are studied to point out the changes of important properties (structure, atomic charge, energy and so on) induced by coupling between the movement of the atoms in the hydrogen bonds and the stacking interaction. The comparison of these results with those for the A-T monomer system explains the role of the stacking interaction in the hydrogen-atom transfer in this biologically important base pair. The results support the idea that this coupling depends on the exact dimer considered and is different for the N-N and N-O hydrogen bonds. In particular, the correlation between the hydrogen transfer and the stacking interaction is more relevant for the N-N bridge than for the N-O one. Also, the two different mechanisms of two-hydrogen transfer (step by step and concerted) can be modified by the stacking interaction between the base pairs.

Structural studies on ligand-DNA systems: A robust approach in drug design

Molecular docking, molecular mechanics, molecular dynamics and relaxation matrix simulation protocols have been extensively used to generate the structural details of ligand-receptor complexes in order to understand the binding interactions between the two entities. Experimental methods like NMR spectroscopy and X-ray crystallography are known to provide structural information about ligand-receptor complexes. In addition, fluorescence spectroscopy, circular dichroism (CD) spectroscopy and molecular docking have also been utilized to decode the phenomenon of the ligand-DNA interactions, with good correlation between experimental and computational results. The DNA binding affinity was demonstrated by analysing fluorescence spectral data. Structural rigidity of DNA upon ligand binding was identified by CD spectroscopy. Docking is carried out using the DNA-Dock program which results in the binding affinity data along with structural information like interatomic distances and H-bonding, etc. The complete structural analyses of various drug-DNA complexes have afforded results that indicate a specific DNA binding pattern of these ligands. It also exhibited that certain structural features of ligands can make a ligand to be AT- or GC-specific. It was also demonstrated that changing specificity from AT base pairs to GC base pairs further improved the DNA topoisomerase inhibiting activity in certain ligands. Thus, a specific molecular recognition signature encrypted in the structure of ligand can be decoded and can be effectively employed in designing more potent antiviral and antitumour agents.

Can We Accurately Describe the Structure of Adenine Tracts in B-DNA? Reference Quantum-Chemical Computations Reveal Overstabilization of Stacking by Molecular Mechanics

Sequence-dependent local variations of helical parameters, structure, and flexibility are crucial for molecular recognition processes involving B-DNA. A-tracts, i.e., stretches of several consecutive adenines in one strand that are in phase with the DNA helical repeat, mediate significant DNA bending. During the past few decades, there have been intense efforts to understand the sequence dependence of helical parameters in DNA. Molecular dynamics (MD) simulations can provide valuable insights into the molecular mechanism behind the relationship between sequence and structure. However, although recent improvements in empirical force fields have helped to capture many sequence-dependent B-DNA properties, several problems remain, such as underestimation of the helical twist and suspected underestimation of the propeller twist in A-tracts. Here, we employ reference quantum mechanical (QM) calculations, explicit solvent MD, and bioinformatics to analyze the underestimation of propeller twisting of A-tracts in simulations. Although we did not identify a straightforward explanation, we discovered two imbalances in the empirical force fields. The first was overestimation of stacking interactions accompanied by underestimation of base-pairing energy, which we attribute to anisotropic polarizabilities that are not reflected by the isotropic force fields. This may lead to overstacking with potentially important consequences for MD simulations of nucleic acids. The second observed imbalance was steric clash between A(N1) and T(N3) nitrogens of AT base pairs in force-field descriptions, resulting in overestimation of the AT pair stretch in MD simulations. We also substantially extend the available set of benchmark estimated CCSD(T)/CBS data for B-DNA base stacking and provide a code that allows the generation of diverse base-stacking geometries suitable for QM computations with predefined intra- and interbase pair parameters.

Determinants of bacteriophage 933W repressor DNA binding specificity

We reported previously that 933W repressor apparently does not cooperatively bind to adjacent sites on DNA and that the relative affinities of 933W repressor for its operators differ significantly from that of any other lambdoid bacteriophage. These findings indicate that the operational details of the lysis-lysogeny switch of bacteriophage 933W are unique among lambdoid bacteriophages. Since the functioning of the lysis-lysogeny switch in 933W bacteriophage uniquely and solely depends on the order of preference of 933W repressor for its operators, we examined the details of how 933W repressor recognizes its DNA sites. To identify the specificity determinants, we first created a molecular model of the 933W repressor-DNA complex and tested the predicted protein-DNA interactions. These results of these studies provide a picture of how 933W repressor recognizes its DNA sites. We also show that, opposite of what is normally observed for lambdoid phages, 933W operator sequences have evolved in such a way that the presence of the most commonly found base sequences at particular operator positions serves to decrease, rather than increase, the affinity of the protein for the site. This finding cautions against assuming that a consensus sequence derived from sequence analysis defines the optimal, highest affinity DNA binding site for a protein.

HYDRATE MODEL OF THE EQUILIBRIUM DNA–WATER SYSTEMS

A hydrate model for the DNA– H2O system is proposed. For the notions developed, available data on a tendency of H2O molecules to form structured hydrates containing atomic groups, molecules, or atoms housed within the structural cavities formed by H-bonded H2O molecules are used. It is shown that the large and small cavities of the hydrate structure II are in close geometric agreement with N-bases and deoxyribose and with other atomic groups of DNA molecules, respectively. On the basis of the model proposed and with the Watson–Crick base-pairing scheme, the number of pairs of N-bases per helix turn (11.25), the density (1.161 g/cm3), and the helix step (0.567 nm) for the quasi-equilibrium DNA (RNA)–water system are computed, the last value nearly coinciding with the well-known alpha-helix step in the protein secondary structure (5.44 nm). Assumptions on the phenomenology of some stages of mitosis are presented.