Dhananjay Bhattacharyya

NUPARM and NUCGEN: software for analysis and generation of sequence dependent nucleic acid structures

Software packages NUPARM and NUCGEN, are described, which can be used to understand sequence directed structural variations in nucleic acids, by analysis and generation of non-uniform structures. A set of local inter basepair parameters (viz. tilt, roll, twist, shift, slide and rise) have been defined, which use geometry and coordinates of two successive basepairs only and can be used to generate polymeric structures with varying geometries for each of the 16 possible dinucleotide steps. Intra basepair parameters, propeller, buckle, opening and the C6...C8 distance can also be varied, if required, while the sugar phosphate backbone atoms are fixed in some standard conformation in each of the nucleotides. NUPARM can be used to analyse both DNA and RNA structures, with single as well as double stranded helices. The NUCGEN software generates double helical models with the backbone fixed in B-form DNA, but with appropriate modifications in the input data, it can also generate A-form DNA and RNA duplex 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)

Structural basis of DNA recognition by anticancer antibiotics, chromomycin A3, and mithramycin: Roles of minor groove width and ligand flexibility

Anticancer antibiotics, chromomycin A(3) (CHR) and mithramycin (MTR), inhibit cellular processes like transcription and replication, by binding reversibly to double-stranded DNA via minor groove, in the presence of bivalent metal ions like Mg(2+) with GC base specificity. Here, we have attempted to assess the roles of two parameters-namely DNA groove dimension and flexibility of the ligand-in the structural recognition between the ligands, (drug)(2)Mg(2+) and DNA. For the purpose we have employed three synthetic oligonucleotides with minor groove width lying between B- and A-type structures as model DNA sequences: d(GCGCGCGC)(2) in B-form, d(CCGGCGCCGG)(2) in B-form with unusual wide minor groove, and (GGGGCCCC)(2) in A-form. Association of the (drug)(2)Mg(2+) with the oligomers have been probed using spectroscopic techniques like absorbance, fluorescence, and CD. The binding and thermodynamic parameters for the different association processes have also been characterized. Major conclusions from the above studies are as follows. Groove size of the oligomers influences the conformation of the bound ligand. A saccharide dependent variation in structural rigidity of the ligands, (MTR)(2)Mg(2+) and (CHR)(2)Mg(2+), has been observed that leads to differences in the energetics of recognition of the same DNA sequence by the two ligands. In contrast to (CHR)(2)Mg(2+), higher flexibility in (MTR)(2)Mg(2+) makes its conformation in the DNA bound form less sensitive to the groove dimension of DNA.

Theoretical analysis of noncanonical base pairing interactions in RNA molecules

Noncanonical base pairs in RNA have strong structural and functional implications but are currently not considered for secondary structure predictions. We present results of comparative ab initio studies of stabilities and interaction energies for the three standard and 24 selected unusual RNA base pairs reported in the literature. Hydrogen added models of isolated base pairs, with heavy atoms frozen in their ‘away from equilibrium’ geometries, built from coordinates extracted from NDB, were geometry optimized using HF/6-31G** basis set, both before and after unfreezing the heavy atoms. Interaction energies, including BSSE and deformation energy corrections, were calculated, compared with respective single point MP2 energies, and correlated with occurrence frequencies and with types and geometries of hydrogen bonding interactions. Systems having two or more N-H...O/N hydrogen bonds had reasonable interaction energies which correlated well with respective occurrence frequencies and highlighted the possibility of some of them playing important roles in improved secondary structure prediction methods. Several of the remaining base pairs with one N-H...O/N and/or one C-H...O/N interactions respectively, had poor interaction energies and negligible occurrences. High geometry variations on optimization of some of these were suggestive of their conformational switch like characteristics.

Groove Width and Depth of B-DNA Structures Depend on Local Variation in Slide

The groove widths of DNA helix, especially minor groove width, are generally believed to be important for recognition of DNA by various types of ligands. It has been postulated earlier that large negative propeller twist, in the AT rich regions compresses the minor groove of duplex DNA. A systematic study has now been carried out by generating models with different values of local doublet and intra-basepair parameters and calculating their minor groove widths. It is found that several local doublet parameters affect the minor groove width but it depends most strongly on the local step parameters roll and slide when each parameter is considered individually. However, a detailed analysis of the various local parameters within the B-DNA family of crystal structures indicates that propeller twist and slide are most strongly correlated with the observed values of minor groove width. The groove depth is also strongly correlated with slide. Thus the local base sequence dependent variations in slide can modify both the groove width and depth and consequently determine the ligand binding properties of DNA.

Conformational specificity of non-canonical base pairs and higher order structures in nucleic acids: crystal structure database analysis

Non-canonical base pairs contribute immensely to the structural and functional variability of RNA, which calls for a detailed characterization of their spatial conformation. Intra-base pair parameters, namely propeller, buckle, open-angle, stagger, shear and stretch describe structure of base pairs indicating planarity and proximity of association between the two bases. In order to study the conformational specificities of non-canonical base pairs occurring in RNA crystal structures, we have upgraded NUPARM software to calculate these intra-base pair parameters using a new base pairing edge specific axis system. Analysis of base pairs and base triples with the new edge specific axis system indicate the presence of specific structural signatures for different classes of non-canonical pairs and triples. Differentiating features could be identified for pairs in cis or trans orientation, as well as those involving sugar edges or C–H-mediated hydrogen bonds. It was seen that propeller for all types of base pairs in cis orientation are generally negative, while those for trans base pairs do not have any preference. Formation of a base triple is seen to reduce propeller of the associated base pair along with reduction of overall flexibility of the pairs. We noticed that base pairs involving sugar edge are generally more non-planar, with large propeller or buckle values, presumably to avoid steric clash between the bulky sugar moieties. These specific conformational signatures often provide an insight into their role in the structural and functional context of RNA.

Structure, stability, and dynamics of canonical and noncanonical base pairs: quantum chemical studies.

The importance of non-Watson-Crick base pairs in the three-dimensional structure of RNA is now well established. The structure and stability of these noncanonical base pairs are, however, poorly understood. We have attempted to understand structural features of 33 frequently occurring base pairs using density functional theory. These are of three types, namely (i) those stabilized by two or more polar hydrogen bonds between the bases, (ii) those having one polar and another C-H...O/N type interactions, and (iii) those having one H-bond between the bases and another involving one of the sugars linked to the bases. We found that the base pairs having two polar H-bonds are very stable as compared to those having one C-H...O/N interaction. Our quantitatively analysis of structures of these optimized base pairs indicates that they possess a different amount of nonplanarity with large propeller or buckle values as also observed in the crystal structures. We further found that geometry optimization does not modify the hydrogen-bonding pattern, as values of shear and open angle of the base pairs remain conserved. The structures of initial crystal geometry and final optimized geometry of some base pairs having only one polar H-bond and a C-H...O/N interaction, however, are significantly different, indicating the weak nature of the nonpolar interaction. The base pair flexibility, as measured from normal-mode analysis, in terms of the intrinsic standard deviations of the base pair structural parameters are in conformity with those calculated from RNA crystal structures. We also noticed that deformation of a base pair along the stretch direction is impossible for all of the base pairs, and movements of the base pairs along shear and open are also quite restricted. The base pair opening mode through alteration of propeller or buckle is considerably less restricted for most of the base pairs.

Unzipping and binding of small interfering RNA with single walled carbon nanotube: a platform for small interfering RNA delivery.

In an effort to design efficient platform for siRNA delivery, we combine all atom classical and quantum simulations to study the binding of small interfering RNA (siRNA) by pristine single wall carbon nanotube (SWCNT). Our results show that siRNA strongly binds to SWCNT surface via unzipping its base-pairs and the propensity of unzipping increases with the increase in the diameter of the SWCNTs. The unzipping and subsequent wrapping events are initiated and driven by van der Waals interactions between the aromatic rings of siRNA nucleobases and the SWCNT surface. However, molecular dynamics (MD) simulations of double strand DNA (dsDNA) of the same sequence show that the dsDNA undergoes much less unzipping and wrapping on the SWCNT in the simulation time scale of 70 ns. This interesting difference is due to smaller interaction energy of thymidine of dsDNA with the SWCNT compared to that of uridine of siRNA, as calculated by dispersion corrected density functional theory (DFT) methods. After the optimal binding of siRNA to SWCNT, the complex is very stable which serves as one of the major mechanisms of siRNA delivery for biomedical applications. Since siRNA has to undergo unwinding process with the effect of RNA-induced silencing complex, our proposed delivery mechanism by SWCNT possesses potential advantages in achieving RNA interference.

Non-Canonical Base Pairs and Higher Order Structures in Nucleic Acids: Crystal Structure Database Analysis

Abstract Non-canonical base pairs, mostly present in the RNA, often play a prominent role towards maintaining their structural diversity. Higher order structures like base triples are also important in defining and stabilizing the tertiary folded structure of RNA. We have developed a new program BPFIND to analyze different types of canonical and non-canonical base pairs and base triples involving at least two direct hydrogen bonds formed between polar atoms of the bases or sugar O2′ only. We considered 104 possible types of base pairs, out of which examples of 87 base pair types are found to occur in the available RNA crystal structures. Analysis indicates that approximately 32.7% base pairs in the functional RNA structures are non-canonical, which include different types of GA and GU Wobble base pairs apart from a wide range of base pair possibilities. We further noticed that more than 10.4% of these base pairs are involved in triplet formation, most of which play important role in maintaining long-range tertiary contacts in the three-dimensional folded structure of RNA. Apart from detection, the program also gives a quantitative estimate of the conformational deformation of detected base pairs in comparison to an ideal planar base pair. This helps us to gain insight into the extent of their structural variations and thus assists in understanding their specific role towards structural and functional diversity.

RNA structure and dynamics: A base pairing perspective

RNA is now known to possess various structural, regulatory and enzymatic functions for survival of cellular organisms. Functional RNA structures are generally created by three-dimensional organization of small structural motifs, formed by base pairing between self-complementary sequences from different parts of the RNA chain. In addition to the canonical Watson-Crick or wobble base pairs, several non-canonical base pairs are found to be crucial to the structural organization of RNA molecules. They appear within different structural motifs and are found to stabilize the molecule through long-range intra-molecular interactions between basic structural motifs like double helices and loops. These base pairs also impart functional variation to the minor groove of A-form RNA helices, thus forming anchoring site for metabolites and ligands. Non-canonical base pairs are formed by edge-to-edge hydrogen bonding interactions between the bases. A large number of theoretical studies have been done to detect and analyze these non-canonical base pairs within crystal or NMR derived structures of different functional RNA. Theoretical studies of these isolated base pairs using ab initio quantum chemical methods as well as molecular dynamics simulations of larger fragments have also established that many of these non-canonical base pairs are as stable as the canonical Watson-Crick base pairs. This review focuses on the various structural aspects of non-canonical base pairs in the organization of RNA molecules and the possible applications of these base pairs in predicting RNA structures with more accuracy.