Under-representation of intrinsic terminators across bacterial genomic islands: Rho as a principal regulator of xenogenic DNA expression

Two transcription termination mechanisms - intrinsic and Rho-dependent - have evolved in bacteria. The Rho factor occurs in most bacterial lineages, and has been hypothesized to play a global regulatory role. Genome-wide studies using microarray, 2D-gel electrophoresis and ChIP-chip provided evidence that Rho serves to silence transcription from horizontally acquired genes and prophages in Escherichia coli K-12, implicating the factor to be a part of the ``cellular immune mechanism'' protecting against deleterious phages and aberrant gene expression from acquired xenogenic DNA. We have investigated this model by adopting an alternate in silico approach and have extended the study to other species. Our analysis shows that several genomic islands across diverse phyla have under-representation of intrinsic terminators, similar to that experimentally observed in E. coli K-12. This implies that Rho-dependent termination is the predominant process operational in these islands and that silencing of foreign DNA is a conserved function of Rho. From the present analysis, it is evident that horizontally acquired islands have lost intrinsic terminators to facilitate Rho-dependent termination. These results underscore the importance of Rho as a conserved, genome-wide sentinel that regulates potentially toxic xenogenic DNA. (C) 2012 Elsevier B.V. All rights reserved.

A Bifunctional Enzyme That Has Both Monoacylglycerol Acyltransferase and Acyl Hydrolase Activities

Monoacylglycerol acyltransferase (MGAT) catalyzes the synthesis of diacylglycerol, the precursor of triacylglycerol biosynthesis and an important signaling molecule. Here, we describe the isolation and characterization of the peanut (Arachis hypogaea) MGAT gene. The soluble enzyme utilizes invariant histidine-62 and aspartate-67 residues of the acyltransferase motif for its MGAT activity. A sequence analysis revealed the presence of a hydrolase (GXSXG) motif, and enzyme assays revealed the presence of monoacylglycerol (MAG) and lysophosphatidylcholine (LPC) hydrolytic activities, indicating the bifunctional nature of the enzyme. The overexpression of the MGAT gene in yeast (Saccharomyces cerevisiae) caused an increase in triacylglycerol accumulation. Similar to the peanut MGAT, the Arabidopsis (Arabidopsis thaliana) homolog (At1g52760) also exhibited both acyltransferase and hydrolase activities. Interestingly, the yeast homolog lacks the conserved HX4D motif, and it is deficient in the acyltransferase function but exhibits MAG and LPC hydrolase activities. This study demonstrates the presence of a soluble MGAT/hydrolase in plants. The predicted three-dimensional homology modeling and substrate docking suggested the presence of two separate substrate (MAG and LPC)-binding sites in a single polypeptide. Our study describes a soluble bifunctional enzyme that has both MGAT and hydrolase functions.

Structural and mechanistic investigations on Salmonella typhimurium acetate kinase (AckA): identification of a putative ligand binding pocket at the dimeric interface

Background: Bacteria such as Escherichia coli and Salmonella typhimurium can utilize acetate as the sole source of carbon and energy. Acetate kinase (AckA) and phosphotransacetylase (Pta), key enzymes of acetate utilization pathway, regulate flux of metabolites in glycolysis, gluconeogenesis, TCA cycle, glyoxylate bypass and fatty acid metabolism. Results: Here we report kinetic characterization of S. typhimurium AckA (StAckA) and structures of its unliganded (Form-I, 2.70 angstrom resolution) and citrate-bound (Form-II, 1.90 angstrom resolution) forms. The enzyme showed broad substrate specificity with k(cat)/K-m in the order of acetate > propionate > formate. Further, the K-m for acetyl-phosphate was significantly lower than for acetate and the enzyme could catalyze the reverse reaction (i.e. ATP synthesis) more efficiently. ATP and Mg2+ could be substituted by other nucleoside 5'-triphosphates (GTP, UTP and CTP) and divalent cations (Mn2+ and Co2+), respectively. Form-I StAckA represents the first structural report of an unliganded AckA. StAckA protomer consists of two domains with characteristic beta beta beta alpha beta alpha beta alpha topology of ASKHA superfamily of proteins. These domains adopt an intermediate conformation compared to that of open and closed forms of ligand-bound Methanosarcina thermophila AckA (MtAckA). Spectroscopic and structural analyses of StAckA further suggested occurrence of inter-domain motion upon ligand-binding. Unexpectedly, Form-II StAckA structure showed a drastic change in the conformation of residues 230-300 compared to that of Form-I. Further investigation revealed electron density corresponding to a citrate molecule in a pocket located at the dimeric interface of Form-II StAckA. Interestingly, a similar dimeric interface pocket lined with largely conserved residues could be identified in Form-I StAckA as well as in other enzymes homologous to AckA suggesting that ligand binding at this pocket may influence the function of these enzymes. Conclusions: The biochemical and structural characterization of StAckA reported here provides insights into the biochemical specificity, overall fold, thermal stability, molecular basis of ligand binding and inter-domain motion in AckA family of enzymes. Dramatic conformational differences observed between unliganded and citrate-bound forms of StAckA led to identification of a putative ligand-binding pocket at the dimeric interface of StAckA with implications for enzymatic function.

HSPIR: a manually annotated heat shock protein information resource

Heat shock protein information resource (HSPIR) is a concerted database of six major heat shock proteins (HSPs), namely, Hsp70, Hsp40, Hsp60, Hsp90, Hsp100 and small HSP. The HSPs are essential for the survival of all living organisms, as they protect the conformations of proteins on exposure to various stress conditions. They are a highly conserved group of proteins involved in diverse physiological functions, including de novo folding, disaggregation and protein trafficking. Moreover, their critical role in the control of disease progression made them a prime target of research. Presently, limited information is available on HSPs in reference to their identification and structural classification across genera. To that extent, HSPIR provides manually curated information on sequence, structure, classification, ontology, domain organization, localization and possible biological functions extracted from UniProt, GenBank, Protein Data Bank and the literature. The database offers interactive search with incorporated tools, which enhances the analysis. HSPIR is a reliable resource for researchers exploring structure, function and evolution of HSPs.

Interaction Signatures Stabilizing the NAD(P)-Binding Rossmann Fold: A Structure Network Approach

The fidelity of the folding pathways being encoded in the amino acid sequence is met with challenge in instances where proteins with no sequence homology, performing different functions and no apparent evolutionary linkage, adopt a similar fold. The problem stated otherwise is that a limited fold space is available to a repertoire of diverse sequences. The key question is what factors lead to the formation of a fold from diverse sequences. Here, with the NAD(P)-binding Rossmann fold domains as a case study and using the concepts of network theory, we have unveiled the consensus structural features that drive the formation of this fold. We have proposed a graph theoretic formalism to capture the structural details in terms of the conserved atomic interactions in global milieu, and hence extract the essential topological features from diverse sequences. A unified mathematical representation of the different structures together with a judicious concoction of several network parameters enabled us to probe into the structural features driving the adoption of the NAD(P)-binding Rossmann fold. The atomic interactions at key positions seem to be better conserved in proteins, as compared to the residues participating in these interactions. We propose a ``spatial motif'' and several ``fold specific hot spots'' that form the signature structural blueprints of the NAD(P)-binding Rossmann fold domain. Excellent agreement of our data with previous experimental and theoretical studies validates the robustness and validity of the approach. Additionally, comparison of our results with statistical coupling analysis (SCA) provides further support. The methodology proposed here is general and can be applied to similar problems of interest.

Structural and molecular basis of interaction of HCV non-structural protein 5A with human casein kinase 1 alpha and PKR

Background: Interaction of non-structural protein 5A (NS5A) of Hepatitis C virus (HCV) with human kinases namely, casein kinase 1 alpha (ck1 alpha) and protein kinase R (PKR) have different functional implications such as regulation of viral replication and evasion of interferon induced immune response respectively. Understanding the structural and molecular basis of interactions of the viral protein with two different human kinases can be useful in developing strategies for treatment against HCV. Results: Serine 232 of NS5A is known to be phosphorylated by human ck1 alpha. A structural model of NS5A peptide containing phosphoacceptor residue Serine 232 bound to ck1 alpha has been generated using the known 3-D structures of kinase-peptide complexes. The substrate interacting residues in ck1 alpha has been identified from the model and these are found to be conserved well in the ck1 family. ck1 alpha - substrate peptide complex has also been used to understand the structural basis of association between ck1 alpha and its other viral stress induced substrate, tumour suppressor p53 transactivation domain which has a crystal structure available. Interaction of NS5A with another human kinase PKR is primarily genotype specific. NS5A from genotype 1b has been shown to interact and inhibit PKR whereas NS5A from genotype 2a/3a are unable to bind and inhibit PKR efficiently. This is one of the main reasons for the varied response to interferon therapy in HCV patients across different genotypes. Using PKR crystal structure, sequence alignment and evolutionary trace analysis some of the critical residues responsible for the interaction of NS5A 1b with PKR have been identified. Conclusions: The substrate interacting residues in ck1 alpha have been identified using the structural model of kinase substrate peptide. The PKR interacting NS5A 1b residues have also been predicted using PKR crystal structure, NS5A sequence analysis along with known experimental results. Functional significance and nature of interaction of interferon sensitivity determining region and variable region 3 of NS5A in different genotypes with PKR which was experimentally shown are also supported by the findings of evolutionary trace analysis. Designing inhibitors to prevent this interaction could enable the HCV genotype 1 infected patients respond well to interferon therapy.

Plasmodium falciparum Prp16 homologue and its role in splicing

Large numbers of Plasmodium genes have been predicted to have introns. However, little information exists on the splicing mechanisms in this organism. Here, we describe the DExD/DExH-box containing Pre-mRNA processing proteins (Prps), PfPrp2p, PfPrp5p, PfPrp16p, PfPrp22p, PfPrp28p, PfPrp43p and PfBrr2p, present in the Plasmodium falciparum genome and characterized the role of one of these factors, PfPrp16p. It is a member of DEAH-box protein family with nine collinear sequence motifs, a characteristic of helicase proteins. Experiments with the recombinantly expressed and purified PfPrp16 helicase domain revealed binding to RNA, hydrolysis of ATP as well as catalytic helicase activities. Expression of helicase domain with the C-terminal helicase-associated domain (HA2) reduced these activities considerably, indicating that the helicase-associated domain may regulate the PfPrp16 function. Localization studies with the PfPrp16 GFP transgenic lines suggested a role of its N-terminal domain (1-80 amino acids) in nuclear targeting. Immunodepletion of PfPrp16p, from nuclear extracts of parasite cultures, blocked the second catalytic step of an in vitro constituted splicing reaction suggesting a role for PfPrp16p in splicing catalysis. Further we show by complementation assay in yeast that a chimeric yeast-Plasmodium Prp16 protein, not the full length PfPrp16, can rescue the yeast prp16 temperature-sensitive mutant. These results suggest that although the role of Prp16p in catalytic step II is highly conserved among Plasmodium, human and yeast, subtle differences exist with regards to its associated factors or its assembly with spliceosomes. (C) 2012 Elsevier B.V. All rights reserved.

Functional characterization, homology modeling and docking studies of beta-glucosidase responsible for bioactivation of cyanogenic hydroxynitrile glucosides from Leucaena leucocephala (subabul)

Glycosyl hydrolase family 1 beta-glucosidases are important enzymes that serve many diverse functions in plants including defense, whereby hydrolyzing the defensive compounds such as hydroxynitrile glucosides. A hydroxynitrile glucoside cleaving beta-glucosidase gene (Llbglu1) was isolated from Leucaena leucocephala, cloned into pET-28a (+) and expressed in E. coli BL21 (DE3) cells. The recombinant enzyme was purified by Ni-NTA affinity chromatography. The optimal temperature and pH for this beta-glucosidase were found to be 45 A degrees C and 4.8, respectively. The purified Llbglu1 enzyme hydrolyzed the synthetic glycosides, pNPGlucoside (pNPGlc) and pNPGalactoside (pNPGal). Also, the enzyme hydrolyzed amygdalin, a hydroxynitrile glycoside and a few of the tested flavonoid and isoflavonoid glucosides. The kinetic parameters K (m) and V (max) were found to be 38.59 mu M and 0.8237 mu M/mg/min for pNPGlc, whereas for pNPGal the values were observed as 1845 mu M and 0.1037 mu M/mg/min. In the present study, a three dimensional (3D) model of the Llbglu1 was built by MODELLER software to find out the substrate binding sites and the quality of the model was examined using the program PROCHEK. Docking studies indicated that conserved active site residues are Glu 199, Glu 413, His 153, Asn 198, Val 270, Asn 340, and Trp 462. Docking of rhodiocyanoside A with the modeled Llbglu1 resulted in a binding with free energy change (Delta G) of -5.52 kcal/mol on which basis rhodiocyanoside A could be considered as a potential substrate.

Structure of an amidohydrolase, SACOL0085, from methicillin-resistant Staphylococcus aureus COL

Staphylococcus aureus is an opportunistic pathogen that rapidly acquires resistance to frontline antibiotics. The characterization of novel protein targets from this bacterium is thus an important step towards future therapeutic strategies. Here, the crystal structure of an amidohydrolase, SACOL0085, from S. aureus COL is described. SACOL0085 is a member of the M20D family of peptidases. Unlike other M20D peptidases, which are either monomers or dimers, SACOL0085 adopts a butterfly-shaped homotetrameric arrangement with extensive intersubunit interactions. Each subunit of SACOL0085 contains two Mn2+ ions at the active site. A conserved cysteine residue at the active site distinguishes M20D peptidases from other M20 family members. This cysteine, Cys103, serves as bidentate ligand coordinating both Mn2+ ions in SACOL0085.

Purification and characterization of prophenoloxidase from cotton bollworm, Helicoverpa armigera

Phenoloxidases are oxidative enzymes, which play an important role in both cell mediated and humoral immunity. Purification and biochemical characterization of prophenoloxidase from cotton bollworm, Helicoverpa armigera (Hubner) were carried out to study its biochemical properties. Prophenoloxidase consists of a single polypeptide chain with a relative molecular weight of 85 kDa as determined by SDSPAGE, MALDITOF MS and LCESI MS. After the final step, the enzyme showed 71.7 fold of purification with a recovery of 49.2%. Purified prophenoloxidase showed high specific activity and homology with phenoloxidase subunit-1 of Bombyx mori and the conserved regions of copper binding (B) site of phenoloxidase. Purified prophenoloxidase has pH optima of 6.8 and has high catalytic efficiency towards the dopamine as a substrate in comparison to catechol and L-Dopa. The PO activity was strongly inhibited by phenylthiourea, thiourea, dithiothreitol and kojic acid.

SuSeFLAV 1.2: Program for supersymmetric mass spectra with seesaw mechanism and rare lepton flavor violating decays

Accurate supersymmetric spectra are required to confront data from direct and indirect searches of supersymmetry. SuSeFLAV is a numerical tool capable of computing supersymmetric spectra precisely for various supersymmetric breaking scenarios applicable even in the presence of flavor violation. The program solves MSSM RGEs with complete 3 x 3 flavor mixing at 2-loop level and one loop finite threshold corrections to all MSSM parameters by incorporating radiative electroweak symmetry breaking conditions. The program also incorporates the Type-I seesaw mechanism with three massive right handed neutrinos at user defined mass scales and mixing. It also computes branching ratios of flavor violating processes such as l(j) -> l(i)gamma, l(j) -> 3 l(i), b -> s gamma and supersymmetric contributions to flavor conserving quantities such as (g(mu) - 2). A large choice of executables suitable for various operations of the program are provided. Program summary Program title: SuSeFLAV Catalogue identifier: AEOD_v1_0 Program summary URL: http://cpc.cs.qub.ac.uk/summaries/AEOD_v1_0.html Program obtainable from: CPC Program Library, Queen's University, Belfast, N. Ireland Licensing provisions: GNU General Public License No. of lines in distributed program, including test data, etc.: 76552 No. of bytes in distributed program, including test data, etc.: 582787 Distribution format: tar.gz Programming language: Fortran 95. Computer: Personal Computer, Work-Station. Operating system: Linux, Unix. Classification: 11.6. Nature of problem: Determination of masses and mixing of supersymmetric particles within the context of MSSM with conserved R-parity with and without the presence of Type-I seesaw. Inter-generational mixing is considered while calculating the mass spectrum. Supersymmetry breaking parameters are taken as inputs at a high scale specified by the mechanism of supersymmetry breaking. RG equations including full inter-generational mixing are then used to evolve these parameters up to the electroweak breaking scale. The low energy supersymmetric spectrum is calculated at the scale where successful radiative electroweak symmetry breaking occurs. At weak scale standard model fermion masses, gauge couplings are determined including the supersymmetric radiative corrections. Once the spectrum is computed, the program proceeds to various lepton flavor violating observables (e.g., BR(mu -> e gamma), BR(tau -> mu gamma) etc.) at the weak scale. Solution method: Two loop RGEs with full 3 x 3 flavor mixing for all supersymmetry breaking parameters are used to compute the low energy supersymmetric mass spectrum. An adaptive step size Runge-Kutta method is used to solve the RGEs numerically between the high scale and the electroweak breaking scale. Iterative procedure is employed to get the consistent radiative electroweak symmetry breaking condition. The masses of the supersymmetric particles are computed at 1-loop order. The third generation SM particles and the gauge couplings are evaluated at the 1-loop order including supersymmetric corrections. A further iteration of the full program is employed such that the SM masses and couplings are consistent with the supersymmetric particle spectrum. Additional comments: Several executables are presented for the user. Running time: 0.2 s on a Intel(R) Core(TM) i5 CPU 650 with 3.20 GHz. (c) 2012 Elsevier B.V. All rights reserved.

Regulation of lipid biosynthesis, sliding motility, and biofilm formation by a membrane-anchored nucleoid-associated protein of mycobacterium tuberculosis

Bacteria use a number of small basic proteins for organization and compaction of their genomes. By their interaction with DNA, these nucleoid-associated proteins (NAPs) also influence gene expression. Rv3852, a NAP of Mycobacterium tuberculosis, is conserved among the pathogenic and slow-growing species of mycobacteria. Here, we show that the protein predominantly localizes in the cell membrane and that the carboxy-terminal region with the propensity to form a transmembrane helix is necessary for its membrane localization. The protein is involved in genome organization, and its ectopic expression in Mycobacterium smegmatis resulted in altered nucleoid morphology, defects in biofilm formation, sliding motility, and change in apolar lipid profile. We demonstrate its crucial role in regulating the expression of KasA, KasB, and GroEL1 proteins, which are in turn involved in controlling the surface phenotypes in mycobacteria.

Mycobacterium tuberculosis RsdA provides a conformational rationale for selective regulation of sigma-factor activity by proteolysis

The relative levels of different sigma factors dictate the expression profile of a bacterium. Extracytoplasmic function sigma factors synchronize the transcriptional profile with environmental conditions. The cellular concentration of free extracytoplasmic function sigma factors is regulated by the localization of this protein in a sigma/anti-sigma complex. Anti-sigma factors are multi-domain proteins with a receptor to sense environmental stimuli and a conserved anti-sigma domain (ASD) that binds a sigma factor. Here we describe the structure of Mycobacterium tuberculosis anti-sigma(D) (RsdA) in complex with the -35 promoter binding domain of sigma(D) (sigma(D)(4)). We note distinct conformational features that enable the release of sigma(D) by the selective proteolysis of the ASD in RsdA. The structural and biochemical features of the sigma(D)/RsdA complex provide a basis to reconcile diverse regulatory mechanisms that govern sigma/anti-sigma interactions despite high overall structural similarity. Multiple regulatory mechanisms embedded in an ASD scaffold thus provide an elegant route to rapidly re-engineer the expression profile of a bacterium in response to an environmental stimulus.

Assay development for identifying inhibitors of the mycobacterial FadD32 activity

FadD32, a fatty acyl-AMP ligase (FAAL32) involved in the biosynthesis of mycolic acids, major and specific lipid components of the mycobacterial cell envelope, is essential for the survival of Mycobacterium tuberculosis, the causative agent of tuberculosis. The protein catalyzes the conversion of fatty acid to acyl-adenylate (acyl-AMP) in the presence of adenosine triphosphate and is conserved in all the mycobacterial species sequenced so far, thus representing a promising target for the development of novel antituberculous drugs. Here, we describe the optimization of the protein purification procedure and the development of a high-throughput screening assay for FadD32 activity. This spectrophotometric assay measuring the release of inorganic phosphate was optimized using the Mycobacterium smegmatis FadD32 as a surrogate enzyme. We describe the use of Tm (melting temperature) shift assay, which measures the modulation of FadD32 thermal stability, as a tool for the identification of potential ligands and for validation of compounds as inhibitors. Screening of a selected library of compounds led to the identification of five novel classes of inhibitors.

Multiple molecular alterations in phosphodegrons 1-3 within AAV2 capsid demonstrates higher hepatic gene transfer efficiency

The success of AAV2 mediated hepatic gene transfer in human trials for diseases such as hemophilia has been hampered by a combination of low transduction efficiency and a robust immune response directed against these vectors. We have previously shown that AAV2 is targeted for destruction in the cytoplasm by the host-cellular kinase/ubiquitination/proteasomal degradation machinery and modification of the serine(S)/threonine(T) kinase and lysine(K) targets on AAV capsid is beneficial. Thus targeted single mutations of S/T>A(S489A, S498A, T251A) and K>R (K532R) improved the efficiency of gene transfer in vivo as compared to wild type (WT)-AAV2 vectors (∼6-14 fold). In the present study, we evaluated if combined alteration of the phosphodegrons (PD), which are the phosphorylation sites recognized as degradation signals by ubiquitin ligases, improves further the gene transfer efficiency. Thus, we generated four multiple mutant vectors (PD: 1+3, S489A+K532R, PD: 1+3, S489A+K532R together with T251 residue which did not lie in any of the phosphodegrons but had shown increased transduction efficiency compared to the WT-AAV2 vector (∼6 fold) and was also conserved in 9 out of 10 AAV serotypes (AAV 1 to 10), PD: 1+3, S489A+K532R+S498A and a fourth combination of PD: 3, K532R+T251. We then evaluated them in vitro and in vivo and compared their gene transfer efficiency with either the WT-AAV2 or the best single mutant S489A-AAV2 vector. The novel multiple mutations on the AAV2 capsid did not affect the overall vector packaging efficiency. All the multiple AAV2 mutants showed superior transduction efficiency in HeLa cells in vitro when compared to either the WT (62-72% Vs 21%) or the single mutant S489A (62-72% Vs 50%) AAV2 vectors as demonstrated by FACS analysis (Fig. 1A). On hepatic gene transfer with 5x10^10 vgs per animal in C57BL/6 mice, all the multiple mutants showed increased transgene expression compared to either the WT-AAV2 (∼15-23 fold) or the S489A single mutant vector (∼2-3 fold) (Fig.1B and C). These novel multiple mutant AAV2 vectors also showed higher vector copy number in murine hepatocytes 4 weeks post transduction, as compared to the WT-AAV2 (∼5-6 Vs 1.4 vector copies/diploid genome) and further higher when compared to the single mutant S489A(∼5-6 fold Vs 3.8 fold) (Fig.1D). Further ongoing studies will demonstrate the therapeutic benefit of one or more of the multiple mutants vectors in preclinical models of hemophilia.