Repertoire of Protein Kinases Encoded in the Genome of Takifugu rubripes

Takifugu rubripes is teleost fish widely used in comparative genomics to understand the human system better due to its similarities both in number of genes and structure of genes. In this work we survey the fugu genome, and, using sensitive computational approaches, we identify the repertoire of putative protein kinases and classify them into groups and subfamilies. The fugu genome encodes 519 protein kinase-like sequences and this number of putative protein kinases is comparable closely to that of human. However, in spite of its similarities to human kinases at the group level, there are differences at the subfamily level as noted in the case of KIS and DYRK subfamilies which contribute to differences which are specific to the adaptation of the organism. Also, certain unique domain combination of galectin domain and YkA domain suggests alternate mechanisms for immune response and binding to lipoproteins. Lastly, an overall similarity with the MAPK pathway of humans suggests its importance to understand signaling mechanisms in humans. Overall the fugu serves as a good model organism to understand roles of human kinases as far as kinases such as LRRK and IRAK and their associated pathways are concerned.

Subdiffusion in hair bundle dynamics: The role of protein conformational fluctuations

The detection of sound signals in vertebrates involves a complex network of different mechano-sensory elements in the inner ear. An especially important element in this network is the hair bundle, an antenna-like array of stereocilia containing gated ion channels that operate under the control of one or more adaptation motors. Deflections of the hair bundle by sound vibrations or thermal fluctuations transiently open the ion channels, allowing the flow of ions through them, and producing an electrical signal in the process, eventually causing the sensation of hearing. Recent high frequency (0.1-10 kHz) measurements by Kozlov et al. Proc. Natl. Acad. Sci. U. S. A. 109, 2896 (2012)] of the power spectrum and the mean square displacement of the thermal fluctuations of the hair bundle suggest that in this regime the dynamics of the hair bundle are subdiffusive. This finding has been explained in terms of the simple Brownian motion of a filament connecting neighboring stereocilia (the tip link), which is modeled as a viscoelastic spring. In the present paper, the diffusive anomalies of the hair bundle are ascribed to tip link fluctuations that evolve by fractional Brownian motion, which originates in fractional Gaussian noise and is characterized by a power law memory. The predictions of this model for the power spectrum of the hair bundle and its mean square displacement are consistent with the experimental data and the known properties of the tip link. (C) 2012 American Institute of Physics. http://dx.doi.org/10.1063/1.4768902]

Distinctive contributions of the ribosomal P-site elements m(2)G966, m(5)C967 and the C-terminal tail of the S9 protein in the fidelity of initiation of translation in Escherichia coli

The accuracy of pairing of the anticodon of the initiator tRNA (tRNA(fMet)) and the initiation codon of an mRNA, in the ribosomal P-site, is crucial for determining the translational reading frame. However, a direct role of any ribosomal element(s) in scrutinizing this pairing is unknown. The P-site elements, m(2)G966 (methylated by RsmD), m(5)C967 (methylated by RsmB) and the C-terminal tail of the protein S9 lie in the vicinity of tRNA(fMet). We investigated the role of these elements in initiation from various codons, namely, AUG, GUG, UUG, CUG, AUA, AUU, AUC and ACG with tRNA(CAU)(fmet) (tRNA(fMet) with CAU anticodon); CAC and CAU with tRNA(GUG)(fme); UAG with tRNA(GAU)(fMet) using in vivo and computational methods. Although RsmB deficiency did not impact initiation from most codons, RsmD deficiency increased initiation from AUA, CAC and CAU (2- to 3.6-fold). Deletion of the S9 C-terminal tail resulted in poorer initiation from UUG, GUG and CUG, but in increased initiation from CAC, CAU and UAC codons (up to 4-fold). Also, the S9 tail suppressed initiation with tRNA(CAU)(fMet)lacking the 3GC base pairs in the anticodon stem. These observations suggest distinctive roles of 966/967 methylations and the S9 tail in initiation.

Human la protein interaction with GCAC near the initiator AUG enhances hepatitis C virus RNA replication by promoting linkage between 5 ` and 3 ` untranslated regions

Human La protein has been implicated in facilitating the internal initiation of translation as well as replication of hepatitis C virus (HCV) RNA. Previously, we demonstrated that La interacts with the HCV internal ribosome entry site (IRES) around the GCAC motif near the initiator AUG within stem-loop IV by its RNA recognition motif (RRM) (residues 112 to 184) and influences HCV translation. In this study, we have deciphered the role of this interaction in HCV replication in a hepatocellular carcinoma cell culture system. We incorporated mutation of the GCAC motif in an HCV monocistronic subgenomic replicon and a pJFH1 construct which altered the binding of La and checked HCV RNA replication by reverse transcriptase PCR (RT-PCR). The mutation drastically affected HCV replication. Furthermore, to address whether the decrease in replication is a consequence of translation inhibition or not, we incorporated the same mutation into a bicistronic replicon and observed a substantial decrease in HCV RNA levels. Interestingly, La overexpression rescued this inhibition of replication. More importantly, we observed that the mutation reduced the association between La and NS5B. The effect of the GCAC mutation on the translation-to-replication switch, which is regulated by the interplay between NS3 and La, was further investigated. Additionally, our analyses of point mutations in the GCAC motif revealed distinct roles of each nucleotide in HCV replication and translation. Finally, we showed that a specific interaction of the GCAC motif with human La protein is crucial for linking 5' and 3' ends of the HCV genome. Taken together, our results demonstrate the mechanism of regulation of HCV replication by interaction of the cis-acting element GCAC within the HCV IRES with human La protein.

Redox regulation of protein tyrosine phosphatase 1B (PTP1B): Importance of steric and electronic effects on the unusual cyclization of the sulfenic acid intermediate to a sulfenyl amide

The redox regulation of protein tyrosine phosphatase 1B (PTP1B) via the unusual transformation of its sulfenic acid (PTP1B-SOH) to a cyclic sulfenyl amide intermediate is studied by using small molecule chemical models. These studies suggest that the sulfenic acids derived from the H2O2-mediated reactions o-amido thiophenols do not efficiently cyclize to sulfenyl amides and the sulfenic acids produced in situ can be trapped by using methyl iodide. Theoretical calculations suggest that the most stable conformer of such sulfenic acids are stabilized by n(O) -> sigma* (S-OH) orbital interactions, which force the -OH group to adopt a position trans to the S center dot center dot center dot O interaction, leading to an almost linear arrangement of the O center dot center dot center dot S-O moiety and this may be the reason for the slow cyclization of such sulfenic acids to their corresponding sulfenyl amides. On the other hand, additional substituents at the 6-position of o-amido phenylsulfenic acids that can induce steric environment and alter the electronic properties around the sulfenic acid moiety by S center dot center dot center dot N or S center dot center dot center dot O nonbonded interactions destabilize the sulfenic acids by inducing strain in the molecule. This may lead to efficient the cyclization of such sulfenic acids. This model study suggests that the amino acid residues in the close proximity of the sulfenic acid moiety in PTP1B may play an important role in the cyclization of PTP1B-SOH to produce the corresponding sulfenyl amide.

An automated approach to network features of protein structure ensembles

Network theory applied to protein structures provides insights into numerous problems of biological relevance. The explosion in structural data available from PDB and simulations establishes a need to introduce a standalone-efficient program that assembles network concepts/parameters under one hood in an automated manner. Herein, we discuss the development/application of an exhaustive, user-friendly, standalone program package named PSN-Ensemble, which can handle structural ensembles generated through molecular dynamics (MD) simulation/NMR studies or from multiple X-ray structures. The novelty in network construction lies in the explicit consideration of side-chain interactions among amino acids. The program evaluates network parameters dealing with topological organization and long-range allosteric communication. The introduction of a flexible weighing scheme in terms of residue pairwise cross-correlation/interaction energy in PSN-Ensemble brings in dynamical/chemical knowledge into the network representation. Also, the results are mapped on a graphical display of the structure, allowing an easy access of network analysis to a general biological community. The potential of PSN-Ensemble toward examining structural ensemble is exemplified using MD trajectories of an ubiquitin-conjugating enzyme (UbcH5b). Furthermore, insights derived from network parameters evaluated using PSN-Ensemble for single-static structures of active/inactive states of 2-adrenergic receptor and the ternary tRNA complexes of tyrosyl tRNA synthetases (from organisms across kingdoms) are discussed. PSN-Ensemble is freely available from http://vishgraph.mbu.iisc.ernet.in/PSN-Ensemble/psn_index.html.

Comparative Study of Protein Unfolding in Aqueous Urea and Dimethyl Sulfoxide Solutions: Surface Polarity, Solvent Specificity, and Sequence of Secondary Structure Melting

Elucidation of possible pathways between folded (native) and unfolded states of a protein is a challenging task, as the intermediates are often hard to detect. Here, we alter the solvent environment in a controlled manner by choosing two different cosolvents of water, urea, and dimethyl sulfoxide (DMSO) and study unfolding of four different proteins to understand the respective sequence of melting by computer simulation methods. We indeed find interesting differences in the sequence of melting of alpha helices and beta sheets in these two solvents. For example, in 8 M urea solution, beta-sheet parts of a protein are found to unfold preferentially, followed by the unfolding of alpha helices. In contrast, 8 M DMSO solution unfolds alpha helices first, followed by the separation of beta sheets for the majority of proteins. Sequence of unfolding events in four different alpha/beta proteins and also in chicken villin head piece (HP-36) both in urea and DMSO solutions demonstrate that the unfolding pathways are determined jointly by relative exposure of polar and nonpolar residues of a protein and the mode of molecular action of a solvent on that protein.

STAT-1 expression is regulated by IGFBP-3 in malignant glioma cells and is a strong predictor of poor survival in patients with glioblastoma Laboratory investigation

Object. Insulin-like growth factor binding proteins (IGEBPs) have been implicated in the pathogenesis of glioma. In a previous study the authors demonstrated that IGFBP-3 is a novel glioblastoma biomarker associated with poor survival. Since signal transducer and activator of transcription 1 (STAT-1) has been shown to be regulated by IGFBP-3 during chondrogenesis and is a prosurvival and radioresistant molecule in different tumors, the aim in the present study was to explore the functional significance of IGFBP-3 in malignant glioma cells, to determine if STAT-1 is indeed regulated by IGFBP-3, and to study the potential of STAT-1 as a biomarker in glioblastoma. Methods. The functional significance of IGFBP-3 was investigated using the short hairpin (sh)RNA gene knockdown approach on U251MG cells. STAT-1 regulation by IGFBP-3 was tested on U251MG and U87MG cells by shRNA gene knockdown and exogenous treatment with recombinant IGFBP-3 protein. Subsequently, the expression of STAT-1 was analyzed with real-time reverse transcription polymerase chain reaction (RT-PCR) and immunohistochemistry (IHC) in glioblastoma and control brain tissues. Survival analyses were done on a uniformly treated prospective cohort of adults with newly diagnosed glioblastoma (136 patients) using Kaplan-Meier and Cox regression models. Results. IGFBP-3 knockdown significantly impaired proliferation, motility, migration, and invasive capacity of U251MG cells in vitro (p < 0.005). Exogenous overexpression of IGFBP-3 in U251MG and U87MG cells demonstrated STAT-1 regulation. The mean transcript levels (by real-time RT-PCR) and the mean labeling index of STAT-1 (by IHC) were significantly higher in glioblastoma than in control brain tissues (p = 0.0239 and p < 0.001, respectively). Multivariate survival analysis revealed that STAT-1 protein expression (HR 1.015, p = 0.033, 95% CI 1.001-1.029) along with patient age (HR 1.025, p = 0.005, 95% CI 1.008-1.042) were significant predictors of shorter survival in patients with glioblastoma. Conclusions. IGFBP-3 influences tumor cell proliferation, migration, and invasion and regulates STAT-1 expression in malignant glioma cells. STAT-1 is overexpressed in human glioblastoma tissues and emerges as a novel prognostic biomarker.

Role of the loop L-4,L-5 in allosteric regulation in mtHsp70s: in vivo significance of domain communication and its implications in protein translocation

Mitochondrial Hsp70 (mtHsp70) is essential for a vast repertoire of functions, including protein import, and requires effective interdomain communication for efficient partner-protein interactions. However, the in vivo functional significance of allosteric regulation in eukaryotes is poorly defined. Using integrated biochemical and yeast genetic approaches, we provide compelling evidence that a conserved substrate-binding domain (SBD) loop, L-4,L-5, plays a critical role in allosteric communication governing mtHsp70 chaperone functions across species. In yeast, a temperature-sensitive L-4,L-5 mutation (E467A) disrupts bidirectional domain communication, leading to compromised protein import and mitochondrial function. Loop L-4,L-5 functions synergistically with the linker in modulating the allosteric interface and conformational transitions between SBD and the nucleotide-binding domain (NBD), thus regulating interdomain communication. Second-site intragenic suppressors of E467A isolated within the SBD suppress domain communication defects by conformationally altering the allosteric interface, thereby restoring import and growth phenotypes. Strikingly, the suppressor mutations highlight that restoration of communication from NBD to SBD alone is the minimum essential requirement for effective in vivo function when primed at higher basal ATPase activity, mimicking the J-protein-bound state. Together these findings provide the first mechanistic insights into critical regions within the SBD of mtHsp70s regulating interdomain communication, thus highlighting its importance in protein translocation and mitochondrial biogenesis.

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.

Inhibition of protein synthesis leading to unfolded protein response is the major event in abrin-mediated apoptosis

Abrin obtained from the plant Abrus precatorius inhibits protein synthesis and also triggers apoptosis in cells. Previous studies from our laboratory suggested a link between these two events. Using an active site mutant of abrin A-chain which exhibits 225-fold lower protein synthesis inhibitory activity than the wild-type abrin A-chain, we demonstrate in this study that inhibition of protein synthesis induced by abrin is the major factor triggering unfolded protein response leading to apoptosis. Since abrin A-chain requires the B-chain for internalization into cells, the wild-type and mutant recombinant abrin A-chains were conjugated to native ricin B-chain to generate hybrid toxins, and the toxic effects of the two conjugates were compared. The rate of inhibition of protein synthesis mediated by the mutant ricin B-rABRA (R167L) conjugate was slower than that of the wild-type ricin B-rABRA conjugate as expected. The mutant conjugate activated p38MAPK and caspase-3 similar to its wild-type counterpart although at later time points. Overall, these results confirm that inhibition of protein synthesis is the major event contributing to abrin-mediated apoptosis.

REPERTOIRE OF PROTEIN KINASES ENCODED IN THE GENOME OF ZEBRAFISH SHOWS REMARKABLY LARGE POPULATION OF PIM KINASES

In recent times, zebrafish has garnered lot of popularity as model organism to study human cancers. Despite high evolutionary divergence from humans, zebrafish develops almost all types of human tumors when induced. However, mechanistic details of tumor formation have remained largely unknown. Present study is aimed at analysis of repertoire of kinases in zebrafish proteome to provide insights into various cellular components. Annotation using highly sensitive remote homology detection methods revealed ``substantial expansion'' of Ser/Thr/Tyr kinase family in zebrafish compared to humans, constituting over 3% of proteome. Subsequent classification of kinases into subfamilies revealed presence of large number of CAMK group of kinases, with massive representation of PIM kinases, important for cell cycle regulation and growth. Extensive sequence comparison between human and zebrafish PIM kinases revealed high conservation of functionally important residues with a few organism specific variations. There are about 300 PIM kinases in zebrafish kinome, while human genome codes for only about 500 kinases altogether. PIM kinases have been implicated in various human cancers and are currently being targeted to explore their therapeutic potentials. Hence, in depth analysis of PIM kinases in zebrafish has opened up new avenues of research to verify the model organism status of zebrafish.

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.

On the Fourier refinement of protein structures

The situation normally encountered in the high-resolution refinement of protein structures is one in which the inaccurate positions of P out of a total of N atoms are known whereas those of the remaining atoms are unknown. Fourier maps with coefficients (FN -- F'P) × exp (i[alpha]'P) and (mFN -- nF'P) exp (i[alpha]'P), where FN is the observed structure factor and F'P and [alpha]'P are the magnitude and the phase angle of the calculated structure factor corresponding to the inaccurate atomic positions, are often used to correct the positions of the P atoms and to determine those of the Q unknown atoms. A general theoretical approach is presented to elucidate the effect of errors in the positions of the known atoms on the corrected positions of the known atoms and the positions of the unknown atoms derived from such maps. The theory also leads to the optimal choice of parameters used in the different syntheses. When the errors in the positions of the input atoms are systematic, their effects are not taken care of automatically by the syntheses.