Two-Angstrom Crystal Structure of Oxidized Chromatium High Potential Iron Protein

Three-dimensional solution structure of the oxidized high potential iron-sulfur protein from Chromatium vinosum through NMR. Comparative analysis with the solution structure of the reduced species

The NMR solution structure of the oxidized HiPIP from Chromatium vinosum has been solved. Despite the fact that the protein is paramagnetic, 85% of the 1H and 80% of the 15N signals have been assigned. Through 1537 NOEs, out of which 1142 were found to be relevant for the structure determination, a family of structures has been obtained by distance geometry calculations. These structures have then been subjected to restrained energy minimization (REM) and restrained molecular dynamics (RMD) calculations in vacuum. Finally, the mean structure of the RMD family has been treated through RMD in water. The RMSD values for the backbone and heavy atoms within the RMD family are 0.57 +/- 0.14 and 1.08 +/- 0.16 A, respectively. These values together with other parameters indicate that the structure is of good quality and as good as the structure of the reduced protein. The RMDw structures of the reduced and oxidized proteins are different beyond the experimental indetermination. The set of constraints for the reduced and oxidized forms have been used to treat the available X-ray structure by RMD in water. The two structures generated in this way are quite similar to their respective solution structures, thus confirming that the experimental constraints are capable of yielding two different structures from the same starting structural model. This is the first time that independently determined solution structures of two redox states of a paramagnetic protein are available. Differences between them and the X-ray structure are discussed.

Feed-forward neural networks for secondary structure prediction

A feed-forward neural network has been employed for protein secondary structure prediction. Attempts were made to improve on previous prediction accuracies using a hierarchical mixture of experts (HME). In this method input data are clustered and used to train a series of different networks. Application of an HME to the prediction of protein secondary structure is shown to provide no advantages over a single network. We have also tried various new input representations, chosen to incorporate the effect of residues a long distance away in the one-dimensional amino acid chain. Prediction accuracy using these methods is comparable to that achieved by other neural networks.

15N resonance assignments of oxidized and reduced Chromatium vinosum high-potential iron protein

The 15N resonances in reduced and oxidized Chromatium vinosum high-potential iron protein have been assigned by use of 1H-1H COSY spectra and 1H-15N HMQC. HMQC-COSY, and HMQC-NOESY spectra. Unambiguous assignment of 70 of 85 backbone 15N resonances in the reduced protein and 62 of 85 resonances in the oxidized protein are made, as are 12 of 21 side-chain 15N resonances.

Sequence-specific assignment of the 1H and 15N nuclear magnetic resonance spectra of the reduced recombinant high-potential iron-sulfur protein I from Ectothiorhodospira halophila

A 1H and 15N NMR investigation through two-dimensional and three-dimensional spectroscopy has been performed on the reduced form ([Fe4S4]2+) of the recombinant high-potential iron-sulfur protein (HiPIP) I from Ectothiorhodospira halophila expressed in Escherichia coli. [Fe4S4]2+ clusters in proteins are paramagnetic with a relatively low mu eff of about 0.8 mu B/iron ion, but the paramagnetic effects on nuclear relaxation are so strong as to yield T1 values of a few milliseconds and linewidths of hundreds of hertz for the nuclei closet to the paramagnetic center. Despite these features, 71 out of 73 residues were identified, most of which were assigned completely as far as proton resonances are concerned; as many as 68 residues could be assigned without any reference to the existing X-ray structure. A total of 88% of all protein protons and 58 out of 69 peptide HN nitrogen signals were assigned. To the best of our knowledge, this is the most extensive 1H assignment of a paramagnetic protein to date. Protons sensitive to the proximity of the cluster were assigned through suitable NOE spectroscopy experiments. Three out of the four coordinated cysteines were assigned, and two residues have been identified whose peptide HN protons give rise to H bonds with coordinated sulfur atoms. The inter-residue NOE cross peaks are in qualitative agreement with the secondary and tertiary structure as obtained from the available X-ray crystallographic analysis of the wild-type protein at 250-pm resolution. It is therefore shown that the expressed protein is properly folded and that it is a reliable model for the wild-type protein. These data are meaningful for the detection of structural differences among mutants in future studies.

The three-dimensional structure in solution of the paramagnetic high-potential iron-sulfur protein I from Ectothiorhodospira halophila through nuclear magnetic resonance

The three-dimensional structure in solution of reduced recombinant high-potential iron-sulfur protein iso-I from Ectothiorhodospira halophila was determined using 948 relevant interproton NOEs out of the 1246 observed NOEs. The determination was accomplished using the XEASY program for spectral analysis and the distance geometry (DG) program DIANA for generation of the structure as described by Wüthrich [Wüthrich, K. (1989) Acc. Chem. Res. 22, 36-44]. The FeS cluster was simulated using an amino acid residue constructed for the present work from a cysteinyl residue with an iron and a sulfur atom attached to the terminal thiol. The family of structures obtained from distance geometry were subjected to energy minimization and molecular dynamics simulations using previously defined force field parameters. The quality of these structures at each stage of the refinement process is discussed with respect to the dihedral angle order parameter and the root-mean-square deviation of the atomic coordinates. The latter values for the backbone atoms vary from 67 pm for the distance-geometry structures to 60 pm for the energy-minimized structures to 51 pm for the structures subjected to restrained molecular dynamics. Finally, the structure in best agreement with the NOE constraints has been further treated with extensive restrained molecular dynamics in water. The solution structure is well defined and is very similar to the available X-ray structure. We do not know of any previous determination of the structure of a paramagnetic protein in solution by NMR. The effect of paramagnetism on the quality of the structure determination is discussed.

Calculation of the redox potentials of iron-sulfur proteins: the 2-/3-couple of [Fe4S*4Cys4] clusters in Peptococcus aerogenes ferredoxin, Azotobacter vinelandii ferredoxin I, and Chromatium vinosum high-potential iron protein

Calculations of the redox potentials of the 2-/3-couples of [Fe4S*4Cys4] clusters in the iron-sulfur proteins Peptococcus aerogenes ferredoxin (PaFd), Azotobacter vinelandii ferredoxin I (AvFdI) and Chromatium vinosum high potential iron protein (CvHiPIP) based on the Protein Dipoles Langevin Dipoles (PDLD) method are reported. The structures of these proteins have been determined by X-ray crystallography; in the case of PaFd the structure has recently been revised due to a change in the sequence close to Cluster II. The large differences between the potentials of the [Fe4S*4Cys4] clusters of PaFd and AvFdI and the potential of the [Fe4S*4Cys4] cluster of CvHiPIP are successfully modeled and originate principally in differences in the configuration of main-chain amide groups near the clusters. The small difference between the potentials of PaFd and AvFdI is also satisfactorily modeled in the case of Cluster I of PaFd. Solvent dipoles close to the cluster in PaFd are an important contributor to its higher potential. The two X-ray structures of PaFd yield similar results for Cluster I of PaFd. In contrast, the results for Cluster II differ substantially; for reasons not yet clear, the recently revised structure leads to results in worse agreement with experiment.

Solution 1H NMR determination of secondary structure for the three-iron form of ferredoxin from the hyperthermophilic archaeon Pyrococcus furiosus

Two-dimensional 1H NMR data have been used to make sequence-specific assignments and define the secondary structure of the three-iron form of the oxidized ferredoxin, Fd, from the hyperthermophilic archaeon Pyrococcus furiosus, Pf. Signals for at least some protons were located for 65 of the 66 amino acids in the sequence, in spite of the paramagnetic (S = 1/2) ground state, but not all could be assigned. Unassigned and missing signals could be qualitatively correlated with the expected proximity of the protons to the paramagnetic cluster. The secondary structure was deduced from qualitative analysis of the 2D nuclear Overhauser effect, which identified two antiparallel beta-sheets, one triple-stranded including Ala1-Ser5, Val39-Glu41, and Thr62-Ala66, and one double-stranded consisting of Glu26-Asn28 and Lys32-Glu34, as well as an alpha-helix involving Glu43-Glu54. Three tight type I turns are located at residues Asp7-Thr10, Pro22-Phe25, and Asp29-Gly31. Comparison with the crystal structure of Desulfovibrio gigas, Dg, Fd (Kissinger et al., 1991) reveals a very similar folding topology, although several secondary structural elements are extended in Pf relative to Dg Fd. Thus the beta-sheet involving the two termini is expanded to include the two terminal residues and incorporates a third strand from the internal loop that is lengthened by several insertions in Pf relative to Dg Fd. The double-stranded beta-sheet in the interior of Pf Fd is lengthened slightly due to a much tighter type I turn between the two strands. The helix near the C-terminus is three residues longer in Pf than in Dg Fd, as well as being shifted toward the N-terminus. The disulfide link between the two nonligating Cys residues (Cys21 and Cys48) is conserved in Pf Fd, but the link near the C-terminus is in the middle of the long alpha-helix in Pf Fd, instead of at the N-terminus of the helix as in Dg Fd. The extensions of the beta-sheets and alpha-helix increase the number of main-chain hydrogen bonds in Pf Fd by approximately 8 relative to those in Dg Fd and likely contribute to its remarkable thermostability (it is unaffected by anaerobic incubation at 95 degrees C for 24 h).(ABSTRACT TRUNCATED AT 400 WORDS)

Effect of replacing conserved proline residues on the EPR and NMR properties of Clostridium pasteurianum 2[4Fe-4S] ferredoxin

Most of [4Fe-4S] proteins bind their metallic center by four cysteine residues, three clustered in a single stretch of seven amino acids and a remote fourth generally followed by a proline residue. Two such prolines in Clostridium pasteurianum 2[4Fe-4S] ferredoxin have been substituted by different amino acids and the resulting molecular variants studied with EPR and NMR spectroscopies. The isolated EPR contributions of the [4Fe-4S]+ clusters do not change much in all variants. The exact positions or the number of features composing the fully reduced EPR spectra built by the two interacting [4Fe-4S]+ S = 1/2 systems vary slightly but, in none of the proteins in which either proline 19 or 48 were substituted, do they indicate a major difference either in the folding of the ferredoxin or in the electronic structure of its clusters. A subset of paramagnetically shifted NMR signals is significantly affected by these replacements at both redox levels. The corresponding protons belong to two cysteines liganding the cluster close to the substitution. These data, combined with the presently available three-dimensional information, form the basis for partial assignments of the most shifted resonances in the NMR spectra of such proteins. The positions of intermediate lines in the NMR spectra of semireduced ferredoxins depend on the difference between the redox potentials of the two clusters; this difference is sensitive to the substitutions of either conserved proline residue by lysine.

Improved calculations of compactness and a reevaluation of continuous compact units

A new method for calculating compactness (Z) that uses look-up table-based algorithms for area and volume computations is introduced. These algorithms can be used in any iterative area or volume calculation, are more than 1000 times faster than the original algorithms, and have equal or better precision. With the faster algorithms it is now possible to calculate the compactness of all continuous units in a protein, and to precisely locate the optimal compact units without the screening functions and limited resolution used previously. These methods have been incorporated into a fully automatic domain finding algorithm, and this method has been applied to the 21 proteins originally analyzed as well as 12 additional proteins. This method is robust, and yields similar units even when applied to coordinates of protein crystals grown under different experimental conditions.

1H-NMR investigation of oxidized and reduced high-potential iron-sulfur protein from Rhodopseudomonas globiformis

1H one-dimensional and two-dimensional NMR spectra have been recorded for the oxidized and reduced forms of the high-potential iron-sulfur protein (HiPIP) from Rhodopseudomonas globiformis which has the highest known reduction potential. The spectrum of the oxidized protein is similar to that of Chromatium vinosum and Rhodocyclus gelatinosus HiPIP but different from that of the HiPIP II from Ectothiorhodospira halophila. Surprisingly, site-specific assignment has shown that in the oxidized protein the distribution of oxidation numbers within the cluster is very similar to that found for E. halophila HiPIP II and different from that of the other two proteins. The spectrum of the reduced species is very similar to that of all other HiPIPs known to date, indicating very similar electronic and geometric structures for the reduced forms. These findings are discussed in terms of cluster structure in HiPIPs and of redox potentials.

Three-dimensional structure of the high-potential iron-sulfur protein isolated from the purple phototrophic bacterium Rhodocyclus tenuis determined and refined at 1.5 A resolution

The molecular structure of the high-potential iron-sulfur protein (HiPIP) isolated from the phototrophic bacterium, Rhodocyclus tenuis, has been solved and refined to a nominal resolution of 1.5 A with a crystallographic R-factor of 17.3% for all measured X-ray data from 30 A to 1.5 A. It is the smallest of the HiPIP structures studied thus far with 62 amino acid residues. Crystals used in the investigation belonged to the space group P2(1) with unit cell dimensions of a = 36.7 A, b = 52.6 A, c = 27.6 A and beta = 90.8 degrees and contained two molecules per asymmetric unit. The structure was solved by a combination of multiple isomorphous replacement with two heavy-atom derivatives, anomalous scattering from the iron-sulfur cluster, symmetry averaging and solvent flattening. The folding motif for this HiPIP is characterized by one small alpha-helix, six Type I turns, an approximate Type II turn and one Type I' turn. As in other HiPIPs, the iron-sulfur cluster is co-ordinated by four cysteinyl ligands and exhibits a cubane-like motif. These cysteinyl ligands are all located in Type I turns. The hydrogen bonding around the metal cluster in the R. tenuis protein is similar to the patterns observed in the Chromatium vinosum and Ectothiorhodospira halophila HiPIPs. Several of the amino acid residues invariant in the previously determined C. vinosum and E. halophila structures are not retained in the R. tenuis molecule. There are 13 solvent molecules structurally conserved between the two R. tenuis HiPIP molecules in the asymmetric unit, some of which are important for stabilizing surface loops. Interestingly, while it is assumed that this HiPIP functions as a monomer in solution, the two molecules in the asymmetric unit pack as a dimer and are related to each other by an approximate twofold rotation axis.

Accurate modeling of protein conformation by automatic segment matching

Segment match modeling uses a data base of highly refined known protein X-ray structures to build an unknown target structure from its amino acid sequence and the atomic coordinates of a few of its atoms (generally only the C alpha atoms). The target structure is first broken into a set of short segments. The data base is then searched for matching segments, which are fitted onto the framework of the target structure. Three criteria are used for choosing a matching data base segment: amino acid sequence similarity, conformational similarity (atomic co-ordinates), and compatibility with the target structure (van der Waals' interactions). The new method works surprisingly well: for eight test proteins ranging in size from 46 to 323 residues, the all-atom root-mean-square deviation of the modeled structures is between 0.93 A and 1.73 A (the average is 1.26 A). Deviations of this magnitude are comparable with those found for protein co-ordinates before and after refinement against X-ray data or for co-ordinates of the same protein in different crystal packings. These results are insensitive to errors in the C alpha positions or to missing C alpha atoms: accurate models can be built with C alpha errors of up to 1 A or by using only half the C alpha atoms. The fit to the X-ray structures is improved significantly by building several independent models based on different random choices and then averaging co-ordinates; this novel concept has general implications for other modeling tasks. The segment match modeling method is fully automatic, yields a complete set of atomic co-ordinates without any human intervention and is efficient (14 s/residue on the Silicon Graphics 4D/25 Personal Iris workstation.

Sequence-specific assignments of the 1H nuclear magnetic resonance spectra of reduced high-potential ferredoxin (HiPIP) from Chromatium vinosum

The 1H resonances of the high-potential [4Fe-4S]2+ ferredoxin from Chromatium vinosum have been assigned through conventional sequential methodology applied to 2D NMR spectra. Almost 80% of the residues were identified using standard 2D COSY, HOHAHA, and NOESY pulse sequences. These residues correspond to four segments of the primary structure that do not interact strongly with the iron-sulfur cluster. A minor correction to the amino acid sequence is strongly suggested by these NMR data. Additional protons more sensitive to the proximity of the cluster were assigned by a combination of NOESY experiments with fast repetition rates and short mixing times and of HOHAHA spectra recorded with reduced spin-lock duration aimed at compensating for the short relaxation rates. Hence, the contributions of 79 residues out of 85 were identified in NMR spectra, among which the assignments of 64 residues were completed. Even the fastest relaxing protons, like those of the cysteine ligands, could be correlated, partly because the strong hyperfine shifts isolate them from the crowded diamagnetic region. However, other protons, in particular those involved in NH-S hydrogen bonds with the iron-sulfur cluster, were more difficult to identify, most probably because their relatively broad signals overlapped with those of protons not or less perturbed by the active site. The availability of the major part of the 1H NMR assignments has enabled the detection and identification of many interresidue NOESY cross peaks. These data are in full agreement with the elements of secondary structure previously revealed by X-ray crystallographic analysis of the protein.(ABSTRACT TRUNCATED AT 250 WORDS)

Evaluation of protein models by atomic solvation preference

Important properties of globular proteins, such as the stability of the folded state, depend sensitively on interactions with solvent molecules. An excluded volume approximation to protein-solvent interaction, the solvent contact model, was used to derive atomic solvation preference parameters from a database of known protein structures. The ability of solvation preference to discriminate between correct and incorrect three-dimensional structures for a given sequence, or to identify the correct sequence placement in a given structure, was tested. Backbone co-ordinates were taken from experimentally known structures or hypothetical models and side-chain conformations (in rotamer space) were optimized by an efficient Monte Carlo algorithm using simulated annealing and simple potential functions. Discrimination by solvation preference was very clear between deliberately misfolded and correct globular models as well as between native-like and non-native-like topologies of combinatorially generated myoglobin models. Due to its statistical nature, the evaluation works best on entire protein models, while the identification of incorrect parts of models is more difficult. In one case locally incorrect chain tracing in a crystal structure was identified. The method is computationally fast compared to methods based on surface area calculations and is recommended for use as a diagnostic tool in model building based on sequence similarity, in folding simulations and in protein design.

Sequential resonance assignments of oxidized high-potential iron-sulfur protein from Chromatium vinosum

2D NMR spectra of the high-potential iron-sulfur protein (HiPIP) from Chromatium vinosum have been used to obtain partial resonance assignments for the oxidized paramagnetic redox state of the protein. Sequence-specific assignments were made using NOESY and COSY spectra in H2O and D2O of the following backbone segments: Asn-5-Arg-33, Glu-39-Asp-45, Gly-55-Cys-63, Gly-68-Ala-78, and Leu-82-Gly-85. NOESY spectra with a spectral width wide enough to include the hyperfine-shifted resonances revealed numerous NOE contacts between these signals and those in the main envelope of the proton spectrum. With the aid of the X-ray crystal structure [Carter, C.W., Kraut, J., Freer, S. T., Xuong, N. H., Alden, R. A., & Bartsch, R. G. (1974) J. Biol. Chem. 249, 4212], these NOEs permitted seven of the nine hyperfine-shifted signals to be assigned to three of the cysteine residues liganded to the metal cluster (Cys-43, Cys-46, and Cys-77). The other two hyperfine-shifted signals produced no detectable NOEs to other resonances in the spectrum and were tentatively assigned to the remaining cysteinyl ligand (Cys-63). These assignments, in conjunction with recent theoretical models of the electronic structure of the Fe4S4 cluster [Noodleman, L. (1988) Inorg. Chem. 27, 3677; Bertini, I., Briganti, F., Luchinat, C., Scozzafava, A., & Sola, M. (1991) J. Am. Chem. Soc. 113, 1237], indicate that the iron atoms coordinated to Cys-63 and Cys-77 are those of the mixed-valence Fe(3+)-Fe2+ pair whereas Cys-43 and Cys-46 are ligands to the Fe(3+)-Fe3+ metal pair.

Generalized protein tertiary structure recognition using associative memory Hamiltonians

In previous papers, a method of protein tertiary structure recognition was described based on the construction of an associative memory Hamiltonian, which encoded the amino acid sequence and the C alpha co-ordinates of a set of database proteins. Using molecular dynamics with simulated annealing, the ability of the Hamiltonian to successfully recall the structure of a protein in the memory database was successfully demonstrated, as long as the total number of database proteins did not exceed a characteristic value, called the capacity of the Hamiltonian, equal to 0.5N to 0.7N, where N is the number of amino acid residues in the protein to be recalled. In this paper, we describe the development of additional methods to increase the capacity of the Hamiltonian, including use of a more complete representation of the protein backbone and the incorporation of contextual information into the Hamiltonian through the use of secondary structure prediction. In addition, we further extend the ability of associative memory models to predict the tertiary structures of proteins not present in the protein data set, by making the Hamiltonian invariant with respect to biological symmetries that represent site mutations and insertions and deletions. The ability of the Hamiltonian to generalize from homologous proteins to an unknown protein in the presence of other unrelated proteins in the data set is demonstrated.

An investigation of Chromatium vinosum high-potential iron-sulfur protein by EPR and Mossbauer spectroscopy; evidence for a freezing-induced dimerization in NaCl solutions

The high-potential iron-sulfur protein (HiPIP) from Chromatium vinosum contains a cubane prosthetic group that shuttles between the [4Fe-4S]3+,2+ states. We find that the EPR spectra from this protein can be explained as a sum of two components, a major one with g = 2.02; 2.04; 2.12, and a minor one with g = 2.04; 2.07; approximately 2.13. In the presence of 0.1-2.0 M NaCl, freezing induces polymerization of the protein (presumably dimers), which is detected as intercluster spin-spin interaction in the EPR. The observed spin-spin interactions are interpreted as being due to two very similar dimeric structures in an approx. 1:2 ratio. Computer simulation of the X- and Q-band EPR spectra shows that the z-components of the g-tensors in each dimer pair must be co-linear, with center-to-center distances between the clusters of approximately 13 A and approximately 16 A. Inspection of possible dimeric structures of C. vinosum HiPIP by standard molecular graphics procedures revealed that the Fe/S cluster is exposed toward a flattened surface and is accessible to solvent. Moreover, the Fe/S clusters in two HiPIP molecules can easily achieve a center-to-center distance of approximately 14 A when approaching along a common 3-fold axis that extends through the S4 sulfur atom of the cubane; the z-component of the EPR g-tensor is co-linear with this symmetry axis.

Modeling of globular proteins. A distance-based data search procedure for the construction of insertion/deletion regions and Pro----non-Pro mutations

A distance-based database search scheme is proposed for modeling Pro----in non-Pro and insertion/deletion regions of homologous globular proteins up to six residues in length. In the first step, geometric descriptors, the number of residues involved and target distances corresponding to the separation of C alpha atom positions adjacent to the "missing" segment, are chosen. In the second step, a database of high-resolution X-ray structures is scanned for segments with similar descriptors and selected segments are binned according to conformational type. In the third and fourth steps, the selected conformations are docked into the protein, and geometric and energetic criteria are used to determine their viability as segment models. The fifth step consists of an interaction scheme in which the geometric descriptors are redefined. This compensates for the use of a limited database and/or for the use of a poor original protein model adjacent to the missing segment. The procedure has been tested on Pro----non-Pro mutations in the homologous proteins penicillopepsin and endothiapepsin, and on the insertion/deletion regions of the homologs penicillopepsin and endothiapepsin, trypsin and gamma-chymotrypsin and hen and human lysozyme. The test cases represent a wide variety of secondary structural elements (helix, sheet, turn and coil) and insertion/deletion lengths (0 to 4 residues). It is shown that 79% of the test cases are accurately modeled (within 0.54 A root-mean-square (r.m.s.) deviation for main-chain atoms) using the proposed scheme. Failure of the scheme (main-chain atom r.m.s. deviations greater than 1.29 A) in 21% of the cases appears to be related to the presence of infrequently observed conformations or locally unique folds of the target proteins with respect to the database (18% of the test cases); the remaining 3% are unexplained. Geometric and energetic criteria are able to discriminate between trial conformations that correspond to the X-ray structures and those that are different in 97% of the conformations generated by the distance-weighted database search scheme. The scheme is shown to be relatively insensitive to uncertainty in the template co-ordinates, since the geometric descriptors were taken from the homologous protein (r.m.s. deviations in the position of descriptors range from 0.18 to 1.35 A for the accurately modeled test cases). It is demonstrated that the scheme can be used to correct local sequence misalignments.

Identification of native protein folds amongst a large number of incorrect models. The calculation of low energy conformations from potentials of mean force

We present an approach that is able to detect native folds amongst a large number of non-native conformations. The method is based on the compilation of potentials of mean force of the interactions of the C beta atoms of all amino acid pairs from a database of known three-dimensional protein structures. These potentials are used to calculate the conformational energy of amino acid sequences in a number of different folds. For a substantial number of proteins we find that the conformational energy of the native state is lowest amongst the alternatives. Exceptions are proteins containing large prosthetic groups, Fe-S clusters or polypeptide chains that do not adopt globular folds. We discuss briefly potential applications in various fields of protein structural research.

Cloning and sequencing of the locus encoding the large and small subunit genes of the periplasmic [NiFe]hydrogenase from Desulfovibrio fructosovorans

The genetic locus encoding the periplasmic [NiFe]hydrogenase (Hyd) from Desulfovibrio fructosovorans was cloned and sequenced. The genes of this two-subunit enzyme have an operon organization in which the 0.94-kb gene encoding the small subunit precedes the 1.69-kb gene encoding the large subunit. A Shine-Dalgarno sequence is centered at -9 bp from the ATG of both subunits. The possible presence of another open reading frame downstream from the large-subunit-encoding gene is considered. The N-terminal sequence of the large 61-kDa subunit deduced from the nucleotide sequence is in perfect agreement with the results of the amino acid (aa) sequence determined by Edman degradation. A 50-aa leader peptide precedes the small 28-kDa subunit. The aa sequence of the enzyme shows nearly 65% homology with the [NiFe]Hyd aa sequence of Desulfovibrio gigas. Comparisons with a large range of Hyds from various bacterial species indicate the presence of highly conserved Cys residues, the implications of which are discussed from the point of view of nickel atom and cluster accommodation.

Improvements in protein secondary structure prediction by an enhanced neural network

Computational neural networks have recently been used to predict the mapping between protein sequence and secondary structure. They have proven adequate for determining the first-order dependence between these two sets, but have, until now, been unable to garner higher-order information that helps determine secondary structure. By adding neural network units that detect periodicities in the input sequence, we have modestly increased the secondary structure prediction accuracy. The use of tertiary structural class causes a marked increase in accuracy. The best case prediction was 79% for the class of all-alpha proteins. A scheme for employing neural networks to validate and refine structural hypotheses is proposed. The operational difficulties of applying a learning algorithm to a dataset where sequence heterogeneity is under-represented and where local and global effects are inadequately partitioned are discussed.

1H NMR studies of oxidized high-potential iron protein from Chromatium vinosum. Nuclear Overhauser effect measurements

1H nuclear Overhauser effect experiments on the isotropically shifted signals of oxidized Chromatium vinosum HiPIP have been used to identify the four beta-CH2 geminal couples of the cysteine ligands. A partial assignment to individual residues has been proposed from a computer graphics analysis of the X-ray structure. Tentative assignments of other resonances are discussed.

Use of techniques derived from graph theory to compare secondary structure motifs in proteins

A substructure matching algorithm is described that can be used for the automatic identification of secondary structural motifs in three-dimensional protein structures from the Protein Data Bank. The proteins and motifs are stored for searching as labelled graphs, with the nodes of a graph corresponding to linear representations of helices and strands and the edges to the inter-line angles and distances. A modification of Ullman's subgraph isomorphism algorithm is described that can be used to search these graph representations. Tests with patterns from the protein structure literature demonstrate both the efficiency and the effectiveness of the search procedure, which has been implemented in FORTRAN 77 on a MicroVAX-II system, coupled to the molecular fitting program FRODO on an Evans and Sutherland PS300 graphics system.

An investigation of oligopeptides linking domains in protein tertiary structures and possible candidates for general gene fusion

Fifty-one examples of oligopeptides linking protein domains were extracted from the Brookhaven database of three-dimensional protein structures. In general, the peptides displayed specific characteristics in composition, conformation, hydrogen bonding, flexibility and the like. The entire database was then searched for pentapeptides that would optimize these natural linker properties. The oligopeptides found are suggested as general candidates to link protein molecules or domains through gene fusion.

Conformation of beta-hairpins in protein structures. A systematic classification with applications to modelling by homology, electron density fitting and protein engineering

A systematic classification of beta-hairpin structures which takes into account the polypeptide chain length and hydrogen bonding between the two antiparallel beta-strands is described. We have used this classification of beta-hairpin structures and their specific sequence pattern to derive rules which demonstrate its usefulness in assisting modelling beta-hairpins. These rules can be applied to comparative model building, modelling into electron density and in the prediction of conformation of beta-hairpins to aid protein engineering.

1H-NMR studies of high-potential iron-sulfur protein from the purple phototrophic bacterium, Rhodospirillum tenue

The high-potential iron-sulfur protein (HiPIP) from Rhodospirillum tenue (strain 3761) shows only a weak (20-25%) sequence similarity to HiPIPs from Chromatium vinosum, Ectothiorhodospira halophila and Ectothiorhodospira vacuolata, including the strict conservation of only two of the twelve residues assumed to be in the 4Fe-4S cluster packing region [Tedro, S. M., Meyer, T. E. and Kamen, M. D. (1979) J. Biol. Chem. 254, 1495-1500]. In spite of these differences, the general range and distribution of hyperfine-shifted 1H-NMR peaks of oxidized and reduced R. tenue HiPIP resemble those of E. halophila HiPIP I [Krishnamoorthi, R., Markley, J. L., Cusanovich, M. A., Pryzycieki, C. T. and Meyer, T. E. (1986) Biochemistry 25, 60-67]. Temperature- and pH-dependence and longitudinal relaxation behavior were determined for hyperfine-shifted peaks of the oxidized protein. Tentative assignments of peaks to ligands and aromatic residues suggest the presence of common apoprotein-active-site interactions in these proteins. Differences occur in the pattern of paramagnetically shifted peaks attributed to hydrogens bonded to the 4Fe-4S cluster. Hyperfine-shifted peaks of R. tenue HiPIP are not perturbed by pH changes in the range 5-9. In contrast, those of the C. vinosum protein exhibit a pH-dependence of chemical shifts that has been attributed to the titration of His42 [Nettesheim, D. G., Meyer, T. E., Feinberg, B. A. and Otvos, J. D. (1983) J. Biol. Chem. 258, 8235-8239]. Since R. tenue HiPIP contains no histidine, the present observation confirms the above hypothesis.