Prediction of protein B-factor profiles

The polypeptide backbones and side chains of proteins are constantly moving due to thermal motion and the kinetic energy of the atoms. The B-factors of protein crystal structures reflect the fluctuation of atoms about their average positions and provide important information about protein dynamics. Computational approaches to predict thermal motion are useful for analyzing the dynamic properties of proteins with unknown structures. In this article, we utilize a novel support vector regression (SVR) approach to predict the B-factor distribution (B-factor profile) of a protein from its sequence. We explore schemes for encoding sequences and various settings for the parameters used in SVR. Based on a large dataset of high-resolution proteins, our method predicts the B-factor distribution with a Pearson correlation coefficient (CC) of 0.53. In addition, our method predicts the B-factor profile with a CC of at least 0.56 for more than half of the proteins. Our method also performs well for classifying residues (rigid vs. flexible). For almost all predicted B-factor thresholds, prediction accuracies (percent of correctly predicted residues) are greater than 70%. These results exceed the best results of other sequence-based prediction methods. (C) 2005 Wiley-Liss, Inc.

Analysis of distinct tartrate-resistant acid phosphatase promoter regions in transgenic mice

The tartrate-resistant acid phosphatase (TRAP) is present in multiple tissues, including kidney, liver, lung, spleen, and bone. Recent study of (TRAP) gene expression has provided evidence for distinct promoters within the (TRAP) gene, suggesting that the gene has alternative, tissue-preferred mRNA transcripts. Examination of endogenous (TRAP) exon 1B and 1C mRNA transcripts revealed tissue-preferred transcript abundance with increased exon 1B transcripts detected in liver and kidney and increased exon 1C transcripts detected in bone and spleen. In this investigation, we have made transgenic mice that express a marker gene driven by two candidate promoters, designated BC and C, within the (TRAP) gene. The BC and C promoters are 2.2 and 1.6 kb, respectively, measured from the translation initiation site. Evaluation of BC transgenic lines demonstrated robust expression in multiple tissues. In contrast, significant transgene expression was not detected in C transgenic lines. Evaluation of transgene mRNAs in BC transgenic lines revealed that virtually all expression was in the form of B transcripts, suggesting that the tissue-preferred pattern of endogenous (TRAP) was not replicated in the BC transgenic line. Likewise, osteoclastogenic cultures from BC, but not C, transgenic bone marrow cells expressed the transgene following receptor activator of NFkappaB ligand/macrophage colony-stimulating factor stimulation. In conclusion, when compared with the 2.2-kb BC portion of the (TRAP) promoter region, the 1.6-kb C portion does not account for significant gene expression in vivo or in vitro; production of the bone- and spleen-preferred (TRAP) C transcript must depend on regulatory elements outside of the 2.2-kb promoter. As the majority of currently investigated transcription factors that influence transcriptional regulation of osteoclast gene expression bind within the 1.6-kb C portion of the (TRAP) promoter, it is likely that transcription binding sites outside of the 2.2-kb region will have profound effects on regulation of the gene in vivo and in vitro.

In vivo analysis of growth hormone receptor signaling domains and their associated transcripts

The growth hormone receptor (GHR) is a critical regulator of postnatal growth and metabolism. However, the GHR signaling domains and pathways that regulate these processes in vivo are not defined. We report the first knock-in mouse models with deletions of specific domains of the receptor that are required for its in vivo actions. Mice expressing truncations at residue m569 (plus Y539/545-F) and at residue m391 displayed a progressive impairment of postnatal growth with receptor truncation. Moreover, after 4 months of age, marked male obesity was observed in both mutant 569 and mutant 391 and was associated with hyperglycemia. Both mutants activated hepatic JAK2 and ERK2, whereas STAT5 phosphorylation was substantially decreased for mutant 569 and absent from mutant 391, correlating with loss of IGF-1 expression and reduction in growth. Microarray analysis of these and GHR(-/-) mice demonstrated that particular signaling domains are responsible for the regulation of different target genes and revealed novel actions of growth hormone. These mice represent the first step in delineating the domains of the GHR regulating body growth and composition and the transcripts associated with these domains.

Functional characterization and genomic organization of the human Na+-sulfate cotransporter hNaS2 gene (SLC13A4)

Sulfate plays an essential role in human growth and development. Here, we characterized the functional properties of the human Na+-sulfate cotransporter (hNaS2), determined its tissue distribution, and identified its gene (SLC13A4) structure. Expression of hNaS2 protein in Xenopus oocytes led to a Na+-dependent transport of sulfate that was inhibited by thiosulfate, phosphate, molybdate. selenate and tungstate, but not by oxalate, citrate, succinate, phenol red or DIDS. Transport kinetics of hNaS2 determined a K, for sulfate of 0.38 mM, suggestive of a high affinity sulfate transporter. Na+ kinetics determined a Hill coefficient of 1.6 +/- 0.6, suggesting a Na: SO42- stoichiometry of 2:1. hNaS2 mRNA was highly expressed in placenta and testis, with intermediate levels in brain and lower levels found in the heart, thymus, and liver. The SLC13A4 gene contains 16 exons, spanning over 47 kb in length. Its 5'-flanking region contains CAAT- and GC-box motifs, and a number of putative transcription factor binding sites, including GATA-1, AP-1, and AP-2 consensus sequences. This is the first study to characterize hNaS2 transport kinetics, define its tissue distribution, and resolve its gene (SLC13A4) structure and 5' flanking region. (C) 2004 Elsevier Inc. All rights reserved.

hnRNP A2, a potential ssDNA/RNA molecular adapter at the telomere

The heterogeneous nuclear ribonucleoprotein (hnRNP) A2 is a multi-tasking protein that acts in the cytoplasm and nucleus. We have explored the possibility that this protein is associated with telomeres and participates in their maintenance. Rat brain hnRNP A2 was shown to have two nucleic acid binding sites. In the presence of heparin one site binds single-stranded oligodeoxyribonucleotides irrespective of sequence but not the corresponding oligoribonucleotides. Both the hnRNP A2-binding cis-acting element for the cytoplasmic RNA trafficking element, A2RE, and the ssDNA telomere repeat match a consensus sequence for binding to a second sequence-specific site identified by mutational analysis. hnRNP A2 protected the telomeric repeat sequence, but not the complementary sequence, against DNase digestion: the glycine-rich domain was found to be necessary, but not sufficient, for protection. The N-terminal RRM (RNA recognition motif) and tandem RRMs of hnRNP A2 also bind the single-stranded, template-containing segment of telomerase RNA. hnRNP A2 colocalizes with telomeric chromatin in the subset of PML bodies that are a hallmark of ALT cells, reinforcing the evidence for hnRNPs having a role in telomere maintenance. Our results support a model in which hnRNP A2 acts as a molecular adapter between single-stranded telomeric repeats, or telomerase RNA, and another segment of ssDNA.

Effect of osmolytes on the interaction of flavin adenine dinucleotide with muscle glycogen phosphorylase b

The effect of three osmolytes, trimethylamine N-oxide (TMAO), betaine and proline, on the interaction of muscle glycogen phosphorylase b with allosteric inhibitor FAD has been examined. In the absence of osmolyte, the interaction is described by a single intrinsic dissociation constant (17.8 muM) for two equivalent and independent binding sites on the dimeric enzyme. However, the addition of osmolytes gives rise to sigmoidal dependencies of fractional enzyme-site saturation upon free inhibitor concentration. The source of this cooperativity has been shown by difference sedimentation velocity to be an osmolyte-mediated isomerization of phosphorylase b to a smaller dimeric state with decreased affinity for FAD. These results thus have substantiated a previous inference that the tendency for osmolyte-enhanced self-association of dimeric glycogen phosphorylase b in the presence of AMP was being countered by the corresponding effect of molecular crowding on an isomerization of dimer to a smaller, nonassociating state. (C) 2004 Elsevier Ltd. Inc. All rights reserved.

Anomalous scattering analysis of Agrobacterium radiobacter phosphotriesterase: the prominent role of iron in the heterobinuclear active site

Bacterial phosphotriesterases are binuclear metalloproteins for which the catalytic mechanism has been studied with a variety of techniques, principally using active sites reconstituted in vitro from apoenzymes. Here, atomic absorption spectroscopy and anomalous X-ray scattering have been used to determine the identity of the metals incorporated into the active site in vivo. We have recombinantly expressed the phosphotriesterase from Agrobacterium radiobacter (OpdA) in Escherichia coli grown in medium supplemented with 1 mM CoCl2 and in unsupplemented medium. Anomalous scattering data, collected from a single crystal at the Fe-K, Co-K and Zn-K edges, indicate that iron and cobalt are the primary constituents of the two metal-binding sites in the catalytic centre (alpha and P) in the protein expressed in E. coli grown in supplemented medium. Comparison with OpdA expressed in unsupplemented medium demonstrates that the cobalt present in the supplemented medium replaced zinc at the beta-position of the active site, which results in an increase in the catalytic efficiency of the enzyme. These results suggest an essential role for iron in the catalytic mechanism of bacterial phosphotriesterases, and that these phosphotriesterases are natively heterobinuclear iron-zinc enzymes.

The mitogenic signaling pathway for fibroblast growth factor-2 involves the tyrosine phosphorylation of cyclin D2 in MCF-7 human breast cancer cells

Fibroblast growth factor-2 (FGF-2) is mitogenic for the human breast cancer cell line MCF-7; here we investigate some of the signaling pathways subserving this activity. FGF-2 stimulation of MCF-7 cells resulted in a global increase of intracellular tyrosine phosphorylation of proteins, particularly FGF receptor substrate-2, the protooncogene product Src and the mitogen-activated protein kinase (MAP kinase) cascade, A major increase in the tyrosine phosphorylation of a 30-kDa protein species was also found. This protein was identified as cyclin D2 by mass spectrometry after trypsin digestion. Immunoprecipitation of cyclin D2 and immunoblotting with anti-phosphotyrosine antibodies confirmed that the tyrosine phosphorylation of cyclin D2 was indeed induced by FGF-2 stimulation. In addition, pharmacological inhibition of Src (with herbimycin A and PP2), and of the MAP kinase cascade (with PD98059), confirmed that Src activity is required for the FGF-2-induced phosphorylation of cyclin D2 whereas MAP kinase activity is not, Thus, tyrosine phosphorylation of cyclin D2 may be a hey regulatory target for FGF-2 signaling. (C) 2000 Federation of European Biochemical Societies. Published by Elsevier Science B.V. All rights reserved.

Comparison of the expression of calcitonin receptor-like receptor (CRLR) and receptor activity modifying proteins (RAMPs) with CGRP and adrenomedullin binding in cell lines

1. The calcitonin receptor-like receptor (CRLR) and specific receptor activity modifying proteins (RAMPs) together form receptors for calcitonin gene-related peptide (CGRP) and/or adrenomedullin in transfected cells. 2. There is less evidence that innate CGRP and adrenomedullin receptors are formed by CRLR/RAMP combinations. We therefore examined whether CGRP and/or adrenomedullin binding correlated with CRLR and RAMP mRNA expression in human and rat cell lines known to express these receptors. Specific human or rat CRLR antibodies were used to examine the presence of CRLR in these cells. 3. We confirmed CGRP subtype 1 receptor (CGRP(1)) pharmacology in SK-N-MC neuroblastoma cells. L6 myoblast cells expressed both CGRP(1) and adrenomedullin receptors whereas Rat-2 fibroblasts expressed only adrenomedullin receptors. In contrast we could not confirm CGRP(2) receptor pharmacology for Col-29 colonic epithelial cells, which, instead were CGRP(1)-like in this study. 4. L6, SK-N-MC and Col-29 cells expressed mRNA for RAMP1 and RAMP2 but Rat-2 fibroblasts had only RAMP2. No cell line had detectable RAMP3 mRNA. 5. SK-N-MC, Col-29 and Rat-2 fibroblast cells expressed CRLR mRNA. By contrast, CRLR mRNA was undetectable by Northern analysis in one source of L6 cells. Conversely, a different source of L6 cells had mRNA for CRLR. All of the cell lines expressed CRLR protein. Thus circumstances where CRLR mRNA is apparently absent by Northern analysis do not exclude the presence of this receptor. 6. These data strongly support CRLR, together with appropriate RAMPs as binding sites for CGRP and adrenomedullin in cultured cells.

Thermodynamics of binding of regulatory ligands to tissue transglutaminase

The transamidating activity of tissue transglutaminase is regulated by the ligands calcium and GTP, via conformational changes which facilitate or interfere with interaction with the peptidyl-glutamine substrate. We have analysed binding of these ligands by calorimetric and computational approaches. In the case of GTP we have detected a single high affinity site (K (D) approximately 1 muM), with moderate thermal effects suggestive that binding GTP involves replacement of GDP, normally bound to the protein. On line with this possibility no significant binding was observed during titration with GDP and computational studies support this view. Titration with calcium at a high cation molar excess yielded a complex binding isotherm with a number of "apparent binding sites" in large excess over those detectable by equilibrium dialysis (6 sites). This binding pattern is ascribed to occurrence of additional thermal contributions, beyond those of binding, due to the occurrence of conformational changes and to catalysis itself (with protein self-crosslinking). In contrast only one site for binding calcium with high affinity (K (D) approximately 0.15 muM) is observed with samples of enzyme inactivated by alkylation at the active site (to prevent enzyme crosslinkage and thermal effects of catalysis). These results indicate an intrinsic ability of tissue transglutaminase to bind calcium with high affinity and the necessity of careful reassessment of the enzyme regulatory pattern in relation to the concentrations of ligands in living cells, taking also in account effects of ligands on protein subcellular compartimentation.

PDNAsite:identification of DNA-binding site from protein sequence by incorporating spatial and sequence context

Protein-DNA interactions are involved in many fundamental biological processes essential for cellular function. Most of the existing computational approaches employed only the sequence context of the target residue for its prediction. In the present study, for each target residue, we applied both the spatial context and the sequence context to construct the feature space. Subsequently, Latent Semantic Analysis (LSA) was applied to remove the redundancies in the feature space. Finally, a predictor (PDNAsite) was developed through the integration of the support vector machines (SVM) classifier and ensemble learning. Results on the PDNA-62 and the PDNA-224 datasets demonstrate that features extracted from spatial context provide more information than those from sequence context and the combination of them gives more performance gain. An analysis of the number of binding sites in the spatial context of the target site indicates that the interactions between binding sites next to each other are important for protein-DNA recognition and their binding ability. The comparison between our proposed PDNAsite method and the existing methods indicate that PDNAsite outperforms most of the existing methods and is a useful tool for DNA-binding site identification. A web-server of our predictor (http://hlt.hitsz.edu.cn:8080/PDNAsite/) is made available for free public accessible to the biological research community.

Specific binding of the mammalian high mobility group protein AT-hook 2 to the minor groove of AT -rich DNAs: Thermodynamic and specificity studies

The mammalian high mobility group protein AT-hook 2 (HMGA2) is a small transcriptional factor involved in cell development and oncogenesis. It contains three "AT-hook" DNA binding domains, which specifically recognize the minor groove of AT-rich DNA sequences. It also has an acidic C-terminal motif. Previous studies showed that HMGA2 mediates all its biological effects through interactions with AT-rich DNA sequences in the promoter regions. In this dissertation, I used a variety of biochemical and biophysical methods to examine the physical properties of HMGA2 and to further investigate HMGA2's interactions with AT-rich DNA sequences. The following are three avenues perused in this study: (1) due to the asymmetrical charge distribution of HMGA2, I have developed a rapid procedure to purify HMGA2 in the milligram range. Preparation of large amounts of HMGA2 makes biophysical studies possible; (2) Since HMGA2 binds to different AT-rich sequences in the promoter regions, I used a combination of isothermal titration calorimetry (ITC) and DNA UV melting experiment to characterize interactions of HMGA2 with poly(dA-dT) 2 and poly(dA)poly(dT). My results demonstrated that (i) each HMGA2 molecule binds to 15 AT bp; (ii) HMGA2 binds to both AT DNAs with very high affinity. However, the binding reaction of HMGA2 to poly(dA-dT) 2 is enthalpy-driven and the binding reaction of HMGA2 with poly(dA)poly(dT) is entropy-driven; (iii) the binding reactions are strongly depended on salt concentrations; (3) Previous studies showed that HMGA2 may have sequence specificity. In this study, I used a PCR-based SELEX procedure to examine the DNA binding specificity of HMGA2. Two consensus sequences for HMGA2 have been identified: 5'-ATATTCGCGAWWATT-3' and 5'-ATATTGCGCAWWATT-3', where W represents A or T. These consensus sequences have a unique feature: the first five base pairs are AT-rich, the middle four to five base pairs are GC-rich, and the last five to six base pairs are AT-rich. All three segments are critical for high affinity binding. Replacing either one of the AT-rich sequences to a non-AT-rich sequence causes at least 100-fold decrease in the binding affinity. Intriguingly, if the GC-segment is substituted by an AT-rich segment, the binding affinity of HMGA2 is reduced approximately 5-fold. Identification of the consensus sequences for HMGA2 represents an important step towards finding its binding sites within the genome.

Characterization of a Full-Length TTP Family Member Association with RNA Sequence Elements

Post-transcriptional regulation of cytoplasmic mRNAs is an efficient mechanism of regulating the amounts of active protein within a eukaryotic cell. RNA sequence elements located in the untranslated regions of mRNAs can influence transcript degradation or translation through associations with RNA-binding proteins. Tristetraprolin (TTP) is the best known member of a family of CCCH zinc finger proteins that targets adenosine-uridine rich element (ARE) binding sites in the 3’ untranslated regions (UTRs) of mRNAs, promoting transcript deadenylation through the recruitment of deadenylases. More specifically, TTP has been shown to bind AREs located in the 3’-UTRs of transcripts with known roles in the inflammatory response. The mRNA-binding region of the protein is the highly conserved CCCH tandem zinc finger (TZF) domain. The synthetic TTP TZF domain has been shown to bind with high affinity to the 13-mer sequence of UUUUAUUUAUUUU. However, the binding affinities of full-length TTP family members to the same sequence and its variants are unknown. Furthermore, the distance needed between two overlapping or neighboring UUAUUUAUU 9-mers for tandem binding events of a full-length TTP family member to a target transcript has not been explored. To address these questions, we recombinantly expressed and purified the full-length C. albicans TTP family member Zfs1. Using full-length Zfs1, tagged at the N-terminus with maltose binding protein (MBP), we determined the binding affinities of the protein to the optimal TTP binding sequence, UUAUUUAUU. Fluorescence anisotropy experiments determined that the binding affinities of MBP-Zfs1 to non-canonical AREs were influenced by ionic buffer strength, suggesting that transcript selectivity may be affected by intracellular conditions. Furthermore, electrophoretic mobility shift assays (EMSAs) revealed that separation of two core AUUUA sequences by two uridines is sufficient for tandem binding of MBP-Zfs1. Finally, we found evidence for tandem binding of MBP-Zfs1 to a 27-base RNA oligonucleotide containing only a single ARE-binding site, and showed that this was concentration and RNA length dependent; this phenomenon had not been seen previously. These data suggest that the association of the TTP TZF domain and the TZF domains of other species, to ARE-binding sites is highly conserved. Domains outside of the TZF domain may mediate transcript selectivity in changing cellular conditions, and promote protein-RNA interactions not associated with the ARE-binding TZF domain.

In summary, the evidence presented here suggests that Zfs1-mediated decay of mRNA targets may require additional interactions, in addition to ARE-TZF domain associations, to promote transcript destabilization and degradation. These studies further our understanding of post-transcriptional steps in gene regulation.

RNA Localization and Translational Regulation on the Endoplasmic Reticulum

mRNA localization is emerging as a critical cellular mechanism for the spatiotemporal regulation of protein expression and serves important roles in oogenesis, embryogenesis, cell fate specification, and synapse formation. Signal sequence-encoding mRNAs are localized to the endoplasmic reticulum (ER) membrane by either of two mechanisms, a canonical mechanism of translation on ER-bound ribosomes (signal recognition particle pathway), or a poorly understood direct ER anchoring mechanism. In this study, we identify that the ER integral membrane proteins function as RNA-binding proteins and play important roles in the direct mRNA anchoring to the ER. We report that one of the ER integral membrane RNA-binding protein, AEG-1 (astrocyte elevated gene-1), functions in the direct ER anchoring and translational regulation of mRNAs encoding endomembrane transmembrane proteins. HITS-CLIP and PAR-CLIP analyses of the AEG-1 mRNA interactome of human hepatocellular carcinoma cells revealed a high enrichment for mRNAs encoding endomembrane organelle proteins, most notably encoding transmembrane proteins. AEG-1 binding sites were highly enriched in the coding sequence and displayed a signature cluster enrichment downstream of encoded transmembrane domains. In overexpression and knockdown models, AEG-1 expression markedly regulates translational efficiency and protein functions of two of its bound transcripts, MDR1 and NPC1. This study reveals a molecular mechanism for the selective localization of mRNAs to the ER and identifies a novel post-transcriptional gene regulation function for AEG-1 in membrane protein expression.