Conformational Polymorphism in G-tetraplex Structures: Strand Reversal by Base Flipover or Sugar Flipover

Guanine rich sequences adopt a variety of four stranded structures, which differ in strand orientation and conformation about the glycosidic bond even though they are all stabilised by Hoogsteen hydrogen bonded guanine tetrads. Detailed model building and molecular mechanics calculations have been carried out to investigate various possible conformations of guanines along a strand and different possible orientations of guanine strands in a G-tetraplex structure. It is found that for an oligo G stretch per se, a parallel four stranded structure with all guanines in anti conformation is favoured over other possible tetraplex structures. Hence an alternating syn-anti arrangement of guanines along a strand is likely to occur only in folded back tetraplex structures with antiparallel G strands. Our study provides a theoretical rationale for the observed alternation of glycosidic conformation and the inverted stacking arrangement arising from base flipover, in antiparallel G-tetraplex structures and also highlights the various structural features arising due to different types of strand orientations. The molecular mechanics calculations help in elucidating the various interactions which stabilize different G-tetraplex structures and indicate that screening of phosphate charge by counterions could have a dramatic effect on groove width in these four stranded structures.

Conformational polymorphism in G-tetraplex structures: strand reversal by base flipover or sugar flipover

Guanlne rich sequences adopt a variety of four stranded structures, which differ in strand orientation and conformation about the glycosldic bond even though they are all stabilised by Hoogsteen hydrogen bonded guanlne tetrads. Detailed model building and molecular mechanics calculations have been carried out to investigate various possible conformations of guanlnes along a strand and different possible orientations of guanlne strands In a G-tetraplex structure. It is found that for an ollgo G stretch per se, a parallel four stranded structure with all guanines In anti conformation is favoured over other possible tetraplex structures. Hence an alternating syn-anti arrangement of guanlnes along a strand is likely to occur only in folded back tetraplex structures with antiparallel G strands. Our study provides a theoretical rationale for the observed alternation of glycosldic conformation and the inverted stacking arrangement arising from base filpover, In antlparallel G-tetraplex structures and also highlights the various structural features arising due to different types of strand orientations. The molecular mechanics calculations help in elucidating the various interactions which stabilize different G-tetraplex structures and indicate that screening of phosphate charge by counterions could have a dramatic effect on groove width in these four stranded structures.

Crucial contribution of the multiple copies of the initiator tRNA genes in the fidelity of tRNA(fMet) selection on the ribosomal P-site in Escherichia coli

The accuracy of the initiator tRNA (tRNA(fMet)) selection in the ribosomal P-site is central to the fidelity of protein synthesis. A highly conserved occurrence of three consecutive G-C base pairs in the anticodon stem of tRNA(fMet) contributes to its preferential selection in the P-site. In a genetic screen, using a plasmid borne copy of an inactive tRNA(fMet) mutant wherein the three G-C base pairs were changed, we isolated Escherichia coli strains that allow efficient initiation with the tRNA(fMet) mutant. Here, extensive characterization of two such strains revealed novel mutations in the metZWV promoter severely compromising tRNA(fMet) levels. Low cellular abundance of the chromosomally encoded tRNA(fMet) allows efficient initiation with the tRNA(fMet) mutant and an elongator tRNA(Gln), revealing that a high abundance of the cellular tRNA(fMet) is crucial for the fidelity of initiator tRNA selection on the ribosomal P-site in E. coli. We discuss possible implications of the changes in the cellular tRNA(fMet) abundance in proteome remodeling.

Formation of a G-quadruplex at the BCL2 major breakpoint region of the t(14;18) translocation in follicular lymphoma

The t(14;18) translocation in follicular lymphoma is one of the most common chromosomal translocations. Most breaks on chromosome 18 are located at the 3'-UTR of the BCL2 gene and are mainly clustered in the major breakpoint region (MBR). Recently, we found that the BCL2 MBR has a non-B DNA character in genomic DNA. Here, we show that single-stranded DNA modeled from the template strand of the BCL2 MBR, forms secondary structures that migrate faster on native PAGE in the presence of potassium, due to the formation of intramolecular G-quadruplexes. Circular dichroism shows evidence for a parallel orientation for G-quadruplex structures in the template strand of the BCL2 MBR. Mutagenesis and the DMS modification assay confirm the presence of three guanine tetrads in the structure. 1H nuclear magnetic resonance studies further confirm the formation of an intramolecular G-quadruplex and a representative model has been built based on all of the experimental evidence. We also provide data consistent with the possible formation of a G-quadruplex structure at the BCL2 MBR within mammalian cells. In summary, these important features could contribute to the single-stranded character at the BCL2 MBR, thereby contributing to chromosomal fragility.

Formation of a G-quadruplex at the BCL2 major breakpoint region of the t(14;18) translocation in follicular lymphoma

The t(14;18) translocation in follicular lymphoma is one of the most common chromosomal translocations. Most breaks on chromosome 18 are located at the 3'-UTR of the BCL2 gene and are mainly clustered in the major breakpoint region (MBR). Recently, we found that the BCL2 MBR has a non-B DNA character in genomic DNA. Here, we show that single-stranded DNA modeled from the template strand of the BCL2 MBR, forms secondary structures that migrate faster on native PAGE in the presence of potassium, due to the formation of intramolecular G-quadruplexes. Circular dichroism shows evidence for a parallel orientation for G-quadruplex structures in the template strand of the BCL2 MBR. Mutagenesis and the DMS modification assay confirm the presence of three guanine tetrads in the structure. 1H nuclear magnetic resonance studies further confirm the formation of an intramolecular G-quadruplex and a representative model has been built based on all of the experimental evidence. We also provide data consistent with the possible formation of a G-quadruplex structure at the BCL2 MBR within mammalian cells. In summary, these important features could contribute to the single-stranded character at the BCL2 MBR, thereby contributing to chromosomal fragility.

WebGeSTer DB-a transcription terminator database

We present WebGeSTer DB, the largest database of intrinsic transcription terminators (http://pallab.serc.iisc.ernet.in/gester). The database comprises of a million terminators identified in 1060 bacterial genome sequences and 798 plasmids. Users can obtain both graphic and tabular results on putative terminators based on default or user-defined parameters. The results are arranged in different tiers to facilitate retrieval, as per the specific requirements. An interactive map has been incorporated to visualize the distribution of terminators across the whole genome. Analysis of the results, both at the whole-genome level and with respect to terminators downstream of specific genes, offers insight into the prevalence of canonical and non-canonical terminators across different phyla. The data in the database reinforce the paradigm that intrinsic termination is a conserved and efficient regulatory mechanism in bacteria. Our database is freely accessible.

Inefficient excision of uracil from loop regions of DNA oligomers by E.coli uracil DNA glycosylase

Kinetic parameters for uracil DNA glycosylase (E.coli)-catalysed excision of uracil from DNA oligomers containing dUMP in different structural contexts were determined. Our results show that single-stranded oligonucleotides (unstructured) are used as somewhat better substrates than the double-stranded oligonucleotides. This is mainly because of the favourable V-max value of the enzyme for single-stranded substrates. More interestingly, however, we found that uracil release from loop regions of DNA hairpins is extremely inefficient. The poor efficiency with which uracil is excised from loop regions is a result of both increased K-m and lowered V-max values. This observation may have significant implications in uracil DNA glycosylase-directed repair of DNA segments that can be extruded as hairpins. In addition, these studies are useful in designing oligonucleotides for various applications in DNA research where the use of uracil DNA glycosylase is sought.

A perspective of biological supramolecular electron transfer

Electron transfer is an essential activity in biological systems. The migrating electron originates from water-oxygen in photosynthesis and reverts to dioxygen in respiration. In this cycle two metal porphyrin complexes possessing circular conjugated system and macrocyclic pi-clouds, chlorophyll and hems, play a decisive role in mobilising electrons for travel over biological structures as extraneous electrons. Transport of electrons within proteins (as in cytochromes) and within DNA (during oxidative damage and repair) is known to occur. Initial evaluations did not favour formation of semiconducting pathways of delocalized electrons of the peptide bonds in proteins and of the bases in nucleic acids. Direct measurement of conductivity of bulk material and quantum chemical calculations of their polymeric structures also did not support electron transfer in both proteins and nucleic acids. New experimental approaches have revived interest in the process of charge transfer through DNA duplex. The fluorescence on photoexcitation of Ru-complex was found to be quenched by Rh-complex, when both were tethered to DNA and intercalated in the base stack. Similar experiments showed that damage to G-bases and repair of T-T dimers in DNA can occur by possible long range electron transfer through the base stack. The novelty of this phenomenon prompted the apt name, chemistry at a distance. Based on experiments with ruthenium modified proteins, intramolecular electron transfer in proteins is now proposed to use pathways that include C-C sigma-bonds and surprisingly hydrogen bonds which remained out of favour for a long time. In support of this, some experimental evidence is now available showing that hydrogen bond-bridges facilitate transfer of electrons between metal-porphyrin complexes. By molecular orbital calculations over 20 years ago. we found that "delocalization of an extraneous electron is pronounced when it enters low-lying virtual orbitals of the electronic structures of peptide units linked by hydrogen bonds". This review focuses on supramolecular electron transfer pathways that can emerge on interlinking by hydrogen bonds and metal coordination of some unnoticed structures with pi-clouds in proteins and nucleic acids, potentially useful in catalysis and energy missions.

Sex-specific organisation of middle repetitive DNA sequences in the mealybug Planococcus lilacinus

Differential organisation of homologous chromosomes is related to both sex determination and genomic imprinting in coccid insects, the mealybugs. We report here the identification of two middle repetitive sequences that are differentially organised between the two sexes and also within the same diploid nucleus. These two sequences form a part of the male-specific nuclease-resistant chromatin (NRC) fraction of a mealybug Planococcus lilacinus. To understand the phenomenon of differential organisation we have analysed the components of NRC by cloning the DNA sequences present, deciphering their primary sequence, nucleosomal organisation, genomic distribution and cytological localisation, Our observations suggest that the middle repetitive sequences within NRC are functionally significant and we discuss their probable involvement in male-specific chromatin organisation.

Structural characterisation of a uracil containing hairpin DNA by NMR and molecular dynamics

Three-dimensional (3D) structure of a hairpin DNA d-CTAGAGGATCCTTTUGGATCCT (22mer; abbreviated as U4-hairpin), which has a uracil nucleotide unit at the fourth position from the 5' end of the tetra-loop has been solved by NMR spectroscopy. The H-1 resonances of this hairpin have been assigned almost completely. NMR restrained molecular dynamics and energy minimisation procedures have been used to describe the 3D structure of the U4 hairpin. This study establishes that the stem of the hairpin adopts a right handed B-DNA conformation while the T-12 and U-15 nucleotide stack upon 3' and 5' ends of the stem, respectively. Further, T-14 stacks upon both T-12 and U-15 while T-13 partially stacks upon T-14. Very weak stacking interaction is observed between T-13 and T-12. All the individual nucleotide bases adopt 'anti' conformation with respect to their sugar moiety. The turning phosphate in the loop is located between T-13 and T-14. The stereochemistry of U-15 mimics the situation where uracil would stack in a B-DNA conformation. This could be the reason as to why the U4-hairpin is found to be the best substrate for its interaction with uracil DNA glycosylase (UDG) compared to the other substrates in which the uracil is at the first, second and third positions of the tetra-loop from its 5' end, as reported previously.

iPBA: a tool for protein structure comparison using sequence alignment strategies

With the immense growth in the number of available protein structures, fast and accurate structure comparison has been essential. We propose an efficient method for structure comparison, based on a structural alphabet. Protein Blocks (PBs) is a widely used structural alphabet with 16 pentapeptide conformations that can fairly approximate a complete protein chain. Thus a 3D structure can be translated into a 1D sequence of PBs. With a simple Needleman-Wunsch approach and a raw PB substitution matrix, PB-based structural alignments were better than many popular methods. iPBA web server presents an improved alignment approach using (i) specialized PB Substitution Matrices (SM) and (ii) anchor-based alignment methodology. With these developments, the quality of similar to 88% of alignments was improved. iPBA alignments were also better than DALI, MUSTANG and GANGSTA(+) in > 80% of the cases. The webserver is designed to for both pairwise comparisons and database searches. Outputs are given as sequence alignment and superposed 3D structures displayed using PyMol and Jmol. A local alignment option for detecting subs-structural similarity is also embedded. As a fast and efficient `sequence-based' structure comparison tool, we believe that it will be quite useful to the scientific community. iPBA can be accessed at http://www.dsimb.inserm.fr/dsimb_tools/ipba/.

PRUNE and PROBE-two modular web services for protein-protein docking

The protein-protein docking programs typically perform four major tasks: (i) generation of docking poses, (ii) selecting a subset of poses, (iii) their structural refinement and (iv) scoring, ranking for the final assessment of the true quaternary structure. Although the tasks can be integrated or performed in a serial order, they are by nature modular, allowing an opportunity to substitute one algorithm with another. We have implemented two modular web services, (i) PRUNE: to select a subset of docking poses generated during sampling search (http://pallab.serc.iisc.ernet.in/prune) and (ii) PROBE: to refine, score and rank them (http://pallab.serc.iisc.ernet.in/probe). The former uses a new interface area based edge-scoring function to eliminate > 95% of the poses generated during docking search. In contrast to other multi-parameter-based screening functions, this single parameter based elimination reduces the computational time significantly, in addition to increasing the chances of selecting native-like models in the top rank list. The PROBE server performs ranking of pruned poses, after structure refinement and scoring using a regression model for geometric compatibility, and normalized interaction energy. While web-service similar to PROBE is infrequent, no web-service akin to PRUNE has been described before. Both the servers are publicly accessible and free for use.

DEPTH: a web server to compute depth and predict small-molecule binding cavities in proteins

Depth measures the extent of atom/residue burial within a protein. It correlates with properties such as protein stability, hydrogen exchange rate, protein-protein interaction hot spots, post-translational modification sites and sequence variability. Our server, DEPTH, accurately computes depth and solvent-accessible surface area (SASA) values. We show that depth can be used to predict small molecule ligand binding cavities in proteins. Often, some of the residues lining a ligand binding cavity are both deep and solvent exposed. Using the depth-SASA pair values for a residue, its likelihood to form part of a small molecule binding cavity is estimated. The parameters of the method were calibrated over a training set of 900 high-resolution X-ray crystal structures of single-domain proteins bound to small molecules (molecular weight < 1.5 KDa). The prediction accuracy of DEPTH is comparable to that of other geometry-based prediction methods including LIGSITE, SURFNET and Pocket-Finder (all with Matthew's correlation coefficient of similar to 0.4) over a testing set of 225 single and multi-chain protein structures. Users have the option of tuning several parameters to detect cavities of different sizes, for example, geometrically flat binding sites. The input to the server is a protein 3D structure in PDB format. The users have the option of tuning the values of four parameters associated with the computation of residue depth and the prediction of binding cavities. The computed depths, SASA and binding cavity predictions are displayed in 2D plots and mapped onto 3D representations of the protein structure using Jmol. Links are provided to download the outputs. Our server is useful for all structural analysis based on residue depth and SASA, such as guiding site-directed mutagenesis experiments and small molecule docking exercises, in the context of protein functional annotation and drug discovery.

Structural basis for poor uracil excision from hairpin DNA

Two-dimensional NMR and molecular dynamics simulations have been used to determine the three-dimensional structures of two hairpin DNA structures: d-CTAGAG GATCCUTTTGGATCCT (abbreviated as U1-hairpin) and d-CTAGAGGATCCTTUTGGATCCT (abbreviated as U3-hairpin). The (1) H resonances of both of these hairpin structures have been assigned almost completely. NMR restrained molecular dynamics and energy minimization procedures have been used to describe the three-dimensional structures of these hairpins. This study and concurrent NMR structural studies on two other d-CTAGAGGA TCCTUTTGGATCCT (abbreviated as U2-hairpin) and d-CTAGAGGATCCTTTUGGATCCT (abbreviated as U4-hairpin) have shed light upon various interactions reported between Echerichia coli uracil DNA glycosylase (UDG) and uracil-containing DNA. The backbone torsion angles, which partially influence the local conformation of U12 and U14 in U1 and U3-hairpins, respectively, are probably locked in the trans conformation as in the case of U-13 in the U2-hairpin. Such a stretched-out backbone conformation in the vicinity of U-12 and U-14 is thought to be the reason why the K-m value is poor for U1- and U3-hairpins as it is for the U2-hairpin. Furthermore, the bases U-12 and U-14 in both U1- and U3-hairpins adopt an anti conformation, in contrast with the base conformation of U-13 in the U2-hairpin, which adopts a syn conformation. The clear discrepancy observed in the U-base orientation with respect to the sugar moieties could explain why the V-max value is 10- to 20-fold higher for the U1- and U3-hairpins compared with the U2-hairpin. Taken together, these observations support our interpretation that the unfavourable backbone results in a poor K-m value, whereas the unfavourable nucleotide conformation results in a poor V-max value. These two parameters therefore make the U1- and U3-hairpins better substrates for UDG compared with the U2-hairpin, as reported earlier [Kumar, N. V. & Varshney, U. (1997) Nucleic Acids Res. 25, 2336-2343.].

How does DNA break during chromosomal translocations?

Chromosomal translocations are one of the most common types of genetic rearrangements and are molecular signatures for many types of cancers. They are considered as primary causes for cancers, especially lymphoma and leukemia. Although many translocations have been reported in the last four decades, the mechanism by which chromosomes break during a translocation remains largely unknown. In this review, we summarize recent advances made in understanding the molecular mechanism of chromosomal translocations.