Detection of native-like models for amino acid sequences of unknown three-dimensional structure in a data base of known protein conformations.

We present an approach which can be used to identify native-like folds in a data base of protein conformations in the absence of any sequence homology to proteins in the data base. The method is based on a knowledge-based force field derived from a set of known protein conformations. A given sequence is mounted on all conformations in the data base and the associated energies are calculated. Using several conformations and sequences from the globin family we show that the native conformation is identified correctly. In fact the resolution of the force field is high enough to discriminate between a native fold and several closely related conformations. We then apply the procedure to several globins of known sequence but unknown three dimensional structure. The homology of these sequences to globins of known structures in the data base ranges from 49 to 17%. With one exception we find that for all globin sequences one of the known globin folds is identified as the most favorable conformation. These results are obtained using a force field derived from a data base devoid of globins of known structure. We briefly discuss useful applications in protein structural research and future development of our approach.     

Application of a directed conformational search for generating 3-D coordinates for protein structures from alpha-carbon coordinates.

A directed conformational search algorithm using the program CONGEN (ref. 3), which samples backbone conformers, is described. The search technique uses information from the partially built structures to direct the search process and is tested on the problem of generating a full set of backbone Cartesian coordinates given only alpha-carbon coordinates. The method has been tested on six proteins of known structure, varying in size and classification, and was able to generate the original backbone coordinates with RMSs ranging from 0.30-0.87A for the alpha-carbons and 0.5-0.99A RMSs for the backbone atoms. Cis peptide linkages were also correctly identified. The procedure was also applied to two proteins available with only alpha-carbon coordinates in the Brookhaven Protein Data Bank; thioredoxin (SRX) and triacylglycerol acylhydrolase (TGL). All-atom models are proposed for the backbone of both these proteins. In addition, the technique was applied to randomized coordinates of flavodoxin to assess the effects of irregularities in the data on the final RMS. This study represents the first time a deterministic conformational search was used on such a large scale.     

Computation of relaxation matrix elements from incomplete NOESY data sets

Summary The structural determination of biological molecules in solution by NMR relies on the determination of a set of interatomic distances obtained by measurement of intramolecular nuclear Overhauser effects (NOE). It is shown in this paper that it is possible to obtain the accurate relaxation rate (and hence the interatomic distance) from the direct measurement of a single NOE signal. The precise analysis of a NOESY peak evolution with respect to the mixing time allows the evaluation of the relaxation parameters for the pair of spins under consideration. This is done without any assumption on the relaxation of unmeasured spins, or on the movement of the molecule. The theoretical basis of this method is presented. In order to evaluate the proposed method, a simulated case on the protein BPTI is studied, which shows that the method performs very well even in the case of noisy data sets.     

Processing of multi-dimensional NMR data with the new software PROSA

Summary The new program PROSA is an efficient implementation of the common data-processing steps for multi-dimensional NMR spectra. PROSA performs linear prediction, digital filtering, Fourier transformation, automatic phase correction, and baseline correction. High efficiency is achieved by avoiding disk storage of intermediate data and by the absence of any graphics display, which enables calculation in the batch mode and facilitates porting PROSA on a variety of different computer systems; including supercomputers. Furthermore, all time-consuming routines are completely vectorized. The elimination of a graphics display was made possible by the use of a new, reliable automatic phase-correction routine. CPU times for complete processing of a typical heteronuclear three-dimensional NMR data set of a protein vary between less than 1 min on a NEC SX3 supercomputer and 40 min on a Sun-4 computer system.     

Environment-specific amino acid substitution tables: tertiary templates and prediction of protein folds.

The local environment of an amino acid in a folded protein determines the acceptability of mutations at that position. In order to characterize and quantify these structural constraints, we have made a comparative analysis of families of homologous proteins. Residues in each structure are classified according to amino acid type, secondary structure, accessibility of the side chain, and existence of hydrogen bonds from the side chains. Analysis of the pattern of observed substitutions as a function of local environment shows that there are distinct patterns, especially for buried polar residues. The substitution data tables are available on diskette with Protein Science. Given the fold of a protein, one is able to predict sequences compatible with the fold (profiles or templates) and potentially to discriminate between a correctly folded and misfolded protein. Conversely, analysis of residue variation across a family of aligned sequences in terms of substitution profiles can allow prediction of secondary structure or tertiary environment.     

The development of the Jamming Avoidance Response (JAR) in Eigenmannia: An innate behavior indeed

Summary In its Jamming Avoidance Response (JAR), the gymnotiform electric fish Eigenmannia shifts its electric organ discharge (EOD) frequency away from similar interfering frequencies. Continual behavioral measurements were carried out in 164 juvenile fish until a correct JAR emerged. Sixty-four of these fish were raised in complete isolation, the remainder in a community of their siblings. A correct JAR emerged in fish of 1.2–1.6 cm in body length, corresponding to a developmental age of 24–32 days. In 6 of 164 fish, the emergence of a correct JAR followed an interim appearance of an incorrect JAR, which involved frequency shifts in the direction opposite to those of a correct JAR. The fish raised in isolation developed the same forms of behavior and showed the same sequence in their appearance as did socially raised fish. This indicates that the JAR and its developmental schedule are innate. The appearance of an incorrect JAR suggests initial errors or incompleteness in the wiring of central nervous connections. A correct JAR ultimately emerged even if a stimulus regimen was offered that rewarded frequency shifts in the direction opposite to those of a correct JAR. This indicates that the development of the JAR is immune to experimental alterations of sensory experience.Abbreviations Df frequency difference between a jamming stimulus and fish's EOD - ELL electrosensory lateral line lobe - EO electric organ - EOD electric organ discharge - JAR Jamming Avoidance Response - nE nucleus electrosensorius - nE subnucleus of nE, causing drop of EOD frequency - nE subnucleus of nE, causing rise of EOD frequency - Pn pacemaker nucleus - PPn prepacemaker nucleus     

The tertiary structure of Aspergillus saitoi minor ribonuclease (Ms) predicted from the structure of RNase T1

Ribonuclease Ms from Aspergillus saitoi is a small acidic protein (11 714 Da) containing 106 amino acids of known sequence. Unlike other enzymes belonging to the RNase T₁ family this ribonuclease is base-unspecific. Using interactive computer graphics and energy minimisation we predicted the structure of RNase Ms on the basis of sequence homology to RNase T₁ of known structure. In this report the predicted structure of this protein is presented and characterised.     

Genotype dependent variation in mycorrhizal colonization and response to inoculation of pearl millet

Summary Genotypes of pearl millet (Pennisetum americanum L. Leeke) were examined for differences in vesicular-arbuscular mycorrhizal (VAM) colonization and response to inoculation. For thirty genotypes tested across three field locations there was a range of mycorrhizal colonization intensity between 25 and 56%. In another experiment with two male-sterile lines, restorer lines and their derived crosses, grown in pots filled with non-sterilized soil there were significant differences between genotypes for colonization by mycorrhiza. This showed hostgenotype dependence for mycorrhizal colonization.Root growth rates, mycorrhizal root length, percentage root colonization and plant growth and P uptake were studied in ten genotypes. A set of 3 genotypes with similar root lengths varied significantly with regard to mycorrhizal root length and the percentage colonization. This supports the suggestion that VAM colonization and spread is dependent on the host genotype. The growth responses differed significantly between the genotypes and they also differed in their responses to P uptake and VAM inoculation. The utility of host-genotype dependent differences in VAM symbiosis in plant breeding is discussed.Journal Article No. 453     

A critique of the utility of the prediction of protein secondary structure

The Chou-Fasman predictive algorithm for determining the secondary structure of proteins from the primary sequence is reviewed. Many examples of its use are presented which illustrate its wide applicability, such as predicting (a) regions with the potential for conformational change, (b) sequences which are capable of assuming several conformations in different environments, (c) effects of single amino acid mutations, (d) amino acid replacements in synthesis of peptides to bring about a change in conformation, (e) guide to the synthesis of polypeptides with definitive secondary structure,e.g. signal sequences, (f) conformational homologues from varying sequences and (g) the amino acid requirements for amphiphilicα-helical peptides.     

Prediction and comparison of the secondary structure of legume lectins

Secondary structure prediction for the 4 legume lectins: Concanavalin A, soybean agglutinin, favabean lectin and lentil lectin, was done by the method of Chou and Fasman. This prediction shows that these four lectins fall into a structurally distinct class of proteins, containing high amounts of β-sheet and β-turns. There is a notable similarity in the gross structure of these proteins; all four of them contain about 40–50% of β-sheet, 35–45 % β-turn and 0–10% of α-helix. When the secondary structure of corresponding residues in each pair of these lectins was compared, there was a striking similarity in the Concanavalin A-soybean agglutinin and favabean lectin-lentil lectin pairs, and considerably less similarity in the other pairs, suggesting that these legume lectins have probably evolved in a divergent manner from a common ancestor. A comparison of the predicted potential β-turn sites also supports the hypothesis of divergent evolution in this class of lectins.