Crystallographic evidence for water-assisted photo-induced peptide cleavage in the stony coral fluorescent protein Kaede

A coral fluorescent protein from Trachyphyllia geoffroyi, Kaede, possesses a tripeptide of His62-Tyr63-Gly64, which forms a chromophore with green fluorescence. This chromophore's fluorescence turns red following UV light irradiation. We have previously shown that such photoconversion is achieved by a formal beta-elimination reaction, which results in a cleavage of the peptide bond found between the amide nitrogen and the alpha-carbon at His62. However, the stereochemical arrangement of the chromophore and the precise structural basis for this reaction mechanism previously remained unknown. Here, we report the crystal structures of the green and red form of Kaede at 1.4 A and 1.6 A resolutions, respectively. Our structures depict the cleaved peptide bond in the red form. The chromophore conformations both in the green and red forms are similar, except a well-defined water molecule in the proximity of the His62 imidazole ring in the green form. We propose a molecular mechanism for green-to-red photoconversion, which is assisted by the water molecule.     

Quantitative Analysis of the Conservation of the Tertiary Structure of Protein Segments

The publication of the crystallographic structure of calmodulin protein has offered an example leading us to believe that it is possible for many protein sequence segments to exhibit multiple 3D structures referred to as multi-structural segments. To this end, this paper presents statistical analysis of uniqueness of the 3D-structure of all possible protein sequence segments stored in the Protein Data Bank (PDB, Jan. of 2003, release 103) that occur at least twice and whose lengths are greater than 10 amino acids (AAs). We refined the set of segments by choosing only those that are not parts of longer segments, which resulted in 9297 segments called a sponge set. By adding 8197 signature segments, which occur uniquely in the PDB, into the sponge set we have generated a benchmark set. Statistical analysis of the sponge set demonstrates that rotating, missing and disarranging operations described in the text, result in the segments becoming multi-structural. It turns out that missing segments do not exhibit a change of shape in the 3D-structure of a multi-structural segment. We use the root mean square distance for unit vector sequence (URMSD) as an improved measure to describe the characteristics of hinge rotations, missing, and disarranging segments. We estimated the rate of occurrence for rotating and disarranging segments in the sponge set and divided it by the number of sequences in the benchmark set which is found to be less than 0.85%. Since two of the structure changing operations concern negligible number of segment and the third one is found not to have impact on the structure, we conclude that the 3D-structure of proteins is conserved statistically for more than 98% of the segments. At the same time, the remaining 2% of the sequences may pose problems for the sequence alignment based structure prediction methods.*Jishou Ruan research was supported by Liuhui Center for Applied Mathematics, China-Canada exchange program administered by MITACS and NSFC (10271061). ^#Ke Chen and Lukasz A. Kurgan research was partially supported by NSERC Canada. ^†Jack A. Tuszynkski research has been supported by MITACS, NSERC Canada and the Allard Foundation.     

Comprehensive evaluation of protein structure alignment methods: scoring by geometric measures

We report the largest and most comprehensive comparison of protein structural alignment methods. Specifically, we evaluate six publicly available structure alignment programs: SSAP, STRUCTAL, DALI, LSQMAN, CE and SSM by aligning all 8,581,970 protein structure pairs in a test set of 2930 protein domains specially selected from CATH v.2.4 to ensure sequence diversity. We consider an alignment good if it matches many residues, and the two substructures are geometrically similar. Even with this definition, evaluating structural alignment methods is not straightforward. At first, we compared the rates of true and false positives using receiver operating characteristic (ROC) curves with the CATH classification taken as a gold standard. This proved unsatisfactory in that the quality of the alignments is not taken into account: sometimes a method that finds less good alignments scores better than a method that finds better alignments. We correct this intrinsic limitation by using four different geometric match measures (SI, MI, SAS, and GSAS) to evaluate the quality of each structural alignment. With this improved analysis we show that there is a wide variation in the performance of different methods; the main reason for this is that it can be difficult to find a good structural alignment between two proteins even when such an alignment exists. We find that STRUCTAL and SSM perform best, followed by LSQMAN and CE. Our focus on the intrinsic quality of each alignment allows us to propose a new method, called "Best-of-All" that combines the best results of all methods. Many commonly used methods miss 10-50% of the good Best-of-All alignments. By putting existing structural alignments into proper perspective, our study allows better comparison of protein structures. By highlighting limitations of existing methods, it will spur the further development of better structural alignment methods. This will have significant biological implications now that structural comparison has come to play a central role in the analysis of experimental work on protein structure, protein function and protein evolution.     

Structural and biochemical study of Bacillus subtilis HmoB in complex with heme

Most bacteria have developed a hemoprotein degradation system to acquire iron from their hosts. Bacillus subtilis HmoB, a heme monooxygenase, is involved in the degradation of heme and subsequent release of iron. HmoB contains a C-terminal ABM domain, which is similar in sequence and structure to other heme monooxygenases. Heme degradation assay showed that highly conserved residues (N70, W128, and H138) near the heme-binding site were critical for activity of HmoB. However, HmoB was shown to be different from other bacterial heme oxygenases due to its longer N-terminal region and formation of a biological monomer instead of a dimer. The degradation product of B. subtilis HmoB was identified as staphylobilin from mass spectrometric analysis of the product and release of formaldehyde during degradation reaction.     

Three-dimensional structure of iminodisuccinate epimerase defines the fold of the MmgE/PrpD protein family

Iminodisuccinate (IDS) epimerase catalyzes the epimerisation of R,R-, S,S- and R,S- iminodisuccinate, one step in the biodegradation of the chelating agent iminodisuccinate by Agrobacterium tumefaciens BY6. The enzyme is a member of the MmgE/PrpD protein family, a diverse and little characterized class of proteins of prokaryotic and eukaryotic origin. IDS epimerase does not show significant overall amino acid sequence similarity to any other protein of known three-dimensional structure. The crystal structure of this novel epimerase has been determined by multi-wavelength diffraction to 1.5 A resolution using selenomethionine-substituted enzyme. In the crystal, the enzyme forms a homo-dimer, and the subunit consists of two domains. The larger domain, not consecutive in sequence and comprising residues Met1-Lys266 and Leu400-Pro446, forms a novel all alpha-helical fold with a central six-helical bundle. The second, smaller domain folds into an alpha+beta domain, related in topology to chorismate mutase by a circular permutation. IDS epimerase is thus not related in three-dimensional structure to other known epimerases. The fold of the IDS epimerase is representative for the whole MmgE/PrpD family. The putative active site is located at the interface between the two domains of the subunit, and is characterized by a positively charged surface, consistent with the binding of a highly negatively charged substrate such as iminodisuccinate. Docking experiments suggest a two-base mechanism for the epimerisation reaction.     

Decoy models for protein structure comparison score normalisation

A method is described to construct sets of decoy models that can be used to generate a background score distribution for protein structure comparison. The models are derived directly from the two proteins being compared and retain all the essential properties of the structures, including length, density, shape and secondary structure composition but have different folds. As each comparison involves a pair of proteins of the same length, no explicit normalisation is required to adjust for the length of the proteins being compared. This allows substructure (or domain) matches to score almost equally to the comparison of isolated domains. A normalised probability measure was derived that allows joint family/family comparison. The method was applied to some of the CASP6 models for targets with new folds.     

The PDB is a covering set of small protein structures

Structure comparisons of all representative proteins have been done. Employing the relative root mean square deviation (RMSD) from native enables the assessment of the statistical significance of structure alignments of different lengths in terms of a Z-score. Two conclusions emerge: first, proteins with their native fold can be distinguished by their Z-score. Second and somewhat surprising, all small proteins up to 100 residues in length have significant structure alignments to other proteins in a different secondary structure and fold class; i.e. 24.0% of them have 60% coverage by a template protein with a RMSD below 3.5 Å and 6.0% have 70% coverage. If the restriction that we align proteins only having different secondary structure types is removed, then in a representative benchmark set of proteins of 200 residues or smaller, 93% can be aligned to a single template structure (with average sequence identity of 9.8%), with a RMSD less than 4 Å, and 79% average coverage. In this sense, the current Protein Data Bank (PDB) is almost a covering set of small protein structures. The length of the aligned region (relative to the whole protein length) does not differ among the top hit proteins, indicating that protein structure space is highly dense. For larger proteins, non-related proteins can cover a significant portion of the structure. Moreover, these top hit proteins are aligned to different parts of the target protein, so that almost the entire molecule can be covered when combined. The number of proteins required to cover a target protein is very small, e.g. the top ten hit proteins can give 90% coverage below a RMSD of 3.5 Å for proteins up to 320 residues long. These results give a new view of the nature of protein structure space, and its implications for protein structure prediction are discussed.     

Acidic pH-induced conformational changes in amyloidogenic mutant transthyretin

Several proteins, including transthyretin (TTR), can generate in tissues extracellular insoluble aggregates, in the form of fibrils, that are associated with pathological states known as amyloidoses. To date, more than 80 different TTR point mutations have been associated with hereditary amyloidosis in humans. In vitro, the formation of amyloid fibrils by human TTR is known to be triggered by acidic pH. We show here that, in vitro, the natural amyloidogenic I84S and the non-natural I84A TTR mutant forms exhibit a propensity to produce fibrils in an acidic medium significantly higher than that of wild-type TTR. The two mutant forms have been crystallized at both neutral and acidic pH. Their neutral pH crystal structures are very similar to that of wild-type TTR, consistent with previous evidence indicating that only minor structural changes are induced by amyloidogenic mutations. On the contrary, their crystal structures at moderately low pH (4.6) show significant conformational differences as compared to their neutral pH structures. Remarkably, such changes are not induced in wild-type TTR crystallized at low pH. The most relevant consist of the unwinding of the TTR short alpha-helix and of the change in conformation of the loop connecting the alpha-helix to beta-strand F. Only one monomer of the crystallographic dimer is affected, causing a disruption of the tetrameric symmetry. This asymmetry and a possible destabilization of the tetrameric quaternary structure of TTR may be responsible for the amyloidogenic potential of the two TTR mutant forms at low pH.     

Clustering of protein structural fragments reveals modular building block approach of nature

Structures of peptide fragments drawn from a protein can potentially occupy a vast conformational continuum. We co-ordinatize this conformational space with the help of geometric invariants and demonstrate that the peptide conformations of the currently available protein structures are heavily biased in favor of a finite number of conformational types or structural building blocks. This is achieved by representing a peptides' backbone structure with geometric invariants and then clustering peptides based on closeness of the geometric invariants. This results in 12,903 clusters, of which 2207 are made up of peptides drawn from functionally and/or structurally related proteins. These are termed "functional" clusters and provide clues about potential functional sites. The rest of the clusters, including the largest few, are made up of peptides drawn from unrelated proteins and are termed "structural" clusters. The largest clusters are of regular secondary structures such as helices and beta strands as well as of beta hairpins. Several categories of helices and strands are discovered based on geometric differences. In addition to the known classes of loops, we discover several new classes, which will be useful in protein structure modeling. Our algorithm does not require assignment of secondary structure and, therefore, overcomes the limitations in loop classification due to ambiguity in secondary structure assignment at loop boundaries.     

The crystal structure of the human (S100A8/S100A9)2 heterotetramer, calprotectin, illustrates how conformational changes of interacting alpha-helices can determine specific association of two EF-hand proteins

The EF-hand proteins S100A8 and S100A9 are important calcium signalling proteins that are involved in wound healing and provide clinically relevant markers of inflammatory processes, such as rheumatoid arthritis and inflammatory bowel disease. Both can form homodimers via distinct modes of association, probably of lesser stability in the case of S100A9, whereas in the presence of calcium S100A8 and S100A9 associate to calprotectin, the physiologically active heterooligomer. Here we describe the crystal structure of the (S100A8/S100A9)(2) heterotetramer at 1.8 A resolution. Its quaternary structure illustrates how specific heteroassociation is energetically driven by a more extensive burial of solvent accessible surface areas in both proteins, most pronounced for S100A9, thus leading to a dimer of heterodimers. A major contribution to tetramer association is made by the canonical calcium binding loops in the C-terminal halves of the two proteins. The mode of heterodimerisation in calprotectin more closely resembles the subunit association previously observed in the S100A8 homodimer and provides trans stabilisation for S100A9, which manifests itself in a significantly elongated C-terminal alpha-helix in the latter. As a consequence, two different putative zinc binding sites emerge at the S100A8/S100A9 subunit interface. One of these corresponds to a high affinity arrangement of three His residues and one Asp side-chain, which is unique to the heterotetramer. This structural feature explains the well known Zn(2+) binding activity of calprotectin, whose overexpression can cause strong dysregulation of zinc homeostasis with severe clinical symptoms.