A self-consistent formulation for analysis and generation of non-uniform DNA structure

A fully self-consistent formulation is described here for the analysis and generation of base-pairs in non-uniform DNA structures, in terms of various local parameters. It is shown that the internal "wedge parameters" are mathematically related to the parameters describing the base-pair orientation with respect to an external helix axis. Hence any one set of three translation and three rotation parameters are necessary and sufficient to completely describe the relative orientation of the base-pairs comprising a step (or doublet). A general procedure is outlined for obtaining an average or global helix axis from the local helix axes for each step. A graphical representation of the local helix axes in the form of a polar plot is also shown and its application for estimating the curvature of oligonucleotide structures is illustrated, with examples of both A and B type structures.

A general procedure for generation of curved DNA molecules.

A general method for generation of base-pairs in a curved DNA structure, for any prescribed values of helical parameters--unit rise (h), unit twist (theta), wedge roll (theta R) and wedge tilt (theta T), propeller twist (theta p) and displacement (D) is described. Its application for generation of uniform as well curved structures is also illustrated with some representative examples. An interesting relationship is observed between helical twist (theta), base-pair parameters theta x, theta y and the wedge parameters theta R, theta T, which has important consequences for the description and estimation of DNA curvature.

Crystal and molecular structure of the ammonium salt of the dinucleoside monophosphate d(CpG)

The crystal and molecular structure of the ammonium salt of deoxycytidylyl-(3'-5')-deoxyguanosine has been determined from 0.85 A resolution single crystal X-ray diffraction data. The crystals obtained by acetone diffusion technique at -20 degrees C, are orthorhombic, P212121, a = 12.880(2), b = 17444(2) and c = 27.642(2) A. The structure was solved by high resolution Patterson and Fourier methods and refined to R = 0.136. There are two d(CpG) molecules in the asymmetric unit forming a mini left handed Z-DNA helix. This is in contrast to the earlier reported forms of d(CpG) where the molecules form self base paired duplexes. There are two ammonium ions in the asymmetric unit. The major groove NH+4 ion interacts with N7 of guanines through water bridges besides making H-bonded interactions directly with the phosphate oxygen atoms. A second NH+4 ion is found in the minor groove interacting directly with the phosphate oxygen atoms. Symmetry related molecules pack in such a way that the cytosine base stacks on cytosine and guanine base on guanine. Our structure demonstrates that alternating d(CpG) sequences have the ability to adopt the left handed Z-DNA structure even at the dimer level i.e., in a sequence which is only two base pairs long.

Stabilization of the bio-membrane by small molecules: interaction of trehalose with the phospholipid bilayer

Anhydrobiotic organisms undergo periods of acute dehydration during their life cycle. It is of interest to understand how the biomembrane remains intact through such stress. A disaccharide, trehalose, which is metabolised during anhydrobiosis is found to prevent disruption of model membrane systems. Molecular modelling techniques are used to investigate the possible mode of interaction of trehalose with a model monolayer. The objective is to maximise hydrogen bonding between the two systems. A phospholipid matrix consisting of 1,2-dimyristoyl-sn-glycero-3-phosphatidylcholine (DMPC) is chosen to represent the monolayer. The crystal structure of DMPC reveals that there are two distinct conformers designated as A and B. An expansion of the monolayer, coplanar with its surface, results in the trehalose molecule being accommodated in a pocket formed by four B conformers. One glucose ring of the sugar rests on the hydrophobic patch provided by the choline methyls of an A conformer. Five hydrogen bonds are formed involving the phosphate oxygens of three of the surrounding B conformers. The model will be discussed with reference to relevant experimental data on the interaction.

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.

Interaction of metal ions with nucleic acids and related compounds. II. Studies on Au(III)-nucleic acid system

The nature of interaction of Au(III) with nucleic acids was studied by using methods such as uv and ir spectrophotometry, viscometry, pH titrations, and melting-temperature measurements. Au(III) is found to interact slowly with nucleic acids over a period of several hours. The uv spectra of native calf-thymus DNA 9pH 5.6 acetate buffer containing (0.01M NaCIO4) showed a shift in λ max to high wavelengths and an increase in optical density at 260 nm. There was a fourfold decrease in viscosity (expressed as ηsp/c). The reaction was faster at pH 4.0 and also with denatured DNA (pH 5.6) and whole yeast RNA (pH 5.6). The order of preference of Au(III) (as deduced from the time of completion of reaction) for the nucleic acids in RNA > denatured DNA > DNA. The reaction was found to be completely reversible with respect KCN. Infrared spectra of DNA-Au(III) complexes showed binding to both the phosphate and bases of DNA. The same conclusions were also arrived at by melting-temperature studies of Au(III)-DNA system. pH titrations showed liberation of two hydroxylions at r = 0.12 [r = moles of HAuCl4 added per mole of DNA-(P)] and one hydrogen ion at r = 0.5. The probable binding sites could be N(1)/N(7) of adenine, N(7) and/or C(6)O of guanine, N(3) of cytosine and N(3) of thymine. DNAs differing in their (G = C)-contents [Clostridium perfingens DNA(G = C, 29%), salmon sperm DNA (G + C, 42%) and Micrococcus lysodeikticus DNA(G + C, 29%), salmon sperm DNA (G = C, 72%)] behaved differently toward Au(III). The hyperchromicity observed for DNAs differing in (G + C)-content and cyanide reversal titrations indicate selectivity toward ( A + T)-rich DNA at lw values of r. Chemical analysis and job's continuous variation studies indicated the existence of possible complexes above and below r = 1. The results indicate that Au(III) ions probably bind to hte phosphate group in the initial stages of the reaction, particularly at low values of r, and participation of the base interaction also increases. Cross-linking of the two strands by Au(III) may take place, but a complete collapse of the doulbe helix is not envisaged. It is probable that tilting of the bases or rotaiton of the bases around the glucosidic bond, resulting in a significant distrotion of the double helix, might take place due to binding of Au(III) to DNA.

Computer modelling studies on the mechanism of action of ribonuclease T1

The mechanism of action of ribonuclease (RNase) T1 is still a matter of considerable debate as the results of x-ray, 2-D nmr and site-directed mutagenesis studies disagree regarding the role of the catalytically important residues. Hence computer modelling studies were carried out by energy minimisation of the complexes of RNase T1 and some of its mutants (His40Ala, His40Lys, and Glu58Ala) with the substrate guanyl cytosine (GpC), and of native RNase T1 with the reaction intermediate guanosine 2',3'-cyclic phosphate (G greater than p). The puckering of the guanosine ribose moiety in the minimum energy conformer of the RNase T1-GpC (substrate) complex was found to be O4'-endo and not C3'-endo as in the RNase T1-3'-guanylic acid (inhibitor/product) complex. A possible scheme for the mechanism of action of RNase T1 has been proposed on the basis of the arrangement of the catalytically important amino acid residues His40, Glu58, Arg77, and His92 around the guanosine ribose and the phosphate moiety in the RNase T1-GpC and RNase T1-G greater than p complexes. In this scheme, Glu58 serves as the general base group and His92 as the general acid group in the transphosphorylation step. His40 may be essential for stabilising the negatively charged phosphate moiety in the enzyme-transition state complex.

Distribution of simple repetitive (TG/CA)n and (CT/AG)n sequences in human and rodent genomes

Sixteen million nucleotide sequence of genome of various organisms have been analysed to detect and study the extent of occurrence of simple repetitive sequences. Two sequence motifs (TG/CA)n and (CT/AG)n capable of adopting unusual DNA structures, left handed Z-conformation and triple-helical conformation respectively, are found to be abundant in rodent and human genomes, but almost completely absent in bacterial genome. (TG/CA)n and (CT/AG)n sequences are present mostly in the intron or 5'/3' flanking regions of the genes. The presence of such repeat motifs in genomic sequence of higher eukaryotes has been correlated with their possible functional significance in nucleosome organization, recombination and gene expression.

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.

Studies on nucleotidases in plants: Fluorescence and kinetic properties of nucleotide pyrophosphatase from mung bean (Phaseolus aureus) seedlings

Nucleotide pyrophosphatase of mung bean seedlings has earlier been isolated in our laboratory in a dimeric form (Mr 65,000) and has been shown to be converted to a tetramer by AMP and to a monomer by p-hydroxymercuribenzoate. All the molecular forms were enzymatically active with different kinetic properties. By a modified procedure using blue-Sepharose affinity chromatography, we have now obtained a dimeric form of the enzyme which is desensitized to AMP interaction. The molecular weight of the desensitized form of the enzyme was found to be the same as that of the native dimeric enzyme. However, the desensitized enzyme functioned with a linear time course, contrary to the biphasic time course exhibited by the native enzyme. In addition, it was not converted to a tetramer on the addition of AMP, had only one binding site for adenine nucleotides, and p-hydroxy-mercuribenzoate had no effect on the time course of the reaction or on the molecular weight of the enzyme. The temperature optimum of the desensitized enzyme was found to be 67 °C in contrast to the optimum of 49 °C for the native dimer. Fifty percent of the tryptophan residues of the desensitized enzyme were not accessible for quenching by iodide. Fluorescence studies gave Kd values of 0.34, 2.2, and 0.8 mImage for AMP, ADP, and ATP, which were close to the Ki values of 0.12, 2.2, and 0.9 mImage , respectively, for these nucleotides. The binding and inhibition studies with AMP and its analogs showed that the 6-amino group and the 5′-phosphate group were essential for the inhibition of the enzyme activity.

Regulation of the activity of mung bean (Phaseolus aureus) glutamine synthetase by amino acids and nucleotides

The activity of glutamine synthetase isolated from the germinated seedlings of Phaseolus aureus was regulated by feedback inhibition by alanine, glycine, histidine, AMP, and ADP. When glutamate was the varied substrate, alanine, histidine, and glycine were partial noncompetitive, competitive, and mixed-type inhibitors, respectively. The type of inhibition by these amino acids was confirmed by fractional inhibition analysis. The adenine nucleotides, AMP and ADP, completely inhibited the enzyme activity and were competitive with respect to ATP. Multiple inhibition analyses revealed the presence of separate and nonexclusive binding sites for the amino acids and mutually exclusive sites for adenine nucleotides. Cumulative inhibition was observed with these end products.

A glossary of DNA structures from A to Z

The right-handed double-helical Watson-Crick model for B-form DNA is the most commonly known DNA structure. In addition to this classic structure, several other forms of DNA have been observed and it is clear that the DNA molecule can assume different structures depending on the base sequence and environment. The various forms of DNA have been identified as A, B, C etc. In fact, a detailed inspection of the literature reveals that only the letters F, Q, U, V and Y are now available to describe any new DNA structure that may appear in the future. It is also apparent that it may be more relevant to talk about the A, B or C type dinucleotide steps, since several recent structures show mixtures of various different geometries and a careful analysis is essential before identifying it as a 'new structure'. This review provides a glossary of currently identified DNA structures and is quite timely as it outlines the present understanding of DNA structure exactly 50 years after the original discovery of DNA structure by Watson and Crick

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

Mechanistic Aspects of Nucleophilic Substitution at Half-Sandwich Metal Complexes

The hydrolysis reactions of organometallic ruthenium(II) piano-stool complexes of the type Ru-II(eta(6)-cymene)(L)Cl](0/+) (1-5, where L = kappa(1)- or kappa(2)-1,1-bis(diphenylphosphino)methane,1,1bis-(diphenylphosphino)methane oxide, kappa(1)-mercaptobenzothiazole) have been studied using density functional theory at the B3LYP level. In addition to considering a syn attack in an associative fashion, where the nucleophile approaches from the same side as the leaving group, we have explored alternative paths such as an anti attack in an associative manner, where the nucleophile attacks from the opposite side of the leaving group. During the anti attack, an intermediate is formed and there is a coordination mode change of the arene ring from eta(6) to eta(2) along with its rotation. When the intermediate goes to the product, the arene ring slips back from eta(2) to eta(6) coordination. This coordinated movement of the arene ring makes the associative anti attack an accessible pathway for the substitution process. Our calculations predict very similar activation barriers for both syn and anti attacks. In the dissociative path, the rate-determining step is the generation of a coordinatively unsaturated 16-electron ruthenium species. This turns out to be viable once solvent effects are included. The large size of the ancillary ligands on Ru makes the dissociative process as favorable as the associative process. Activation energy calculations reveal that although the dissociative path is favorable for kappa(1) complexes, both dissociative and associative processes can have significant contribution to the hydrolysis reaction in kappa(2) complexes. Once activated by hydrolysis, these complexes react with guanine and adenine bases of DNA. The thermodynamic stabilities of complexes formed with the nucleobases are also presented.

A nanosecond molecular dynamics study of antiparallel d(G)(7) quadruplex structures: Effect of the coordinated cations

Nanosecond scale molecular dynamics simulations have been performed on antiparallel Greek key type d(G(7)) quadruplex structures with different coordinated ions, namely Na+ and K+ ion, water and Na+ counter ions, using the AMBER force field and Particle Mesh Ewald technique for electrostatic interactions. Antiparallel structures are stable during the simulation, with root mean square deviation values of similar to1.5 Angstrom from the initial structures. Hydrogen bonding patterns within the G-tetrads depend on the nature of the coordinated ion, with the G-tetrad undergoing local structural variation to accommodate different cations. However, alternating syn-anti arrangement of bases along a chain as well as in a quartet is maintained through out the MD simulation. Coordinated Na+ ions, within the quadruplex cavity are quite mobile within the central channel and can even enter or exit from the quadruplex core, whereas coordinated K+ ions are quite immobile. MD studies at 400 K indicate that K+ ion cannot come out from the quadruplex core without breaking the terminal G-tetrads. Smaller grooves in antiparallel structures are better binding sites for hydrated counter ions, while a string of hydrogen bonded water molecules are observed within both the small and large grooves. The hydration free energy for the K+ ion coordinated structure is more favourable than that for the Na+ ion coordinated antiparallel quadruplex structure.