Narrow Minor Groove Research Articles

Unusual bending patterns of spermidine3+ bound to DNA double helix

Natural polyamines play a fundamental role in the cell cycle. Despite being recognized as the most abundant organic counterions of DNA in the cell nucleus, their interactions with DNA have not been fully characterized. In a recent work [S. Perepelytsya, T. Vasiliu, A. Laaksonen, L. Engelbrecht, G. Brancato, and F. Mocci, J. Molec. Liq.389, 122828 (2023)], we have shown how the interactions between spermidine3+ and the DNA double helix induce significant conformational variations in the polyamine molecule. Specifically, we found that DNA induces conformations that are not observed in solution. Following that study, we present here a detailed investigation of the most compact conformation of the polyamine, analyzing its connection to the interaction with the DNA duplex. The analysis reveals that anomalous bent conformations of the spermidine3+ molecule result from the interaction of all three amino groups of the polyamine with the DNA phosphate groups on the minor groove side of the double helix. The changes in dihedral angles of the bent spermidine3+ molecule can be explained in terms of conformational transformations of six- and seven-membered rings, analogous to cyclohexane and cycloheptane. The analysis of the position of spermidine3+ molecule along the DNA surface reveals a sequence specificity of this binding mode with a marked preference for the narrow minor groove of A-tracts. The formation of the anomalous bent conformations of spermidine3+ in the complex with the DNA double helix is expected to be of paramount importance in understanding the mechanisms underlying DNA’s biological function.

Conformational Morphing by a DNA Analogue Featuring 7-Deazapurines and 5-Halogenpyrimidines and the Origins of Adenine-Tract Geometry.

Several efforts are currently directed at the creation and cellular implementation of alternative genetic systems composed of pairing components that are orthogonal to the natural dA/dT and dG/dC base pairs. In an alternative approach, Watson-Crick-type pairing is conserved, but one or all of the four letters of the A, C, G, and T alphabet are substituted by modified components. Thus, all four nucleobases were altered to create halogenated deazanucleic acid (DZA): dA was replaced by 7-deaza-2'-deoxyadenosine (dzA), dG by 7-deaza-2'-deoxyguanosine (dzG), dC by 5-fluoro-2'-deoxycytidine (FdC), and dT by 5-chloro-2'-deoxyuridine (CldU). This base-pairing system was previously shown to retain function in Escherichia coli. Here, we analyze the stability, hydration, structure, and dynamics of a DZA Dickerson-Drew Dodecamer (DDD) of sequence 5'-FdC-dzG-FdC-dzG-dzA-dzA-CldU-CldU-FdC-dzG-FdC-dzG-3'. Contrary to similar stabilities of DDD and DZA-DDD, osmotic stressing revealed a dramatic loss of hydration for the DZA-DDD relative to that for the DDD. The parent DDD 5'-d(CGCGAATTCGCG)-3' features an A-tract, a run of adenosines uninterrupted by a TpA step, and exhibits a hallmark narrow minor groove. Crystal structures─in the presence of RNase H─and MD simulations show increased conformational plasticity ("morphing") of DZA-DDD relative to that of the DDD. The narrow dzA-tract minor groove in one structure widens to resemble that in canonical B-DNA in a second structure. These changes reflect an indirect consequence of altered DZA major groove electrostatics (less negatively polarized compared to that in DNA) and hydration (reduced compared to that in DNA). Therefore, chemical modifications outside the minor groove that lead to collapse of major groove electrostatics and hydration can modulate A-tract geometry.

Modulation of the tetrameric I-motif folding of C-rich Tetrahymena telomeric sequences by hexitol nucleic acid (HNA) modifications

I-motifs are non-canonical DNA structures consisting of two parallel strands held together by hemiprotonated cytosine-cytosine+ base pairs, which intercalate to form a ordered column of stacked base pairs. This unique structure covers potential relevance in various fields, including gene regulation and biotechnological applications. A unique structural feature of I-motifs (iM), is the presence of sugar-sugar interactions through their extremely narrow minor grooves. Consistently, oligonucleotides containing pentose derivatives such as ribose, 2′-deoxyribose, arabinose, and 2′-deoxy-2′-fluoroarabinose highlighted a very different attitude to fold into iM. On the other hand, there is significant attention focused on exploring sugar-modifications that can increase nucleic acids resistance to nuclease degradation, a crucial requirement for therapeutic applications.An interesting example, not addressed in the iM field yet, is represented by hexitol nucleic acid (HNA), a metabolically stable six-membered ring analogue compatible with A-like double helix formation. Herein, we selected two DNA C-rich Tetrahymena telomeric sequences whose tetrameric iMs were already resolved by NMR and we investigated the iM folding of related HNA and RNA oligonucleotides by circular dichroism, differential scanning calorimetry and NMR. The comparison of their behaviours vs the DNA counterparts provided interesting insights into the influence of the sugar on iM folding. In particular, ribose and hexitol prevented iM formation. However, by clustering the hexitol-containing residues at the 3′-end, it was possible to modulate the distribution of the different topological species described for the DNA iMs. These data open new avenues for the exploitation of sugar modifications for I-motif characterization and applications.

DNA structural properties of DNA binding sites for 21 transcription factors in the mycobacterial genome.

Mycobacterium tuberculosis, the causative agent of tuberculosis, has evolved over time into a multidrug resistance strain that poses a serious global pandemic health threat. The ability to survive and remain dormant within the host macrophage relies on multiple transcription factors contributing to virulence. To date, very limited structural insights from crystallographic and NMR studies are available for TFs and TF-DNA binding events. Understanding the role of DNA structure in TF binding is critical to deciphering MTB pathogenicity and has yet to be resolved at the genome scale. In this work, we analyzed the compositional and conformational preference of 21 mycobacterial TFs, evident at their DNA binding sites, in local and global scales. Results suggest that most TFs prefer binding to genomic regions characterized by unique DNA structural signatures, namely, high electrostatic potential, narrow minor grooves, high propeller twist, helical twist, intrinsic curvature, and DNA rigidity compared to the flanking sequences. Additionally, preference for specific trinucleotide motifs, with clear periodic signals of tetranucleotide motifs, are observed in the vicinity of the TF-DNA interactions. Altogether, our study reports nuanced DNA shape and structural preferences of 21 TFs.

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.

C-B-A Test of DNA Force Fields.

The DNA duplex may be locally strongly bent in complexes with proteins, for example, with polymerases or in a nucleosome. At such bends, the DNA helix is locally in the noncanonical forms A (with a narrow major groove and a large amount of north sugars) or C (with a narrow minor groove and a large share of BII phosphates). To model the formation of such complexes by molecular dynamics methods, the force field is required to reproduce these conformational transitions for a naked DNA. We analyzed the available experimental data on the B-C and B-A transitions under the conditions easily implemented in modeling: in an aqueous NaCl solution. We selected six DNA duplexes which conformations at different salt concentrations are known reliably enough. At low salt concentrations, poly(GC) and poly(A) are in the B-form, classical and slightly shifted to the A-form, respectively. The duplexes ATAT and GGTATACC have a strong and salt concentration dependent bias toward the A-form. The polymers poly(AC) and poly(G) take the C- and A-forms, respectively, at high salt concentrations. The reproduction of the behavior of these oligomers can serve as a test for the balance of interactions between the base stacking and the conformational flexibility of the sugar-phosphate backbone in a DNA force field. We tested the AMBER bsc1 and CHARMM36 force fields and their hybrids, and we failed to reproduce the experiment. In all the force fields, the salt concentration dependence is very weak. The known B-philicity of the AMBER force field proved to result from the B-philicity of its excessively strong base stacking. In the CHARMM force field, the B-form is a result of a fragile balance between the A-philic base stacking (especially for G:C pairs) and the C-philic backbone. Finally, we analyzed some recent simulations of the LacI-, SOX-4-, and Sac7d-DNA complex formation in the framework of the AMBER force field.

Role of Water in Defining the Structure and Properties of B-Form DNA

DNA in the cell is rarely naked but normally protein-bound in nucleosomes. Of special interest is the DNA bound to other factors that control its key functions of transcription, replication, and repair. For these several transactions of DNA, the state of hydration plays an important role in its function, and therefore needs to be defined in as much detail as possible. High-resolution crystallography of short B-form duplexes shows that the mixed polar and apolar surface of the major groove binds water molecules over the broad polar floor of the groove in a sequence-dependent varied manner. In contrast, the narrower minor groove, particularly at AT-rich segments, binds water molecules to the polar groups of the bases in a regular double layer reminiscent of the structure of ice. This review is largely devoted to measurements made in solution, principally calorimetric, that are fully consistent with the location of water molecules seen in crystals, thereby emphasizing the substantial difference between the hydration patterns of the two grooves.

Dynamic effects of the spine of hydrated magnesium on viral RNA pseudoknot structure.

In the cellular environment, a viral RNA Pseudoknot (PK) structure is responsive to its prevailing ion atmosphere created by a mixture of monovalent and divalent cations. We investigate the influence of such a mixed-salt environment on RNA-PK structure at an atomic resolution through three sets of 1.5 μs explicit solvent molecular dynamics (MD) simulations and also by building a dynamic counterion-condensation (DCC) model at varying divalent Mg2+ concentrations. The DCC model includes explicit interaction of the ligand and adjacent chelated Mg2+ by extending the recently developed generalized Manning condensation model. Its implementation within an all-atom structure-based molecular dynamics framework bolsters its opportunity to explore large-length scale and long-timescale phenomena associated with complex RNA systems immersed in its dynamic ion environment. In the present case of RNA-PK, both explicit MD and DCC simulations reveal a spine of hydrated-Mg2+ to induce stem-I and stem-II closure where the minor groove between these stems is akin to breathing. Mg2+ mediated minor-groove narrowing is coupled with local base-flipping dynamics of a base triple and quadruple, changing the stem structure of such RNA. Contrary to the conversational view of the indispensable need for Mg2+ for the tertiary structure of RNA, the study warns about the plausible detrimental effect of specific Mg2+-phosphate interactions on the RNA-PK structure beyond a certain concentration of Mg2+.

Analysis of nucleoid-associated protein-binding regions reveals DNA structural features influencing genome organization in Mycobacteriumtuberculosis.

Nucleoid-associated proteins (NAPs) maintain bacterial nucleoid configuration through their architectural properties of DNA bending, wrapping, and bridging. However, the contribution of DNA structural alterations to DNA-NAP recognition at the genomic scale remains unresolved. Present work dissects the DNA sequence, shape and altered structural preferences at a genomic scale for six NAPs in Mycobacteriumtuberculosis. Results suggest narrower minor groove width (MGW) and higher DNA rigidity are marked for the binding sites of EspR and Lsr2, while mIHF, MtHU and NapM haveheterogeneous DNA structural predilections. In contrast, WhiB4-DNA-binding sites were characterized by wider MGW, highly deformable and less curved DNA. This work provides systematic insight into NAP-mediated genome organization as a function of DNA structural features.

DNA sequence-dependent positioning of the linker histone in a nucleosome: A single-pair FRET study

Linker histones (LHs) bind to nucleosomes with their globular domain (gH) positioned in either an on- or an off-dyad binding mode. Here, we study the effect of the linker DNA (L-DNA) sequence on the binding of a full-length LH, Xenopus laevis H1.0b, to a Widom 601 nucleosome core particle (NCP) flanked by two 40 bp long L-DNA arms, by single-pair FRET spectroscopy. We varied the sequence of the 11 bp of L-DNA adjoining the NCP on either side, making the sequence either A-tract, purely GC, or mixed with 64% AT. The labeled gH consistently exhibited higher FRET efficiency with the labeled L-DNA containing the A-tract than that with the pure-GC stretch, even when the stretches were swapped. However, it did not exhibit higher FRET efficiency with the L-DNA containing 64% AT-rich mixed DNA when compared to the pure-GC stretch. We explain our observations with a model that shows that the gH binds on dyad and that two arginines mediate recognition of the A-tract via its characteristically narrow minor groove. To investigate whether this on-dyad minor groove-based recognition was distinct from previously identified off-dyad major groove-based recognition, a nucleosome was designed with A-tracts on both the L-DNA arms. One A-tract was complementary to thymine and the other to deoxyuridine. The major groove of the thymine-tract was lined with methyl groups that were absent from the major groove of the deoxyuridine tract. The gH exhibited similar FRET for both these A-tracts, suggesting that it does not interact with the thymine methyl groups exposed on the major groove. Our observations thus complement previous studies that suggest that different LH isoforms may employ different ways of recognizing AT-rich DNA and A-tracts. This adaptability may enable the LH to universally compact scaffold-associated regions and constitutive heterochromatin, which are rich in such sequences.

Understanding sequence effect in DNA bending elasticity by molecular dynamic simulations

Structural elasticity of double-strand DNAs is very important for their biological functions such as DNA-ligand binding and DNA-protein recognition. By all-atom molecular dynamics simulations, we investigated the bending elasticity of DNA with three typical sequences including poly(A)-poly(T) (AA-TT), poly(AT)-poly(TA) (AT-TA), and a generic sequence (GENE). Our calculations indicate that, AA-TT has an apparently larger bending persistence length (P ∼63 nm) than GENE (P ∼49 nm) and AT-TA (P ∼48 nm) while the persistence length of AT-TA is only very slightly smaller than that of GENE, which agrees well with those from existing works. Moreover, through extensive electrostatic calculations, we found that the sequence-dependent bending elasticity is attributed to the sequence-dependent electrostatic bending energy for AA-TT, AT-TA and GENE, which is coupled to their backbone structures. Particularly, the apparently stronger bending stiffness of AA-TT is attributed to its narrower minor groove. Interestingly, for the three DNAs, we predicted the non-electrostatic persistence length of ∼17 nm, thus electrostatic interaction makes the major contribution to DNA bending elasticity. The mechanism of electrostatic energy dominating sequence effect in DNA bending elasticity is furtherly illustrated through the electrostatic calculations for a grooved coarse-grained DNA model where minor groove width and other microscopic structural parameters can be artificially adjusted.

Dynamics of K+ counterions around DNA double helix in the external electric field: A molecular dynamics study.

The structure of DNA double helix is stabilized by metal counterions condensed to a diffuse layer around the macromolecule. The dynamics of counterions in real conditions is governed by the electric fields from DNA and other biological macromolecules. In the present work the molecular dynamics study was performed for the system of DNA double helix with neutralizing K+ counterions and for the system of KCl salt solution in an external electric field of different strength (up to 32mV/Å). The analysis of ionic conductivities of these systems has shown that the counterions around the DNA double helix are slowed down compared with the KCl salt solution. The calculated values of ion mobility are within (0.05-0.4)mS/cm depending on the orientation of the external electric field relatively to the double helix. Under the electric field parallel to the macromolecule K+ counterions move along the grooves of the double helix staying longer in the places with narrower minor groove. Under the electric field perpendicular to the macromolecule the dynamics of counterions is less affected by DNA atoms, and starting with the electric field values about 30mV/Å the double helix undergoes a phase transition from a double-stranded to a single-strand state.

Molecular docking of potential cytotoxic alkylating carmustine derivatives 2-chloroethylnitrososulfamides analogues of 2-chloroethylnitrosoureas

The objective of this work was the molecular modelling by bio-isostery of 2-chloroethylnitrosoureas CENU into 2-chloroethylnitrososulfamides CENS derived from Carmustine. We evaluated the pharmacodynamic profile of the new chemical class by studying molecular docking using innovative software. Good molecular docking scores were obtained through Auto-dock vina of the PyRx 0.8 software, the energy of the complexes formed (Target-Ligand) during the interaction varies from - 5,400 to -5,700 Kcal/mol, the total average between the 45 conformers is −5,213 Kcal/mol. The results were validated by Auto-dock vina 1.5.6 in collaboration with the Molecular Chemistry and Natural Substances Laboratory at the Meknes Faculty of Science - Morocco, a range of -4,900 to -5,100 Kcal/mol was noted for CENS complexes derived from Carmustine with the 2DND target, reflecting a better CENS chemical affinity to the biological target and the stability of the ligand-DNA complex, compared with the analogue reference Carmustine with a score of - 4,700 Kcal/mol. By superimposing the results of molecular docking, analysis of data from the study of electrophilia based on load transfer ECT and publications on CENS, we can predict that inter-strand crosslink is likely to occur between the Guanine dG22 of strand B and the Cytosine dC3 of strand A, located in the poly dA-poly dT segment end within the narrow minor groove of the DNA target (2DND). The molecular docking study was a preliminary approach to understand the therapeutic mode of action of CENS. Communicated by Ramaswamy H. Sarma

Delving into Eukaryotic Origins of Replication Using DNA Structural Features.

DNA replication in eukaryotes is an intricate process, which is precisely synchronized by a set of regulatory proteins, and the replication fork emanates from discrete sites on chromatin called origins of replication (Oris). These spots are considered as the gateway to chromosomal replication and are stereotyped by sequence motifs. The cognate sequences are noticeable in a small group of entire origin regions or totally absent across different metazoans. Alternatively, the use of DNA secondary structural features can provide additional information compared to the primary sequence. In this article, we report the trends in DNA sequence-based structural properties of origin sequences in nine eukaryotic systems representing different families of life. Biologically relevant DNA secondary structural properties, namely, stability, propeller twist, flexibility, and minor groove shape were studied in the sequences flanking replication start sites. Results indicate that Oris in yeasts show lower stability, more rigidity, and narrow minor groove preferences compared to genomic sequences surrounding them. Yeast Oris also show preference for A-tracts and the promoter element TATA box in the vicinity of replication start sites. On the contrary, Drosophila melanogaster, humans, and Arabidopsis thaliana do not have such features in their Oris, and instead, they show high preponderance of G-rich sequence motifs such as putative G-quadruplexes or i-motifs and CpG islands. Our extensive study applies the DNA structural feature computation to delve into origins of replication across organisms ranging from yeasts to mammals and including a plant. Insights from this study would be significant in understanding origin architecture and help in designing new algorithms for predicting DNA trans-acting factor recognition events.

HIV chromatin is a preferred target for drugs that bind in the DNA minor groove.

The HIV genome is rich in A but not G or U and deficient in C. This nucleotide bias controls HIV phenotype by determining the highly unusual composition of all major HIV proteins. The bias is also responsible for the high frequency of narrow DNA minor groove sites in the double-stranded HIV genome as compared to cellular protein coding sequences and the bulk of the human genome. Since drugs that bind in the DNA minor groove disrupt nucleosomes on sequences that contain closely spaced oligo-A tracts which are prevalent in HIV DNA because of its bias, it was of interest to determine if these drugs exert this selective inhibitory effect on HIV chromatin. To test this possibility, nucleosomes were reconstituted onto five double-stranded DNA fragments from the HIV-1 pol gene in the presence and in the absence of several minor groove binding drugs (MGBDs). The results demonstrated that the MGBDs inhibited the assembly of nucleosomes onto all of the HIV-1 segments in a manner that was proportional to the A-bias, but had no detectable effect on the formation of nucleosomes on control cloned fragments or genomic DNA from chicken and human. Nucleosomes preassembled onto HIV DNA were also preferentially destabilized by the drugs as evidenced by enhanced nuclease accessibility in physiological ionic strength and by the preferential loss of the histone octamer in hyper-physiological salt solutions. The drugs also selectively disrupted HIV-containing nucleosomes in yeast as revealed by enhanced nuclease accessibility of the in vivo assembled HIV chromatin and reductions in superhelical densities of plasmid chromatin containing HIV sequences. A comparison of these results to the density of A-tracts in the HIV genome indicates that a large fraction of the nucleosomes that make up HIV chromatin should be preferred in vitro targets for the MGBDs. These results show that the MGBDs preferentially disrupt HIV-1 chromatin in vitro and in vivo and raise the possibility that non-toxic derivatives of certain MGBDs might serve as a novel class of anti-HIV agents.

The HRP3 PWWP domain recognizes the minor groove of double-stranded DNA and recruits HRP3 to chromatin.

HDGF-related protein 3 (HRP3, also known as HDGFL3) belongs to the family of HDGF-related proteins (HRPs) and plays an essential role in hepatocellular carcinoma pathogenesis. All HRPs have a PWWP domain at the N-terminus that binds both histone and DNA substrates. Despite previous advances in PWWP domains, the molecular basis by which HRP3 interacts with chromatin is unclear. In this study, we solved the crystal structures of the HRP3 PWWP domain in complex with various double-stranded DNAs with/without bound histone peptides. We found that HRP3 PWWP bound to the phosphate backbone of the DNA minor groove and showed a preference for DNA molecules bearing a narrow minor groove width. In addition, HRP3 PWWP preferentially bound to histone peptides bearing the H3K36me3/2 modification. HRP3 PWWP uses two adjacent surfaces to bind both DNA and histone substrates simultaneously, enabling us to generate a model illustrating the recruitment of PWWP to H3K36me3-containing nucleosomes. Cell-based analysis indicated that both DNA and histone binding by the HRP3 PWWP domain is important for HRP3 recruitment to chromatin in vivo. Our work establishes that HRP3 PWWP is a new family of minor groove-specific DNA-binding proteins, which improves our understanding of HRP3 and other PWWP domain-containing proteins.

Cryo-neutron crystallographic data collection and preliminary refinement of left-handed Z-DNA d(CGCGCG).

Crystals of left-handed Z-DNA [d(CGCGCG)]2 diffract X-rays to beyond 1 Å resolution, feature a small unit cell (∼18 × 31 × 44 Å) and are well hydrated, with around 90 water molecules surrounding the duplex in the asymmetric unit. The duplex shows regular hydration patterns in the narrow minor groove, on the convex surface and around sugar-phosphate backbones. Therefore, Z-DNA offers an ideal case to test the benefits of low-temperature neutron diffraction data collection to potentially determine the donor-acceptor patterns of first- and second-shell water molecules. Nucleic acid fragments pose challenges for neutron crystallography because water molecules are located on the surface rather than inside sequestered spaces such as protein active sites or channels. Water molecules can be expected to display dynamic behavior, particularly in cases where water is not part of an inner shell and directly coordinated to DNA atoms. Thus, nuclear density maps based on room-temperature diffraction data with a resolution of 1.6 Å did not allow an unequivocal determination of the orientations of water molecules. Here, cryo-neutron diffraction data collection for a Z-DNA crystal on the Macromolecular Neutron Diffractometer at the Spallation Neutron Source at Oak Ridge National Laboratory and the outcome of an initial refinement of the structure are reported. A total of 12 diffraction images were recorded with an exposure time of 3.5 h per image, whereby the crystal was static for each diffraction image with a 10° ϕ rotation between images. Initial refinements using these neutron data indicated the positions and orientations of 30 water molecules within the first hydration shell of the DNA molecule. This experiment constitutes a state-of-the-art approach and is the first attempt to our knowledge to determine the low-temperature neutron structure of a DNA crystal.

The floral homeotic protein SEPALLATA3 recognizes target DNA sequences by shape readout involving a conserved arginine residue in the MADS-domain.

SEPALLATA3 of Arabidopsis thaliana is a MADS-domain transcription factor (TF) and a key regulator of flower development. MADS-domain proteins bind to sequences termed 'CArG-boxes' [consensus 5'-CC(A/T)6 GG-3']. Because only a fraction of the CArG-boxes in the Arabidopsis genome are bound by SEPALLATA3, more elaborate principles have to be discovered to better understand which features turn CArG-boxes into genuine recognition sites. Here, we investigate to what extent the shape of the DNA is involved in a 'shape readout' that contributes to the binding of SEPALLATA3. We determined invitro binding affinities of SEPALLATA3 to DNA probes that all contain the CArG-box motif, but differ in their predicted DNA shape. We found that binding affinity correlates well with a narrow minor groove of the DNA. Substitution of canonical bases with non-standard bases supports the hypothesis of minor groove shape readout by SEPALLATA3. Analysis of mutant SEPALLATA3 proteins further revealed that a highly conserved arginine residue, which is expected to contact the DNA minor groove, contributes significantly to the shape readout. Our studies show that the specific recognition of cis-regulatory elements by a plant MADS-domain TF, and by inference probably also of other TFs of this type, heavily depends on shape readout mechanisms.

The sequence preference of DNA cleavage by T4 endonuclease VII

The enzyme T4 endonuclease VII is a resolvase that acts on branched DNA intermediates during genetic recombination, by cleaving DNA with staggered cuts approximately 3–6 bp apart. In this paper, we investigated the sequence preference of this cleavage reaction utilising two different DNA sequences. For the first time, the DNA sequence preference of T4 endonuclease VII cleavage sites has been examined without the presence of a known DNA substrate to mask any inherent nucleotide preference. The use of the ABI3730 platform enables the cleavage site to be determined at nucleotide resolution. We found that T4 endonuclease VII cleaves DNA with a sequence preference. We calculated the frequency of nucleotides surrounding the cleavage sites and found that following nucleotides had the highest incidence: AWTAN*STC, where N* indicates the cleavage site between positions 0 and 1, N is any base, W is A or T, and S is G or C. An A at position −1 and T at position +2 were the most predominant nucleotides at the cleavage site. Using a Sequence Logo method, the sequence TATTAN*CT was derived at the cleavage site. Note that A and T nucleotides were highly preferred 5′ to the cleavage sites in both methods of analysis. It was proposed that the enzyme recognises the narrower minor groove of these consecutive AT base pairs and cleaves DNA 3′ to this feature.

Comparative analysis of inosine-substituted duplex DNA by circular dichroism and X-ray crystallography

Leveraging structural biology tools, we report the results of experiments seeking to determine if the different mechanical properties of DNA polymers with base analog substitutions can be attributed, at least in part, to induced changes from classical B-form DNA. The underlying hypothesis is that different inherent bending and twisting flexibilities may characterize non-canonical B-DNA, so that it is inappropriate to interpret mechanical changes caused by base analog substitution as resulting simply from ‘electrostatic’ or ‘base stacking’ influences without considering the larger context of altered helical geometry. Circular dichroism spectra of inosine-substituted oligonucleotides and longer base-substituted DNAs in solution indicated non-canonical helical conformations, with the degree of deviation from a standard B-form geometry depending on the number of I⋅C pairs. X-ray diffraction of a highly inosine-substituted DNA decamer crystal (eight I⋅C and two A⋅T pairs) revealed an A-tract-like conformation with a uniformly narrow minor groove, reduced helical rise, and the majority of sugars adopting a C1′-exo (southeastern) conformation. This contrasts with the standard B-DNA geometry with C2′-endo sugar puckers (south conformation). In contrast, the crystal structure of a decamer with only four I⋅C pairs has a geometry similar to that of the reference duplex with eight G⋅C and two A⋅T pairs. The unique crystal geometry of the inosine-rich duplex is noteworthy given its unusual CD signature in solution and the altered mechanical properties of some inosine-containing DNAs.