Crystal structure of the 2[4Fe-4S] ferredoxin from Chromatium vinosum: evolutionary and mechanistic inferences for [3/4Fe-4S] ferredoxins.

The crystal structure of the 2[4Fe-4S] ferredoxin from Chromatium vinosum has been solved by molecular replacement using data recorded with synchrotron radiation. The crystals were hexagonal prisms that showed a strong tendency to develop into long tubes. The hexagonal prisms diffracted to 2.1 A resolution at best, and a structural model for C. vinosum ferredoxin has been built with a final R of 19.2%. The N-terminal domain coordinates the two [4Fe-4S] clusters in a fold that is almost identical to that of other known ferredoxins. However, the structure has two unique features. One is a six-residue insertion between two ligands of one cluster forming a two-turn external loop; this short loop changes the conformation of the Cys 40 ligand compared to other ferredoxins and hampers the building of one NH...S H-bond to one of the inorganic sulfurs. The other remarkable structural element is a 3.5-turn alpha-helix at the C-terminus that covers one side of the same cluster and is linked to the cluster-binding domain by a six-residue external chain segment. The charge distribution is highly asymmetric over the molecule. The structure of C. vinosum ferredoxin strongly suggests divergent evolution for bacterial [3/4Fe-4S] ferredoxins from a common ancestral cluster-binding core. The unexpected slow intramolecular electron transfer rate between the clusters in C. vinosum ferredoxin, compared to other similar proteins, may be attributed to the unusual electronic properties of one of the clusters arising from localized changes in its vicinity rather than to a global structural rearrangement.     

普通烟草LBD基因家族的全基因组序列鉴定与表达分析

LBD是一类具有LOB(lateral organ boundaries)结构域的基因家族,在植物发育过程中起到非常重要的作用。采用生物信息学方法,根据拟南芥LBD基因序列鉴定了普通烟草基因组中的LBD基因,并对家族成员进行了序列特征、系统发育和表达谱分析。结果表明:普通烟草基因组中共有98个LBD基因成员,其基因结构相对简单,一般含有1~3个外显子。LBD基因家族可分成I和II两大类,两类均含有CX_2CX_6CX_3C保守结构域,但II类不含有LX_6LX_3LX_6L形成的"卷曲螺旋"二级结构,根据与拟南芥LBD蛋白构建的系统发育树则可细分成5个亚家族(Ia、Ib、Ic、Id和II)。将LBD基因与表达序列标签(EST)比对,发现36个基因有EST证据;EST、芯片数据和转录组数据分析表明:LBD基因具有不同的组织表达模式,部分基因表现出组织特异性。这些研究结果为普通烟草LBD基因家族功能的深入研究奠定了基础。     

Database of homology-derived protein structures and the structural meaning of sequence alignment

The database of known protein three-dimensional structures can be significantly increased by the use of sequence homology, based on the following observations. (1) The database of known sequences, currently at more than 12,000 proteins, is two orders of magnitude larger than the database of known structures. (2) The currently most powerful method of predicting protein structures is model building by homology. (3) Structural homology can be inferred from the level of sequence similarity. (4) The threshold of sequence similarity sufficient for structural homology depends strongly on the length of the alignment. Here, we first quantify the relation between sequence similarity, structure similarity, and alignment length by an exhaustive survey of alignments between proteins of known structure and report a homology threshold curve as a function of alignment length. We then produce a database of homology-derived secondary structure of proteins (HSSP) by aligning to each protein of known structure all sequences deemed homologous on the basis of the threshold curve. For each known protein structure, the derived database contains the aligned sequences, secondary structure, sequence variability, and sequence profile. Tertiary structures of the aligned sequences are implied, but not modeled explicitly. The database effectively increases the number of known protein structures by a factor of five to more than 1800. The results may be useful in assessing the structural significance of matches in sequence database searches, in deriving preferences and patterns for structure prediction, in elucidating the structural role of conserved residues, and in modeling three-dimensional detail by homology.     

Trematode and Monogenean rRNA ITS2 Secondary Structures Support a Four-Domain Model

The secondary structure of rRNA internal transcribed spacer 2 is important in the process of ribosomal biogenesis. Trematode ITS sequences are poorly conserved and difficult to align for phylogenetic comparisons above a family level. If a conserved secondary structure can be identified, it can be used to guide primary sequence alignments. ITS2 sequences from 39 species were compared. These species span four orders of trematodes (Echinostomiformes, Plagiorchiformes, Strigeiformes, and Paramphistomiformes) and one monogenean (Gyrodactyliformes). The sequences vary in length from 251 to 431 bases, with an average GC content of 48%. The monogenean sequence could not be aligned with confidence to the trematodes. Above the family level trematode sequences were alignable from the 5′ end for 139 bases. Secondary structure foldings predicted a four-domain model. Three folding patterns were required for the apex of domain B. The folding pattern of domains C and D varies for each family. The structures display a high GC content within stems. Bases A and U are favored in unpaired regions and variable sites cluster. This produces a mosaic of conserved and variable regions with a structural conformation resistant to change. Two conserved strings were identified, one in domain B and the other in domain C. The first site can be aligned to a processing site identified in yeast and rat. The second site has been found in plants, and structural location appears to be important. A phylogenetic tree of the trematode sequences, aligned with the aid of secondary structures, distinguishes the four recognized orders. Received: 21 November 1997 / Accepted: 9 February 1998     

Sequence patterns derived from the automated prediction of functional residues in structurally-aligned homologous protein families

MOTIVATION: Most proteins have evolved to perform specific functions that are dependent on the adoption of well-defined three-dimensional (3D) structures. Specific patterns of conserved residues in amino acid sequences of divergently evolved proteins are frequently observed; these may reflect evolutionary restraints arising both from the need to maintain tertiary structure and the requirement to conserve residues more directly involved in function. Databases of such sequence patterns are valuable in identifying distant homologues, in predicting function and in the study of evolution. RESULTS: A fully automated database of protein sequence patterns, Functional Protein Sequence Pattern Database (FPSPD), has been derived from the analysis of the conserved residues that are predicted to be functional in structurally aligned homologous families in the HOMSTRAD database. Environment-dependent substitution tables, evolutionary trace analysis, solvent accessibility calculations and 3D-structures were used to obtain the FPSPD. The method yielded 3584 patterns that are considered functional and 3049 patterns that are probably functional. FPSPD could be useful for assigning a protein to a homologous superfamily and thereby providing clues about function. AVAILABILITY: FPSPD is available at http://www-cryst.bioc.cam.ac.uk/~fpspd/     

Conservation analysis and structure prediction of the SH2 family of phosphotyrosine binding domains.

Src homology 2 (SH2) regions are short (approximately 100 amino acids), non-catalytic domains conserved among a wide variety of proteins involved in cytoplasmic signaling induced by growth factors. It is thought that SH2 domains play an important role in the intracellular response to growth factor stimulation by binding to phosphotyrosine containing proteins. In this paper we apply the techniques of multiple sequence alignment, secondary structure prediction and conservation analysis to 67 SH2 domain amino acid sequences. This combined approach predicts seven core secondary structure regions with the pattern beta-alpha-beta-beta-beta-beta-alpha, identifies those residues most likely to be buried in the hydrophobic core of the native SH2 domain, and highlights patterns of conservation indicative of secondary structural elements. Residues likely to be involved in phosphotyrosine binding are shown and orientations of the predicted secondary structures suggested which could enable such residues to cooperate in phosphate binding. We propose a consensus pattern that encapsulates the principal conserved features of the SH2 domains. Comparison of the proposed SH2 domain of akt to this pattern shows only 12/40 matches, suggesting that this domain may not exhibit SH2-like properties.     

Relationship between sequence conservation and three-dimensional structure in a large family of esterases, lipases, and related proteins.

Based on the recently determined X-ray structures of Torpedo californica acetylcholinesterase and Geotrichum candidum lipase and on their three-dimensional superposition, an improved alignment of a collection of 32 related amino acid sequences of other esterases, lipases, and related proteins was obtained. On the basis of this alignment, 24 residues are found to be invariant in 29 sequences of hydrolytic enzymes, and an additional 49 are well conserved. The conservation in the three remaining sequences is somewhat lower. The conserved residues include the active site, disulfide bridges, salt bridges, and residues in the core of the proteins. Most invariant residues are located at the edges of secondary structural elements. A clear structural basis for the preservation of many of these residues can be determined from comparison of the two X-ray structures.     

The identification of conserved interactions within the SH3 domain by alignment of sequences and structures

The SH3 domain, comprised of approximately 60 residues, is found within a wide variety of proteins, and is a mediator of protein-protein interactions. Due to the large number of SH3 domain sequences and structures in the databases, this domain provides one of the best available systems for the examination of sequence and structural conservation within a protein family. In this study, a large and diverse alignment of SH3 domain sequences was constructed, and the pattern of conservation within this alignment was compared to conserved structural features, as deduced from analysis of eighteen different SH3 domain structures. Seventeen SH3 domain structures solved in the presence of bound peptide were also examined to identify positions that are consistently most important in mediating the peptide-binding function of this domain. Although residues at the two most conserved positions in the alignment are directly involved in peptide binding, residues at most other conserved positions play structural roles, such as stabilizing turns or comprising the hydrophobic core. Surprisingly, several highly conserved side-chain to main-chain hydrogen bonds were observed in the functionally crucial RT-Src loop between residues with little direct involvement in peptide binding. These hydrogen bonds may be important for maintaining this region in the precise conformation necessary for specific peptide recognition. In addition, a previously unrecognized yet highly conserved beta-bulge was identified in the second beta-strand of the domain, which appears to provide a necessary kink in this strand, allowing it to hydrogen bond to both sheets comprising the fold.     

Amino acid sequences of Nostoc strain MAC ferredoxins I and II

The amino acid sequences of ferredoxins I and II from a blue-green alga, Nostoc strain MAC were determined. This alga is able to grow autotrophically in the light or heterotrophically in the dark. Analyses of tryptic peptides of Cm-proteins by conventional methods including solid-phase Edman degradation gave the complete amino acid sequences. Both molecules consisted of 98 amino acid residues and 34 amino acid differences including two deletions were found between the two. Comparing these sequences with those of ferredoxins from Chlorogloeopsis fritschii and Synechocystis 6714, which are also capable of growing under both conditions, showed that Nostoc strain MAC ferredoxin II had unique amino acids around the [2Fe-2S] cluster. This finding provides a structural basis for explaining the different chemical and functional properties of Nostoc strain MAC ferredoxin II reported in a previous paper (Hutson et al. (1978) Biochem. J. 172, 465-477).     

Identification and mutational analysis of Arabidopsis FLS2 leucine-rich repeat domain residues that contribute to flagellin perception

Mutational, phylogenetic, and structural modeling approaches were combined to develop a general method to study leucine-rich repeat (LRR) domains and were used to identify residues within the Arabidopsis thaliana FLAGELLIN-SENSING2 (FLS2) LRR that contribute to flagellin perception. FLS2 is a transmembrane receptor kinase that binds bacterial flagellin or a flagellin-based flg22 peptide through a presumed physical interaction within the FLS2 extracellular domain. Double-Ala scanning mutagenesis of solvent-exposed beta-strand/beta-turn residues across the FLS2 LRR domain identified LRRs 9 to 15 as contributors to flagellin responsiveness. FLS2 LRR-encoding domains from 15 Arabidopsis ecotypes and 20 diverse Brassicaceae accessions were isolated and sequenced. FLS2 is highly conserved across most Arabidopsis ecotypes, whereas more diversified functional FLS2 homologs were found in many but not all Brassicaceae accessions. flg22 responsiveness was correlated with conserved LRR regions using Conserved Functional Group software to analyze structural models of the LRR for diverse FLS2 proteins. This identified conserved spatial clusters of residues across the beta-strand/beta-turn residues of LRRs 12 to 14, the same area identified by the Ala scan, as well as other conserved sites. Site-directed randomizing mutagenesis of solvent-exposed beta-strand/beta-turn residues across LRRs 9 to 15 identified mutations that disrupt flg22 binding and showed that flagellin perception is dependent on a limited number of tightly constrained residues of LRRs 9 to 15 that make quantitative contributions to the overall phenotypic response.     

Refinement of 3D models of horseradish peroxidase isoenzyme C: Predictions of 2D NMR assignments and substrate binding sites

In this study, two alternative three-dimensional (3D) models of horseradish peroxidase (HRP-C)—differing mainly in the structure of a long untemplated insertion—were refined, systematically assessed, and used to make predictions that can both guide and be tested by future experimental studies. A key first step in the model-building process was a procedure for multiple sequence alignment based on structurally conserved regions and key conserved residues, including those side chains providing ligands to the two Ca²⁺ binding sites. The model refinements reported here include (1) optimization of side-chain conformations; (3) addition of structural waters using a template-independent procedure; (2) structural refinement of the untemplated 34 amino acid insertion located between the F and G helices, using both energy criteria and NMR data; (4) unconstrained energy optimization of the refined models. Using these procedures, two refined structures of HRP-C were obtained, differing mainly in the conformation of this long insertion. The presence of residues in this insertion that could potentially interact with bound substrates suggests a functional role that may be related to the general ability of class III peroxidases to form stable 1:1 complexes with a variety of substrates. The structural validity of the models was systematically assessed by a variety of criteria. Most notably, the ProsaII z scores and Profiles 3D scores of the two HRP-C models indicated that they are significantly better than would be obtained by simple amino acid replacement, using any of the known structures as a template. These two 3D HRP-C models, were then used to predict candidate residues for the assignment of NOESY cross-peaks previously noted in 2D-NMR studies. Specifically, the residues known as Ile X, Phe A, Phe B, aliphatic residue Q, and Ile T. Candidate substrate binding sites were also identified and compared with experimentally based predictions. This work is timely because new X-ray structures are anticipated that will facilitate the validation of these procedures. © 1996 Wiley-Liss, Inc.     

Structure-function studies of [2Fe-2S] ferredoxins

The ability to overexpress [2Fe-2S] ferredoxins inEscherichia coli has opened up exciting research opportunities. High-resolution x-ray structures have been determined for the wild-type ferredoxins produced by the vegetative and heterocyst forms ofAnabaena strain 7120 (in their oxidized states), and these have been compared to structural information derived from multidimensional, multinuclear NMR spectroscopy. The electron delocalization in these proteins in their oxidized and reduced states has been studied by¹H,²H,¹³C, and¹⁵N NMR spectroscopy. Site-directed mutagenesis has been used to prepare variants of these ferredoxins. Mutants (over 50) of the vegetative ferredoxin have been designed to explore questions about cluster assembly and stabilization and to determine which residues are important for recognition and electron transfer to the redox partnerAnabaena ferredoxin reductase. The results have shown that serine can replace cysteine at each of the four cluster attachment sites and still support cluster assembly. Electron transfer has been demonstrated with three of the four mutants. Although these mutants are less stable than the wild-type ferredoxin, it has been possible to determine the x-ray structure of one (C49S) and to characterize all four by EPR and NMR. Mutagenesis has identified residues 65 and 94 of the vegetative ferredoxin as crucial to interaction with the reductase. Three-dimensional models have been obtained by x-ray diffraction analysis for several additional mutants: T48S, A50V, E94K (four orders of magnitude less active than wild type in functional assays), and A43S/A45S/T48S/A50N (quadruple mutant).     

Protein structural similarities predicted by a sequence-structure compatibility method.

A method for protein structure prediction has been developed, which evaluates the compatibility of an amino acid sequence with known 3-dimensional structures and identifies the most likely structure. The method was applied to a large number of sequences in a database, and the structures of the following proteins were predicted: (1) shikimate kinase (SKase), (2) the hydrophilic subunit of mannose permease (IIABMan), (3) rat tyrosine aminotransferase (Tyr AT), and (4) threonine dehydratase (TDH). The functional and evolutionary implications of the predictions are discussed. (1) The structural similarity between SKase and adenylate kinase was predicted. Alignment of their sequences reveals that the ATP-binding type A sequence motif and 2 ATP-binding arginine residues are conserved. The prediction suggests a similarity in their functional mechanisms as well as an evolutionary relationship. (2) The structural similarity between IIABMan and galactose/glucose-binding protein (GGBP) was predicted. The IIA and IIB domains are aligned with the N- and C-terminal domains of GGBP, respectively. The 2 phosphorylated residues, His 10 and His 175, of IIABMan are threaded onto loops located in the substrate-binding cleft of GGBP. The prediction accounts for the phosphoryl transfer from His 10 to His 175, and to the sugar substrate. (3) The structural similarity between rat Tyr AT and Escherichia coli aspartate AT was predicted, as well as (4) the structural similarity between TDH and the tryptophan synthase beta subunit. Predictions (3) and (4) support the previous predictions based on observations of the functional similarities between the proteins.     

Sequence conservation in families whose members have little or no sequence similarity: the four-helical cytokines and cytochromes

Proteins for which there are good structural, functional and genetic similarities that imply a common evolutionary origin, can have sequences whose similarities are low or undetectable by conventional sequence comparison procedures. Do these proteins have sequence conservation beyond the simple conservation of hydrophobic and hydrophilic character at specific sites and if they do what is its nature? To answer these questions we have analysed the structures and sequences of two superfamilies: the four-helical cytokines and cytochromes c'-b(562). Members of these superfamilies have sequence similarities that are either very low or not detectable. The cytokine superfamily has within it a long chain family and a short chain family. The sequences of known representative structures of the two families were aligned using structural information. From these alignments we identified the regions that conserve the same main-chain conformation: the common core (CC). For members of the same family, the CC comprises some 50% of the individual structures; for the combination of both families it is 30%. We added homologous sequences to the structural alignment. Analysis of the residues occurring at sites within the CCs showed that 30% have little or no conservation, whereas about 40% conserve the polar/neutral or hydrophobic/neutral character of their residues. The remaining 30% conserve hydrophobic residues with strong or medium limitations on their volume variations. Almost all of these residues are found at sites that form the "buried spine" of each helix (at sites i, i+3, i+7, i+10, etc., or i, i+4, i+7, i+11, etc.) and they pack together at the centre of each structure to give a pattern of residue-residue contacts that is almost absolutely conserved. These CC conserved hydrophobic residues form only 10-15% of all the residues in the individual structures.A similar analysis of the cytochromes c'-b(562), which bind haem and have a very different function to that of the cytokines, gave very similar results. Again some 30% of the CC residues have hydrophobic residues with strong or medium conservation. Most of these form the buried spine of each helix and play the same role as those in the cytokines. The others, and some spine residues bind the haem co-factor.     

Structural and functional relations among thioredoxins of different species

Three-dimensional models have been constructed of homologous thioredoxins and protein disulfide isomerases based on the high resolution x-ray crystallographic structure of the oxidized form of Escherichia coli thioredoxin. The thioredoxins, from archebacteria to humans, have 27-69% sequence identity to E. coli thioredoxin. The models indicate that all the proteins have similar three-dimensional structures despite the large variation in amino acid sequences. As expected, residues in the active site region of thioredoxins are highly conserved. These include Asp-26, Ala-29, Trp-31, Cys-32, Gly-33, Pro-34, Cys-35, Asp-61, Pro-76, and Gly-92. Similar residues occur in most protein disulfide isomerase sequences. Most of these residues form the surface around the active site that appears to facilitate interactions with other enzymes. Other structurally important residues are also conserved. A proline at position 40 causes a kink in the alpha-2 helix and thus provides the proper position of the active site residues at the amino end of this helix. Pro-76 is important in maintaining the native structure of the molecule. In addition, residues forming the internal contact surfaces between the secondary structural elements are generally unchanged such as Phe-12, Val-25, and Phe-27.     

Combining protein sequence,structure, and dynamics: A novel approach for functional evolution analysis of PAS domain superfamily

PAS domains are widespread in archaea, bacteria, and eukaryota, and play important roles in various functions. In this study, we aim to explore functional evolutionary relationship among proteins in the PAS domain superfamily in view of the sequence‐structure‐dynamics‐function relationship. We collected protein sequences and crystal structure data from RCSB Protein Data Bank of the PAS domain superfamily belonging to three biological functions (nucleotide binding, photoreceptor activity, and transferase activity). Protein sequences were aligned and then used to select sequence‐conserved residues and build phylogenetic tree. Three‐dimensional structure alignment was also applied to obtain structure‐conserved residues. The protein dynamics were analyzed using elastic network model (ENM) and validated by molecular dynamics (MD) simulation. The result showed that the proteins with same function could be grouped by sequence similarity, and proteins in different functional groups displayed statistically significant difference in their vibrational patterns. Interestingly, in all three functional groups, conserved amino acid residues identified by sequence and structure conservation analysis generally have a lower fluctuation than other residues. In addition, the fluctuation of conserved residues in each biological function group was strongly correlated with the corresponding biological function. This research suggested a direct connection in which the protein sequences were related to various functions through structural dynamics. This is a new attempt to delineate functional evolution of proteins using the integrated information of sequence, structure, and dynamics.     

SPLASH: structural pattern localization analysis by sequential histograms

MOTIVATION: The discovery of sparse amino acid patterns that match repeatedly in a set of protein sequences is an important problem in computational biology. Statistically significant patterns, that is patterns that occur more frequently than expected, may identify regions that have been preserved by evolution and which may therefore play a key functional or structural role. Sparseness can be important because a handful of non-contiguous residues may play a key role, while others, in between, may be changed without significant loss of function or structure. Similar arguments may be applied to conserved DNA patterns. Available sparse pattern discovery algorithms are either inefficient or impose limitations on the type of patterns that can be discovered. RESULTS: This paper introduces a deterministic pattern discovery algorithm, called Splash, which can find sparse amino or nucleic acid patterns matching identically or similarly in a set of protein or DNA sequences. Sparse patterns of any length, up to the size of the input sequence, can be discovered without significant loss in performances. Splash is extremely efficient and embarrassingly parallel by nature. Large databases, such as a complete genome or the non-redundant SWISS-PROT database can be processed in a few hours on a typical workstation. Alternatively, a protein family or superfamily, with low overall homology, can be analyzed to discover common functional or structural signatures. Some examples of biologically interesting motifs discovered by Splash are reported for the histone I and for the G-Protein Coupled Receptor families. Due to its efficiency, Splash can be used to systematically and exhaustively identify conserved regions in protein family sets. These can then be used to build accurate and sensitive PSSM or HMM models for sequence analysis. AVAILABILITY: Splash is available to non-commercial research centers upon request, conditional on the signing of a test field agreement. CONTACT: acal@us.ibm.com, Splash main page http://www.research.ibm.com/splash     

Binary patterning of polar and nonpolar amino acids in the sequences and structures of native proteins

Protein sequences can be represented as binary patterns of polar (○) and nonpolar (?) amino acids. These binary sequence patterns are categorized into two classes: Class A patterns match the structural repeat of an idealized amphiphilic α-helix (3.6 residues per turn), and class B patterns match the structural repeat of an idealized amphiphilic β-strand (2 residues per turn). The difference between these two classes of sequence patterns has led to a strategy for de novo protein design based on binary patterning of polar and nonpolar amino acids. Here we ask whether similar binary patterning is incorporated in the sequences and structures of natural proteins. Analysis of the Protein Data Bank demonstrates the following. (1) Class A sequence patterns occur considerably more frequently in the sequences of natural proteins than would be expected at random, but class B patterns occur less often than expected. (2) Each pattern is found predominantly in the secondary structure expected from the binary strategy for protein design. Thus, class A patterns are found more frequently in α-helices than in β-strands, and class B patterns are found more frequently in β-strands than in α-helices. (3) Among the α-helices of natural proteins, the most commonly used binary patterns are indeed the class A patterns. (4) Among all β-strands in the database, the most commonly used binary patterns are not the expected class B patterns. (5) However, for solvent-exposed β-strands, the correlation is striking: All β-strands in the database that contain the class B patterns are exposed to solvent. (6) The bias of class A patterns for α-structure over β-structure and the bias of class B patterns for β-structure over α-structure are significant, not merely when compared to other binary patterns of polar (○) and nonpolar (?) amino acids, but also when compared to the full range of sequences in the database. The implications for the design of novel proteins are discussed.     

Nature disfavors sequences of alternating polar and non-polar amino acids: implications for amyloidogenesis

Recent experiments with combinatorial libraries of de novo proteins have demonstrated that sequences designed to contain polar and non-polar amino acid residues arranged in an alternating pattern form fibrillar structures resembling beta-amyloid. This finding prompted us to probe the distribution of alternating patterns in the sequences of natural proteins. Analysis of a database of 250,514 protein sequences (79,708,024 residues) for all possible binary patterns of polar and non-polar amino acid residues revealed that alternating patterns occur significantly less often than other patterns with similar compositions. The under-representation of alternating binary patterns in natural protein sequences, coupled with the observation that such patterns promote amyloid-like structures in de novo proteins, suggests that sequences of alternating polar and non-polar amino acids are inherently amyloidogenic and consequently have been disfavored by evolutionary selection.