Structural and Biochemical Bases for the Redox Sensitivity of Mycobacterium tuberculosis RslA

An effective transcriptional response to redox stimuli is of particular importance for Mycobacterium tuberculosis, as it adapts to the environment of host alveoli and macrophages. The M. tuberculosis a factor sigma(L) regulates the expression of genes involved in cell-wall and polyketide syntheses. sigma(L) interacts with the cytosolic anti-sigma domain of a membrane-associated protein, RslA. Here we demonstrate that RslA binds Zn2+ and can sequester sigma(L) in a reducing environment. In response to an oxidative stimulus, proximal cysteines in the CXXC motif of RslA form a disulfide bond, releasing bound Zn2+. This results in a substantial rearrangement of the sigma(L)/RslA complex, leading to an 8-fold decrease in the affinity of RslA for sigma(L). The crystal structure of the -35-element recognition domain of sigma(L), sigma(L)(4), bound to RslA reveals that RslA inactivates sigma(L) by sterically occluding promoter DNA and RNpolymerase binding sites. The crystal structure further reveals that the cysteine residues that coordinate Zn2+ in RslA are solvent exposed in the complex, thus providing a structural basis for the redox sensitivity of RslA. The biophysical parameters of sigma(L)/RslA interactions provide a template for understanding how variations in the rate of Zn2+ release and associated conformational changes could regulate the activity of a Zn2+-associated anti-sigma factor. (C) 2010 Elsevier Ltd. All rights reserved.

A structural basis for the role of nucleotide specifying residues in regulating the Rv1625c adenylyl cyclase oligomerization of the from M-tuberculosis

The Rv1625c Class III adenylyl cyclase from Mycobacterium tuberculosis is a homodimeric enzyme with two catalytic centers at the dimer interface, and shows sequence similarity with the mammalian adenylyl and guanylyl cyclases. Mutation of the substrate-specifying residues in the catalytic domain of Rv1625c, either independently or together, to those present in guanylyl cyclases not only failed to confer guanylyl cyclase activity to the protein, but also severely abrogated the adenylyl cyclase activity of the enzyme. Biochemical analysis revealed alterations in the behavior of the mutants on ion-exchange chromatography, indicating differences in the surface-exposed charge upon mutation of substrate-specifying residues. The mutant proteins showed alterations in oligomeric status as compared to the wild-type enzyme, and differing abilities to heterodimerize with the wild-type protein. The crystal structure of a mutant has been solved to a resolution of 2.7 angstrom. On the basis of the structure, and additional biochemical studies, we provide possible reasons for the altered properties of the mutant proteins, as well as highlight unique structural features of the Rv1625c adenylyl cyclase. (c) 2005 Elsevier Ltd. All rights reserved.

5-Fluorouracil induces structural changes in rat liver tRNA

5-fluorouracil (FUra) has been shown to modulate the aminoacylation function of rat liver tRNA. The present study was aimed at studying the structure-function relationship of FUra-substituted tRNA. Male Wistar rats (2-3 month old) were given a single i.p. injection of FUra at 50, 250, or 500 mg/kg body wt. and FUra-substituted total liver tRNA, i.e. tRNA(FUra50, 250, and 500, respectively, were isolated 3 h later. Normal tRNA (tRNA(N)) was isolated from saline-treated control rats. Thermal denaturation studies showed higher melting temperatures for tRNA(FUra) compared to tRNA(N). Heat denaturation followed by renaturation of total tRNA did not affect the activity of tRNA(N) and tRNA(FUra50), where as tRNA(FUra250 and 500) lost 35% and 72% of activity, respectively, compared to the corresponding group of non-denatured tRNA. Antibodies specific to rat liver tRNA recognized normal and FUra-substituted tRNA in the order of tRNA(N) > tRNA(FUra50) > or = tRNA(FUra250) > tRNA(FUra500) in an avidin-biotin micro-enzyme linked immunosorbant assay. tRNA(N) or tRNA(FUra50) preincubated with tRNA antiserum showed 74% and 59% of aminoacylation activity, respectively, compared to that of corresponding tRNA preincubated with normal rabbit IgG. However, activities of similarly treated tRNA(FUra250 and 500) were not affected. The observations of possible changes in the secondary structure of rat liver tRNA upon incorporation of FUra are discussed.

Structural proteins of mycobacteriophage I3: cloning, expression and sequence analysis of a gene encoding a 70-kDa structural protein

The structural proteins of mycobacteriophage I3 have been analysed by sodium dodecyl sulfate-polyacrylamide-gel electrophoresis (SDS-PAGE), radioiodination and immunoblotting. Based on their abundance the 34- and 70-kDa bands appeared to represent the major structural proteins. Successful cloning and expression of the 70-kDa protein-encoding gene of phage I3 in Escherichia coli and its complete nucleotide sequence determination have been accomplished, A second (partial) open reading frame following the stop codon for the 70-kDa protein was also identified within the cloned fragment. The deduced amino-acid sequence of the 70-kDa protein and the codon usage patterns indicated the preponderance of codons, as predicted from the high G+C content of the genomic DNA of phage I3.

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.

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

Structural Annotation of Mycobacterium tuberculosis Proteome

Of the similar to 4000 ORFs identified through the genome sequence of Mycobacterium tuberculosis (TB) H37Rv, experimentally determined structures are available for 312. Since knowledge of protein structures is essential to obtain a high-resolution understanding of the underlying biology, we seek to obtain a structural annotation for the genome, using computational methods. Structural models were obtained and validated for similar to 2877 ORFs, covering similar to 70% of the genome. Functional annotation of each protein was based on fold-based functional assignments and a novel binding site based ligand association. New algorithms for binding site detection and genome scale binding site comparison at the structural level, recently reported from the laboratory, were utilized. Besides these, the annotation covers detection of various sequence and sub-structural motifs and quaternary structure predictions based on the corresponding templates. The study provides an opportunity to obtain a global perspective of the fold distribution in the genome. The annotation indicates that cellular metabolism can be achieved with only 219 folds. New insights about the folds that predominate in the genome, as well as the fold-combinations that make up multi-domain proteins are also obtained. 1728 binding pockets have been associated with ligands through binding site identification and sub-structure similarity analyses. The resource (http://proline.physics.iisc.ernet.in/Tbstructuralannotation), being one of the first to be based on structure-derived functional annotations at a genome scale, is expected to be useful for better understanding of TB and for application in drug discovery. The reported annotation pipeline is fairly generic and can be applied to other genomes as well.

A sequential dual method for structural SVMs

In many real world prediction problems the output is a structured object like a sequence or a tree or a graph. Such problems range from natural language processing to compu- tational biology or computer vision and have been tackled using algorithms, referred to as structured output learning algorithms. We consider the problem of structured classifi- cation. In the last few years, large margin classifiers like sup-port vector machines (SVMs) have shown much promise for structured output learning. The related optimization prob -lem is a convex quadratic program (QP) with a large num-ber of constraints, which makes the problem intractable for large data sets. This paper proposes a fast sequential dual method (SDM) for structural SVMs. The method makes re-peated passes over the training set and optimizes the dual variables associated with one example at a time. The use of additional heuristics makes the proposed method more efficient. We present an extensive empirical evaluation of the proposed method on several sequence learning problems.Our experiments on large data sets demonstrate that the proposed method is an order of magnitude faster than state of the art methods like cutting-plane method and stochastic gradient descent method (SGD). Further, SDM reaches steady state generalization performance faster than the SGD method. The proposed SDM is thus a useful alternative for large scale structured output learning.

DNA STRUCTURAL FEATURES AND ARCHITECTURE OF PROMOTER REGIONS PLAY A ROLE IN GENE RESPONSIVENESS OF S. CEREVISIAE

Gene expression is the most fundamental biological process, which is essential for phenotypic variation. It is regulated by various external (environment and evolution) and internal (genetic) factors. The level of gene expression depends on promoter architecture, along with other external factors. Presence of sequence motifs, such as transcription factor binding sites (TFBSs) and TATA-box, or DNA methylation in vertebrates has been implicated in the regulation of expression of some genes in eukaryotes, but a large number of genes lack these sequences. On the other hand, several experimental and computational studies have shown that promoter sequences possess some special structural properties, such as low stability, less bendability, low nucleosome occupancy, and more curvature, which are prevalent across all organisms. These structural features may play role in transcription initiation and regulation of gene expression. We have studied the relationship between the structural features of promoter DNA, promoter directionality and gene expression variability in S. cerevisiae. This relationship has been analyzed for seven different measures of gene expression variability, along with two different regulatory effect measures. We find that a few of the variability measures of gene expression are linked to DNA structural properties, nucleosome occupancy, TATA-box presence, and bidirectionality of promoter regions. Interestingly, gene responsiveness is most intimately correlated with DNA structural features and promoter architecture.

Re-analysis of cryoEM data on HCV IRES bound to 40S subunit of human ribosome integrated with recent structural information suggests new contact regions between ribosomal proteins and HCV RNA

In this study, we combine available high resolution structural information on eukaryotic ribosomes with low resolution cryo-EM data on the Hepatitis C Viral RNA (IRES) human ribosome complex. Aided further by the prediction of RNA-protein interactions and restrained docking studies, we gain insights on their interaction at the residue level. We identified the components involved at the major and minor contact regions, and propose that there are energetically favorable local interactions between 40S ribosomal proteins and IRES domains. Domain II of the IRES interacts with ribosomal proteins S5 and S25 while the pseudoknot and the downstream domain IV region bind to ribosomal proteins S26, S28 and S5. We also provide support using UV cross-linking studies to validate our proposition of interaction between the S5 and IRES domains II and IV. We found that domain IIIe makes contact with the ribosomal protein S3a (S1e). Our model also suggests that the ribosomal protein S27 interacts with domain IIIc while S7 has a weak contact with a single base RNA bulge between junction IIIabc and IIId. The interacting residues are highly conserved among mammalian homologs while IRES RNA bases involved in contact do not show strict conservation. IRES RNA binding sites for S25 and S3a show the best conservation among related viral IRESs. The new contacts identified between ribosomal proteins and RNA are consistent with previous independent studies on RNA-binding properties of ribosomal proteins reported in literature, though information at the residue level is not available in previous studies.

An overview of recent advances in structural bioinformatics of protein-protein interactions and a guide to their principles

Rich data bearing on the structural and evolutionary principles of protein protein interactions are paving the way to a better understanding of the regulation of function in the cell. This is particularly the case when these interactions are considered in the framework of key pathways. Knowledge of the interactions may provide insights into the mechanisms of crucial `driver' mutations in oncogenesis. They also provide the foundation toward the design of protein protein interfaces and inhibitors that can abrogate their formation or enhance them. The main features to learn from known 3-D structures of protein protein complexes and the extensive literature which analyzes them computationally and experimentally include the interaction details which permit undertaking structure-based drug discovery, the evolution of complexes and their interactions, the consequences of alterations such as post-translational modifications, ligand binding, disease causing mutations, host pathogen interactions, oligomerization, aggregation and the roles of disorder, dynamics, allostery and more to the protein and the cell. This review highlights some of the recent advances in these areas, including design, inhibition and prediction of protein protein complexes. The field is broad, and much work has been carried out in these areas, making it challenging to cover it in its entirety. Much of this is due to the fast increase in the number of molecules whose structures have been determined experimentally and the vast increase in computational power. Here we provide a concise overview. (C) 2014 Elsevier Ltd. All rights reserved.

A Gaussian network model study suggests that structural fluctuations are higher for inactive states than active states of protein kinases

We performed Gaussian network model based normal mode analysis of 3-dimensional structures of multiple active and inactive forms of protein kinases. In 14 different kinases, a more number of residues (1095) show higher structural fluctuations in inactive states than those in active states (525), suggesting that, in general, mobility of inactive states is higher than active states. This statistically significant difference is consistent with higher crystallographic B-factors and conformational energies for inactive than active states, suggesting lower stability of inactive forms. Only a small number of inactive conformations with the DFG motif in the ``in'' state were found to have fluctuation magnitudes comparable to the active conformation. Therefore our study reports for the first time, intrinsic higher structural fluctuation for almost all inactive conformations compared to the active forms. Regions with higher fluctuations in the inactive states are often localized to the aC-helix, aG-helix and activation loop which are involved in the regulation and/or in structural transitions between active and inactive states. Further analysis of 476 kinase structures involved in interactions with another domain/protein showed that many of the regions with higher inactive-state fluctuation correspond to contact interfaces. We also performed extensive GNM analysis of (i) insulin receptor kinase bound to another protein and (ii) holo and apo forms of active and inactive conformations followed by multi-factor analysis of variance. We conclude that binding of small molecules or other domains/proteins reduce the extent of fluctuation irrespective of active or inactive forms. Finally, we show that the perceived fluctuations serve as a useful input to predict the functional state of a kinase.