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.

Amino acid sequence of the winged bean (Psophocarpus tetragonolobus) basic lectin. Adenine binding and identification of the active-site tryptophan residue

The complete amino acid sequence of winged bean basic agglutinin (WBA I) was obtained by a combination of manual and gas-phase sequencing methods. Peptide fragments for sequence analyses were obtained by enzymatic cleavages using trypsin and Staphylococcus aureus V8 endoproteinase and by chemical cleavages using iodosobenzoic acid, hydroxylamine, and formic acid. COOH-terminal sequence analysis of WBA I and other peptides was performed using carboxypeptidase Y. The primary structure of WBA I was homologous to those of other legume lectins and more so to Erythrina corallodendron. Interestingly, the sequence shows remarkable identities in the regions involved in the association of the two monomers of E. corallodendron lectin. Other conserved regions are the double metal-binding site and residues contributing to the formation of the hydrophobic cavity and the carbohydrate-binding site. Chemical modification studies both in the presence and absence of N-acetylgalactosamine together with sequence analyses of tryptophan-containing tryptic peptides demonstrate that tryptophan 133 is involved in the binding of carbohydrate ligands by the lectin. The location of tryptophan 133 at the active center of WBA I for the first time subserves to explain a role for one of the most conserved residues in legume lectins.

Analysis of the effects of amino-acid-sequence on the structure of proteins

The conformation of amino acid side chains as observed in well-determined structures of globular proteins has earlier been extensively investigated. In contrast, the structural features of the polypeptide backbone that result from the occurrence of specific amino acids along the polypeptide have not been analysed. In this article, we present the statistically significant features in the backbone geometry that appear to be a consequence of the occurrence of rotamers of different amino acid side chains by analysing 102 well-refined structures that form a random collection of proteins. It is found that the persistence of helical segments around each residue is influenced by the residue type. Several residues exert asymmetrical influence between the carboxyl and amino terminal polypeptide segments. The degree to which secondary structures depart from an average geometry also appears to depend on residue type. These departures are correlated to the corresponding Chou and Fasman parameters of amino acid residues. The frequency distribution of the side chain rotamers is influenced by polypeptide secondary structure. In turn, the rotamer conformation of side chain affects the extension of the secondary structure of the backbone. The strongest correlation is found between the occurrence of g+ conformation and helix propagation on the carboxyl side of many residues.

Comparative Nucleotide and Amino Acid Sequence Analysis of the Sequence-Specific RNA-Binding Rotavirus Nonstructural Protein NSP3

NSP3, an acidic nonstructural protein, encoded by gene 7 has been implicated as the key player in the assembly of the 11 viral plus-strand RNAs into the early replication intermediates during rotavirus morphogenesis. To date, the sequence or NSP3 from only three animal rotaviruses (SA11, SA114F, and bovine UK) has been determined and that from a human strain has not been reported. To determine the genetic diversity among gene 7 alleles from group A rotaviruses, the nucleotide sequence of the NSP3 gene from 13 strains belonging to nine different G serotypes, from both humans and animals, has been determined. Based on the amino acid sequence identity as well as phylogenetic analysis, NSP3 from group A rotaviruses falls into three evolutionarily related groups, i.e., the SA11 group, the Wa group, and the S2 group. The SA 11/SA114F gene appears to have a distant ancestral origin from that of the others and codes for a polypeptide of 315 amino acids (aa) in length. NSP3 from all other group A rotaviruses is only 313 aa in length because of a 2-amino-acid deletion near the carboxy-terminus, While the SA114F gene has the longest 3' untranslated region (UTR) of 132 nucleotides, that from other strains suffered deletions of varying lengths at two positions downstream of the translational termination codon. In spite of the divergence of the nucleotide (nt) sequence in the protein coding region, a stretch of about 80 nt in the 3' UTR is highly conserved in the NSP3 gene from all the strains. This conserved sequence in the 3' UTR might play an important role in the regulation of expression of the NSP3 gene. (C) 1995 Academic Press, Inc.

The complete amino acid sequence of growth hormone and partial amino acid sequence of prolactin and somatolactin from sea perch (Lateolabrax japonicus)

Growth hormone (GH), prolactin (PRL) and somatolactin (SL) were purified simultaneously under alkaline condition (pH 9.0) from pituitary glands of sea perch (Lateolabrax japonicas) by a two-step procedure involving gel filtration on Sephadex G-100 and reverse-phase high-performance liquid chromatography (rpHPLC). At each step of purification, fractions were monitored by sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE) and by immunoblotting with chum salmon GH. PRL and SL antisera. The yields of sea perch GH, PRL and SL were 4.2, 1.0 and 0.28 mg/g wet tissue, respectively. The molecular weights of 19,200 and 20,370 Da were estimated by SDS-PAGE for sea perch GH and PRL, respectively. Two forms of sea perch SL were found: one (28,400 Da) is probably glycosylated, while the other one (23,200 Da) is believed to be deglycosylated. GH bioactivity was examined by an in vivo assay. Intraperitoneal injection of sea perch GH at a dose of 0.01 and 0.1 mug/g body weight at 7-day intervals resulted in a significant increase in body weight and length of juvenile rainbow trout. The complete sea-perch GH amino acid sequence of 187 residues was determined by sequencing fragments cleaved by chemicals and enzymes. Alignment of sea-perch GH with those of other fish GHs revealed that sea-perch GH is most similar to advanced marine fish, such as tuna, gilthead sea bream, yellowfin porgy, red sea bream, bonito and yellow tail with 98.4, 96.2%, 95.7%, 95.2%, 94.1% and 91% sequence identity, respectively. Sea-perch GH has low identity to Atlantic cod (76.5%), hardtail (73.3%), flounder (68.4%), chum salmon (66.3%), carp (54%) and blue shark (38%). Partial amino-acid sequences of 127 of sea-perch PRL and the N-terminal of 16 amino-acid sequence of sea-perch SL have been determined. The data show that sea-perch PRL has a slightly higher sequence identity with tilapia PRL( 73.2%) than with chum salmon PRL(70%) in this 127 amino-acid sequence. (C) 2001 Elsevier Science B.V. All rights reserved.

The natriuretic peptide/helokinestatin precursor from Mexican beaded lizard (Heloderma horridum) venom: Amino acid sequence deduced from cloned cDNA and identification of two novel encoded helokinestatins

Natriuretic peptides are common components of reptile venoms and molecular cloning of their biosynthetic precursors has revealed that in snakes, they co-encode bradykinin-potentiating peptides and in venomous lizards, some co-encode bradykinin inhibitory peptides such as the helokinestatins. The common natriuretic peptide/helokinestatin precursor of the Gila Monster, Heloderma suspectum, encodes five helokinestatins of differing primary structures. Here we report the molecular cloning of a natriuretic peptide/helokinestatin precursor cDNA from a venom-derived cDNA library of the Mexican beaded lizard (Heloderma horridum). Deduction of the primary structure of the encoded precursor protein from this cloned cDNA template revealed that it consisted of 196 amino acid residues encoding a single natriuretic peptide and five helokinestatins. While the natriuretic peptide was of identical primary structure to its Gila Monster (H. suspectum) homolog, the encoded helokinestatins were not, with this region of the common precursor displaying some significant differences to its H. suspectum homolog. The helokinestatin-encoding region contained a single copy of helokinestatin-1, 2 copies of helokinestatin-3 and single copies of 2 novel peptides, (Phe)(5)-helokinestatin-2 (VPPAFVPLVPR) and helokinestatin-6 (GPPFNPPPFVDYEPR). All predicted peptides were found in reverse phase HPLC fractions of the same venom. Synthetic replicates of both novel helokinestatins were found to antagonize the relaxing effect of bradykinin on rat tail artery smooth muscle. Thus lizard venom continues to provide a source of novel biologically active peptides. (C) 2011 Published by Elsevier Inc.

Rapid protein domain assignment from amino acid sequence using predicted secondary structure

The elucidation of the domain content of a given protein sequence in the absence of determined structure or significant sequence homology to known domains is an important problem in structural biology. Here we address how successfully the delineation of continuous domains can be accomplished in the absence of sequence homology using simple baseline methods, an existing prediction algorithm (Domain Guess by Size), and a newly developed method (DomSSEA). The study was undertaken with a view to measuring the usefulness of these prediction methods in terms of their application to fully automatic domain assignment. Thus, the sensitivity of each domain assignment method was measured by calculating the number of correctly assigned top scoring predictions. We have implemented a new continuous domain identification method using the alignment of predicted secondary structures of target sequences against observed secondary structures of chains with known domain boundaries as assigned by Class Architecture Topology Homology (CATH). Taking top predictions only, the success rate of the method in correctly assigning domain number to the representative chain set is 73.3%. The top prediction for domain number and location of domain boundaries was correct for 24% of the multidomain set (±20 residues). These results have been put into context in relation to the results obtained from the other prediction methods assessed

Amino acid sequence of a myotoxic Lys49-phospholipase A(2) homologue from the venom of Cerrophidion (Bothrops) godmani

The complete amino acid sequence of myotoxin II (godMT-II), a myotoxic phospholipase A( 2 )(PLA(2)) homologue from the venom of the Central American crotaline snake Cerrophidion (Bothrops) godmani, was determined by direct protein sequencing methods. GodMT-II is a class II PLA, showing a Lys instead of Asp at position 49. An additional substitution in the calcium binding loop region (Asn instead of Tyr at position 28) suggests the lack of enzymatic activity observed in this toxin is due to loss of its ability to bind the co-factor Ca2+, since the residues involved in forming the catalytic network of PLA(2)s (His-48, Tyr-52 and Asp-99) an conserved in godMT-II. This myotoxin shows highest sequence homology with other Lys-49 PLA(2)s from Bothrops, Agkistrodon and Trimeresurus species, suggesting that they constitute a conserved family of proteins, yet in contrast presents lower homology with Bothrops asper myotoxin III, a catalytically-active PLA(2). The C-terminal region of godMT-II, which is rich in cationic and hydrophobic residues, shares high sequence homology to the corresponding region in the myotoxin II from B. asper, which has been proposed to play an important role in the Ca2+-independent membrane damaging activity. (C) 1998 Elsevier B.V. B.V. All rights reserved.

Amino acid sequence and crystal structure of BaP1, a metalloproteinase from Bothrops asper snake venom that exerts multiple tissue-damaging activities

BaP1 is a 22.7-kD P-I-type zinc-dependent metalloproteinase isolated from the venom of the snake Bothrops asper, a medically relevant species in Central America. This enzyme exerts multiple tissue-damaging activities, including hemorrhage, myonecrosis, dermonecrosis, blistering, and edema. BaP1 is a single chain of 202 amino acids that shows highest sequence identity with metalloproteinases isolated front the venoms of snakes of the subfamily Crotalinae. It has six Cys residues involved in three disulfide bridges (Cys 117-Cys 197, Cys 159-Cys 181, Cys 157-Cys 164). It has the consensus sequence H(142)E(143)XXH(146)XXGXXH(152), as well as the sequence C164I165M166, which characterize the metzincin superfamily of metalloproteinases. The active-site cleft separates a major subdomain (residues 1-152), comprising four a-helices and a five-stranded beta-sheet, from the minor subdomain, which is formed by a single a-helix and several loops. The catalytic zinc ion is coordinated by the N-epsilon2 nitrogen atoms of His 142, His 146, and His 152, in addition to a solvent water molecule, which in turn is bound to Glu 143. Several conserved residues contribute to the formation of the hydrophobic pocket, and Met 166 serves as a hydrophobic base for the active-site groups. Sequence and structural comparisons of hemorrhagic and nonhemorrhagic P-I metalloproteinases from snake venoms revealed differences in several regions. In particular, the loop comprising residues 153 to 176 has marked structural differences between metalloproteinases with very different hemorrhagic activities. Because this region lies in close proximity to the active-site microenvironment, it may influence the interaction of these enzymes with physiologically relevant substrates in the extracellular matrix.

AMINO-ACID-SEQUENCE OF TSTX-V, AN ALPHA-TOXIN FROM TITYUS-SERRULATUS SCORPION-VENOM, AND ITS EFFECT ON K+ PERMEABILITY OF BETA-CELLS FROM ISOLATED RAT ISLETS OF LANGERHANS

Highly purified Tityustoxin V (TsTX-V), an alpha-toxin isolated from the venom of the Brazilian scorpion Tityus serrulatus, was obtained by ion exchange chromatography on carboxymethylcellulose-52. It was shown to be homogeneous by reverse phase high performance liquid chromatography, N-terminal sequencing (first 39 residues) of the reduced and alkylated protein and by polyacrylamide gel electrophoresis in the presence of sodium dodecylsulfate and tricine. Following enzymatic digestion, the complete amino acid sequence (64 residues) was determined. The sequence showed higher homology with the toxins from the venoms of the North African than with those of the North and South American scorpions. Using the rate of Rb-86(+) release from depolarized rat pancreatic beta-cells as a measure of K+ permeability changes, TsTX-V (5.6 mu g/ml) was found to increase by 2.0-2.4-fold the rate of marker outflow in the presence of 8.3 mM glucose. This effect was persistent and slowly reversible, showing similarity to that induced by 100 mu-M veratridine, an agent that increases the open period of Na+ channels, delaying their inactivation. It is suggested that, by extending the depolarized period, TsTX-V indirectly affects beta-cell voltage-dependent K+ channels, thus increasing K+ permeability.

Identities among actin-encoding cDNAs of the Nile tilapia (Oreochromis niloticus) and other eukaryote species revealed by nucleotide and amino acid sequence analyses

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Amino acid sequence of the alpha subunit and computer modelling of the alpha and beta subunits of echicetin from the venom of Echis carinatus (saw-scaled viper)

Echicetin, a heterodimeric protein from the venom of Echis carinatus, binds to platelet glycoprotein Ib (GPIb) and so inhibits platelet aggregation or agglutination induced by various platelet agonists acting via GPIb. The amino acid sequence of the beta subunit of echicetin has been reported and found to belong to the recently identified snake venom subclass of the C-type lectin protein family. Echicetin alpha and beta subunits were purified. N-terminal sequence analysis provided direct evidence that the protein purified was echicetin. The paper presents the complete amino acid sequence of the alpha subunit and computer models of the alpha and beta subunits. The sequence of alpha echicetin is highly similar to the alpha and beta chains of various heterodimeric and homodimeric C-type lectins. Neither of the fully reduced and alkylated alpha or beta subunits of echicetin inhibited the platelet agglutination induced by von Willebrand factor-ristocetin or alpha-thrombin. Earlier reports about the inhibitory activity of reduced and alkylated echicetin beta subunit might have been due to partial reduction of the protein.

Amino Acid Sequence Requirements of the Transmembrane and Cytoplasmic Domains of Influenza Virus Hemagglutinin for Viable Membrane Fusion

The amino acid sequence requirements of the transmembrane (TM) domain and cytoplasmic tail (CT) of the hemagglutinin (HA) of influenza virus in membrane fusion have been investigated. Fusion properties of wild-type HA were compared with those of chimeras consisting of the ectodomain of HA and the TM domain and/or CT of polyimmunoglobulin receptor, a nonviral integral membrane protein. The presence of a CT was not required for fusion. But when a TM domain and CT were present, fusion activity was greater when they were derived from the same protein than derived from different proteins. In fact, the chimera with a TM domain of HA and truncated CT of polyimmunoglobulin receptor did not support full fusion, indicating that the two regions are not functionally independent. Despite the fact that there is wide latitude in the sequence of the TM domain that supports fusion, a point mutation of a semiconserved residue within the TM domain of HA inhibited fusion. The ability of a foreign TM domain to support fusion contradicts the hypothesis that a pore is composed solely of fusion proteins and supports the theory that the TM domain creates fusion pores after a stage of hemifusion has been achieved.

Anti-peptide aptamers recognize amino acid sequence and bind a protein epitope.

In vitro selection of nucleic acid binding species (aptamers) is superficially similar to the immune response. Both processes produce biopolymers that can recognize targets with high affinity and specificity. While antibodies are known to recognize the sequence and conformation of protein surface features (epitopes), very little is known about the precise interactions between aptamers and their epitopes. Therefore, aptamers that could recognize a particular epitope, a peptide fragment of human immunodeficiency virus type I Rev, were selected from a random sequence RNA pool. Several of the selected RNAs could bind the free peptide more tightly than a natural RNA ligand, the Rev-binding element. In accord with the hypothesis that protein and nucleic acid binding cusps are functionally similar, interactions between aptamers and the peptide target could be disrupted by sequence substitutions. Moreover, the aptamers appeared to be able to bind peptides with different solution conformations, implying an induced fit mechanism for binding. Just as anti-peptide antibodies can sometimes recognize the corresponding epitope when presented in a protein, the anti-peptide aptamers were found to specifically bind to Rev.

Global properties of the mapping between local amino acid sequence and local structure in proteins.

Local protein structure prediction efforts have consistently failed to exceed approximately 70% accuracy. We characterize the degeneracy of the mapping from local sequence to local structure responsible for this failure by investigating the extent to which similar sequence segments found in different proteins adopt similar three-dimensional structures. Sequence segments 3-15 residues in length from 154 different protein families are partitioned into neighborhoods containing segments with similar sequences using cluster analysis. The consistency of the sequence-to-structure mapping is assessed by comparing the local structures adopted by sequence segments in the same neighborhood in proteins of known structure. In the 154 families, 45% and 28% of the positions occur in neighborhoods in which one and two local structures predominate, respectively. The sequence patterns that characterize the neighborhoods in the first class probably include virtually all of the short sequence motifs in proteins that consistently occur in a particular local structure. These patterns, many of which occur in transitions between secondary structural elements, are an interesting combination of previously studied and novel motifs. The identification of sequence patterns that consistently occur in one or a small number of local structures in proteins should contribute to the prediction of protein structure from sequence.