CLAP: A web-server for automatic classification of proteins with special reference to multi-domain proteins

Background: The function of a protein can be deciphered with higher accuracy from its structure than from its amino acid sequence. Due to the huge gap in the available protein sequence and structural space, tools that can generate functionally homogeneous clusters using only the sequence information, hold great importance. For this, traditional alignment-based tools work well in most cases and clustering is performed on the basis of sequence similarity. But, in the case of multi-domain proteins, the alignment quality might be poor due to varied lengths of the proteins, domain shuffling or circular permutations. Multi-domain proteins are ubiquitous in nature, hence alignment-free tools, which overcome the shortcomings of alignment-based protein comparison methods, are required. Further, existing tools classify proteins using only domain-level information and hence miss out on the information encoded in the tethered regions or accessory domains. Our method, on the other hand, takes into account the full-length sequence of a protein, consolidating the complete sequence information to understand a given protein better. Results: Our web-server, CLAP (Classification of Proteins), is one such alignment-free software for automatic classification of protein sequences. It utilizes a pattern-matching algorithm that assigns local matching scores (LMS) to residues that are a part of the matched patterns between two sequences being compared. CLAP works on full-length sequences and does not require prior domain definitions. Pilot studies undertaken previously on protein kinases and immunoglobulins have shown that CLAP yields clusters, which have high functional and domain architectural similarity. Moreover, parsing at a statistically determined cut-off resulted in clusters that corroborated with the sub-family level classification of that particular domain family. Conclusions: CLAP is a useful protein-clustering tool, independent of domain assignment, domain order, sequence length and domain diversity. Our method can be used for any set of protein sequences, yielding functionally relevant clusters with high domain architectural homogeneity. The CLAP web server is freely available for academic use at http://nslab.mbu.iisc.ernet.in/clap/.

Direct regulation of topoisomerase activity by a nucleoid-associated protein

The topological homeostasis of bacterial chromosomes is maintained by the balance between compaction and the topological organization of genomes. Two classes of proteins play major roles in chromosome organization: the nucleoid-associated proteins (NAPs) and topoisomerases. The NAPs bind DNA to compact the chromosome, whereas topoisomerases catalytically remove or introduce supercoils into the genome. We demonstrate that HU, a major NAP of Mycobacterium tuberculosis specifically stimulates the DNA relaxation ability of mycobacterial topoisomerase I (TopoI) at lower concentrations but interferes at higher concentrations. A direct physical interaction between M. tuberculosis HU (MtHU) and TopoI is necessary for enhancing enzyme activity both in vitro and in vivo. The interaction is between the amino terminal domain of MtHU and the carboxyl terminal domain of TopoI. Binding of MtHU did not affect the two catalytic trans-esterification steps but enhanced the DNA strand passage, requisite for the completion of DNA relaxation, a new mechanism for the regulation of topoisomerase activity. An interaction-deficient mutant of MtHU was compromised in enhancing the strand passage activity. The species-specific physical and functional cooperation between MtHU and TopoI may be the key to achieve the DNA relaxation levels needed to maintain the optimal superhelical density of mycobacterial genomes.

Sensitivity of polarization fluctuations to the nature of protein-water interactions: Study of biological water in four different protein-water systems

Since the time of Kirkwood, observed deviations in magnitude of the dielectric constant of aqueous protein solution from that of neat water (similar to 80) and slower decay of polarization have been subjects of enormous interest, controversy, and debate. Most of the common proteins have large permanent dipole moments (often more than 100 D) that can influence structure and dynamics of even distant water molecules, thereby affecting collective polarization fluctuation of the solution, which in turn can significantly alter solution's dielectric constant. Therefore, distance dependence of polarization fluctuation can provide important insight into the nature of biological water. We explore these aspects by studying aqueous solutions of four different proteins of different characteristics and varying sizes, chicken villin headpiece subdomain (HP-36), immunoglobulin binding domain protein G (GB1), hen-egg white lysozyme (LYS), and Myoglobin (MYO). We simulate fairly large systems consisting of single protein molecule and 20000-30000 water molecules (varied according to the protein size), providing a concentration in the range of similar to 2-3 mM. We find that the calculated dielectric constant of the system shows a noticeable increment in all the cases compared to that of neat water. Total dipole moment auto time correlation function of water < dM(W) (0)delta M-W (t) > is found to be sensitive to the nature of the protein. Surprisingly, dipole moment of the protein and total dipole moment of the water molecules are found to be only weakly coupled. Shellwise decomposition of water molecules around protein reveals higher density of first layer compared to the succeeding ones. We also calculate heuristic effective dielectric constant of successive layers and find that the layer adjacent to protein has much lower value (similar to 50). However, progressive layers exhibit successive increment of dielectric constant, finally reaching a value close to that of bulk 4-5 layers away. We also calculate shellwise orientational correlation function and tetrahedral order parameter to understand the local dynamics and structural re-arrangement of water. Theoretical analysis providing simple method for calculation of shellwise local dielectric constant and implication of these findings are elaborately discussed in the present work. (C) 2014 AIP Publishing LLC.

Chaperone protein HYPK interacts with the first 17 amino acid region of Huntingtin and modulates mutant HTT-mediated aggregation and cytotoxicity

Huntington's disease is a polyglutamine expansion disorder, characterized by mutant HTT-mediated aggregate formation and cytotoxicity. Many reports suggests roles of N-terminal 17 amino acid domain of HTT (HTT-N17) towards subcellular localization, aggregate formation and subsequent pathogenicity induced by N-terminal HTT harboring polyQ stretch in pathogenic range. HYPK is a HTT-interacting chaperone which can reduce N-terminal mutant HTT-mediated aggregate formation and cytotoxicity in neuronal cell lines. However, how HYPK interacts with N-terminal fragment of HTT remained unknown. Here we report that specific interaction of HYPK with HTT-N17 is crucial for the chaperone activity of HYPK. Deletion of HTT-N17 leads to formation of tinier, SDS-soluble nuclear aggregates formed by N-terminal mutant HTT. The increased cytotoxicity imparted by these tiny aggregates might be contributed due to loss of interaction with HYPK. (C) 2014 Elsevier Inc. All rights reserved.

NrichD database: sequence databases enriched with computationally designed protein-like sequences aid in remote homology detection

NrichD ( ext-link-type=''uri'' xlink:href=''http://proline.biochem.iisc.ernet.in/NRICHD/'' xlink:type=''simple''>http://proline.biochem.iisc.ernet.in/NRICHD/)< /named-content> is a database of computationally designed protein-like sequences, augmented into natural sequence databases that can perform hops in protein sequence space to assist in the detection of remote relationships. Establishing protein relationships in the absence of structural evidence or natural `intermediately related sequences' is a challenging task. Recently, we have demonstrated that the computational design of artificial intermediary sequences/linkers is an effective approach to fill naturally occurring voids in protein sequence space. Through a large-scale assessment we have demonstrated that such sequences can be plugged into commonly employed search databases to improve the performance of routinely used sequence search methods in detecting remote relationships. Since it is anticipated that such data sets will be employed to establish protein relationships, two databases that have already captured these relationships at the structural and functional domain level, namely, the SCOP database and the Pfam database, have been `enriched' with these artificial intermediary sequences. NrichD database currently contains 3 611 010 artificial sequences that have been generated between 27 882 pairs of families from 374 SCOP folds. The data sets are freely available for download. Additional features include the design of artificial sequences between any two protein families of interest to the user.

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.

A Universal Stress Protein (USP) in Mycobacteria Binds cAMP

Mycobacteria are endowed with rich and diverse machinery for the synthesis, utilization, and degradation of cAMP. The actions of cyclic nucleotides are generally mediated by binding of cAMP to conserved and well characterized cyclic nucleotide binding domains or structurally distinct cGMP-specific and -regulated cyclic nucleotide phosphodiesterase, adenylyl cyclase, and E. coli transcription factor FhlA (GAF) domain-containing proteins. Proteins with cyclic nucleotide binding and GAF domains can be identified in the genome of mycobacterial species, and some of them have been characterized. Here, we show that a significant fraction of intracellular cAMP is bound to protein in mycobacterial species, and by using affinity chromatography techniques, we identify specific universal stress proteins (USP) as abundantly expressed cAMP-binding proteins in slow growing as well as fast growing mycobacteria. We have characterized the biochemical and thermodynamic parameters for binding of cAMP, and we show that these USPs bind cAMP with a higher affinity than ATP, an established ligand for other USPs. We determined the structure of the USP MSMEG_3811 bound to cAMP, and we confirmed through structure-guided mutagenesis, the residues important for cAMP binding. This family of USPs is conserved in all mycobacteria, and we suggest that they serve as ``sinks'' for cAMP, making this second messenger available for downstream effectors as and when ATP levels are altered in the cell.

The nuclease activities of both the Smr domain and an additional LDLK motif are required for an efficient anti-recombination function of Helicobacter pyloriMutS2

Helicobacter pylori, a human pathogen, is a naturally and constitutively competent bacteria, displaying a high rate of intergenomic recombination. While recombination events are essential for evolution and adaptation of H.pylori to dynamic gastric niches and new hosts, such events should be regulated tightly to maintain genomic integrity. Here, we analyze the role of the nuclease activity of MutS2, a protein that limits recombination during transformation in H.pylori. In previously studied MutS2 proteins, the C-terminal Smr domain was mapped as the region responsible for its nuclease activity. We report here that deletion of Smr domain does not completely abolish the nuclease activity of HpMutS2. Using bioinformatics analysis and mutagenesis, we identified an additional and novel nuclease motif (LDLK) at the N-terminus of HpMutS2 unique to Helicobacter and related epsilon-proteobacterial species. A single point mutation (D30A) in the LDLK motif and the deletion of Smr domain resulted in approximate to 5-10-fold loss of DNA cleavage ability of HpMutS2. Interestingly, the mutant forms of HpMutS2 wherein the LDLK motif was mutated or the Smr domain was deleted were unable to complement the hyper-recombination phenotype of a mutS2(-) strain, suggesting that both nuclease sites are indispensable for an efficient anti-recombinase activity of HpMutS2.

Do we see what we should see? Describing non-covalent interactions in protein structures including precision

The power of X-ray crystal structure analysis as a technique is to `see where the atoms are'. The results are extensively used by a wide variety of research communities. However, this `seeing where the atoms are' can give a false sense of security unless the precision of the placement of the atoms has been taken into account. Indeed, the presentation of bond distances and angles to a false precision (i.e. to too many decimal places) is commonplace. This article has three themes. Firstly, a basis for a proper representation of protein crystal structure results is detailed and demonstrated with respect to analyses of Protein Data Bank entries. The basis for establishing the precision of placement of each atom in a protein crystal structure is non-trivial. Secondly, a knowledge base harnessing such a descriptor of precision is presented. It is applied here to the case of salt bridges, i.e. ion pairs, in protein structures; this is the most fundamental place to start with such structure-precision representations since salt bridges are one of the tenets of protein structure stability. Ion pairs also play a central role in protein oligomerization, molecular recognition of ligands and substrates, allosteric regulation, domain motion and alpha-helix capping. A new knowledge base, SBPS (Salt Bridges in Protein Structures), takes these structural precisions into account and is the first of its kind. The third theme of the article is to indicate natural extensions of the need for such a description of precision, such as those involving metalloproteins and the determination of the protonation states of ionizable amino acids. Overall, it is also noted that this work and these examples are also relevant to protein three-dimensional structure molecular graphics software.

Evidence for N-7 guanine methyl transferase activity encoded within the modular domain of RNA-dependent RNA polymerase L of a Morbillivirus

Post-transcriptional modification of viral mRNA is essential for the translation of viral proteins by cellular translation machinery. Due to the cytoplasmic replication of Paramyxoviruses, the viral-encoded RNA-dependent RNA polymerase (RdRP) is thought to possess all activities required for mRNA capping and methylation. In the present work, using partially purified recombinant RNA polymerase complex of rinderpest virus expressed in insect cells, we demonstrate the in vitro methylation of capped mRNA. Further, we show that a recombinant C-terminal fragment (1717-2183 aa) of L protein is capable of methylating capped mRNA, suggesting that the various post-transcriptional activities of the L protein are located in independently folding domains.

Ac2PIM-responsive miR-150 and miR-143 Target Receptor-interacting Protein Kinase 2 and Transforming Growth Factor Beta-activated Kinase 1 to Suppress NOD2-induced Immunomodulators

Specific and coordinated regulation of innate immune receptor-driven signaling networks often determines the net outcome of the immune responses. Here, we investigated the cross-regulation of toll-like receptor (TLR)2 and nucleotide-binding oligomerization domain (NOD)2 pathways mediated by Ac2PIM, a tetra-acylated form of mycobacterial cell wall component and muramyl dipeptide (MDP), a peptidoglycan derivative respectively. While Ac2PIM treatment of macrophages compromised their ability to induce NOD2-dependent immunomodulators like cyclooxygenase (COX)-2, suppressor of cytokine signaling (SOCS)-3, and matrix metalloproteinase (MMP)-9, no change in the NOD2-responsive NO, TNF-alpha, VEGF-A, and IL-12 levels was observed. Further, genome-wide microRNA expression profiling identified Ac2PIM-responsive miR-150 and miR-143 to target NOD2 signaling adaptors, RIP2 and TAK1, respectively. Interestingly, Ac2PIM was found to activate the SRC-FAK-PYK2-CREB cascade via TLR2 to recruit CBP/P300 at the promoters of miR-150 and miR-143 and epigenetically induce their expression. Loss-of-function studies utilizing specific miRNA inhibitors establish that Ac2PIM, via the miRNAs, abrogate NOD2-induced PI3K-PKC delta-MAPK pathway to suppress beta-catenin-mediated expression of COX-2, SOCS-3, and MMP-9. Our investigation has thus underscored the negative regulatory role of Ac2PIM-TLR2 signaling on NOD2 pathway which could broaden our understanding on vaccine potential or adjuvant utilities of Ac2PIM and/or MDP.

Negative Cooperativity and High Affinity in Chitooligosaccharide Binding by a Mycobacterium smegmatis Protein Containing LysM and Lectin Domains

LysM domains have been recognized in bacteria and eukaryotes as carbohydrate-binding protein modules, but the mechanism of their binding to chitooligosaccharides has been underexplored. Binding of a Mycobacterium smegmatis protein containing a lectin (MSL) and one LysM domain to chitooligosaccharides has been studied using isothermal titration calorimetry and fluorescence titration that demonstrate the presence of two binding sites of nonidentical affinities per dimeric MSL-LysM molecule. The affinity of the molecule for chitooligosaccharides correlates with the length of the carbohydrate chain. Its binding to chitooligosaccharides is characterized by negative cooperativity in the interactions of the two domains. Apparently, the flexibility of the long linker that connects the LysM and MSL domains plays a facilitating role in this recognition. The LysM domain in the MSL-LysM molecule, like other bacterial domains but unlike plant LysM domains, recognizes equally well peptidoglycan fragments as well as chitin polymers. Interestingly, in the case presented here, two LysM domains are enough for binding to peptidoglycan in contrast to the three reportedly required by the LysM domains of Bacillus subtilis and Lactococcus lactis. Also, the affinity of the MSL-LysM molecule for chitooligosaccharides is higher than that of LysM-chitooligosaccharide interactions reported so far.

Structural studies on chimeric Sesbania mosaic virus coat protein: Revisiting SeMV assembly

The capsid protein (CP) of Sesbania mosaic virus (SeMV, a T=3 plant virus) consists of a disordered N-terminal R-domain and an ordered S-domain. Removal of the R-domain results in the formation of T=1 particles. In the current study, the R-domain was replaced with unrelated polypeptides of similar lengths: the B-domain of Staphylococcus aureus SpA, and SeMV encoded polypeptides P8 and P10. The chimeric proteins contained T=3 or larger virus-like particles (VLPs) and could not be crystallized. The presence of metal ions during purification resulted in a large number of heterogeneous nucleoprotein complexes. N Delta 65-B (R domain replaced with B domain) could also be purified in a dimeric form. Its crystal structure revealed T=1 particles devoid of metal ions and the B-domain was disordered. However, the B-domain was functional in N Delta 65-B VLPs, suggesting possible biotechnological applications. These studies illustrate the importance of N-terminal residues, metal ions and robustness of the assembly process. (C) 2015 Elsevier Inc. All rights reserved.

Folding of Protein L with Implications for Collapse in the Denatured State Ensemble

A fundamental question in protein folding is whether the coil to globule collapse transition occurs during the initial stages of folding (burst phase) or simultaneously with the protein folding transition. Single molecule fluorescence resonance energy transfer (FRET) and small-angle X-ray scattering (SAXS) experiments disagree on whether Protein L collapse transition occurs during the burst phase of folding. We study Protein L folding using a coarse-grained model and molecular dynamics simulations. The collapse transition in Protein L is found to be concomitant with the folding transition. In the burst phase of folding, we find that FRET experiments overestimate radius of gyration, R-g, of the protein due to the application of Gaussian polymer chain end-to-end distribution to extract R-g from the FRET efficiency. FRET experiments estimate approximate to 6 angstrom decrease in R-g when the actual decrease is approximate to 3 angstrom on guanidinium chloride denaturant dilution from 7.5 to 1 M, thereby suggesting pronounced compaction in the protein dimensions in the burst phase. The approximate to 3 angstrom decrease is close to the statistical uncertainties of the R-g data measured from SAXS experiments, which suggest no compaction, leading to a disagreement with the FRET experiments. The transition-state ensemble (TSE) structures in Protein L folding are globular and extensive in agreement with the Psi-analysis experiments. The results support the hypothesis that the TSE of single domain proteins depends on protein topology and is not stabilized by local interactions alone.