Classification of Protein Kinases on the Basis of Both Kinase and Non-Kinase Regions

Background: Protein phosphorylation is a generic way to regulate signal transduction pathways in all kingdoms of life. In many organisms, it is achieved by the large family of Ser/Thr/Tyr protein kinases which are traditionally classified into groups and subfamilies on the basis of the amino acid sequence of their catalytic domains. Many protein kinases are multidomain in nature but the diversity of the accessory domains and their organization are usually not taken into account while classifying kinases into groups or subfamilies. Methodology: Here, we present an approach which considers amino acid sequences of complete gene products, in order to suggest refinements in sets of pre-classified sequences. The strategy is based on alignment-free similarity scores and iterative Area Under the Curve (AUC) computation. Similarity scores are computed by detecting common patterns between two sequences and scoring them using a substitution matrix, with a consistent normalization scheme. This allows us to handle full-length sequences, and implicitly takes into account domain diversity and domain shuffling. We quantitatively validate our approach on a subset of 212 human protein kinases. We then employ it on the complete repertoire of human protein kinases and suggest few qualitative refinements in the subfamily assignment stored in the KinG database, which is based on catalytic domains only. Based on our new measure, we delineate 37 cases of potential hybrid kinases: sequences for which classical classification based entirely on catalytic domains is inconsistent with the full-length similarity scores computed here, which implicitly consider multi-domain nature and regions outside the catalytic kinase domain. We also provide some examples of hybrid kinases of the protozoan parasite Entamoeba histolytica. Conclusions: The implicit consideration of multi-domain architectures is a valuable inclusion to complement other classification schemes. The proposed algorithm may also be employed to classify other families of enzymes with multidomain architecture.

An algorithm to find distant repeats in a pair of protein sequences

Distant repeats between a pair of protein sequences can be exploited to study the various aspects of proteins such as structure-function relationship, disorders due to protein malfunction, evolutionary analysis, etc. An in-depth analysis of the distant repeats would facilitate to establish a stable evolutionary relation of the repeats with respect to their three-dimensional structure. To this effect, an algorithm has been devised to identify the distant repeats in a pair of protein sequences by essentially using the scores of PAM (Percent Accepted Mutation) matrices. The proposed algorithm will be of much use to researchers involved in the comparative study of various organisms based on the amino-acid repeats in protein sequences. (C) 2010 Elsevier B.V. All rights reserved.

The Multifunctional PE_PGRS11 Protein from Mycobacterium tuberculosis Plays a Role in Regulating Resistance to Oxidative Stress

Mycobacterium tuberculosis utilizes unique strategies to survive amid the hostile environment of infected host cells. Infection-specific expression of a unique mycobacterial cell surface antigen that could modulate key signaling cascades can act as a key survival strategy in curtailing host effector responses like oxidative stress. We demonstrate here that hypothetical PE_PGRS11 ORF encodes a functional phosphoglycerate mutase. The transcriptional analysis revealed that PE_PGRS11 is a hypoxia-responsive gene, and enforced expression of PE_PGRS11 by recombinant adenovirus or Mycobacterium smegmatis imparted resistance to alveolar epithelial cells against oxidative stress. PE_PGRS11-induced resistance to oxidative stress necessitated the modulation of genetic signatures like induced expression of Bcl2 or COX-2. This modulation of specific antiapoptotic molecular signatures involved recognition of PE_PGRS11 by TLR2 and subsequent activation of the PI3K-ERK1/ 2-NF-kappa B signaling axis. Furthermore, PE_PGRS11 markedly diminished H2O2-induced p38 MAPK activation. Interestingly, PE_PGRS11 protein was exposed at the mycobacterial cell surface and was involved in survival of mycobacteria under oxidative stress. Furthermore, PE_PGRS11 displayed differential B cell responses during tuberculosis infection. Taken together, our investigation identified PE_PGRS11 as an in vivo expressed immunodominant antigen that plays a crucial role in modulating cellular life span restrictions imposed during oxidative stress by triggering TLR2-dependent expression of COX-2 and Bcl2. These observations clearly provide a mechanistic basis for the rescue of pathogenic Mycobacterium-infected lung epithelial cells from oxidative stress.

Prediction of protein-protein interactions in dengue virus coat proteins guided by low resolution cryoEM structures

Background: Dengue virus along with the other members of the flaviviridae family has reemerged as deadly human pathogens. Understanding the mechanistic details of these infections can be highly rewarding in developing effective antivirals. During maturation of the virus inside the host cell, the coat proteins E and M undergo conformational changes, altering the morphology of the viral coat. However, due to low resolution nature of the available 3-D structures of viral assemblies, the atomic details of these changes are still elusive. Results: In the present analysis, starting from C alpha positions of low resolution cryo electron microscopic structures the residue level details of protein-protein interaction interfaces of dengue virus coat proteins have been predicted. By comparing the preexisting structures of virus in different phases of life cycle, the changes taking place in these predicted protein-protein interaction interfaces were followed as a function of maturation process of the virus. Besides changing the current notion about the presence of only homodimers in the mature viral coat, the present analysis indicated presence of a proline-rich motif at the protein-protein interaction interface of the coat protein. Investigating the conservation status of these seemingly functionally crucial residues across other members of flaviviridae family enabled dissecting common mechanisms used for infections by these viruses. Conclusions: Thus, using computational approach the present analysis has provided better insights into the preexisting low resolution structures of virus assemblies, the findings of which can be made use of in designing effective antivirals against these deadly human pathogens.

Anti-apoptotic Protein BCL2 Down-regulates DNA End Joining in Cancer Cells

Cancer cells are often associated with secondary chromosomal rearrangements, such as deletions, inversions, and translocations, which could be the consequence of unrepaired/misrepaired DNA double strand breaks (DSBs). Nonhomologous DNA end joining is one of the most common pathways to repair DSBs in higher eukaryotes. By using oligomeric DNA substrates mimicking various endogenous DSBs in a cell-free system, we studied end joining (EJ) in different cancer cell lines. We found that the efficiency of EJ varies among cancer cells; however, there was no remarkable difference in the mechanism and expression of EJ proteins. Interestingly, cancer cells with lower levels of EJ possessed elevated expression of BCL2 and vice versa. Removal of BCL2 by immunoprecipitation or protein fractionation led to elevated EJ. More importantly, we show that overexpression of BCL2 or the addition of purified BCL2 led to the down-regulation of EJ. Further, we found that BCL2 interacts with KU proteins both in vitro and in vivo. Hence, our results suggest that EJ in cancer cells could be negatively regulated by the anti-apoptotic protein, BCL2, and this may contribute toward increased chromosomal abnormalities in cancer.

Inhibition of peroxidase-catalyzed protein tyrosine nitration by antithyroid drugs and their analogues

In this paper, we describe the effect of some commonly used thiourea-based antithyroid drugs and their analogues on the peroxidase-catalyzed nitration reactions. The nitration of bovine serum albumin (BSA) and cytochrome c was studied using the antibody against 3-nitro-L-tyrosine. This study reveals that the thione-based antithyroid drugs effectively inhibit lactoperoxidase (LPO)-catalyzed nitration of BSA. These compounds show very weak inhibition towards the nitration of cytochrome c. Some of these compounds also inhibit myeloperoxidase (MPO)-catalyzed nitration of L-tyrosine. A structure-activity correlation study on the peroxidase-catalyzed nitration of L-tyrosine reveals that the presence of thione/selone moiety is important for the inhibition. Although the presence of a free N-H group adjacent to C=S moiety is necessary for most of the thiones to inhibit the LPO-catalyzed nitration, the corresponding selones do not require the presence of any free N-H group for their activity. Furthermore, experiments with different concentrations of H2O2 suggest that the antithyroid drugs and related thiones inhibit the nitration reaction mainly by coordinating to the Fe(III)-center of the enzyme active site as previously proposed for the inhibition of peroxidase-catalyzed iodination. On the other hand, the selenium compounds inhibit the nitration by scavenging H2O2 without interacting with the enzyme active site. This assumption is based on the observations that catalase effectively inhibits tyrosine nitration by scavenging H2O2, which is one of the substrates for the nitration. In contrast, superoxide dismutase (SOD) does not alter the nitration reactions, indicating the absence of superoxide radical anion (O-2 center dot(-)) during the peroxidase-catalyzed nitration reactions. (C) 2010 Elsevier B.V. All rights reserved.

Empirical torsional potential functions from protein structure data. Phi- and psi-potentials for non-glycyl amino acid residues.

The torsional potential functions Vt(phi) and Vt(psi) around single bonds N--C alpha and C alpha--C, which can be used in conformational studies of oligopeptides, polypeptides and proteins, have been derived, using crystal structure data of 22 globular proteins, fitting the observed distribution in the (phi, psi)-plane with the value of Vtot(phi, psi), using the Boltzmann distribution. The averaged torsional potential functions, obtained from various amino acid residues in L-configuration, are Vt(phi) = 1.0 cos (phi + 60 degrees); Vt(psi) = 0.5 cos (psi + 60 degrees) - 1.0 cos (2 psi + 30 degrees) - 0.5 cos (3 psi + 30 degrees). The dipeptide energy maps Vtot(phi, psi) obtained using these functions, instead of the normally accepted torsional functions, were found to explain various observations, such as the absence of the left-handed alpha helix and the C7 conformation, and the relatively high density of points near the line psi = 0 degrees. These functions derived from observational data on protein structures, will, it is hoped, explain various previously unexplained facts in polypeptide conformation.

A monoclonal antibody recognizing the C-terminal region of chicken egg white riboflavin carrier protein terminates early pregnancy in mice

In order to identify the functionally relevant epitopes on chicken riboflavin carrier protein, we have raised monoclonal antibodies to the vitamin carrier. One of these, 6B2C12, was found to interact specifically with a synthetic oligopeptide corresponding to the C-terminal 17 amino acid residues of the chicken egg white riboflavin carrier protein, which is missing in part in the egg yolk riboflavin carrier protein. This epitope is conserved through evolution in mammals including humans. Administration of the ascites fluid of 6B2C12 to pregnant mice intraperitoneally, resulted in the termination of pregnancy indicating that this epitope is involved in or closely associated with the transplacental transport of the vitamin from the maternal circulation to the growing fetus.

Over-expression and purification strategies for recombinant multi-protein oligomers: A case study of Mycobacterium tuberculosis sigma/anti-sigma factor protein complexes

The function of a protein in a cell often involves coordinated interactions with one or several regulatory partners. It is thus imperative to characterize a protein both in isolation as well as in the context of its complex with an interacting partner. High resolution structural information determined by X-ray crystallography and Nuclear Magnetic Resonance offer the best route to characterize protein complexes. These techniques, however, require highly purified and homogenous protein samples at high concentration. This requirement often presents a major hurdle for structural studies. Here we present a strategy based on co-expression and co-purification to obtain recombinant multi-protein complexes in the quantity and concentration range that can enable hitherto intractable structural projects. The feasibility of this strategy was examined using the sigma factor/anti-sigma factor protein complexes from Mycobacterium tuberculosis. The approach was successful across a wide range of sigma factors and their cognate interacting partners. It thus appears likely that the analysis of these complexes based on variations in expression constructs and procedures for the purification and characterization of these recombinant protein samples would be widely applicable for other multi-protein systems. (C) 2010 Elsevier Inc. All rights reserved.

High-resolution methyl edited GFT NMR experiments for protein resonance assignments and structure determination

Three-dimensional (3D) structure determination of proteins is benefitted by long-range distance constraints comprising the methyl groups, which constitute the hydrophobic core of proteins. However, in methyl groups (of Ala, Ile, Leu, Met, Thr and Val) there is a significant overlap of C-13 and H-1 chemical shifts. Such overlap can be resolved using the recently proposed (3,2)D HCCH-COSY, a G-matrix Fourier transform (GFT) NMR based experiment, which facilitates editing of methyl groups into distinct spectral regions by combining their C-13 chemical shifts with that of the neighboring, directly attached, C-13 nucleus. Using this principle, we present three GFT experiments: (a) (4,3)D NOESY-HCCH, (b) (4,3)D H-1-TOCSY-HCCH and (c) (4,3)D C-13-TOCSY-HCCH. These experiments provide unique 4D spectral information rapidly with high sensitivity and resolution for side-chain resonance assignments and NOE analysis of methyl groups. This is exemplified by (4,3)D NOESY-HCCH data acquired for 17.9 kDa non-deuterated cytosolic human J-protein co-chaperone, which provided crucial long-range distance constraints for its 3D structure determination.

PeptideMine - A webserver for the design of peptides for protein-peptide binding studies derived from protein-protein interactomes

Background: Signal transduction events often involve transient, yet specific, interactions between structurally conserved protein domains and polypeptide sequences in target proteins. The identification and validation of these associating domains is crucial to understand signal transduction pathways that modulate different cellular or developmental processes. Bioinformatics strategies to extract and integrate information from diverse sources have been shown to facilitate the experimental design to understand complex biological events. These methods, primarily based on information from high-throughput experiments, have also led to the identification of new connections thus providing hypothetical models for cellular events. Such models, in turn, provide a framework for directing experimental efforts for validating the predicted molecular rationale for complex cellular processes. In this context, it is envisaged that the rational design of peptides for protein-peptide binding studies could substantially facilitate the experimental strategies to evaluate a predicted interaction. This rational design procedure involves the integration of protein-protein interaction data, gene ontology, physico-chemical calculations, domain-domain interaction data and information on functional sites or critical residues. Results: Here we describe an integrated approach called ``PeptideMine'' for the identification of peptides based on specific functional patterns present in the sequence of an interacting protein. This approach based on sequence searches in the interacting sequence space has been developed into a webserver, which can be used for the identification and analysis of peptides, peptide homologues or functional patterns from the interacting sequence space of a protein. To further facilitate experimental validation, the PeptideMine webserver also provides a list of physico-chemical parameters corresponding to the peptide to determine the feasibility of using the peptide for in vitro biochemical or biophysical studies. Conclusions: The strategy described here involves the integration of data and tools to identify potential interacting partners for a protein and design criteria for peptides based on desired biochemical properties. Alongside the search for interacting protein sequences using three different search programs, the server also provides the biochemical characteristics of candidate peptides to prune peptide sequences based on features that are most suited for a given experiment. The PeptideMine server is available at the URL: http://caps.ncbs.res.in/peptidemine

X-ray and molecular-dynamics studies on Mycobacterium leprae single-stranded DNA-binding protein and comparison with other eubacterial SSB structures

The crystal structures of two forms of Mycobacterium leprae single-stranded DNA-binding protein (SSB) have been determined at 2.05 and 2.8 A resolution. Comparison of these structures with the structures of other eubacterial SSBs indicates considerable variation in their quaternary association, although the DNA-binding domains in all of them exhibit the same OB-fold. This variation has no linear correlation with sequence variation, but could be related to variation in protein stability. Molecular-dynamics simulations have been carried out on tetrameric molecules derived from the two forms and the prototype Escherichia coli SSB and the individual subunits of both proteins. Together, the X-ray studies and molecular-dynamics simulations yield information on the relatively rigid and flexible regions of the molecule and on the effect of oligomerization on flexibility. The simulations provide insight into the changes in subunit structure on oligomerization. They also provide insight into the stability and time evolution of the hydrogen bonds/water bridges that connect the two pairs of monomers in the tetramer.