Transglutaminase type 2 affects cell migration through post-translational modification of platelet-derived growth factor-BB.

Migration is a key cellular function with important implications in cell physiology. Impairment of such function is observed in angiogenesis, cancer, central nervous system development, and many other physiological and pathological events. Serum is considered among the most potent physiological chemotactic stimuli. Transglutaminase 2 (TG2) is involved in most of the mentioned processes, suggesting the hypothesis that TG2 may modulate cell movement and chemotaxis by acting on serum factors. Cell biology and biochemistry studies confirmed this hypothesis, showing that human serum contains potent chemotactic signals significantly impaired by activated TG2. Bioinformatics studies indicated that one potent serum factor potential substrate of TG2-dependent transamidation is platelet-derived growth factor-BB (PDGF-BB). Cell biology and immunometric experiments carried out with U87MG human glioma cell line showed that human recombinant PDGF-BB pre-incubated with calcium-activated TG2 lost about 70 % of its chemotactic activity and antigenicity. These data indicate that PDGF-BB is a substrate of TG2-transamidating activity, and such modification may play a key role in the modulation of PDGF's chemotactic features. Further, these findings suggest a novel point of view to study the extracellular functions of TG2 and to understand how protein signals, such as growth factors and cytokines, act in the extracellular space to reach their specific targets.

Tissue transglutaminase activity protects from cutaneous melanoma metastatic dissemination: an in vivo study.

The role of tissue transglutaminase (TG-2, TGase-2) in cancer development is still a fascinating field of research. The available reports do not elucidate fully its mechanism of action, due to the limitations of in vitro approaches. Therefore, to understand TG-2 role in cancer, we carried out an in vivo study with a more direct approach. TG-2 was in vivo overexpressed in a murine model of melanoma (intravenous injection of B16 melanoma cells in C57BL/6N mice) by means of a plasmid carrying the TG-2 cDNA. The evaluation of the frequency and size of the metastases indicated that the number of melanoma lung foci was more markedly reduced by TG-2 overexpression than the metastatic size. Then, TG-2 overexpressing mice showed a prolonged survival with respect to control mice. Further analyses were carried by means of proteomic analysis of melanoma cell lysates and meta-analysis of published transcriptomic datasets. Proteomic analysis of cell lysates from a human melanoma cell line compared to human keratinocytes showed significant differences in the expression of TG-2 substrates known to be involved in proliferation/differentiation and cancer progression. Taken together, these findings indicate a protective role of TG-2 enzymatic activity in melanoma progression in vivo.

Functional hypothesis on miraculin' sweetness by a molecular dynamics approach.

Miraculin differs from other sweet-tasting proteins because it is a taste-modifier having the unusual property of modifying sourness into sweetness. Its dimer is covalently linked by an inter-chain disulphide bond, and shows its taste-modifying activity at acidic pH, with maximum at pH 3.0, while it is flat at neutral pH. Previous studies suggested the importance of two histidine residues for the taste-modifying activity of miraculin. In this work, we have conducted molecular dynamics simulations on wild type miraculin and on three mutated dimers (H29A, H59A and H29A/H59A) both at neutral and acidic pH to investigate the structural and functional role of these two His residues. Our results suggested that at acidic pH the presence of two charged His at the interface induced a structural rearrangement of the two monomers, thus leading to their relative opening and the following adaptation of their conformation to the receptor surface. On the other hand the simulations on three mutants showed that the mutated dimers had a closed form, and highlighted the important role of H29 in stabilizing/destabilizing the dimer arrangement and also a cooperative effect of the two histidines.

Molecular characterisation and structural analysis of an interferon homologue in sea bass (Dicentrarchus labrax L.).

The interferons (IFNs) are a large family of soluble cytokines involved in the immune response against viral pathogens. Three families of IFNs have been identified in mammals (type I, type II and type III) and, recently, homologues of type I and type II genes have been found in various teleost fish species. In this paper we report the cloning of a cDNA encoding an type I IFN molecule from sea bass (Dicentrarchus labrax L.), its expression analysis and gene structure and, finally, its 3D structure obtained by template-based modelling. The sea bass IFN cDNA consists of 1047bp that translates in one reading frame to give the entire molecule containing 185 amino acids. The analysis of the sequence revealed the presence of a putative 22 amino acid signal peptide, two cysteine residues and three potential N-glycosylation sites. The sea bass IFN gene contains four introns as with other type I IFN teleost genes, except medaka that contains three introns. Real time PCR was performed after poly I:C stimulation of DLEC cell line to investigate the expression of sea bass IFN and Mx and an induction was observed for both genes. The predicted 3D structure of sea bass IFN is characterized by an "all-alpha" domain that shows an "up-down bundle" architecture made of six helices (ABB'CDE). The two cysteine residues present in the sequence (i.e. Cys(23) and Cys(126)) are in a position and at a distance that suggest the possible formation of a disulfide bridge that may stabilize the structure. Our results will give the opportunity to investigate more in detail antiviral immune responses in sea bass and add to studies on the evolution of the IFN system in teleosts and vertebrates more generally.

A CD4 homologue in sea bass (Dicentrarchus labrax): molecular characterization and structural analysis.

CD4 is a transmembrane glycoprotein fundamental for cell-mediated immunity. Its action as a T cell co-receptor increases the avidity of association between a T cell and an antigen-presenting cell by interacting with portions of the complex between MHC class II and TR molecules. In this paper we report the cDNA cloning, expression and structural analysis of a CD4 homologue from sea bass (Dicentrarchus labrax). The sea bass CD4 cDNA consists of 2071 bp that translates in one reading frame to give the entire molecule containing 480 amino acids. The analysis of the sequence shows the presence of four putative Ig-like domains and that some fundamental structural features, like a disulphide bond in domain D2 and the CXC signalling motif in the cytoplasmic tail, are conserved from sea bass to mammals. Real-time PCR analysis showed that very high levels of CD4 mRNA transcripts are present in thymus, followed by gut and gills. In vitro stimulation of head kidney leukocytes with LPS and PHA-L gave an increase of CD4 mRNA levels after 4h and a decrease after 24h. Homology modelling has been applied to create a 3D model of sea bass CD4 and to investigate its interaction with sea bass MHC-II. The analysis of the 3D complex between sea bass CD4 and sea bass MHC-II suggests that the absence of a disulfide bond in the CD4 D1 domain could make this molecule more flexible, inducing a different conformation and affecting the binding and the way of interaction between CD4 and MHC-II. Our results will add new insights into the sea bass T cell immune responses and will help in the identification of T cell subsets in teleost fishes to better understand the evolution of cell-mediated immunity from fish to mammals.

Time-course analysis of genome-wide gene expression data from hormone-responsive human breast cancer cells.

BACKGROUND: Microarray experiments enable simultaneous measurement of the expression levels of virtually all transcripts present in cells, thereby providing a 'molecular picture' of the cell state. On the other hand, the genomic responses to a pharmacological or hormonal stimulus are dynamic molecular processes, where time influences gene activity and expression. The potential use of the statistical analysis of microarray data in time series has not been fully exploited so far, due to the fact that only few methods are available which take into proper account temporal relationships between samples. RESULTS: We compared here four different methods to analyze data derived from a time course mRNA expression profiling experiment which consisted in the study of the effects of estrogen on hormone-responsive human breast cancer cells. Gene expression was monitored with the innovative Illumina BeadArray platform, which includes an average of 30-40 replicates for each probe sequence randomly distributed on the chip surface. We present and discuss the results obtained by applying to these datasets different statistical methods for serial gene expression analysis. The influence of the normalization algorithm applied on data and of different parameter or threshold choices for the selection of differentially expressed transcripts has also been evaluated. In most cases, the selection was found fairly robust with respect to changes in parameters and type of normalization. We then identified which genes showed an expression profile significantly affected by the hormonal treatment over time. The final list of differentially expressed genes underwent cluster analysis of functional type, to identify groups of genes with similar regulation dynamics. CONCLUSIONS: Several methods for processing time series gene expression data are presented, including evaluation of benefits and drawbacks of the different methods applied. The resulting protocol for data analysis was applied to characterization of the gene expression changes induced by estrogen in human breast cancer ZR-75.1 cells over an entire cell cycle.

Molecular modelling of miraculin: Structural analyses and functional hypotheses.

Miraculin is a plant protein that displays the peculiar property of modifying taste by swiching sour into a sweet taste. Its monomer is flavourless at all pH as well as at high concentration; the dimer form elicits its taste-modifying activity at acidic pH; a tetrameric form is also reported as active. Two histidine residues, located in exposed regions, are the main responsible of miraculin activity, as demonstrated by mutagenesis studies. Since structural data of miraculin are not available, we have predicted its three-dimensional structure and simulated both its dimer and tetramer forms by comparative modelling and molecular docking techniques. Finally, molecular dynamics simulations at different pH conditions have indicated that at acidic pH the dimer assumes a widely open conformation, in agreement with the hypotheses coming from other studies.

The N-terminal 1-16 peptide derived in vivo from protein seminal vesicle protein IV modulates alpha-thrombin activity: potential clinical implications.

We have previously shown that seminal vesicle protein IV (SV-IV) and its 1-70 N-terminal fragment have anti-inflammatory activity and modulate anti-thrombin III (AT) activity. Moreover, mass spectrometry analysis of purified SV-IV has shown that the protein was found to be highly heterogeneous and 14% of the total SV-IV molecules are truncated forms, of particular interest the 1-16, 1-17, and 1-18 peptides. In this work we report experimental data which demonstrate that the 1-16 peptide (P1-16) possesses a marked effect on the AT activity by preventing the formation of the thrombin-AT complex. We found that the formation of thrombin-AT complex is markedly decreased in the presence of P1-16 used at equimolar concentration with thrombin as evaluated with SDS-PAGE. We also monitored the conformational changes of thrombin in the presence of different P1-16 concentrations, and calculated the K(d) of thrombin/P1-16 system by circular dichroism technique. The probable interaction sites of P1-16 with thrombin have been also evaluated by molecular graphics and computational analyses. These results have potential implications in the treatment of sterility and thrombotic diseases.

Evaluation of the structural quality of modeled proteins by using globularity criteria.

BACKGROUND: The knowledge of the three-dimensional structure of globular proteins is fundamental for a detailed investigation of their functional properties. Experimental methods are too slow for structure investigation on a large scale, while computational prediction methods offer alternatives that are continuously being improved. The international Comparative Assessment of Structure Prediction (CASP), an "a posteriori" evaluation of the quality of theoretical models when the experimental structure becomes available, demonstrates that predictions can be successful as well as unsuccessful, and this suggests the necessity for evaluations able to discard "a priori" the wrong models. RESULTS: We analyzed different structural properties of globular proteins for experimentally solved proteins belonging to the four different structural classes: "mainly alpha", "mainly beta", "alpha/beta" and "alpha+beta". The properties were found to be linearly correlated to protein molecular weight, but with some differences among the four classes. These results were applied to develop an evaluation test of theoretical models based on the expected globular properties of proteins. To verify the success of our test, we applied it to several protein models submitted to the sixth edition of CASP. The best theoretical models, as judged by CASP assessors, were in agreement with the expected properties, while most of the low-quality models had not passed our evaluations. CONCLUSION: This study supports the need for careful checks to avoid the diffusion of incorrect structural models. Our test allows the evaluation of models in the absence of experimental reference structures, thereby preventing the diffusion of incorrect structural models and the formulation of incorrect functional hypotheses. It can be used to check the globularity of predicted models, and to supplement other methods already used to evaluate their quality.

PreSSAPro: a software for the prediction of secondary structure by amino acid properties.

PreSSAPro is a software, available to the scientific community as a free web service designed to provide predictions of secondary structures starting from the amino acid sequence of a given protein. Predictions are based on our recently published work on the amino acid propensities for secondary structures in either large but not homogeneous protein data sets, as well as in smaller but homogeneous data sets corresponding to protein structural classes, i.e. all-alpha, all-beta, or alpha-beta proteins. Predictions result improved by the use of propensities evaluated for the right protein class. PreSSAPro predicts the secondary structure according to the right protein class, if known, or gives a multiple prediction with reference to the different structural classes. The comparison of these predictions represents a novel tool to evaluate what sequence regions can assume different secondary structures depending on the structural class assignment, in the perspective of identifying proteins able to fold in different conformations. The service is available at the URL http://bioinformatica.isa.cnr.it/PRESSAPRO/.

FASMA: a service to format and analyze sequences in multiple alignments.

Multiple sequence alignments are successfully applied in many studies for under- standing the structural and functional relations among single nucleic acids and protein sequences as well as whole families. Because of the rapid growth of sequence databases, multiple sequence alignments can often be very large and difficult to visualize and analyze. We offer a new service aimed to visualize and analyze the multiple alignments obtained with different external algorithms, with new features useful for the comparison of the aligned sequences as well as for the creation of a final image of the alignment. The service is named FASMA and is available at http://bioinformatica.isa.cnr.it/FASMA/.

The role of transglutaminase-2 and its substrates in human diseases.

The most characteristic enzymatic function of the class of enzymes known as transglutaminases (TG, EC 2.3.2.13) is the formation of covalent bonds between epsilon-amino groups of primary amines (from lysines or others) and the gamma-carboxamine group of glutamine residues of proteins. In the last years, a growing body of evidence indicate that the most interesting member of the TG family, namely the tissue TG (tTG, also called transglutaminase type 2, TG2), possesses more than one catalytic function. In fact, TG2 is able to catalyze a crosslinking reaction, a deamidation reaction and also shows GTP-binding/hydrolyzing and isopeptidase activities. Therefore, it can act on several classes of substrates, ranging from proteins to peptides, small reactive molecules like mono- and polyamines, and nucleotides. Given the broad spectrum of potentially different activities, elucidating the role of TG2 and its substrates in cellular functions and human diseases is a difficult task. In this study we focus our attention on substrates of TG2 and report a number of interesting considerations about their possible interplay in biological processes and involvement in human diseases, including genetic disorders. A significant improvement in understanding this complex scenario may come from a "multi-interfaced" approach, by exploiting different bioinformatic tools. Starting from a database of known TG2 substrates and using bioinformatic cross-search among other databases, we generated relational tables from which an involvement of TG2 in several genetic disorders can be hypothesized. Developing new bioinformatic tools and strategies to investigate the role of TG2 in molecular mechanisms underlying human diseases will add new light to this fascinating field of research.

Active Sequences Collection (ASC) database: a new tool to assign functions to protein sequences.

Active Sequences Collection (ASC) is a collection of amino acid sequences, with an unique feature: only short sequences are collected, with a demonstrated biological activity. The current version of ASC consists of three sections: DORRS, a collection of active RGD-containing peptides; TRANSIT, a collection of protein regions active as substrates of transglutaminase enzyme (TGase), and BAC, a collection of short peptides with demonstrated biological activity. Literature references for each entry are reported, as well as cross references to other databases, when available. The current version of ASC includes more than 800 different entries. The main scope of this collection is to offer a new tool to investigate the structural features of protein active sites, additionally to similarity searches against large protein databases or searching for known functional patterns. ASC database is available at the web address http://crisceb.unina2.it/ASC/ which also offers a dedicated query interface to compare user-defined protein sequences with the database, as well as an updating interface to allow contribution of new referenced active sequences.

PDGFR-alpha inhibits melanoma growth via CXCL10/IP-10: a multi-omics approach.

Melanoma is the most aggressive skin-cancer, showing high mortality at advanced stages. Platelet Derived Growth Factor Receptor-alpha (PDGFR-alpha) potently inhibits melanoma- and endothelium-proliferation and its expression is significantly reduced in melanoma-biopsies, suggesting that melanoma progression eliminates cells expressing PDGFR-alpha. In the present study transient overexpression of PDGFR-alpha in endothelial (HUVEC) and melanoma (SKMel-28, A375, Preyer) human-cells shows strong anti-proliferative effects, with profound transcriptome and miRNome deregulation. PDGFR-alpha overexpression strongly affects expression of 82 genes in HUVEC (41 up-, 41 down-regulated), and 52 genes in SKMel-28 (43 up-, 9 down-regulated). CXCL10/IP-10 transcript showed up to 20 fold-increase, with similar changes detectable at the protein level. miRNA expression profiling in cells overexpressing PDGFR-alpha identified 14 miRNAs up- and 40 down-regulated, with miR-503 being the most down-regulated (6.4 fold-reduction). miR-503, miR-630 and miR-424 deregulation was confirmed by qRT-PCR. Interestingly, the most upregulated transcript (i.e., CXCL10/IP-10) was a validated miR-503 target and CXCL10/IP-10 neutralization significantly reverted the anti-proliferative action of PDGFR-alpha, and PDGFR-alpha inhibition by Dasatinb totally reverted the CXCL10/IP10 induction, further supporting a functional interplay of these factors. Finally, integration of transcriptomics and miRNomics data highlighted several pathways affected by PDGFR-alpha.This study demonstrates for the first time that PDGFR-alpha strongly inhibits endothelial and melanoma cells proliferation in a CXCL10/IP-10 dependent way, via miR-503 down-regulation.

Homology modeling studies on human galactose-1-phosphate uridylyltransferase and on its galactosemia-related mutant Q188R provide an explanation of molecular effects of the mutation on homo- and heterodimers.

We have created theoretical models of the three-dimensional dimeric structure of human galactose-1-phosphate uridylyltransferase as well as of homo- and heterodimers carrying the Q188R mutation by using comparative modeling procedures. These mutants are associated to the most frequent form of the genetic disease galactosemia. We have analyzed the impact of this mutation both on enzyme-substrate interactions as well as on interchain interactions in the heterodimers and in the homodimer. We suggest a molecular explanation for the altered function, caused by different enzyme-substrate interactions, and for the partial dominant negative effect of the mutant allele that is present in heterozygotes for this gene, related to a substantial loss of interchain hydrogen bonds. These results can be considered a starting point for a more extensive characterization at the molecular level of the other mutations linked to this genetic disease.

Modelling of HLA-DQ2 and its interaction with gluten peptides to explain molecular recognition in celiac disease.

Celiac disease (CD) is sustained by abnormal intestinal mucosal T-cell response to gluten and it is strongly associated with HLA class II molecules encoded by DQA1*0501/DQB1*02 (DQ2) or DQA1*03/DQB1*0302 (DQ8). The in vitro stimulatory activity of gliadin increases after treatment with tissue transglutaminase (tTG) which catalyses the deamidation of specific residues of glutamine to glutamate that can serve as anchors for binding to DQ2 as well as to DQ8 molecules. We modelled the three-dimensional structure of the DQ2 dimer protein, the most frequent in celiac patients, by using a homology modelling strategy, and deposited the model in the Protein Data Bank (PDB). Then, we simulated the interactions of DQ2 with different gluten peptides and the deamidation of specific peptide glutamines in the known p4, p6, p7 and p9 anchor positions, as well as in p1 and p5 positions, and other substitutions for which experimental effects on binding are available by previous experimental studies. By evaluating the energy of interaction and the H-bond interactions, we were able to distinguish what substitutions improve the interaction peptide-DQ2, in agreement with previously published experimental data. By analysing the peptide-DQ2 complex at the atom level, we observed that these glutamate side chains can interact with specific positively charged amino acids of DQ2, absent in other HLA alleles not related to celiac disease. The simulation was also extended to other peptides, related to celiac disease but for which no experimental data exists about the effects of glutamine deamidation. Our results give an interpretation at the molecular level of previously reported binding experimental data and open a new window to gain further insights about peptide recognition in celiac disease.

Theoretical model of the three-dimensional structure of a sugar-binding protein from Pyrococcus horikoshii: structural analysis and sugar-binding simulations.

The three-dimensional structure of a sugar-binding protein from the thermophilic archaea Pyrococcus horikoshii has been predicted by a homology modelling procedure and investigated for its stability and its ability to bind different sugars. The model was created by using as templates the three-dimensional structures of a maltodextrin-binding protein from Pyrococcus furiosus, a trehalose-maltose-binding protein from Thermococcus litoralis and a maltodextrin-binding protein from Escherichia coli. According to the suggestions from the CASP (Critical Assessment of Structure Prediction) meetings, the homology modelling strategy was applied by assessing an accurate multiple sequence alignment, based on the high structural conservation in the family of ATP-binding cassette transporters to which all these proteins belong. The model has been deposited in the Protein Data Bank with the code 1R25. According to the origin of the protein, several characteristics in the organization of the secondary-structure elements and in the distribution of polar and non-polar amino acids are very similar to those of thermophilic proteins, compared with proteins from mesophilic organisms, and are analysed in detail. Finally, a simulation of the binding of several sugars in the binding site of this protein is presented, and interactions with amino acids are highlighted in detail.

Integrating mutation data and structural analysis of the TP53 tumor-suppressor protein.

TP53 encodes p53, which is a nuclear phosphoprotein with cancer-inhibiting properties. In response to DNA damage, p53 is activated and mediates a set of antiproliferative responses including cell-cycle arrest and apoptosis. Mutations in the TP53 gene are associated with more than 50% of human cancers, and 90% of these affect p53-DNA interactions, resulting in a partial or complete loss of transactivation functions. These mutations affect the structural integrity and/or p53-DNA interactions, leading to the partial or complete loss of the protein's function. We report here the results of a systematic automated analysis of the effects of p53 mutations on the structure of the core domain of the protein. We found that 304 of the 882 (34.4%) distinct mutations reported in the core domain can be explained in structural terms by their predicted effects on protein folding or on protein-DNA contacts. The proportion of "explained" mutations increased to 55.6% when substitutions of evolutionary conserved amino acids were included. The automated method of structural analysis developed here may be applied to other frequently mutated gene mutations such as dystrophin, BRCA1, and G6PD.

A thermostable sugar-binding protein from the Archaeon Pyrococcus horikoshii as a probe for the development of a stable fluorescence biosensor for diabetic patients.

In this work is presented the first attempt to develop an innovative ultrastable protein-based biosensor for blood glucose detections. The gene of a putative thermostable sugar-binding protein has been cloned and expressed in E. coli. The recombinant protein has been purified to homogeneity by thermoprecipitation and affinity chromatography steps. The recombinant protein is a monomer with an apparent molecular weight of 55,000 as judged by gel filtration and sodium dodecyl sulfate polyacrylamide gel eletrophoresis. Circular dichroism experiments showed that the protein possesses a secondary structure content rich in alpha-helices and beta-structures and that the protein is highly stable as investigated in the range of temperature between 20 and 95 degrees C. Fluorescence spectroscopy experiments demonstrated that the recombinant protein binds glucose with a dissociation constant of about 10 mM, a concentration of sugar very close to the concentration of glucose present in the human blood. A docking simulation on the modeled structure of the protein confirms its ability to bind glucose and proposes possible modifications to improve the affinity for glucose and/or its detection. The obtained results suggest the use of the protein as a probe for a stable glucose biosensor.