DNA-induced conformational changes in RecA protein. Evidence for structural heterogeneity among nucleoprotein filaments and implications for homologous pairing

We have used circular dichroism as a probe to characterize the solution conformational changes in RecA protein upon binding to DNA. This approach revealed that RecA protein acquires significant amounts of alpha-helix upon interaction with DNA. These observations, consistent with the data from crystal structure (Story, R. M., Weber, I., and Steitz, T. (1992) Nature 355, 318-325), support the notion that some basic domains including the DNA binding motifs of RecA protein are unstructured and might contribute to the formation of alpha-helix. A comparison of nucleoprotein filaments comprised of RecA protein and a variety of DNA substrates revealed important structural heterogeneity. The most significant difference was observed with poly(dG). poly(dC) and related polymers, rich in GC sequences, which induced minimal amounts of alpha-helix in RecA protein. The magnitude of induction of alpha-helix in RecA protein, which occurred concomitant with the production of ternary complexes, was 2-fold higher with homologous than heterologous duplex DNA. Most importantly, the stimulation of ATP hydrolysis by high salt coincided with that of the induction of alpha-helix in RecA protein. These conformational differences provide a basis for thinking about the biochemical and structural transitions that RecA protein experiences during the formal steps of presynapsis, recognition, and alignment of homologous sequences.

Isolation and visualization of active presynaptic filaments of recA protein and single-stranded DNA

In the presence of ATP, recA protein forms a presynaptic complex with single-stranded DNA that is an obligatory intermediate in homologous pairing. Presynaptic complexes of recA protein and circular single strands that are active in forming joint molecules can be isolated by gel filtration. These isolated active complexes are nucleoprotein filaments with the following characteristics: (i) a contour length that is at least 1.5 times that of the corresponding duplex DNA molecule, (ii) an ordered structure visualized by negative staining as a striated filament with a repeat distance of 9.0 nm and a width of 9.3 nm, (iii) approximately 8 molecules of recA protein and 20 nucleotide residues per striation. The widened spacing between bases in the nucleoprotein filament means that the initial matching of complementary sequences must involve intertwining of the filament and duplex DNA, unwinding of the latter, or some combination of both to equalize the spacing between nascent base pairs. These experiments support the concept that recA protein first forms a filament with single-stranded DNA, which in turn binds to duplex DNA to mediate both homologous pairing and subsequent strand exchange.

The complete mitochondrial genome of the eastern grey kangaroo (Macropus giganteus)

We present the complete mitochondrial genome (accession number: LK995454) of an iconic Australian species, the eastern grey kangaroo (Macropus giganteus). The mitogenomic organization is consistent with other marsupials, encoding 13 protein-coding genes, 22 tRNA genes, 2 ribosomal RNA genes, an origin of light strand replication and a control region or Dloop. No repetitive sequences were detected in the control region. The M. giganteus mitogenome exemplifies a combination of tRNA gene order and structural peculiarities that appear to be unique to marsupials. We present a maximum likelihood phylogeny based on complete mitochondrial protein and RNA coding sequences that confirms the phylogenetic position of the grey kangaroo among macropodids.

Enhanced Tetrahedral Ordering of Water Molecules in Minor Grooves of DNA: Relative Role of DNA Rigidity, Nanoconfinement, and Surface Specific Interactions

Confinement and Surface specific interactions call induce Structures otherwise unstable at that temperature and pressure. Here we Study the groove specific water dynamics ill the nucleic acid sequences, poly-AT and poly-GC, in long B-DNA duplex chains by large scale atomistic molecular dynamics simulations, accompanied by thermodynamic analysis. While water dynamics in the major groove remains insensitive to the sequence differences, exactly the opposite is true for the minor groove water. Much slower water dynamics observed in the minor grooves (especially in the AT minor) call be attributed to all enhanced tetrahedral ordering (< t(h)>) of water. The largest value of < t(h)> in the AT minor groove is related to the spine of hydration found in X-ray Structure. The calculated configurational entropy (S-C) of the water molecules is found to be correlated with the self-diffusion coefficient of water in different region via Adam-Gibbs relation D = A exp(-B/TSC), and also with < t(h)>.

Polypeptide Helices in Hybrid Peptide Sequences

A new class of polypeptide helices in hybrid sequences containing alpha-, beta-, and gamma-residues is described. The molecular conformations in crystals determined for the synthetic peptides Boc-Leu-Phe-Val-Aib-beta Phe-Leu-Phe-Val-OMe 1 (beta Phe: (S)-beta(3)-homophenylalanine) and Boc-Aib-Gpn-AibGpn-OM2(Gpn:1-(aminomethyl)cycl hexaneacetic acid) reveal expanded helical turns in the hybrid sequences (alpha alpha beta)(n) and (ay), In 1, a repetitive helical structure composed Of C-14 hydrogen-bonded units is observed, whereas 2 provides an example of a repetitive C-12 hydrogen-bonded structure. Using experimentally determined backbone torsion angles for the hydrogen-bonded units formed by hybrid sequences, we have generated energetically favorable hybrid helices. Conformational parameters are provided for C-11, C-12, C-13, C-14, and C-15 helices in hybrid sequences.

DNA duplex with the potential to change handedness after every half a turn

Polymorphic forms of the DNA duplex with long stretches of structural monotony are known. Several alternating purine-pyrimidine sequences have been shown to adopt left-handed Z-conformation. We report a DNA sequence d(CGCGCGATCGAT)n exhibiting alternating right-handed B and left-handed Z helical conformation after every half a turn. Further, this unusual conformation with change in handedness after every six base pairs was induced at physiological superhelical density.

Structure of uracil-DNA glycosylase from Mycobacterium tuberculosis: insights into interactions with ligands

Uracil N-glycosylase (Ung) is the most thoroughly studied of the group of uracil DNA-glycosylase (UDG) enzymes that catalyse the first step in the uracil excision-repair pathway. The overall structure of the enzyme from Mycobacterium tuberculosis is essentially the same as that of the enzyme from other sources. However, differences exist in the N- and C-terminal stretches and some catalytic loops. Comparison with appropriate structures indicate that the two-domain enzyme closes slightly when binding to DNA, while it opens slightly when binding to the proteinaceous inhibitor Ugi. The structural changes in the catalytic loops on complexation reflect the special features of their structure in the mycobacterial protein. A comparative analysis of available sequences of the enzyme from different sources indicates high conservation of amino-acid residues in the catalytic loops. The uracil-binding pocket in the structure is occupied by a citrate ion. The interactions of the citrate ion with the protein mimic those of uracil, in addition to providing insights into other possible interactions that inhibitors could be involved in.

Sequence specificity of Z-DNA formation in oligonucleotides

The sequence specific requirement for B----Z transition in solution was examined in d(CGTGCGCACG), d(CGTACGTACG), d(ACGTACGT) in presence of various Z-inducing factors. Conformational studies show that inspite of the alternating nature of purines and pyrimidines, the aforementioned sequences do not undergo B----Z transition under the influence of NaCl, hexamine cobalt chloride and ethanol. A comparison with the crystal structures of an assorted array of purine and pyrimidine sequences show that the sequence requirement for B----Z transition is much more stringent in solution as compared to the solid state. The disruptive influence of AT base pairs in B to Z transition is discussed.

Theoretical studies on alpha-helix--DNA interactions

The interaction energies between (Ala)10 and alpha-helix fragment and different nucleotide sequences in right-handed B-form have been optimized using semi-empirical potential energy functions. The energies are calculated for two different orientations of the alpha-helix, viz., when the alpha-helix axis taken in the N----C direction is (i) parallel and (ii) antiparallel to the 5'-3' ascending strand of DNA, proximal to it. When both the DNA molecule as well as the alpha-helix are treated as rigid molecules it is found that a polyalanine alpha-helix has slightly more favourable contacts when it is in the proximity of a four nucleotide sequence of 5'-(N-A-T-N)-3' type, where N is either a purine or a pyrimidine. However, when the two interacting molecules are allowed to undergo local structural variations then the interaction energy appears to be independent of the base sequence confirming the non-specific nature of these interactions.

DNA polymorphism and local variation in base-pair orientation: a theoretical rationale

Basepair stacking calculations have been carried out to understand the conformational polymorphism of DNA and its sequence dependence. The recently developed self-consistent parameter set, which is specially suitable for describing irregular DNA structures, has been used to describe the geometry of a basepair doublet. While for basepairs without any propeller, the favourable stacking patterns do not appear to have very strong features, much more noticeable sequence dependent stacking patterns emerge once a propeller is applied to the basepairs. The absolute minima for most sequences occurs for a doublet geometry close to the B-DNA fibre models. Hence in the B-DNA region, no strong sequence dependent features are found, but the range of doublet geometries observed in the crystal structures generally lie within the low energy contours, obtained from stacking energy calculations. The doublet geometry corresponding to the A-DNA fibre model is not energetically favourable for the purine-pyrimidine sequences, which prefer small roll angle values when the slide has a large negative value as in A-DNA. However positive roll with large negative slide is allowed for GG, GA, AG and the pyrimidine-purine steps. This is consistent with the observed geometries of various steps in A-DNA crystals. Thus the general features of the basepair doublets predicted from these theoretical studies agree very well with the results from crystal structure analysis. However, since most sequences show an overall preference for B-type doublet geometry, the B --> A transition for random sequence DNA cannot be explained on the basis of basepair stacking interactions.

Modelling studies on neurodegenerative disease causing triplet repeat sequences d(CGG/GCC) and d(CAG/CTG): Duplexes vs Quadruplexes

Model building and molecular mechanics studies have been carried out to examine the potential structures for d(GGC/GCC)5 and d(CAG/CTG)5 that might relate to their biological function and association with triplet repeat expansion diseases. Model building studies suggested that hairpin and quadruplex structures could be formed with these repeat sequences. Molecular mechanics studies have demonstrated that the hairpin and hairpin dimmer structures of triplet repeat sequences formed by looping out of the two strands are as favourable as the corresponding B-DNA type hetero duplex structures. Further, at high salt condition, Greek key type quadruplex structures are energetically comparable with hairpin dimer and B-DNA type duplex structures. All tetrads in the quadruplex structures are well stacked and provide favourable stacking energy values. Interestingly, in the energy minimized hairpin dimer and Greek key type quadruplex structures, all the bases even in the non-G tetrads are cyclically hydrogen bonded, even though the A, C and T-tetrads were not hydrogen bonded in the starting structures.

Structural features of B-DNA dodecamer crystal structures: Influence of crystal packing versus base sequence

We have analyzed the set of inter and intra base pair parameters for each dinucleotide step in single crystal structures of dodecamers, solved at high and medium resolution and all crystallized in P2(1)2(1)2(1) space group. The objective was to identify whether all the structures which have either the Drew-Dickerson (DD) sequence d[CGCGAATTCGCG] with some base modification or related sequence (non-DD), would display the same sequence dependent structural variability about its palindromic sequence, despite the molecule being bent at one end because of similar crystal lattice packing effect. Most of the local doublet parameters for base pairs steps G2-C3 and G10-C11 positions, symmetrically situated about the lateral twofold, were significantly correlated between themselves. In non-DD sequences, significant correlations between these positional parameters were absent. The different range of local step parameter values at each sequence position contributed to the gross feature of smooth helix axis bending in all structures. The base pair parameters in some of the positions, for medium resolution DD sequence, were quite unlike the high-resolution set and encompassed a higher range of values. Twist and slide are the two main parameters that show wider conformational range for the middle region of non-DD sequence structures in comparison to DD sequence structures. On the contrary, the minor and major groove features bear good resemblance between DD and non-DD sequence crystal structure datasets. The sugar-phosphate backbone torsion angles are similar in all structures, in sharp contrast to base pair parameter variation for high and low resolution DD and non-DD sequence structures, consisting of unusual (epsilon =g(-), xi =t) B-II conformation at the 10(th) position of the dodecamer sequence. Thus examining DD and non-DD sequence structures packed in the same crystal lattice arrangement, we infer that inter and intra base pair parameters are as symmetrically equivalent in its value as the symmetry related step for the palindromic DD sequence about lateral two-fold axis. This feature would lead us to agree with the conclusion that DNA conformation is not substantially affected by end-to-end or lateral inter-molecular interaction due to crystal lattice packing effect. Non-DD sequence structures acquire step parameter values which reflect the altered sequence at each of the dodecamer sequence position in the orthorhombic lattice while showing similar gross features of DD sequence structures

Local variability and base sequence effects in DNA crystal structures.

The importance and usefulness of local doublet parameters in understanding sequence dependent effects has been described for A- and B-DNA oligonucleotide crystal structures. Each of the two sets of local parameters described by us in the NUPARM algorithm, namely the local doublet parameters, calculated with reference to the mean z-axis, and the local helical parameters, calculated with reference to the local helix axis, is sufficient to describe the oligonucleotide structures, with the local helical parameters giving a slightly magnified picture of the variations in the structures. The values of local doublet parameters calculated by NUPARM algorithm are similar to those calculated by NEWHELIX90 program, only if the oligonucleotide fragment is not too distorted. The mean values obtained using all the available data for B-DNA crystals are not significantly different from those obtained when a limited data set is used, consisting only of structures with a data resolution of better than 2.4 A and without any bound drug molecule. Thus the variation observed in the oligonucleotide crystals appears to be independent of the quality of their crystallinity. No strong correlation is seen between any pair of local doublet parameters but the local helical parameters are interrelated by geometric relationships. An interesting feature that emerges from this analysis is that the local rise along the z-axis is highly correlated with the difference in the buckle values of the two basepairs in the doublet, as suggested earlier for the dodecamer structures (Bansal and Bhattacharyya, in Structure & Methods: DNA & RNA, Vol. 3 (Eds., R.H. Sarma and M.H. Sarma), pp. 139-153 (1990)). In fact the local rise values become almost constant for both A- and B-forms, if a correction is applied for the buckling of the basepairs. In B-DNA the AA, AT, TA and GA basepair sequences generally have a smaller local rise (3.25 A) compared to the other sequences (3.4 A) and this seems to be an intrinsic feature of basepair stacking interaction and not related to any other local doublet parameter. The roll angles in B-DNA oligonucleotides have small values (less than +/- 8 degrees), while mean local twist varies from 24 degrees to 45 degrees. The CA/TG doublet sequences show two types of preferred geometries, one with positive roll, small positive slide and reduced twist and another with negative roll, large positive slide and increased twist.(ABSTRACT TRUNCATED AT 400 WORDS)

CA/TG sequence at the 5' end of oligo(A)-tracts strongly modulates DNA curvature

An analysis of the base pair doublet geometries in available crystal structures indicates that the often reported intrinsic curvature of DNA containing oligo-(d(A).d(T)) tracts may also depend on the nature of the flanking sequences. The presence of CA/TG doublet in particular at the 5' end of these tracts is expected to enhance their intrinsic bending property. To test this proposition, three oligonucleotides, d(GAAAAACCCCCC), d(CCCCCCAAAAAG), d(GAAAAATTTTTC), and their complementary sequences were synthesized to study the effect of various flanking sequences, at the 5' and 3' ends of the A-tracts, on the curvature of DNA in solution. An analysis of the polyacrylamide gel electrophoretic mobilities of these sequences under different conditions of salts and temperatures (below their melting points) clearly showed that the oligomer with CA/TG sequence in the center was always more retarded than the oligomer with AC/GT sequence, as well as the oligomer with AT/AT sequence. Hydroxyl radical probing of the sequences with AC/GT and CA/TG doublet junctions gives a similar cutting pattern in the A-tracts, which is quite different from that in the C-tracts, indicating that the oligo(A)-tracts have similar structures in the two oligomers. KMnO4 probing shows that the oligomer with a CA/TG doublet junction forms a kink that is responsible for its inherent curvature and unusual electrophoretic mobility. UV melting shows a reduced thermal stability of the duplex with CA/TG doublet junction, and circular dichroism (CD) studies indicate that a premelting transition occurs in the oligomer with CA/TG doublet step before global melting but not in the oligomer with AC/GT doublet step, which may correspond to thermally induced unbending of the oligomer. These observations indicate that the CA/TG doublet junction at the 5' end of the oligo(A)-tract has a crucial role in modulating the overall curvature in DNA.

CA\TG Sequence at the 5'end of Oligo(A)-tracts Strongly Modulates DNA Curvature

An analysis of the base pair doublet geometries in available crystal structures indicates that the often reported intrinsic curvature of DNA containing oligo-(d(A).d(T)) tracts may also depend on the nature of the flanking sequences. The presence of CA/TG doublet in particular at the 5' end of these tracts is expected to enhance their intrinsic bending property. To test this proposition, three oligonucleotides, d(GAAAAACCCCCC), d(CCCCCCAAAAAG), d(GAAAAATTTTTC), and their complementary sequences were synthesized to study the effect of various flanking sequences, at the 5' and 3' ends of the A-tracts, on the curvature of DNA in solution. An analysis of the polyacrylamide gel electrophoretic mobilities of these sequences under different conditions of salts and temperatures (below their melting points) clearly showed that the oligomer with CA/TG sequence in the center was always more retarded than the oligomer with AC/GT sequence, as well as the oligomer with AT/AT sequence. Hydroxyl radical probing of the sequences with AC/GT and CA/TG doublet junctions gives a similar cutting pattern in the A-tracts, which is quite different from that in the C-tracts, indicating that the oligo(A)-tracts have similar structures in the two oligomers. KMnO4 probing shows that the oligomer with a CA/TG doublet junction forms a kink that is responsible for its inherent curvature and unusual electrophoretic mobility. UV melting shows a reduced thermal stability of the duplex with CA/TG doublet junction, and circular dichroism (CD) studies indicate that a premelting transition occurs in the oligomer with CA/TG doublet step before global melting but not in the oligomer with AC/GT doublet step, which may correspond to thermally induced unbending of the oligomer. These observations indicate that the CA/TG doublet junction at the 5' end of the oligo(A)-tract has a crucial role in modulating the overall curvature in DNA.