Better prediction of protein contact number using a support vector regression analysis of amino acid sequence

Background: Protein tertiary structure can be partly characterized via each amino acid's contact number measuring how residues are spatially arranged. The contact number of a residue in a folded protein is a measure of its exposure to the local environment, and is defined as the number of C-beta atoms in other residues within a sphere around the C-beta atom of the residue of interest. Contact number is partly conserved between protein folds and thus is useful for protein fold and structure prediction. In turn, each residue's contact number can be partially predicted from primary amino acid sequence, assisting tertiary fold analysis from sequence data. In this study, we provide a more accurate contact number prediction method from protein primary sequence. Results: We predict contact number from protein sequence using a novel support vector regression algorithm. Using protein local sequences with multiple sequence alignments (PSI-BLAST profiles), we demonstrate a correlation coefficient between predicted and observed contact numbers of 0.70, which outperforms previously achieved accuracies. Including additional information about sequence weight and amino acid composition further improves prediction accuracies significantly with the correlation coefficient reaching 0.73. If residues are classified as being either contacted or non-contacted, the prediction accuracies are all greater than 77%, regardless of the choice of classification thresholds. Conclusion: The successful application of support vector regression to the prediction of protein contact number reported here, together with previous applications of this approach to the prediction of protein accessible surface area and B-factor profile, suggests that a support vector regression approach may be very useful for determining the structure-function relation between primary sequence and higher order consecutive protein structural and functional properties.

STAR: predicting recombination sites from amino acid sequence

Background: Designing novel proteins with site-directed recombination has enormous prospects. By locating effective recombination sites for swapping sequence parts, the probability that hybrid sequences have the desired properties is increased dramatically. The prohibitive requirements for applying current tools led us to investigate machine learning to assist in finding useful recombination sites from amino acid sequence alone. Results: We present STAR, Site Targeted Amino acid Recombination predictor, which produces a score indicating the structural disruption caused by recombination, for each position in an amino acid sequence. Example predictions contrasted with those of alternative tools, illustrate STAR'S utility to assist in determining useful recombination sites. Overall, the correlation coefficient between the output of the experimentally validated protein design algorithm SCHEMA and the prediction of STAR is very high (0.89). Conclusion: STAR allows the user to explore useful recombination sites in amino acid sequences with unknown structure and unknown evolutionary origin. The predictor service is available from http://pprowler.itee.uq.edu.au/star.

Phylogenetic analysis of evolutionary relationships of the planctomycete division of the domain bacteria based on amino acid sequences of elongation factor Tu

Sequences from the tuf gene coding for the elongation factor EF-Tu were amplified and sequenced from the genomic DNA of Pirellula marina and Isosphaera pallida, two species of bacteria within the order Planctomycetales. A near-complete (1140-bp) sequence was obtained from Pi. marina and a partial (759-bp) sequence was obtained for I. pallida. Alignment of the deduced Pi. marina EF-Tu amino acid sequence against reference sequences demonstrated the presence of a unique Il-amino acid sequence motif not present in any other division of the domain Bacteria. Pi. marina shared the highest percentage amino acid sequence identity with I. pallida but showed only a low percentage identity with other members of the domain Bacteria. This is consistent with the concept of the planctomycetes as a unique division of the Bacteria. Neither primary sequence comparison of EF-Tu nor phylogenetic analysis supports any close relationship between planctomycetes and the chlamydiae, which has previously been postulated on the basis of 16S rRNA. Phylogenetic analysis of aligned EF-Tu amino acid sequences performed using distance, maximum-parsimony, and maximum likelihood approaches yielded contradictory results with respect to the position of planctomycetes relative to other bacteria, It is hypothesized that long-branch attraction effects due to unequal evolutionary rates and mutational saturation effects may account for some of the contradictions.

D-amino acid residue in a defensin-like peptide from platypus venom: effect on structure and chromatographic properties

The recent discovery that the natriuretic peptide OvCNPb (Ornithorhynchus venom C-type natriuretic peptide B) from platypus (Ornithorynchus anatinus) venom contains a D-amino acid residue suggested that other D-amino-acid-containing peptides might be present in the venom. In the present study, we show that DLP-2 (defensin-like peptide-2), a 42-amino-acid residue polypeptide in the platypus venom, also contains a D-amino acid residue, D-methionine, at position 2, while DLP-4, which has an identical amino acid sequence, has all amino acids in the L-form. These findings were supported further by the detection of isomerase activity in the platypus gland venom extract that converts DLP-4 into DLP-2. In the light of this new information, the tertiary structure of DLP-2 was recalculated using a new structural template with D-Met(2). The structure of DLP-4 was also determined in order to evaluate the effect of a D-amino acid at position 2 on the structure and possibly to explain the large retention time difference observed for the two molecules in reverse-phase HPLC. The solution structures of the DLP-2 and DLP-4 are very similar to each other and to the earlier reported structure of DLP-2, which assumed that all amino acids were in the L-form. Our results suggest that the incorporation of the D-amino acid at position 2 has minimal effect on the overall fold in solution.

Deduced amino acid sequences surrounding the fusion glycoprotein cleavage site and of the carboxyl-terminus of haemagglutinin-neuraminidase protein of the avirulent thermostable vaccine strain I-2 of Newcastle disease virus

A single-tube RT-PCR technique generated a 387 bp or 300 bp cDNA amplicon covering the F-0 cleavage site or the carboxyl (C)-terminus of the HN gene, respectively, of Newcastle disease virus (NDV) strain 1-2. Sequence analysis was used to deduce the amino acid sequences of the cleavage site of F protein and the C-terminus of HN protein, which were then compared with sequences for other NDV strains. The cleavage site of NDV strain 1-2 had a sequence Motif of (112)RKQGRLIG(119), consistent with an avirulent phenotype. Nucleotide sequencing and deduction of amino acids at the C-terminus of HN revealed that strain 1-2 had a 7-amino-acid extension (VEILKDGVREARSSR). This differs from the virulent viruses that caused outbreaks of Newcastle disease in Australia in the 1930s and 1990s, which have HN extensions of 0 and 9 amino acids, respectively. Amino acid sequence analyses of the F and HN genes of strain 1-2 confirmed its avirulent nature and its Australian origin.

Sequence and tissue distribution of chicken cellular nucleic acid binding protein cDNA

We sequenced cDNAs coding for chicken cellular nucleic acid binding protein (CNBP). Two slightly different variations of the open reading frame were found, each of which translates into a protein with seven zinc finger domains. The longest transcript contains an in-frame insert of 3 bp. The sequence conservation between chick CNBP cDNAs with human, rat and mouse CNBP cDNAs is extreme, especially in the coding region, where the deduced amino acid sequence identity with human, rat and mouse CNBP is 99%. CNBP-like transcripts were also found in various tissues from insect, shrimp, fish and lizard. Regions with remarkable nucleotide conservation were also found in the 3' untranslated region, indicating important functions for these regions. Quantitative reverse transcription polymerase chain reaction (RT-PCR) indicated that in the chick, CNBP is present in all tissues examined in approximately equal ratios to total RNA. RT-PCR of total RNA isolated from different phyla indicate CNBP-like proteins art widespread throughout the animal kingdom. The extraordinary level of conservation suggests an important physiological role for CNBP. (C) 1997 Elsevier Science Inc.

Encephalitogenicity of murine, but not bovine, DM20 in SJL mice is due to a single amino acid difference in the immunodominant encephalitogenic epitope

Myelin proteolipid protein (PLP) contains 2 immunodominant encephalitogenic epitopes in SJL mice, namely PLP residues 139-151 and 178-191. DM20, a minor isoform of PLP, lacks residues 116-150 and consequently contains only the single major encephalitogenic epitope 178-191. However, it has been found previously that bovine DM20 is not encephalitogenic in SJL mice. Since residue 188 within peptide 178-191 is phenylalanine (F) in murine DM20 and alanine (A) in bovine DM20, we tested the effect of this difference on the immune responses and induction of EAE. SJL mice were immunized with either highly purified murine or bovine DM20. Residues 178-191 were found to be immunodominant for each, but only murine and not bovine DM20 was encephalitogenic. A synthetic peptide corresponding to the murine 178-191 sequence (F188) was also encephalitogenic, whereas the peptide corresponding to the bovine sequence (A188) was not. Both F188 and A188 bind with high affinity to I-A(s) and both are recognized by the SJL T cell repertoire. A188-specific T cell lines reacted to both A188 and F188, but F188-specific T cell lines were not stimulated by A188. F188-specific T cell lines produced mRNA for the Th1 cytokines IL2 and IFN gamma and, in passive transfer experiments, were encephalitogenic upon stimulation with F188, but not A188. In contrast, A188-specific T cell lines produced mRNA for IL4, IL5 and IL10, in addition to IL2 and IFN gamma, and were not encephalitogenic after stimulation with either F188 or A188. Cotransfer of A188-specific T cell lines with F188-specific T cell lines resulted in protection from EAE. Thus, A188 induces a functionally different phenotype of T cells from that induced by F188. Taken together these data suggest that the failure of bovine DM20 to induce EAE may be attributable to induction of protective rather than pathogenic T cells by the immunodominant epitope.

Sulfadoxine resistance in Plasmodium vivax is associated with a specific amino acid in dihydropteroate synthase at the putative sulfadoxine-binding site

Sulfadoxine is predominantly used in combination with pyrimethamine, commonly known as Fansidar, for the treatment of Plasmodium falciparum. This combination is usually less effective against Plasmodium vivax, probably due to the innate refractoriness of parasites to the sulfadoxine component. To investigate this mechanism of resistance by P. vivax to sulfadoxine, we cloned and sequenced the P. vivax dhps (pvdhps) gene. The protein sequence was determined, and three-dimensional homology models of dihydropteroate synthase (DHPS) from P. vivax as well as P. falciparum were created. The docking of sulfadoxine to the two DHPS models allowed us to compare contact residues in the putative sulfadoxine-binding site in both species. The predicted sulfadoxine-binding sites between the species differ by one residue, V585 in P. vivax, equivalent to A613 in P. falciparum. V585 in P. vivax is predicted by energy minimization to cause a reduction in binding of sulfadoxine to DHPS in P. vivax compared to P. falciparum. Sequencing dhps genes from a limited set of geographically different P. vivax isolates revealed that V585 was present in all of the samples, suggesting that V585 may be responsible for innate resistance of P. vivax to sulfadoxine. Additionally, amino acid mutations were observed in some P. vivax isolates in positions known to cause resistance in P. falciparum, suggesting that, as in P. falciparum, these mutations are responsible for acquired increases in resistance of P. vivax to sulfadoxine.

Isolation and characterisation of conomap-Vt, a D-amino acid containing excitatory peptide from the venom of a vermivorous cone snail

Cone snail venom is a rich source of bioactives, in particular small disulfide rich peptides that disrupt synaptic transmission. Here, we report the discovery of conomap-Vt (Conp-Vt), an unusual linear tetradecapeptide isolated from Conus vitulinus venom. The sequence displays no homology to known conopeptides, but displays significant homology to peptides of the MATP (myoactive tetradecapeptide) family, which are important endogenous neuromodulators in molluscs, annelids and insects. Conp-Vt showed potent excitatory activity in several snail isolated tissue preparations. Similar to ACh, repeated doses of Conp-Vt were tachyphylactic. Since nicotinic and muscarinic antagonists failed to block its effect and Conp-Vt desensitised tissue remained responsive to ACh, it appears that Conp-Vt contractions were non-cholinergic in origin. Finally, biochemical studies revealed that Conp-Vt is the first member of the MATP family with a D-amino acid. Interestingly, the isomerization of L-Phe to D-Phe enhanced biological activity, suggesting that this post-translational modified conopeptide may have evolved for prey capture. (c) 2006 Published by Elsevier B.V. on behalf of the Federation of European Biochemical Societies.

Molecular cloning and sequence of the ovine gastrin gene

A clone encoding ovine preprogastrin was isolated from a sheep genomic library. The deduced 104 amino acid sequence of ovine preprogastrin was 92% and 68% identical to the sequences of bovine and human preprogastrin, respectively. While the similarity was greatest in the gastrin-17 sequence, an unexpected similarity was also observed in the N-terminus of mature progastrin.

cis-3-hydroxy-N-methyl-L-proline: a new amino acid from a southern Australian marine sponge, Dendrilla sp.

A specimen of the sponge Dendrilla sp. collected during commercial trawling operations in the Great Australian Eight, Australia, analyses for a very high natural abundance of the new amino acid cis-3-hydroxy-N-methyl-L-proline (1). The complete stereostructure for (1) was determined by spectroscopic analysis and chemical derivatization.

Amino acid synthesis in the symbiotic sea anemone Aiptasia pulchella

Symbiotic Aiptasia pulchella and freshly isolated zooxanthellae were incubated in (NaHCO3)-C-14 and NH4Cl for 1 to 240 min, and samples were analysed by reverse-phase high-performance liquid chromatography (HPLC) and an online radiochemical detector. NH4+ was first assimilated into C-14-glutamate and C-14-glutamine in the zooxanthellae residing in A. pulchella. The specific activities (dpm nmol(-1)) of C-14-glutamate and C-14-glutamine in vivo, were far greater in the zooxanthellae than in the host tissue, indicating that NH4+ was principally incorporated into the glutamate and glutamine pools of the zooxanthellae. C-14-alpha-ketoglutarate was taken up from the medium by intact A. pulchella and assimilated into a small amount of C-14-glutamate in the host tissue, but no C-14-glutamine was detected in the host fraction. The C-14-glutamate that was synthesized was most likely produced from transamination reactions as opposed to the direct assimilation of NH4+. The free aminoacid composition of the host tissue and zooxanthellae of A. pulchella was also measured. The results presented here demonstrate that NH4+ was initially assimilated by the zooxanthellae of A. pulchella.

Protein sequence threading, the alignment problem, and a two-step strategy

Conventionally, protein structure prediction via threading relies on some nonoptimal method to align a protein sequence to each member of a library of known structures. We show how a score function (force field) can be modified so as to allow the direct application of a dynamic programming algorithm to the problem. This involves an approximation whose damage can be minimized by an optimization process during score function parameter determination. The method is compared to sequence to structure alignments using a more conventional pair-wise score function and the frozen approximation. The new method produces results comparable to the frozen approximation, but is faster and has fewer adjustable parameters. It is also free of memory of the template's original amino acid sequence, and does not suffer from a problem of nonconvergence, which can be shown to occur with the frozen approximation. Alignments generated by the simplified score function can then be ranked using a second score function with the approximations removed. (C) 1999 John Wiley & Sons, Inc.

Single amino acid substitutions in alpha-conotoxin PnIA shift selectivity for subtypes of the mammalian neuronal nicotinic acetylcholine receptor

The alpha-conotoxins, a class of nicotinic acetylcholine receptor (nAChR) antagonists, are emerging as important probes of the role played by different nAChR subtypes in cell function and communication, In this study, the native alpha-conotoxins PnIA and PnIB were found to cause concentration-dependent inhibition of the ACh-induced current in all rat parasympathetic neurons examined, with IC50 values of 14 and 33 nM, and a maximal reduction in current amplitude of 87% and 71%, respectively. The modified alpha-conotoxin [N11S]PnIA reduced the ACh-induced current with an IC50 value of 375 nM and a maximally effective concentration caused 91% block, [A10L]PnIA was the most potent inhibitor, reducing the ACh-induced current in similar to 80% of neurons, with an IC50 value of 1.4 nM and 46% maximal block of the total current, The residual current was not inhibited further by alpha-bungarotoxin, but was further reduced by the cu-conotoxins PnIA or PnIB, and by mecamylamine. H-1 NMR studies indicate that PnIA, PnIB, and the analogues, [A10L]PnIA and [N11S]PnIA, have identical backbone structures. We propose that positions 10 and II of PnIA and PnIB influence potency and determine selectivity among alpha 7 and other nAChR subtypes, including alpha 3 beta 2 and alpha 3 beta 4, Four distinct components of the nicotinic ACh-induced current in mammalian parasympathetic neurons have been dissected with these conopeptides.