High-quality annotation of promoter regions for 913 bacterial genomes

Motivation: The number of bacterial genomes being sequenced is increasing very rapidly and hence, it is crucial to have procedures for rapid and reliable annotation of their functional elements such as promoter regions, which control the expression of each gene or each transcription unit of the genome. The present work addresses this requirement and presents a generic method applicable across organisms. Results: Relative stability of the DNA double helical sequences has been used to discriminate promoter regions from non-promoter regions. Based on the difference in stability between neighboring regions, an algorithm has been implemented to predict promoter regions on a large scale over 913 microbial genome sequences. The average free energy values for the promoter regions as well as their downstream regions are found to differ, depending on their GC content. Threshold values to identify promoter regions have been derived using sequences flanking a subset of translation start sites from all microbial genomes and then used to predict promoters over the complete genome sequences. An average recall value of 72% (which indicates the percentage of protein and RNA coding genes with predicted promoter regions assigned to them) and precision of 56% is achieved over the 913 microbial genome dataset.

X-ray studies on crystalline complexes involving amino acids and peptides XIX. Effects of change in chirality in the complexes of succinic acid with dl- and l-arginine

Crystals of dl-arginine hemisuccinate dihydrate (I)(monoclinic; P21/c; a = 5.292, b = 16.296, c = 15.203 Å; α= 92.89°; Z = 4) and l-arginine hemisuccinate hemisuccinic acid monohydrate (II) (triclinic; P1; a = 5.099; b = 10.222, c = 14.626 Å; α= 77.31, β= 89.46, γ= 78.42°; Z = 2) were grown under identical conditions from aqueous solutions of the components in molar proportions. The structures were solved by direct methods and refined to R = 0.068 for 2585 observed reflections in the case of (I) and R = 0.036 for 2154 observed reflections in the case of (11). Two of the three crystallographically independent arginine molecules in the complexes have conformations different from those observed so far in the crystal structures containing arginine. The succinic acid molecules and the succinate ions in the structures are centrosymmetric and planar. The crystal structure of (II) is highly pseudosymmetric. Arginine-succinate interactions in both the complexes involve specific guanidyl-carboxylate interactions. The basic elements of aggregation in both the structures are ribbons made up of alternating arginine dimers and succinate ions. However, the ribbons pack in different ways in the two structures. (II) presents an interesting case in which two ionisation states of the same molecule coexist in a crystal. The two complexes provide a good example of the effect of change in chirality on stoichiometry, conformation, aggregation, and ionisation state in the solid state.

X-ray studies on crystalline complexes involving amino acids and peptides. XXI. Structure of a (1:1) complex between L-phenylalanine and D-valine

CsHllNO2.C9HilNO2, Mr = 282.3, P1, a = 5.245 (1), b = 5.424 (1), c = 14.414 (2) A, a = 97.86 (1), fl = 93-69 (2), y = 70-48 (2) °, V= 356 A 3, Z = 1, O m = 1-32 (2), Dx = 1.32 g cm-3, h(Mo Ka) = 0-7107 A, g = 5-9 cm-1, F(000) = 158, T= 298 K, R=0.035 for 1518 observed reflections with I>2tr(I). The molecules aggregate in double layers, one ayer made up of L-phenylalanine molecules and the other of D-valine molecules. Each double layer is stabilized by interactions involving main-chain atoms of both types of molecules. The interactions include hydrogen bonds which give rise to two head-to-tail sequences. The arrangement of molecules in the complex is almost the same as that in the structure of DL-valine (and DL-leucine and DL-isoleucine) except for the change in the side chain of L molecules. The molecules in crystals containing an equal number of L and O hydrophobic amino-acid molecules thus appear to aggregate in a similar fashion, irrespective of the precise details of the side chain.

Protein Hydration and Water Structure: X-ray Analysis of a Closely Packed Protein Crystal with Very Low Solvent Content

Low-humidity monoclinic lysozyme, resulting from a water-mediated transformation, has one of the lowest solvent contents (22% by volume) observed in a protein crystal. Its structure has been solved by the molecular replacement method and refined to an R value of 0.175 for 7684 observed reflections in the 10–1.75 Å resolution shell. 90% of the solvent in the well ordered crystals could be located. Favourable sites of hydration on the protein surface include side chains with multiple hydrogen-bonding centres, and regions between short hydrophilic side chains and the main-chain CO or NH groups of the same or nearby residues. Major secondary structural features are not disrupted by hydration. However, the free CO groups at the C terminii and, to a lesser extent, the NH groups at the N terminii of helices provide favourable sites for water interactions, as do reverse turns and regions which connect β-structure and helices. The hydration shell consists of discontinuous networks of water molecules, the maximum number of molecules in a network being ten. The substrate-binding cleft is heavily hydrated, as is the main loop region which is stabilized by water interactions. The protein molecules are close packed in the crystals with a molecular coordination number of 14. Arginyl residues are extensively involved in intermolecular hydrogen bonds and water bridges. The water molecules in the crystal are organized into discrete clusters. A distinctive feature of the clusters is the frequent occurrence of three-membered rings. The protein molecules undergo substantial rearrangement during the transformation from the native to the low-humidity form. The main-chain conformations in the two forms are nearly the same, but differences exist in the side-chain conformation. The differences are particularly pronounced in relation to Trp 62 and Trp 63. The shift in Trp 62 is especially interesting as it is also known to move during inhibitor binding.

X-ray studies on crystalline complexes involving amino acids and peptides. XXV. Structures of DL-proline hemisuccinic acid and glycyl-L-histidinium semisuccinate monohydrate and a comparative study of amino-acid and peptide complexes of succinic acid

DL-Proline hemisuccinic acid, C5H9NO2.1/2C4H6O4, M(r) = 174.2, P2(1/c) a = 5.254 (1), b = 17.480 (1), c = 10.230 (i) angstrom, beta = 119.60 (6)-degrees Z = 4, D(m) = 1.41 (4), D(x) = 1.42 g cm-3, R = 0.045 for 973 observed reflections. Glycyl-L-histidinium semisuccinate monohydrate, C8H13N4O3+.C4H5O4-.H2O, M(r) = 348.4, P2(1), a = 4.864 (1), b = 17.071 (2), c = 9.397 (1) angstrom, beta = 90.58-degrees, Z = 2, D(m) = 1.45 (1), D(x) = 1.48 g cm-3, R = 0.027 for 1610 observed reflections. Normal amino-acid and dipeptide aggregation patterns are preserved in the structures in spite of the presence of succinic acid/semisuccinate ions. In both the structures, the amino-acid/dipeptide layers stack in such a way that the succinic acid molecules/semisuccinate ions are enclosed in voids created during stacking. Substantial variability in the ionization state and the stoichiometry is observed in amino-acid and peptide complexes of succinic acid. Succinic acid molecules and succinate ions appear to prefer a planar centro-symmetric conformation with the two carboxyl (carboxylate) groups trans with respect to the central C=C bond. Considerable variation is seen in the departure from and modification of normal amino-acid aggregation patterns produced by the presence of succinic acid. Some of the complexes can be described as inclusion compounds with the amino acid/dipeptide as the 'host' and succinic acid/semisuccinate/succinate as the 'guest'. The effects of change in chirality, though very substantial, are not the same in different pairs of complexes involving DL and L isomers of the same amino acid.

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

Stereochemical studies on cyclic peptides. Part XI. Conformation of cyclic pentapeptides having intramolecular 3 leads to 1 hydrogen bonds.

Conformational analysis of cyclic pentapeptides having two intra-ring 3 leads to 1 hydrogen bonds has been carried out. It is found that the structure can easily be formed with trans planar peptide units without causing significant angular strain at the alpha-carbon atoms. Four different types of conformations designated Types I--IV are possible for the backbone structure. Details of these four types of conformations and also the accommodating possibility of these types for allglycyl and all-alanyl residues are presented. Three of the four types have relatively low energies for glycyl residues whereas the other one has a slightly higher energy. When alanyl residues are introduced at the five alpha-carbon atoms, the types that are energetically favourable depend upon the sequence of isomers. Energy calculations have also been carried out for the combinations of glycyl, L- and D-alanyl residues. The theoretical results are compared with available experimental observations both from solution and solid state studies.

Cloning, expression, purification, crystallization and preliminary X-ray studies of a secreted lectin (Rv1419) from Mycobacterium tuberculosis

A secreted lectin, Rv1419, from Mycobacterium tuberculosis has been cloned, expressed, purified and crystallized and the crystals have been characterized. This represents the first X-ray investigation of a lectin or lectin-like molecule from the pathogen. The cubic crystals contain one molecule in the asymmetric unit. Sequence comparisons indicate that the lectin has a beta-trefoil fold and belongs to a well characterized family of carbohydrate-binding modules. Structural analysis of the crystals is in progress.

A standard reference frame for the description of Nucleic Acid Base-Pair Geometry

A common point of reference is needed to describe the three-dimensional arrangements of bases and base-pairs in nucleic acid structures. The different standards used in computer programs created for this purpose give rise to con¯icting interpretations of the same structure.1 For example, parts of a structure that appear ``normal'' according to one computational scheme may be highly unusual according to another and vice versa. It is thus dif®cult to carry out comprehensive comparisons of nucleic acid structures and to pinpoint unique conformational features in individual structures

HELANAL: A program to characterise helix geometry in proteins,

A detailed analysis of structural and position dependent characteristic features of helices will give a better understanding of the secondary structure formation in globular proteins. Here we describe an algorithm that quantifies the geometry of helices in proteins on the basis of their C-alpha atoms alone. The Fortran program HELANAL can extract the helices from the PDB files and then characterises the overall geometry of each helix as being linear, curved or kinked, in terms of its local structural features, viz. local helical twist and rise, virtual torsion angle, local helix origins and bending angles between successive local helix axes. Even helices with large radius of curvature are unambiguously identified as being linear or curved. The program can also be used to differentiate a kinked helix and other motifs, such as helix-loop-helix or a helix-turn-helix (with a single residue linker) with the help of local bending angles. In addition to these, the program can also be used to characterise the helix start and end as well as other types of secondary structures.