Detection of protein fold similarity based on correlation of amino acid properties

An increasing number of proteins with weak sequence similarity have been found to assume similar three-dimensional fold and often have similar or related biochemical or biophysical functions. We propose a method for detecting the fold similarity between two proteins with low sequence similarity based on their amino acid properties alone. The method, the proximity correlation matrix (PCM) method, is built on the observation that the physical properties of neighboring amino acid residues in sequence at structurally equivalent positions of two proteins of similar fold are often correlated even when amino acid sequences are different. The hydrophobicity is shown to be the most strongly correlated property for all protein fold classes. The PCM method was tested on 420 proteins belonging to 64 different known folds, each having at least three proteins with little sequence similarity. The method was able to detect fold similarities for 40% of the 420 sequences. Compared with sequence comparison and several fold-recognition methods, the method demonstrates good performance in detecting fold similarities among the proteins with low sequence identity. Applied to the complete genome of Methanococcus jannaschii, the method recognized the folds for 22 hypothetical proteins.

Herbicide sensitivity determinant of wheat plastid acetyl-CoA carboxylase is located in a 400-amino acid fragment of the carboxyltransferase domain

A series of chimeral genes, consisting of the yeast GAL10 promoter, yeast ACC1 leader, wheat acetyl-CoA carboxylase (ACCase; EC 6.4.1.2) cDNA, and yeast ACC1 3′-tail, was used to complement a yeast ACC1 mutation. These genes encode a full-length plastid enzyme, with and without the putative chloroplast transit peptide, as well as five chimeric cytosolic/plastid proteins. Four of the genes, all containing at least half of the wheat cytosolic ACCase coding region at the 5′-end, complement the yeast mutation. Aryloxyphenoxypropionate and cyclohexanedione herbicides, at concentrations below 10 μM, inhibit the growth of haploid yeast strains that express two of the chimeric ACCases. This inhibition resembles the inhibition of wheat plastid ACCase observed in vitro and in vivo. The differential response to herbicides localizes the sensitivity determinant to the third quarter of the multidomain plastid ACCase. Sequence comparisons of different multidomain and multisubunit ACCases suggest that this region includes part of the carboxyltransferase domain, and therefore that the carboxyltransferase activity of ACCase (second half-reaction) is the target of the inhibitors. The highly sensitive yeast gene-replacement strains described here provide a convenient system to study herbicide interaction with the enzyme and a powerful screening system for new inhibitors.

CKAAPs DB: a conserved key amino acid positions database

The Conserved Key Amino Acid Positions DataBase (CKAAPs DB) provides access to an analysis of structurally similar proteins with dissimilar sequences where key residues within a common fold are identified. The derivation and significance of CKAAPs starting from pairwise structure alignments is described fully in Reddy et al. [Reddy,B.V.B., Li,W.W., Shindyalov,I.N. and Bourne,P.E. (2000) Proteins, in press]. The CKAAPs identified from this theoretical analysis are provided to experimentalists and theoreticians for potential use in protein engineering and modeling. It has been suggested that CKAAPs may be crucial features for protein folding, structural stability and function. Over 170 substructures, as defined by the Combinatorial Extension (CE) database, which are found in approximately 3000 representative polypeptide chains have been analyzed and are available in the CKAAPs DB. CKAAPs DB also provides CKAAPs of the representative set of proteins derived from the CE and FSSP databases. Thus the database contains over 5000 representative poly­peptide chains, covering all known structures in the PDB. A web interface to a relational database permits fast retrieval of structure-sequence alignments, CKAAPs and associated statistics. Users may query by PDB ID, protein name, function and Enzyme Classification number. Users may also submit protein alignments of their own to obtain CKAAPs. An interface to display CKAAPs on each structure from a web browser is also being implemented. CKAAPs DB is maintained by the San Diego Supercomputer Center and accessible at the URL http://ckaaps.sdsc.edu.

Interactions between tRNA identity nucleotides and their recognition sites in glutaminyl-tRNA synthetase determine the cognate amino acid affinity of the enzyme.

Sequence-specific interactions between aminoacyl-tRNA synthetases and their cognate tRNAs both ensure accurate RNA recognition and prevent the binding of noncognate substrates. Here we show for Escherichia coli glutaminyl-tRNA synthetase (GlnRS; EC 6.1.1.18) that the accuracy of tRNA recognition also determines the efficiency of cognate amino acid recognition. Steady-state kinetics revealed that interactions between tRNA identity nucleotides and their recognition sites in the enzyme modulate the amino acid affinity of GlnRS. Perturbation of any of the protein-RNA interactions through mutation of either component led to considerable changes in glutamine affinity with the most marked effects seen at the discriminator base, the 10:25 base pair, and the anticodon. Reexamination of the identity set of tRNA(Gln) in the light of these results indicates that its constituents can be differentiated based upon biochemical function and their contribution to the apparent Gibbs' free energy of tRNA binding. Interactions with the acceptor stem act as strong determinants of tRNA specificity, with the discriminator base positioning the 3' end. The 10:25 base pair and U35 are apparently the major binding sites to GlnRS, with G36 contributing both to binding and recognition. Furthermore, we show that E. coli tryptophanyl-tRNA synthetase also displays tRNA-dependent changes in tryptophan affinity when charging a noncognate tRNA. The ability of tRNA to optimize amino acid recognition reveals a novel mechanism for maintaining translational fidelity and also provides a strong basis for the coevolution of tRNAs and their cognate synthetases.

Posttranslational amino acid epimerization: enzyme-catalyzed isomerization of amino acid residues in peptide chains.

Since ribosomally mediated protein biosynthesis is confined to the L-amino acid pool, the presence of D-amino acids in peptides was considered for many years to be restricted to proteins of prokaryotic origin. Unicellular microorganisms have been responsible for the generation of a host of D-amino acid-containing peptide antibiotics (gramicidin, actinomycin, bacitracin, polymyxins). Recently, a series of mu and delta opioid receptor agonists [dermorphins and deltorphins] and neuroactive tetrapeptides containing a D-amino acid residue have been isolated from amphibian (frog) skin and mollusks. Amino acid sequences obtained from the cDNA libraries coincide with the observed dermorphin and deltorphin sequences, suggesting a stereospecific posttranslational amino acid isomerization of unknown mechanism. A cofactor-independent serine isomerase found in the venom of the Agelenopsis aperta spider provides the first major clue to explain how multicellular organisms are capable of incorporating single D-amino acid residues into these and other eukaryotic peptides. The enzyme is capable of isomerizing serine, cysteine, O-methylserine, and alanine residues in the middle of peptide chains, thereby providing a biochemical capability that, until now, had not been observed. Both D- and L-amino acid residues are susceptible to isomerization. The substrates share a common Leu-Xaa-Phe-Ala recognition site. Early in the reaction sequence, solvent-derived deuterium resides solely with the epimerized product (not substrate) in isomerizations carried out in 2H2O. Significant deuterium isotope effects are obtained in these reactions in addition to isomerizations of isotopically labeled substrates (2H at the epimerizeable serine alpha-carbon atom). The combined kinetic and structural data suggests a two-base mechanism in which abstraction of a proton from one face is concomitant with delivery from the opposite face by the conjugate acid of the second enzymic base.

Mass spectrometric amino acid sequencing of a mixture of seed storage proteins (napin) from Brassica napus, products of a multigene family.

The amino acid sequences of a number of closely related proteins ("napin") isolated from Brassica napus were determined by mass spectrometry without prior separation into individual components. Some of these proteins correspond to those previously deduced (napA, BngNAP1, and gNa), chiefly from DNA sequences. Others were found to differ to a varying extent (BngNAP1', BngNAP1A, BngNAP1B, BngNAP1C, gNa', and gNaA). The short chains of gNa and gNa' and of BngNAP1 and BngNAP1' differ by the replacement of N-terminal proline by pyroglutamic acid; the long chains of gNaA and BngNAP1B contain a six amino acid stretch, MQGQQM, which is present in gNa (according to its DNA sequence) but absent from BngNAP1 and BngNAP1C. These alternations of sequences between napin isoforms are most likely due to homologous recombination of the genetic material, but some of the changes may also be due to RNA editing. The amino acids that follow the untruncated C termini of those napin chains for which the DNA sequences are known (napA, BngNAP1, and gNa) are aromatic amino acids. This suggests that the processing of the proprotein leading to the C termini of the two chains is due to the action of a protease that specifically cleaves a G/S-F/Y/W bond.

Single-amino acid substitutions eliminate lysine inhibition of maize dihydrodipicolinate synthase.

Dihydrodipicolinate synthase (DHPS; EC 4.2.1.52) catalyzes the first step in biosynthesis of lysine in plants and bacteria. DHPS in plants is highly sensitive to end-product inhibition by lysine and, therefore, has an important role in regulating metabolite flux into lysine. To better understand the feedback inhibition properties of the plant enzyme, we transformed a maize cDNA for lysine-sensitive DHPS into an Escherichia coli strain lacking DHPS activity. Cells were mutagenized with ethylmethanesulfonate, and potential DHPS mutants were selected by growth on minimal medium containing the inhibitory lysine analogue S-2-aminoethyl-L-cysteine. DHPS assays identified surviving colonies expressing lysine-insensitive DHPS activity. Ten single-base-pair mutations were identified in the maize DHPS cDNA sequence; these mutations were specific to one of three amino acid residues (amino acids 157, 162, and 166) localized within a short region of the polypeptide. No other mutations were present in the remaining DHPS cDNA sequence, indicating that altering only one of the three residues suffices to eliminate lysine inhibition of maize DHPS. Identification of these specific mutations that change the highly sensitive maize DHPS to a lysine-insensitive isoform will help resolve the lysine-binding mechanism and the resultant conformational changes involved in inhibition of DHPS activity. The plant-derived mutant DHPS genes may also be used to improve nutritional quality of maize or other cereal grains that have inadequate lysine content when fed to animals such as poultry, swine, or humans.

Conformation, energy, and folding ability of selected amino acid sequences.

Evolutionary selection of sequences is studied with a knowledge-based Hamiltonian to find the design principle for folding to a model protein structure. With sequences selected by naive energy minimization, the model structure tends to be unstable and the folding ability is low. Sequences with high folding ability have only the low-lying energy minimum but also an energy landscape which is similar to that found for the native sequence over a wide region of the conformation space. Though there is a large fluctuation in foldable sequences, the hydrophobicity pattern and the glycine locations are preserved among them. Implications of the design principle for the molecular mechanism of folding are discussed.

On the distribution of amino acid residues in transmembrane alpha-helix bundles.

The periodic distribution of residues in the sequence of 469 putative transmembrane alpha-helices from eukaryotic plasma membrane polytopic proteins has been analyzed with correlation matrices. The method does not involve any a priori assumption about the secondary structure of the segments or about the physicochemical properties of individual amino acid residues. Maximal correlation is observed at 3.6 residues per period, characteristic of alpha-helices. A scale extracted from the data describes the propensity of the various residues to lie on the same or on opposite helix faces. The most polar face of transmembrane helices, presumably that buried in the protein core, shows a strong enrichment in aromatic residues, while residues likely to face the fatty acyl chains of lipids are largely aliphatic.

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.

On the utility of alternative amino acid scripts

In this work we propose the hypothesis that replacing the current system of representing the chemical entities known as amino acids using Latin letters with one of several possible alternative symbolic representations will bring significant benefits to the human construction, modification, and analysis of multiple protein sequence alignments. We propose ways in which this might be done without prescribing the choice of actual scripts used. Specifically we propose and explore three ways to encode amino acid texts using novel symbolic alphabets free from precedents. Primary orthographic encoding is the direct substitution of a new alphabet for the standard, Latin-based amino acid code. Secondary encoding imposes static residue groupings onto the orthography of the alphabet by manipulating the shape and/or orientation of amino acid symbols. Tertiary encoding renders each residue as a composite symbol; each such symbol thus representing several alternative amino acid groupings simultaneously. We also propose that the use of a new group-focussed alphabet will free the colouring of amino acid residues often used as a tool to facilitate the representation or construction of multiple alignments for other purposes, possibly to indicate dynamic properties of an alignment such as position-wise residue conservation.

Identifying candidate subunit vaccines using an alignment-independent method based on principal amino acid properties

Subunit vaccine discovery is an accepted clinical priority. The empirical approach is time- and labor-consuming and can often end in failure. Rational information-driven approaches can overcome these limitations in a fast and efficient manner. However, informatics solutions require reliable algorithms for antigen identification. All known algorithms use sequence similarity to identify antigens. However, antigenicity may be encoded subtly in a sequence and may not be directly identifiable by sequence alignment. We propose a new alignment-independent method for antigen recognition based on the principal chemical properties of protein amino acid sequences. The method is tested by cross-validation on a training set of bacterial antigens and external validation on a test set of known antigens. The prediction accuracy is 83% for the cross-validation and 80% for the external test set. Our approach is accurate and robust, and provides a potent tool for the in silico discovery of medically relevant subunit vaccines.

Identification of DNA-binding proteins by incorporating evolutionary information into pseudo amino acid composition via the top-n-gram approach

DNA-binding proteins are crucial for various cellular processes and hence have become an important target for both basic research and drug development. With the avalanche of protein sequences generated in the postgenomic age, it is highly desired to establish an automated method for rapidly and accurately identifying DNA-binding proteins based on their sequence information alone. Owing to the fact that all biological species have developed beginning from a very limited number of ancestral species, it is important to take into account the evolutionary information in developing such a high-throughput tool. In view of this, a new predictor was proposed by incorporating the evolutionary information into the general form of pseudo amino acid composition via the top-n-gram approach. It was observed by comparing the new predictor with the existing methods via both jackknife test and independent data-set test that the new predictor outperformed its counterparts. It is anticipated that the new predictor may become a useful vehicle for identifying DNA-binding proteins. It has not escaped our notice that the novel approach to extract evolutionary information into the formulation of statistical samples can be used to identify many other protein attributes as well.