Prediction of Golgi Type II membrane proteins based on their transmembrane domains

Motivation: A major issue in cell biology today is how distinct intracellular regions of the cell, like the Golgi Apparatus, maintain their unique composition of proteins and lipids. The cell differentially separates Golgi resident proteins from proteins that move through the organelle to other subcellular destinations. We set out to determine if we could distinguish these two types of transmembrane proteins using computational approaches. Results: A new method has been developed to predict Golgi membrane proteins based on their transmembrane domains. To establish the prediction procedure, we took the hydrophobicity values and frequencies of different residues within the transmembrane domains into consideration. A simple linear discriminant function was developed with a small number of parameters derived from a dataset of Type II transmembrane proteins of known localization. This can discriminate between proteins destined for Golgi apparatus or other locations (post-Golgi) with a success rate of 89.3% or 85.2%, respectively on our redundancy-reduced data sets.

Secondary structure of the third extracellular loop responsible for ligand selectivity of a mammalian gonadotropin-releasing hormone receptor

The extracellular loop 3 (ECL3) of the mammalian gonadotropin-releasing hormone receptor (GnRH-R) contains an acidic amino acid (Glu(301) in the mouse GnRH-R,) that confers agonist selectivity for Are in mammalian GnRH. It is proposed that a specific conformation of ECL3 is necessary to orientate the carboxyl side chain of the acidic residue for interaction with Arg(8) of GnRH, which is supported by decreased affinity for Arg(8) GnRH but not Gln(8) GnRH when an adjacent Pro is mutated to Ala. To probe the structural contribution of the loop domain to the proposed presentation of the carboxyl side chain, we synthesized a model peptide (CGPEMLNRVSEPGC) representing residues 293-302 of mouse ECL3, where Cys and Gly residues are added symmetrically at the N and C termini, respectively, allowing the introduction of a disulfide bridge to simulate the distances at which the ECL3 is tethered to the transmembrane domains 6 and 7 of the receptor. The ability of the ECL3 peptide to bind GnRH with low affinity was demonstrated by its inhibition of GnRH stimulation of inositol phosphate production in cells expressing the GnRH-R. The CD bands of the ECL3 peptides exhibited a superposition of predominantly unordered structure and partial contributions from beta-sheet structure. Likewise, the analysis of the amide I and amide III bands from micro-Raman and FT Raman experiments revealed mainly unordered conformations of the cyclic and of the linear peptide. NMR data demonstrated the presence of a beta-hairpin among an ensemble of largely disordered structures in the cyclic peptide. The location of the turn linking the two strands of the hairpin was assigned to the three central residues L-296, N-297, and R-298. A small population of structured species among an ensemble of predominantly random coil conformation suggests that the unliganded receptor represents a variety of structural conformers, some of which have the potential to make contacts with the ligand. We propose a mechanism of receptor activation whereby binding of the agonist to the inactive receptor state induces and stabilizes a particular structural state of the loop domain, leading to further conformational rearrangements across the transmembrane domain and signal propagating interaction with G proteins. Interaction of the Glu(301) of the receptor with Arg(8) of GnRH induces a folded configuration of the ligand. Our proposal thus suggests that conformational changes of both ligand and receptor result from this interaction.

Macropodid herpesvirus 1 encodes genes for both thymidylate synthase and ICP34.5

Macropodid herpesvirus 1 (MaHV-1) is an unclassified alphaherpesvirus linked with the fatal infections of kangaroos and other marsupials. During the characterisation of the internal repeat region of MaHV-1, an open reading frame (ORF) encoding for thymidylate synthase (TS) gene was identified and completely sequenced. Southern blot analysis confirmed the presence of two copies of the TS gene in the MaHV-1 genome as expected. Computer analysis of the TS ORF showed it was 948 nucleotides in length. A putative polyadenylation signal was identified 17-22 bp inside the ORF implying a minimal or absent 3' untranslated region. The predicted polypeptide was 316 amino acid residues in length and contained the highly conserved motifs for folate binding and F-dUMP binding, typical of all TS enzymes. Interestingly, MaHV-1 TS polypeptide had highest similarity to the human TS polypeptide (81%) compared to the TS polypeptides of other herpesviruses (72-75%). Immediately upstream of the TS gene, a second ORF of 510 bp, encoding a polypeptide with 170 amino acid residues, was identified. The carboxyl domain of this MaHV-1 polypeptide shared 68% similarity to a 59 amino acid motif of human herpesvirus 1 ICP34.5, identifying it as the MaHV-1 ICP34.5 homologue. This is the first report of a herpesvirus that encodes for both TS and ICP34.5.

Isolation and characterisation of a novel rabbit sulfotransferase isoform belonging to the SULT1A subfamily

Sulfotransferases (SULTs) catalyse the sulfonation of both endogenous and exogenous compounds including hormones, catecholamines. drugs and xenobiotics. While in most occasions, sulfonation is a detoxication pathway. in the case of certain drugs and carcinogens. it leads to metabolic activation. Since, the rabbit has been extensively used for both pharmacological and toxicological studies, the purpose of this study was to further characterise the sulfotransferase system of this animal. In the present study, a novel sulfotransferase isoform (GenBank Accession no. AF360872) was isolated from a rabbit liver cDNA lambdaZAP 11 library. The full-length sequence of the clone was 1138 bp long and contained a coding region of 888 bp encoding a cytosolic protein of 295 amino acids (deduced molecular weight 34,193 Da). The amino acid sequence of this novel SULT isoform showed >70% identity with members of the SULT1A subfamily of sulfotransferases from other species. Upon expression of the encoded rabbit sulfotransferase in Escherchia coli (E. coli), it was shown that the enzyme was capable of sulfonating both p-nitrophenot (K-m and V-max values of 0.15 muM and 897.5 nmol/min/mg protein. respectively) and dopamine (K-m and V-max values of 175.3 muM and 151.1 nmol/min/mg protein, respectively). Based on the sequence data obtained and substrate specificity, this new rabbit sulfotransferase was named rabSULT1A1. Immunoblotting was used to demonstrate that rabSULT1A1 protein is expressed in liver, duodenum, jejunum, ileum, colon and recturm. (C) 2002 Elsevier Science Ltd. All rights reserved.

Association of the SULT1A1 R213H polymorphism with colorectal cancer

1. Sulphotransferases are a superfamily of enzymes involved in both detoxification and bioactivation of endogenous and exogenous compounds. The arylsulphotransferase SULT1A1 has been implicated in a decreased activity and thermostability when the wild-type arginine at position 213 of the coding sequence is substituted by a histidine. SULT1A1 is the isoform primarily associated with the conversion of dietary N -OH arylamines to DNA binding adducts and is therefore of interest to determine whether this polymorphism is linked to colorectal cancer. 2. Genotyping, using polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) analysis, was performed using DNA samples of healthy control subjects (n = 402) and patients with histologically proven colorectal cancer (n = 383). Both control and test populations possessed similar frequencies for the mutant allele (32.1 and 31%, respectively; P = 0.935). Results were not altered when age and gender were considered as potential confounders in a logistic regression analysis. 3. Examination of the sulphonating ability of the two allozymes with respect to the substrates p -nitrophenol and paracetamol showed that the affinity and rate of sulphonation was unaffected by substitution of arginine to histidine at position 213 of the amino acid sequence. 4. From this study, we conclude that the SULT1A1 R213H polymorphism is not linked with colorectal cancer in this elderly Australian population.

Molecular and functional analyses of Kunjin virus infectious cDNA clones demonstrate the essential roles for NS2A in virus assembly and for a nonconservative residue in NS3 in RNA replication

A number of full-length cDNA clones of Kunjin virus (KUN) were previously prepared; it was shown that two of them, pAKUN and FLSDX, differed in specific infectivities of corresponding in vitro transcribed RNAs by similar to100,000-fold (A. A. Khromykh et al., J. Virol. 72:7270-7279, 1998). In this study, we analyzed a possible genetic determinant(s) of the observed differences in infectivity initially by sequencing the entire cDNAs of both clones and comparing them with the published sequence of the parental KUN strain MRM61C. We found six common amino acid residues in both cDNA clones that were different from those in the published MRM61C sequence but were similar to those in the published sequences of other flaviviruses from the same subgroup. pAKUN clone had four additional codon changes, i.e., Ile59 to Asn and Arg175 to Lys in NS2A and Tyr518 to His and Ser557 to Pro in NS3. Three of these substitutions except the previously shown marker mutation, Arg175 to Lys in NS2A, reverted to the wild-type sequence in the virus eventually recovered from pAKUN RNA-transfected BHK cells, demonstrating the functional importance of these residues in viral replication and/or viral assembly. Exchange of corresponding DNA fragments between pAKUN and FLSDX clones and site-directed mutagenesis revealed that the Tyr518-to-His mutation in NS3 was responsible for an similar to5-fold decrease in specific infectivity of transcribed RNA, while the Ile59-to-Asn mutation in NS2A completely blocked virus production. Correction of the Asn59 in pAKUN NS2A to the wild-type lie residue resulted in complete restoration of RNA infectivity. Replication of KUN replicon RNA with an Ile59-to-Asn substitution in NS2A and with a Ser557-to-Pro substitution in NS3 was not affected, while the Tyr518-to-His substitution in NS3 led to severe inhibition of RNA replication. The impaired function of the mutated NS2A in production of infectious virus was complemented in trans by the helper wild-type NS2A produced from the KUN replicon RNA. However, replicon RNA with mutated NS2A could not be packaged in trans by the KUN structural proteins. The data demonstrated essential roles for the KUN nonstructural protein NS2A in virus assembly and for NS3 in RNA replication and identified specific single-amino-acid residues involved in these functions.

Identification of residues and domains of Raf important for function in vivo and in vitro

Random mutagenesis and genetic screens for impaired Raf function in Caenorhabditis elegans were used to identify six loss-of-function alleles of lin-45 raf that result in a substitution of a single amino acid. The mutations were classified as weak, intermediate, and strong based on phenotypic severity. We engineered these mutations into the homologous residues of vertebrate Raf-1 and analyzed the mutant proteins for their underlying biochemical defects. Surprisingly, phenotype strength did not correlate with the catalytic activity of the mutant proteins. Amino acid substitutions Val-589 and Ser-619 severely compromised Raf kinase activity, yet these mutants displayed weak phenotypes in the genetic screen. Interestingly, this is because these mutant Raf proteins efficiently activate the MAPK (mitogen-activated protein kinase) cascade in living cells, a result that may inform the analysis of knockout mice. Equally intriguing was the observation that mutant proteins with non-functional Ras-binding domains, and thereby deficient in Ras-mediated membrane recruitment, displayed only intermediate strength phenotypes. This confirms that secondary mechanisms exist to couple Ras to Raf in vivo. The strongest phenotype in the genetic screens was displayed by a S508N mutation that again did not correlate with a significant loss of kinase activity or membrane recruitment by oncogenic Ras in biochemical assays. Ser-508 lies within the Raf-1 activation loop, and mutation of this residue in Raf-1 and the equivalent Ser-615 in B-Raf revealed that this residue regulates Raf binding to MEK. Further characterization revealed that in response to activation by epidermal growth factor, the Raf-S508N mutant protein displayed both reduced catalytic activity and aberrant activation kinetics: characteristics that may explain the C. elegans phenotype.

Quantitation of alternatively spliced NMDA receptor NR1 isoform mRNA transcripts in human brain by competitive RT-PCR

The N-methyl-D-aspartate (NMDA)-selective subtype of ionotropic glutamate receptor is of importance in neuronal differentiation and synapse consolidation, activity-dependent forms of synaptic plasticity, and excitatory amino acid-mediated neuronal toxicity [Neurosci. Res. Program, Bull. 19 (1981) 1; Lab. Invest. 68 (1993) 372]. NMDA receptors exist in vivo as tetrameric or pentameric complexes comprising proteins from two families of homologous subunits, designated NR1 and NR2(A-D) [Biochem. Biophys. Res. Commun. 185 (1992) 826]. The gene coding for the human NR1 subunit (hNR1) is composed of 21 exons, three of which (4, 20 and 21) can be differentially spliced to generate a total of eight distinct subunit variants. We detail here a competitive RT-PCR (cRT-PCR) protocol to quantify endogenous levels of hNR1 splice variants in autopsied human brain. Quantitation of each hNR1 splice variant is performed using standard curve methodology in which a known amount of synthetic ribonucleic acid competitor (internal standard) is co-amplified against total RNA. This method can be used for the quantitation of hNR1 mRNA levels in response to acute or chronic disease states, in particular in the glutamatergic-associated neuronal loss observed in Alzheimer's disease [J. Neurochem. 78 (2001) 175]. Furthermore, alterations in hNR1 mRNA expression may be reflected at the translational level, resulting in functional changes in the NMDA receptor. (C) 2003 Elsevier Science B.V. All rights reserved.

A naturally selected dimorphism within the HLA-B44 supertype alters class I structure, peptide repertoire, and T cell recognition

HLA-B*4402 and B*4403 are naturally occurring MHC class I alleles that are both found at a hi,,h frequency in all human populations, and vet they only differ by one residue on the alpha2 helix (B*4402 Aspl56-->B*4403 Leu156) CTLs discriminate between HLA-B*4402 and B*4403, and these allotypes stimulate strong mutual allogeneic responses reflecting their known barrier to hemopoeitic stem cell transplantation. Although HLA-B*4402 and B*4403 share >95% of their peptide repertoire, B*4403 presents more unique peptides than B*4402, consistent with the stronger T cell alloreactivity observed toward B*4403 compared with B*4402. Crystal structures of B*4402 and B*4403 show how the polymorphisin at position 156 is completely buried and yet alters both the peptide and the heavy chain conformation, relaxing ligand selection by B*4403 compared with B*4402. Thus, the polymorphism between HLA-B*4402 and B 4403 modifies both peptide repertoire and T cell recognition, and is reflected lit the paradoxically powerful alloreactivity that occurs across this minimal mismatch. The findings suggest that these closely related class I genes are maintained lit diverse human populations through their differential impact on the selection of peptide ligands and the T cell repertoire.

eMelanoBase: An online locus-specific variant database for familial melanoma

A proportion of melanoma,prone individuals in both familial and non,familial contexts has been shown to carry inactivating mutations in either CDKN2A or, rarely, CDK4. CDKN2A is a complex locus that encodes two unrelated proteins from alternately spliced transcripts that are read in different frames. The alpha transcript (exons 1a, 2, and 3) produces the p16INK4A cyclin-dependent kinase inhibitor, while the beta transcript (exons 1beta and 2) is translated as p14ARF, a stabilizing factor of p53 levels through binding to MDM2. Mutations in exon 2 can impair both polypeptides and insertions and deletions in exons 1alpha, 1beta, and 2, which can theoretically generate p16INK4A,p14ARF fusion proteins. No online database currently takes into account all the consequences of these genotypes, a situation compounded by some problematic previous annotations of CDKN2A related sequences and descriptions of their mutations. As an initiative of the international Melanoma Genetics Consortium, we have therefore established a database of germline variants observed in all loci implicated in familial melanoma susceptibility. Such a comprehensive, publicly accessible database is an essential foundation for research on melanoma susceptibility and its clinical application. Our database serves two types of data as defined by HUGO. The core dataset includes the nucleotide variants on the genomic and transcript levels, amino acid variants, and citation. The ancillary dataset includes keyword description of events at the transcription and translation levels and epidemiological data. The application that handles users' queries was designed in the model,view. controller architecture and was implemented in Java. The object-relational database schema was deduced using functional dependency analysis. We hereby present our first functional prototype of eMelanoBase. The service is accessible via the URL www.wmi.usyd.e, du.au:8080/melanoma.html.

High frequency of BRAF mutations in nevi

To evaluate the timing of mutations in BRAF (v-raf murine sarcoma viral oncogene homolog B1) during melanocytic neoplasia, we carried out mutation analysis on microdissected melanoma and nevi samples. We observed mutations resulting in the V599E amino-acid substitution in 41 of 60 (68%) melanoma metastases, 4 of 5 (80%) primary melanomas and, unexpectedly, in 63 of 77 (82%) nevi. These data suggest that mutational activation of the RAS/RAF/MAPK pathway in nevi is a critical step in the initiation of melanocytic neoplasia but alone is insufficient for melanoma tumorigenesis.

Linearization of a naturally occurring circular protein maintains structure but eliminates hemolytic activity

Cyclotides are a recently discovered family of disulfide rich proteins from plants that contain a circular protein backbone. They are exceptionally stable, as exemplified by their use in native medicine of the prototypic cyclotide kalata B1. The peptide retains uterotonic activity after the plant from which it is derived is boiled to make a medicinal tea. The circular backbone is thought to be in part responsible for the stability of the cyclotides, and to investigate its role in determining structure and biological activity, an acyclic derivative, des-(24-28)-kalata B1, was chemically synthesized and purified. This derivative has five residues removed from the 29-amino acid circular backbone of kalata B1 in a loop region corresponding to a processing site in the biosynthetic precursor protein. Two-dimensional NMR spectra of the peptide were recorded, assigned, and used to identify a series of distance, angle, and hydrogen bonding restraints. These were in turn used to determine a representative family of solution structures. Of particular interest was a determination of the structural similarities and differences between des-(2428)-kalata B1 and native kalata B1. Although the overall three-dimensional fold remains very similar to that of the native circular protein, removal of residues 24-28 of kalata B1 causes disruption of some structural features that are important to the overall stability. Furthermore, loss of hemolytic activity is associated with backbone truncation and linearization.

Computational differentiation of N-terminal signal peptides and transmembrane helices

Signal peptides and transmembrane helices both contain a stretch of hydrophobic amino acids. This common feature makes it difficult for signal peptide and transmembrane helix predictors to correctly assign identity to stretches of hydrophobic residues near the N-terminal methionine of a protein sequence. The inability to reliably distinguish between N-terminal transmembrane helix and signal peptide is an error with serious consequences for the prediction of protein secretory status or transmembrane topology. In this study, we report a new method for differentiating protein N-terminal signal peptides and transmembrane helices. Based on the sequence features extracted from hydrophobic regions (amino acid frequency, hydrophobicity, and the start position), we set up discriminant functions and examined them on non-redundant datasets with jackknife tests. This method can incorporate other signal peptide prediction methods and achieve higher prediction accuracy. For Gram-negative bacterial proteins, 95.7% of N-terminal signal peptides and transmembrane helices can be correctly predicted (coefficient 0.90). Given a sensitivity of 90%, transmembrane helices can be identified from signal peptides with a precision of 99% (coefficient 0.92). For eukaryotic proteins, 94.2% of N-terminal signal peptides and transmembrane helices can be correctly predicted with coefficient 0.83. Given a sensitivity of 90%, transmembrane helices can be identified from signal peptides with a precision of 87% (coefficient 0.85). The method can be used to complement current transmembrane protein prediction and signal peptide prediction methods to improve their prediction accuracies. (C) 2003 Elsevier Inc. All rights reserved.

Screening of human primary melanocytes of defined melanocortin-1 receptor genotype: Pigmentation marker, ultrastructural and UV-survival studies

Recent population studies have demonstrated an association with the red-hair and fair-skin phenotype with variant alleles of the melanocortin-1 receptor (MC1R) which result in amino acid substitutions within the coding region leading to an altered receptor activity. In particular, Arg151Cys, Arg160Trp and Asp294His were the most commonly associated variants seen in the south-east Queensland population with at least one of these alleles found in 93% of those with red hair. In order to study the individual effects of these variants on melanocyte biology and melanocytic pigmentation, we established a series of human melanocyte strains genotyped for the MC1R receptor which included wild-type consensus, variant heterozygotes, compound heterozygotes and homozygotes for Arg151Cys, Arg160Trp, Val60Leu and Val92Met alleles. These strains ranged from darkly pigmented to amelanotic, with all strains of consensus sequence having dark pigmentation. UV sensitivity was found not to be associated with either MC1R genotype or the level of pigmentation with a range of sensitivities seen across all genotypes. Ultrastructural analysis demonstrated that while consensus strains contained stage IV melanosomes in their terminal dendrites, Arg151Cys and Arg160Trp homozygote strains contained only stage II melanosomes. This was despite being able to show expression of tyrosinase and tyrosinase-related protein-1 markers, although at reduced levels and an ability to convert exogenous 3,4-dihydroxyphenyl-alanine (DOPA) to melanin in these strains.