Genome analysis: A new approach for visualization of sequence organization in genomes

In this article we describe and demonstrate the versatility of a computer program, GENOME MAPPING, that uses interactive graphics and runs on an IRIS workstation. The program helps to visualize as well as analyse global and local patterns of genomic DNA sequences. It was developed keeping in mind the requirements of the human genome sequencing programme, which requires rapid analysis of the data. Using GENOME MAPPING one can discern signature patterns of different kinds of sequences and analyse such patterns for repetitive as well as rare sequence strings. Further, one can visualize the extent of global homology between different genomic sequences. An application of our method to the published yeast mitochondrial genome data shows similar sequence organizations in the entire sequence and in smaller subsequences.

Geometry of proline-containing alpha-helices in proteins

Crystal structure analysis of proline-containing alpha-helices in proteins has been carried out. High resolution crystal structures were selected from the Protein Data Bank. Apart from the standard internal parameters, some parameters which are specifically related to the bend in the helix due to proline have been developed and analyzed. Finally the position and nature of these helices and their interactions with the rest of the protein have been analyzed.

Unveiling the Role of Dps in the Organization of Mycobacterial Nucleoid

In order to preserve genetic information in stress conditions, bacterial DNA is organized into higher order nucleoid structure. In this paper, with the help of Atomic Force Microscopy, we show the different structural changes in mycobacterial nucleoid at different points of growth in the presence of different concentrations of glucose in the medium. We also observe that in Mycobacterium smegmatis, two different Dps proteins (Dps1 and Dps2) promote two types of nucleoid organizations. At the late stationary phase, under low glucose availability, Dps1 binds to DNA to form a very stable toroid structure. On the other hand, under the same condition, Dps2-DNA complex forms an incompletely condensed toroid and finally forms a further stable coral reef structure in the presence of RNA. This coral reef structure is stable in high concentration of bivalent ion like Mg2+.

Diversity in Functional Organization of Class I and Class II Biotin Protein Ligase

The cell envelope of Mycobacterium tuberculosis (M. tuberculosis) is composed of a variety of lipids including mycolic acids, sulpholipids, lipoarabinomannans, etc., which impart rigidity crucial for its survival and pathogenesis. Acyl CoA carboxylase (ACC) provides malonyl-CoA and methylmalonyl-CoA, committed precursors for fatty acid and essential for mycolic acid synthesis respectively. Biotin Protein Ligase (BPL/BirA) activates apo-biotin carboxyl carrier protein (BCCP) by biotinylating it to an active holo-BCCP. A minimal peptide (Schatz), an efficient substrate for Escherichia coli BirA, failed to serve as substrate for M. tuberculosis Biotin Protein Ligase (MtBPL). MtBPL specifically biotinylates homologous BCCP domain, MtBCCP87, but not EcBCCP87. This is a unique feature of MtBPL as EcBirA lacks such a stringent substrate specificity. This feature is also reflected in the lack of self/promiscuous biotinylation by MtBPL. The N-terminus/HTH domain of EcBirA has the selfbiotinable lysine residue that is inhibited in the presence of Schatz peptide, a peptide designed to act as a universal acceptor for EcBirA. This suggests that when biotin is limiting, EcBirA preferentially catalyzes, biotinylation of BCCP over selfbiotinylation. R118G mutant of EcBirA showed enhanced self and promiscuous biotinylation but its homologue, R69A MtBPL did not exhibit these properties. The catalytic domain of MtBPL was characterized further by limited proteolysis. Holo-MtBPL is protected from proteolysis by biotinyl-59 AMP, an intermediate of MtBPL catalyzed reaction. In contrast, apo-MtBPL is completely digested by trypsin within 20 min of co-incubation. Substrate selectivity and inability to promote self biotinylation are exquisite features of MtBPL and are a consequence of the unique molecular mechanism of an enzyme adapted for the high turnover of fatty acid biosynthesis.

Effect of Substrate Binding Loop Mutations on the Structure, Kinetics, and Inhibition of Enoyl Acyl Carrier Protein Reductase from Plasmodium falciparum

Enoyl acyl carrier protein reductase (ENR), which catalyzes the final and rate limiting step of fatty acid elongation, has been validated as a potential drug target. Triclosan is known to be an effective inhibitor for this enzyme. We mutated the substrate binding site residue Ala372 of the ENR of Plasmodium falciparum (PfENR) to Methionine and Valine which increased the affinity of the enzyme towards triclosan to almost double, close to that of Escherichia coli ENR (EcENR) which has a Methionine at the structurally similar position of Ala372 of PfENR. Kinetic studies of the mutants of PfENR and the crystal structure analysis of the A372M mutant revealed that a more hydrophobic environment enhances the affinity of the enzyme for the inhibitor. A triclosan derivative showed a threefold increase in the affinity towards the mutants compared to the wild type, due to additional interactions with the A372M mutant as revealed by the crystal structure. The enzyme has a conserved salt bridge which stabilizes the substrate binding loop and appears to be important for the active conformation of the enzyme. We generated a second set of mutants to check this hypothesis. These mutants showed loss of function, except in one case, where the crystal structure showed that the substrate binding loop is stabilized by a water bridge network. (C) 2011 IUBMB mum Life, 63(1): 30-41,2011

Modification of the sugar specificity of a plant lectin: structural studies on a point mutant of Erythrina corallodendron lectin

A mutant of Erythrina corallodendron lectin was generated with the aim of enhancing its affinity for N-acetylgalactosamine. A tyrosine residue close to the binding site of the lectin was mutated to a glycine in order to facilitate stronger interactions between the acetamido group of the sugar and the lectin which were prevented by the side chain of the tyrosine in the wild-type lectin. The crystal structures of this Y106G mutant lectin in complex with galactose and N-acetylgalactosamine have been determined. A structural rationale has been provided for the differences in the relative binding affinities of the wild-type and mutant lectins towards the two sugars based on the structures. A hydrogen bond between the O6 atom of the sugars and the variable loop of the carbohydrate-binding site of the lectin is lost in the mutant complexes owing to a conformational change in the loop. This loss is compensated by an additional hydrogen bond that is formed between the acetamido group of the sugar and the mutant lectin in the complex with N-acetylgalactosamine, resulting in a higher affinity of the mutant lectin for N-acetylgalactosamine compared with that for galactose, in contrast to the almost equal affinity of the wild-type lectin for the two sugars. The structure of a complex of the mutant with a citrate ion bound at the carbohydrate-binding site that was obtained while attempting to crystallize the complexes with sugars is also presented.

CAITG 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.

Characterization of substrate UpA binding to RNase A--computer modelling and energetics approach.

In the past two decades RNase A has been the focus of diverse investigations in order to understand the nature of substrate binding and to know the mechanism of enzyme action. Although this system is reasonably well characterized from the view point of some of the binding sites, the details of interactions in the second base binding (B2) site is insufficient. Further, the nature of ligand-protein interaction is elucidated generally by studies on RNase A-substrate analog complexes (mainly with the help of X-ray crystallography). Hence, the details of interactions at atomic level arising due to substrates are inferred indirectly. In the present paper, the dinucleotide substrate UpA is fitted into the active site of RNase A Several possible substrate conformations are investigated and the binding modes have been selected based on Contact Criteria. Thus identified RNase A-UpA complexes are energy minimized in coordinate space and are analysed in terms of conformations, energetics and interactions. The best possible ligand conformations for binding to RNase A are identified by experimentally known interactions and by the energetics. Upon binding of UpA to RNase A the changes associated,with protein back bone, Side chains in general and at the binding sites in particular are described. Further, the detailed interactions between UpA and RNase A are characterized in terms of hydrogen bonds and energetics. An extensive study has helped in interpreting the diverse results obtained from a number of experiments and also in evaluating the extent of changes the protein and the substrate undergo in order to maximize their interactions.

Conformation of polymerizable diacetylenic lecithin, DC8,9PC, by 1H-NMR spectroscopy

Diacetylenic phospholipid, 1,2 bis-(tricosa-10,12-diynoyl)-sn-glycero-3-phosphocholine (DC8,9PC), forms helices and tubules in addition to liposomes. The diacetylenic moiety responsible for the transformation is probed by 2-D NMR correlated spectroscopy. Chemical shift assignments and the analysis of 2D-COSY measurements were done on the lipid in chloroform-d solution. Based on this analysis, a model for the lipid is proposed. The geometry of the headgroup, glycerol backbone and acyl chains up to three methylenes from glycerol backbone [-(CH2)(3)-] is similar to that of dipalmitoyl phosphatidylcholine. The estimated torsional angle for methylene groups adjacent to diacetylenic moieties suggested an overall tilt of the diacetylenic lipid molecule from the bilayer axis of 25-30 degrees. This tilt could be negative or positive depending on the handedness of the resultant microstructures.

High level stable expression of rat brain type IIA sodium channel ?-subunit in CHO cells

The neuronal sodium channels are responsible for the rising phase of action potential and are composed of three subunits, of which the alpha-subunit has been shown to be adequate for most of its functional properties. We have stably expressed the rat brain type IIA sodium channel alpha-subunit in CHO cell tine using a CMV promoter-based vector. The expression was confirmed by detecting a 6.5 kb RNA corresponding to sodium channel alpha-subunit using Northern hybridization. The cells stably expressing the alpha-subunit, yield isolated sodium currents of amplitudes greater than 4nA when studied in whole-cell configuration of the patch-clamp technique. The sodium currents are characterized by activation and inactivation properties similar to neuronal sodium channels, and are blocked by the voltage gated sodium channel blocker tetrodotoxin (TTX).

Identification of Local Conformational Similarity in Structurally Variable Regions of Homologous Proteins Using Protein Blocks

Structure comparison tools can be used to align related protein structures to identify structurally conserved and variable regions and to infer functional and evolutionary relationships. While the conserved regions often superimpose well, the variable regions appear non superimposable. Differences in homologous protein structures are thought to be due to evolutionary plasticity to accommodate diverged sequences during evolution. One of the kinds of differences between 3-D structures of homologous proteins is rigid body displacement. A glaring example is not well superimposed equivalent regions of homologous proteins corresponding to a-helical conformation with different spatial orientations. In a rigid body superimposition, these regions would appear variable although they may contain local similarity. Also, due to high spatial deviation in the variable region, one-to-one correspondence at the residue level cannot be determined accurately. Another kind of difference is conformational variability and the most common example is topologically equivalent loops of two homologues but with different conformations. In the current study, we present a refined view of the ``structurally variable'' regions which may contain local similarity obscured in global alignment of homologous protein structures. As structural alphabet is able to describe local structures of proteins precisely through Protein Blocks approach, conformational similarity has been identified in a substantial number of `variable' regions in a large data set of protein structural alignments; optimal residue-residue equivalences could be achieved on the basis of Protein Blocks which led to improved local alignments. Also, through an example, we have demonstrated how the additional information on local backbone structures through protein blocks can aid in comparative modeling of a loop region. In addition, understanding on sequence-structure relationships can be enhanced through our approach. This has been illustrated through examples where the equivalent regions in homologous protein structures share sequence similarity to varied extent but do not preserve local structure.