Multiple domains of the mouse p19ARF tumor suppressor are involved in p53-independent apoptosis

The ARF (p19ARF for the mouse ARF consisting of 169 amino acids and p14ARF for the human ARF consisting of 132 amino acids) genes upregulate p53 activities to induce cell cycle arrest and sensitize cells to apoptosis by inhibiting Mdm2 activity. p53-independent apoptosis also is induced by ectopic expression of p19ARF. We constructed various deletion mutants of p19ARF with a cre/loxP-regulated adenoviral vector to determine the regions of p19ARF which are responsible for p53-independent apoptosis. Ectopic expression of the C-terminal region (named C40) of p19ARF whose primary sequence is unique to the rodent ARF induced prominent apoptosis in p53-deficient mouse embryo fibroblasts. Relatively low-grade but significant apoptosis also was induced in p53-deficient mouse embryo fibroblasts by ectopic expression of p19ARF1-129, a p19ARF deletion mutant deficient in the C40 region. In contrast, ectopic expression of the wild-type p14ARF did not induce significant apoptosis in human cells. Taken together, we concluded that p53-independent apoptosis was mediated through multiple regions of the mouse ARF including C40, and the ability of the ARF gene to mediate p53-independent apoptosis has been not well conserved during mammalian evolution.     

The evolution of proteins from random amino acid sequences. I. Evidence from the lengthwise distribution of amino acids in modern protein sequences

Summary We examine in this paper one of the expected consequences of the hypothesis that modern proteins evolved from random heteropeptide sequences. Specifically, we investigate the lengthwise distributions of amino acids in a set of 1,789 protein sequences with little sequence identity using the run test statistic (r _o) of Mood (1940,Ann. Math. Stat. 11, 367–392). The probability density ofr _o for a collection of random sequences has mean=0 and variance=1 [the N(0,1) distribution] and can be used to measure the tendency of amino acids of a given type to cluster together in a sequence relative to that of a random sequence. We implement the run test using binary representations of protein sequences in which the amino acids of interest are assigned a value of 1 and all others a value of 0. We consider individual amino acids and sets of various combinations of them based upon hydrophobicity (4 sets), charge (3 sets), volume (4 sets), and secondary structure propensity (3 sets). We find that any sequence chosen randomly has a 90% or greater chance of having a lengthwise distribution of amino acids that is indistinguishable from the random expectation regardless of amino acid type. We regard this as strong support for the random-origin hypothesis. However, we do observe significant deviations from the random expectation as might be expected after billions years of evolution. Two important global trends are found: (1) Amino acids with a strong α-helix propensity show a strong tendency to cluster whereas those with β-sheet or reverse-turn propensity do not. (2) Clustered rather than evenly distributed patterns tend to be preferred by the individual amino acids and this is particularly so for methionine. Finally, we consider the problem of reconciling the random nature of protein sequences with structurally meaningful periodic “patterns” that can be detected by sliding-window, autocorrelation, and Fourier analyses. Two examples, rhodopsin and bacteriorhodopsin, show that such patterns are a natural feature of random sequences.     

Information-theoretic analysis of protein sequences shows that amino acids self-cluster

We analyse for each of 20 amino acids X the statistics of spacings between consecutive occurrences of X within the well-characterized Saccharomyces cerevisiae genome. The occurrences of amino acids may exhibit near random, clustered or smoothed out behaviour, like one-dimensional stochastic processes along the protein chain. If amino acids are distributed randomly within a sequence, then they follow a Poisson process, and a histogram of the number of observations of each gap size would asymptotically follow a negative exponential distribution. The novelty of the present approach lies in the use of differential geometric methods to quantify information on sequencing of amino acids and groups of amino acids, via the sequences of intervals between their occurrences. The differential geometry arises from an information-theoretic distance function on the two-dimensional space of stochastic processes subordinate to gamma distributions-which latter include the random process as a special case. We find that maximum-likelihood estimates of parametric statistics show that all 20 amino acids tend to cluster, some substantially. In other words, the frequencies of short gap lengths tend to be higher and the variance of the gap lengths is greater than expected by chance. This may be because localizing amino acids with the same properties may favour secondary structure formation or transmembrane domains. Gap sizes of 1 or 2 are generally disfavoured, 1 strongly so. The only exceptions to this are Gln and Ser, as a result of poly(Gln) or poly(Ser) sequences. There are preferences for gaps of 4 and 7 that can be attributed to alpha -helices. In particular, a favoured gap of 7 for Leu is found in coiled coils. Our method contributes to the characterization of whole sequences by extracting and quantifying stable stochastic features.     

Two distinct regions of the murine p53 primary amino acid sequence are implicated in stable complex formation with simian virus 40 T antigen.

We mapped regions of the mouse p53 primary amino acid sequence implicated in stable complex formation with simian virus 40 T antigen. A number of mutant p53 proteins failed to complex stably with T antigen in vivo but formed stable complexes with T antigen in in vitro association assays. In contrast to an earlier report (T.-H. Tan, H. Wallis, and A. J. Levine, J. Virol. 59:574-583, 1986), our study showed that two distinct regions of p53 primary amino acid sequence, highly conserved between mouse and Xenopus laevis, were implicated in stable complex formation. Our data support the proposal that, when in complex, T antigen may occupy a site on p53 that is implicated in the normal function of the protein.     

Structural comparisons between mouse and human prealbumin

In an attempt to construct model systems for familial amyloidotic polyneuropathy, prealbumin cDNA was cloned from a mouse liver cDNA library, using previously cloned human prealbumin cDNA as a hybridization probe. The primary structure of mouse prealbumin deduced from the cDNA sequence shows that it consists of 147 amino acids, including a whole prealbumin sequence (127 amino acids) and a putative signal sequence (20 amino acids). These numbers are in complete agreement with those determined for the human prealbumin. Among the 127 amino acid residues of the mature human prealbumin, 25 are replaced by different amino acids in the mouse prealbumin. Interestingly, 24 out of the 25 substituted amino acids are located at the outer surface of the protein, and the regions corresponding to the core and central channel of the protein are almost completely conserved. The cloned cDNA provided essential information for manipulating amyloidosis in mice.     

Purification and partial amino acid sequence analysis of the cellular tumour antigen, p53, from mouse SV40-transformed cells

The cellular tumour antigen p53 is implicated in the transformation process. To compare p53 from transformed cells and their normal counterparts in detail, and so to identify any structural differences that might alter p53 function, requires information on the primary structure of the protein. By making use of immunochemical techniques we have been able to purify nanomole amounts of p53. This was sufficient, using high sensitivity automated gas-phase sequencing techniques to determine the amino acid sequence of two tryptic peptides from p53. Their sequences agree completely with the predicted polypeptide sequence derived from a cloned cDNA for p53 mRNA and provide the first data on the amino acid sequence of p53. A combination of the high sensitivity amino acid sequencing procedures used here and cDNA sequence analysis should provide the complete amino acid sequence of p53.     

An accurate,sensitive, and scalable method to identify functional sites in protein structures

Functional sites determine the activity and interactions of proteins and as such constitute the targets of most drugs. However, the exponential growth of sequence and structure data far exceeds the ability of experimental techniques to identify their locations and key amino acids. To fill this gap we developed a computational Evolutionary Trace method that ranks the evolutionary importance of amino acids in protein sequences. Studies show that the best-ranked residues form fewer and larger structural clusters than expected by chance and overlap with functional sites, but until now the significance of this overlap has remained qualitative. Here, we use 86 diverse protein structures, including 20 determined by the structural genomics initiative, to show that this overlap is a recurrent and statistically significant feature. An automated ET correctly identifies seven of ten functional sites by the least favorable statistical measure, and nine of ten by the most favorable one. These results quantitatively demonstrate that a large fraction of functional sites in the proteome may be accurately identified from sequence and structure. This should help focus structure-function studies, rational drug design, protein engineering, and functional annotation to the relevant regions of a protein.     

Characterization and expression of the gene encoding membralin, an evolutionary conserved protein expressed in the central nervous system

We have identified a novel gene that encodes an evolutionary conserved protein that we name membralin since it contains multiple transmembrane regions. The human gene C19orf6 localizes to chromosome 19p13.3. Splice variant membralin-1 is encoded by 11 exons, translating into 620 amino acids. In addition, we found evidence for two additional splice variants in the human. The mouse gene ORF61 localizes to chromosome 10. We cloned two splice variants in mouse: membralin-1, which is encoded by 12 exons, translating into 574 amino acids, and membralin-2, which translates into 598 amino acids. The existence of rat membralin-1 (574 amino acids long) is, so far, only supported by in situ hybridization result, whereas the existence of rat membralin-2 (598 amino acids long) is strongly supported by overlapping ESTs. Gene homologues were also identified in fruit-fly (CG8405, chromosome 2R 52; two splice variants), nematode (chromosome III), and Arabidopsis thaliana (chromosome 1). Sequence analysis revealed no closely related genes, suggesting that membralin represents the sole member of a unique protein family.     

Activating mutations for transformation by p53 produce a gene product that forms an hsc70-p53 complex with an altered half-life.

The 11-4 p53 cDNA clone failed to transform primary rat fibroblasts when cotransfected with the ras oncogene. Two linker insertion mutations at amino acid 158 or 215 (of 390 amino acids) activated this p53 cDNA for transformation with ras. These mutant cDNAs produced a p53 protein that lacked an epitope, recognized by monoclonal antibody PAb246 (localized at amino acids 88 to 110 in the protein) and preferentially bound to a heat shock protein, hsc70. In rat cells transformed by a genomic p53 clone plus ras, two populations of p53 proteins were detected, PAb246+ and PAb246-, which did or did not bind to this monoclonal antibody, respectively. The PAb246- p53 preferentially associated with hsc70, and this protein had a half-life 4- to 20-fold longer than free p53 (PAb246+). These data suggest a possible functional role for hsc70 in the transformation process. cDNAs for p53 derived from methylcholanthrene-transformed cells transform rat cells in cooperation with the ras oncogene and produce a protein that bound with the heat shock proteins. Recombinant clones produced between a Meth A cDNA and 11-4 were tested for the ability to transform rat cells. A single amino acid substitution at residue 132 was sufficient to activate the 11-4 p53 cDNA for transformation. These studies have identified a region between amino acids 132 and 215 in the p53 protein which, when mutated, can activate the p53 cDNA. These results also call into question what the correct p53 wild-type sequence is and whether a wild-type p53 gene can transform cells in culture.     

Protein disorder prediction by condensed PSSM considering propensity for order or disorder

Background  

More and more disordered regions have been discovered in protein sequences, and many of them are found to be functionally significant. Previous studies reveal that disordered regions of a protein can be predicted by its primary structure, the amino acid sequence. One observation that has been widely accepted is that ordered regions usually have compositional bias toward hydrophobic amino acids, and disordered regions are toward charged amino acids. Recent studies further show that employing evolutionary information such as position specific scoring matrices (PSSMs) improves the prediction accuracy of protein disorder. As more and more machine learning techniques have been introduced to protein disorder detection, extracting more useful features with biological insights attracts more attention.     

Distribution of amino acid residues in the primary protein structure

The amino acid composition and sequence in primary structure of 180 proteins have been studied. It is shown that the distribution of amino acid residues is near to a random one, i.e. it is determined by the amino acid composition. The ratio between statistical and unique character of protein primary structures has been discussed. The amino acid sequence is suggested to be unique in fibrous proteins. In contrast the amino acid sequence in globular proteins is a statistical one. The statistical character of amino acids distribution in globular proteins explains the possibility of sensible text generation under the frame shift mutations, deletions and insertions.     

Simian virus 40 large T antigen contains two independent activities that cooperate with a ras oncogene to transform rat embryo fibroblasts.

The simian virus 40 large T antigen immortalizes growing primary cells in culture. In addition, this viral oncoprotein cooperates with an activated ras protein to produce dense foci on monolayers of rat embryo fibroblasts (REF). The relationship between independent immortalization and cooperative transformation with ras has not been defined. Previously, two regions of T antigen were shown to contain immortalization activities. An N-terminal fragment consisting of amino acids 1 to 147 immortalizes rodent cells (L. Sompayrac and K. J. Danna, Virology 181:412-415, 1991). Loss-of-function analysis indicated that immortalization depended on integrity of the T-antigen segments containing amino acids 351 to 450 and 533 to 626 (T. D. Kierstead and M. J. Tevethia, J. Virol. 67:1817-1829, 1993). The experiments described here were directed toward determining whether these same T-antigen regions were sufficient for cooperation with ras. Initially, constructs that produce T antigens containing amino acids 176 to 708 (T176-708) or 1 to 147 were tested in a ras cooperation assay. Both polypeptides cooperated with ras to produce dense foci on monolayers of primary REF. These results showed that T antigen contains two separate ras cooperation activities. In order to determine the N-terminal limit of the ras cooperation activity contained within the T176-708 polypeptide, a series of constructs designed to produce fusion proteins containing T-antigen segments beginning at residues 251, 301, 337, 351, 371, 401, 451, 501, 551, 601, and 651 was generated. Each of these constructs was tested for the capacity to cooperate with ras to produce dense foci on REF monolayers. The results indicated that a polypeptide containing T-antigen amino acids 251 to 708 (T251-708) was sufficient to cooperate with ras, whereas the more extensively truncated products were not. The abilities of the N-terminally truncated T antigens to bind p53 were examined in p53-deficient cells infected with a recombinant vaccinia virus expressing a phenotypically wild-type mouse p53. The results showed that polypeptides containing T-antigen amino acids 251 to 708, 301 to 708, 337 to 708, or 351 to 708 retained p53-binding capacity. The introduction into the T251-708 polypeptide of deletions that either prevented p53 binding (dl434-444) or did not prevent p53 binding (dl400) abrogated ras cooperation. These results indicated that although p53 binding may be necessary for ras cooperation, an additional, as-yet-undefined activity contained within the T251-708 polypeptide is needed.     

Maximum likelihood analysis of mammalian p53 indicates the presence of positively selected sites and higher tumorigenic mutations in purifying sites

The tumor suppressor gene TP53 (p53) maintains genome stability. Mutation or loss of p53 is found in most cancers. Analysis of evolutionary constrains and p53 mutations reveal important sites for concomitant functional studies. In this study, phylogenetic analyses of the coding sequences of p53 from 26 mammals were carried out by applying a maximum likelihood method. The results display two branches under adaptive evolution in mammals. Moreover, each codon of p53 was analyzed by the PAML method for presence of positively selected sites. PAML identified several statistically significant amino acids that undergo positive selection. The data indicates that amino acids responsible for the core functions of p53 are highly conserved, while positively selected sites are predominantly located in the N- and C-terminus of p53. Further analysis of evolutionary pressure and mutations showed the occurrence of more frequent tumorigenic mutations in purifying sites of p53.     

Molecular evolution of the thermosensitive PAb1620 epitope of human p53 by DNA shuffling.

Conformational stability of the p53 protein is an absolute necessity for its physiological function as a tumor suppressor. Recent in vitro studies have shown that wild-type p53 is a highly temperature-sensitive protein at the structural and functional levels. Upon heat treatment at 37 degrees C, p53 loses its wild-type (PAb1620(+)) conformation and its ability to bind DNA, but can be stabilized by different classes of ligands. To further investigate the thermal instability of p53, we isolated p53 mutants resistant to heat denaturation. For this purpose, we applied a recently developed random mutagenesis technique called DNA shuffling and screened for p53 variants that could retain reactivity to the native conformation-specific anti-p53 antibody PAb1620 upon thermal treatment. After three rounds of mutagenesis and screening, mutants were isolated with the desired phenotype. The isolated mutants were translated in vitro in either Escherichia coli or rabbit reticulocyte lysate and characterized biochemically. Mutational analysis identified 20 amino acid residues in the core domain of p53 (amino acids 101-120) responsible for the thermostable phenotype. Furthermore, the thermostable mutants could partially protect the PAb1620(+) conformation of tumor-derived p53 mutants from thermal unfolding, providing a novel approach for restoration of wild-type structure and possibly function to a subset of p53 mutants in tumor cells.     

What limits the primary sequence space of natural proteins?

Abstract

Number of naturally occurring primary sequences of proteins is an infinitesimally small subset of the possible number of primary sequences that can be synthesized using 20 amino acids. Prevailing views ascribe this to slow and incremental mutational/selection evolutionary mechanisms. However, considering the large number of avenues available in form of diversity of emerging/evolving and/or disappearing living systems for exploring the primary sequence space over the evolutionary time scale of ～3.5 billion years, this remains a conjecture. Therefore, to investigate primary sequence space limitations, we carried out a systematic study for finding primary sequences absent in nature. We report the discovery of the smallest peptide sequence “Cysteine-Glutamine-Tryptophan-Tryptophan” that is not found in over half-a-million curated protein sequences in the Uniprot (Swiss-Prot) database. Additionally, we report a library of 83605 pentapeptides that are not found in any of the known protein sequences. Compositional analyses of these absent primary sequences yield a remarkably strong power relationship between the percentage occurrence of individual amino acids in all known protein sequences and their respective frequency of occurrence in the absent peptides, regardless of their specific position in the sequences. If random evolutionary mechanisms were responsible for limitations to the primary sequence space, then one would not expect any relationship between compositions of available and absent primary sequences. Thus, we conclusively show that stoichiometric constraints on amino acids limit the primary sequence space of proteins in nature. We discuss the possibly profound implications of our findings in both evolutionary and synthetic biology.

Communicated by Ramaswamy H. Sarma     

Analysis of trk A and p53 association

trk A tyrosine kinase (the high affinity receptor for nerve growth factor) binds to the p53 tumour suppressor protein in vitro and in vivo. Our aim was to determine which regions of p53 are involved in trk A association. In vitro binding experiments using baculovirus expressed trk A and in vitro transcribed and translated C-terminus p53 deletion mutants show amino acids 327-338 critical for association. Also, analysis with mutants at the N-terminus, conserved regions II, III, IV and V or amino acid positions 173, 175, 181, 248 and 249 (which are amino acids frequently mutated in a variety of neoplasms and transformed cell lines), show that these sites are not involved in trk A binding. Importantly, similar results are obtained after immunoprecipitation of lysates from p53 negative fibroblasts expressing trk A and the above p53 mutant proteins. These data suggest that the amino-terminus of the oligomerisation domain of p53 is involved in p53/trk A association.     

Loss of oncogenic H-ras-induced cell cycle arrest and p38 mitogen-activated protein kinase activation by disruption of Gadd45a

The activation of p53 is a guardian mechanism to protect primary cells from malignant transformation; however, the details of the activation of p53 by oncogenic stress are still incomplete. In this report we show that in Gadd45a(-/-) mouse embryo fibroblasts (MEF), overexpression of H-ras activates extracellular signal-regulated kinase (ERK) and c-Jun N-terminal kinase (JNK) but not p38 kinase, and this correlates with the loss of H-ras-induced cell cycle arrest (premature senescence). Inhibition of p38 mitogen-activated protein kinase (MAPK) activation correlated with the deregulation of p53 activation, and both a p38 MAPK chemical inhibitor and the expression of a dominant-negative p38alpha inhibited p53 activation in the presence of H-ras in wild-type MEF. p38, but not ERK or JNK, was found in a complex with Gadd45 proteins. The region of interaction was mapped to amino acids 71 to 96, and the central portion (amino acids 71 to 124) of Gadd45a was required for p38 MAPK activation in the presence of H-ras. Our results indicate that this Gadd45/p38 pathway plays an important role in preventing oncogene-induced growth at least in part by regulating the p53 tumor suppressor.     

Intrinsically disordered regions of p53 family are highly diversified in evolution

Proteins of the p53 family are expressed in vertebrates and in some invertebrate species. The main function of these proteins is to control and regulate cell cycle in response to various cellular signals, and therefore to control the organism's development. The regulatory functions of the p53 family members originate mostly from their highly-conserved and well-structured DNA-binding domains. Many human diseases (including various types of cancer) are related to the missense mutations within this domain. The ordered DNA-binding domains of the p53 family members are surrounded by functionally important intrinsically disordered regions. In this study, substitution rates and propensities in different regions of p53 were analyzed. The analyses revealed that the ordered DNA-binding domain is conserved, whereas disordered regions are characterized by high sequence diversity. This diversity was reflected both in the number of substitutions and in the types of substitutions to which each amino acid was prone. These results support the existence of a positive correlation between protein intrinsic disorder and sequence divergence during the evolutionary process. This higher sequence divergence provides strong support for the existence of disordered regions in p53 in vivo for if they were structured, they would evolve at similar rates as the rest of the protein.