Designing amino acid sequences to fold with good hydrophobic cores

We present two methods for designing amino acid sequences of proteins that will fold to have good hydrophobic cores. Given the coordinates of the desired target protein or polymer structure, the methods generate sequences of hydrophobic (H) and polar (P) monomers that are intended to fold to these structures. One method designs hydrophobic inside, polar outside; the other minimizes an energy function in a sequence evolution process. The sequences generated by these methods agree at the level of 60-80% of the sequence positions in 20 proteins in the Protein Data Bank. A major challenge in protein design is to create sequences that can fold uniquely, i.e. to a single conformation rather than to many. While an earlier lattice-based sequence evolution method was shown not to design unique folders, our method generates unique folders in lattice model tests. These methods may also be useful in designing other types of foldable polymer not based on amino acids.     

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

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 identify using the run test statistic (ro) of Mood (1940, Ann. Math. Stat. 11, 367-392). The probability density of ro 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 alpha-helix propensity show a strong tendency to cluster whereas those with beta-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.     

Nucleotide and deduced amino acid sequences of rat myosin binding protein H (MyBP-H)

The complete nucleotide sequence of the cDNA clone encoding rat skeletal muscle myosin-binding protein H (MyBP-H) was determined and amino acid sequence was deduced from the nucleotide sequence (GenBank accession number AF077338). The full-length cDNA of 1782 base pairs(bp) contains a single open reading frame of 1454 bp encoding a rat MyBP-H protein of the predicted molecular mass 52.7 kDa and includes the common consensus 'CA__TG' protein binding motif. The cDNA sequence of rat MyBP-H show 92%, 84% and 41% homology with those of mouse, human and chicken, respectively. The protein contains tandem internal motifs array (-FN III-Ig C2-FN III-Ig C2-) in the C-terminal region which resembles to the immunoglobulin superfamily C2 and fibronectin type III motifs. The amino acid sequence of the C-terminal Ig C2 was highly conserved among MyBPs family and other thick filament binding proteins, suggesting that the C-terminal Ig C2 might play an important role in its function. All proteins belonging to MyBP-H member contains 'RKPS' sequence which is assumed to be cAMP- and cGMP-dependent protein kinase A phosphorylation site. Computer analysis of the primary sequence of rat MyBP-H predicted 11 protein kinase C (PKC) phosphorylation site, 7 casein kinase II (CK2) phosphorylation site and 4 N-myristoylation site.     

Representation of amino acid sequences as two-dimensional point patterns

An algorithm for the representation of amino acid sequences as two-dimensional point patterns (2-D plot) is described. The algorithm is based on chaos game representation (CGR) for DNA sequences and was extended for amino acid sequences. The 2-D plot depicts the sequentiality of amino acids and the amino acid composition of a protein. Changes in a protein sequence as insertion, deletion and repeats of amino acids are characterized by specific geometrical properties and changes in the 2-D plots. The 2-D plot may be considered as a two-dimensional "fingerprint" of a protein. The properties of the algorithm are explained by user-defined amino acid sequences. As an example the 2-D plots of two selected heart proteins are generated. The sequences of these proteins are obtained from the protein sequence database SWISS-PROT.     

Analysis of amino acid sequences of Rhodobacter sphaeroides glutamine synthetase

A comparative analysis of the amino acid sequence of glutamine synthetase (GS) of the photosynthetic purple bacterium Rhodobacter sphaeroides revealed that the enzyme is typical for first type procaryotic GSs and structurally resembles GSs of enteric bacteria. The data obtained indicate that the complex phenotype of purple bacterial mutants at the glnA gene coding for GS may be conditioned by specific regulation of nitrogen metabolism in bacterial cells rather than by structural-and-functional peculiarities of GS.     

Evolution of DNA or amino acid sequences with dependent sites

A framework is outlined to study the evolution of DNA or amino acid sequences, if sequence sites do not evolve independently. The units of evolution are nonoverlapping subsequences of length l. Each subsequence evolves independently of the others, but within a subsequence the sequences show a Markov order one dependency. We describe an algorithm to mimic the evolution of such sequences. The influence of dependencies between sites on distance estimates and the reliability of tree reconstruction methods is investigated. We show that an inappropriate model of sequence evolution in the tree reconstruction process will lead to a nonempty Felsenstein zone. Finally, we describe a method to infer l from sequence data. Examples from the evolution of DNA sequences as well as from amino acids are given.     

Amino acid sequence of human muscle glyceraldehyde-3-phosphate dehydrogenase. Isolation and amino acid sequences of tryptic peptides

1. The amino acid compostion, N- and C-terminal amino acid sequences, and the subunit molecular weight of glyceraldehyde phosphate dehydrogenase from human muscle, were determined. The obtained results and the maps of tryptic peptides suggest that the enzyme is composed of four identical or very similar polypeptide chains. 2. From the tryptic digest of performic acid-oxidized enzyme, 32 peptides were isolated. The amino acid sequence analysis showed a high degree of homology with the corresponding tryptic peptides of the dehydrogenase from pig muscle, with 9 replacements and probably two additional amino acids in the examined sequences of the human muscle enzyme.     

Purification and partial amino acid sequences of the binding protein from Bombyx mori for CryIAa delta-endotoxin of Bacillus thuringiensis

The binding protein for Bacillus thuringiensis delta-endotoxin, CryIAa, from the brush border membrane of the midgut of Bombyx mori was purified by the dot blot method and delta-endotoxin affinity chromatography. The binding protein was purified to 235-fold enrichment from cholic acid extracts of brush border membranes from B. mori midgut by activated CryIAa-affinity chromatography and DEAE ion-exchange chromatography. The purified binding protein showed a single band of 180 kDa by sodium dodecyl sulfate polyacrylamide gel electrophoresis and this band specifically reacted to 125I-labeled CryIAa on Immobilon membrane. The affinity of the binding protein for CryIAa was equivalent to that of the brush border membrane vesicles and solubilized membrane proteins. Partial amino acid sequences of the binding protein showed sequence similarity to the cadherin-like binding protein for CryIAb from Manduca sexta, but not for CryIAc binding protein from M. sexta and Heliothis virescens.     

Estimating the number of protein folds

A number of fundamental questions in structural biology concern the diversity of protein architectures (or folds). Here, we address two of them, the size of the universe of folds, and the distribution of sequence families among them, using an analysis based on a new and rigorous statistical sampling method. In particular we show that the number of known non-transmembrane protein folds is approximately one half of the total that exist, and that certain superfolds should exist, which accommodate dozens of non-homologous sequence families.     

A new method of inference of ancestral nucleotide and amino acid sequences

A statistical method was developed for reconstructing the nucleotide or amino acid sequences of extinct ancestors, given the phylogeny and sequences of the extant species. A model of nucleotide or amino acid substitution was employed to analyze data of the present-day sequences, and maximum likelihood estimates of parameters such as branch lengths were used to compare the posterior probabilities of assignments of character states (nucleotides or amino acids) to interior nodes of the tree; the assignment having the highest probability was the best reconstruction at the site. The lysozyme c sequences of six mammals were analyzed by using the likelihood and parsimony methods. The new likelihood-based method was found to be superior to the parsimony method. The probability that the amino acids for all interior nodes at a site reconstructed by the new method are correct was calculated to be 0.91, 0.86, and 0.73 for all, variable, and parsimony-informative sites, respectively, whereas the corresponding probabilities for the parsimony method were 0.84, 0.76, and 0.51, respectively. The probability that an amino acid in an ancestral sequence is correctly reconstructed by the likelihood analysis ranged from 91.3 to 98.7% for the four ancestral sequences.     

A tobacco soluble protoporphyrinogen-oxidizing enzyme similar to plant peroxidases in their amino acid sequences and immunochemical reactivity

Partial amino acid sequences of a soluble protoporphyrinogen-oxidizing enzyme (PPO) purified from tobacco cultured cells were identified. The sequences of two regions of the soluble PPO corresponded to the acid/base catalysis and heme binding regions of plant peroxidases. Anti-soluble PPO IgG cross-reacted with these plant peroxidases. Thus, the soluble PPO seems to be a kind of peroxidase.     

Rice embryo globulins: amino-terminal amino acid sequences, cDNA cloning and expression

A globulin fraction prepared from rice embryos contained polypeptides or polypeptide groups of 49 kDa (designated REG1), 46 kDa (designated REG2), about 35 kDa, 32 kDa and 25 kDa. The amino-terminal sequences of REG1 and the major polypeptide in the 35-kDa group were identical, suggesting that the REG1 polypeptide undergoes partial proteolytic processing that removes a carboxy-terminal region. A cDNA clone, designated pcREG2, encoding REG2 was isolated, and its nucleotide sequence was determined. The deduced amino acid sequence of REG2 was found to be 68% identical to that of the maize GLB2 globulin. Reg2 mRNA was present at high levels during embryo development for up to 14 days after flowering (DAF). Lower levels were found 20 DAF when the maturation of embryos was almost completed, and at the dry mature stage. Reg2 mRNA almost disappeared upon imbibition of isolated dry mature embryos but it was re-induced at a low level by further treatment with ABA. The expression of Reg2 was not induced by ABA in suspension-cultured cells, unlike that of Osem, one of the late embryogenesis abundant protein (LEA) genes.     

Recombinant protein sequences can trigger methylation of N-terminal amino acids in Escherichia coli

Recombinant human hemoglobin rHb1.1 has been genetically engineered with the replacement of the wild-type valine residues at all N-termini with methionine, an Asn 108 Lys substitution on the beta globins, and a fusion of the two alpha globins with a glycine linker. When rHb1.1 was expressed in Escherichia coli, methylation of the N-terminal methionine of the alpha globin was discovered. Another mutant has been engineered with the alpha globin gene coding for N-terminal methionine followed by an insertion of alanine. Characterization of expressed hemoglobin from this variant revealed a methylated N-terminal alanine that occurred through two posttranslational events: initial excision of the N-terminal methionine, followed by methylation of alanine as the newly generated N-terminus. No methylation was observed for variants expressed with wild-type valine at the N-terminus of the alpha globin. The methylation of N-terminal amino acids was attributed to a specific protein sequence that can trigger methylation of proteins expressed in E. coli. Here we demonstrate that proline at position 4 in the protein sequence of alpha globin seems an essential part of that signaling. Although N-terminal methylation has been observed previously for native E. coli proteins with similar N-terminal sequences, methylation of the recombinant globins has allowed further delineation of the recognition sequence, and indicates that methylation of heterologous proteins can occur in E. coli.     

Patterns of amino acid variability in NSI-like and SI-like V3 sequences and a linked change in the CD4-binding domain of the HIV-1 Env protein

We used intracellular recording techniques to investigate the actions of clonidine on hypoglossal motoneurons (HMs) in rat brainstem slices. Clonidine (10-100 microM) produced a small (2-6 mV), dose-dependent hyperpolarization in HMs, accompanied by an increase in peak input resistance (RN). It also slowed the time course of the depolarizing 'sag' of the voltage response to constant hyperpolarizing current steps. These effects were mimicked by the alpha2-adrenoceptor (alpha2-AR) agonist guanabenz, but not by the Ih-imidazoline receptor agonists moxonidine or rilmenidine. Recorded in single-electrode voltage clamp mode, clonidine decreased input conductance of HMs and reduced the amplitude of a hyperpolarization-activated inward current (Ih). Clonidine's effect on Ih was three-fold: it shifted the half-activation voltage (V1/2) in the hyperpolarizing direction (by 4.4 +/- 0.7 mV at a dose of 10 microM), decreased the maximal current (by approximately 20%), and slowed the time course of Ih activation at all voltage steps. At the most hyperpolarized potential steps, clonidine slowed activation of Ih dramatically, yielding a striking increase in the activation time constant. The alpha2-AR antagonists yohimbine and idazoxan reduced clonidine's effect on V1/2 and on the Ih activation time course, but neither blocked clonidine's reduction of the maximal current, nor its strong slowing of Ih activation at the most hyperpolarized steps. We were unable to mimic or occlude clonidine's actions with the adenylate cyclase inhibitor SQ 22536 nor with the non-specific protein kinase inhibitor H-7. We conclude that clonidine hyperpolarizes HMs via a reduction of the amount of Ih that is active at rest, and that the response is mediated in part by alpha2-ARs. Some effects of clonidine on these neurons do not appear to be receptor-mediated, and may be due to physical block by clonidine of Ih channels.     

Disallowed Ramachandran conformations of amino acid residues in protein structures

An analysis of the nature and distribution of disallowed Ramachandran conformations of amino acid residues observed in high resolution protein crystal structures has been carried out. A data set consisting of 110 high resolution, non-homologous, protein crystal structures from the Brookhaven Protein Data Bank was examined. The data set consisted of a total of 18,708 non-Gly residues, which were characterized on the basis of their backbone dihedral angles (phi, psi). Residues falling outside the defined "broad allowed limits" on the Ramachandran map were chosen and the reported B-factor value of the alpha-carbon atom was used to further select well defined disallowed conformations. The conformations of the selected 66 disallowed residues clustered in distinct regions of the Ramachandran map indicating that specific phi, psi angle distortions are preferred under compulsions imposed by local constraints. The distribution of various amino acid residues in the disallowed residue data set showed a predominance of small polar/charged residues, with bulky hydrophobic residues being infrequent. As a further check, for all the 66 cases non-hydrogen van der Waals short contacts in the protein structures were evaluated and compared with the ideal "Ala-dipeptide" constructed using disallowed dihedral angle (phi, psi) values. The analysis reveals that short contacts are eliminated in most cases by local distortions of bond angles. An analysis of the conformation of the identified disallowed residues in related protein structures reveals instances of conservation of unusual stereochemistry.     

Barley beta-glucosidase: expression during seed germination and maturation and partial amino acid sequences

Confusion regarding proper use of the terms rate and risk persists in the literature. This has implications for the proper modeling of prognosis and transition between health states in decision analysis and related techniques. The issue is complicated by the plethora of terms related to rate and risk. Although the suggestion to use the terms force and probability as substitutes for rate and risk has some appeal, the change in terminology by itself is unlikely to solve all the confusion or misuse of terms. This paper clarifies the proper definitions and estimations of rates and risks and suggests critical factors for the decision analyst to remember when using, modeling, or interpreting transition rates and risks.