# Deciphering structure and topology of conserved COG2042 orphan proteins

Jean Armengaud, Alain Dedieu, Olivier Solques, Jean-Luc Pellequer, Eric Quemeneur

► To cite this version:

Jean Armengaud, Alain Dedieu, Olivier Solques, Jean-Luc Pellequer, Eric Quemeneur. Deciphering structure and topology of conserved COG2042 orphan proteins. BMC Structural Biology, 2005, 5 (3), pp.1-16. 10.1186/1472-6807-5-3. hal-04269608

HAL Id: hal-04269608

https://hal.science/hal-04269608v1

Submitted on 3 Nov 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons CC BY 4.0 - Attribution - International License
Research article

Open Access

# Deciphering structure and topology of conserved COG2042 orphan proteins

Jean Armengaud*, Alain Dedieu, Olivier Solques, Jean-Luc Pellequer and Eric Quemeneur

Address: CEA-VALRHO, DSV-DIEP-SBTN, Service de Biochimie post-génomique & Toxicologie Nucléaire, Bagnols-sur-Cèze, France

Email: Jean Armengaud* - armengaud@cea.fr; Alain Dedieu - alain.dedieu@cea.fr; Olivier Solques - solques_olivier@yahoo.fr; Jean-Luc Pellequer - jlpellequer@cea.fr; Eric Quemeneur - eric.quemeneur@cea.fr

Published: 08 February 2005

Received: 13 October 2004

BMC Structural Biology 2005, 5:3 doi:10.1186/1472-6807-5-3

Accepted: 08 February 2005

This article is available from: http://www.biomedcentral.com/1472-6807/5/3

© 2005 Armengaud et al; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

## Abstract

**Background:** The cluster of orthologous group COG2042 has members in all sequenced Eukaryota as well as in many Archaea. The cellular function of these proteins of ancient origin remains unknown. PSI-BLAST analysis does not indicate a possible link with even remotely-related proteins that have been functionally or structurally characterized. As a prototype among COG2042 orthologs, SSO0551 protein from the hyperthermophilic archaeon *Sulfolobus solfataricus* was purified to homogeneity for biophysical characterization.

**Results:** The untagged protein is thermostable and behaves as a monomeric protein in gel filtration experiment. Several mass spectrometry-based strategies were combined to obtain a set of low resolution structural information. Kinetic data from limited proteolysis with various endoproteases are concordant in pointing out that region Glu<sup>73</sup>-Arg<sup>78</sup> is hyper-sensitive, and thus accessible and flexible. Lysine labeling with NHS-biotin and cross-linking with DTSSP revealed that the 35 amino acid RLI motif at the N terminus is solvent exposed. Cross-links between Lys<sup>10</sup>-Lys<sup>14</sup> and Lys<sup>23</sup>-Lys<sup>25</sup> indicate that these residues are spatially close and in adequate conformation to be cross-linked. These experimental data have been used to rank multiple three-dimensional models generated by a *de novo* procedure.

**Conclusion:** Our data indicate that COG2042 proteins may share a novel fold. Combining biophysical, mass-spectrometry data and molecular model is a useful strategy to obtain structural information and to help in prioritizing targets in structural genomics programs.

## Background

Genomic comparative studies on entirely sequenced genomes from the three domains of life, i.e. Bacteria, Archaea and Eukaryota [1], evidenced that proteins involved in the organization or processing of genetic information (structures of ribosome and chromatin, translation, transcription, replication and DNA repair)

display a closer relationship between Archaea and Eukaryota than between Bacteria and Eukaryota [2-4]. To identify new proteins involved in such important cellular mechanisms, an exhaustive inventory of proteins of unknown function common to only Eukaryota and Archaea but not in Bacteria has been devised [5-7]. Among such proteins, the Cluster of Orthologous Group COG2042 comprises

<footnote>* Corresponding author</footnote>
proteins ubiquitously present in Eukaryota and present in many, but not all, Archaea; a hallmark of their ancient origin. The corresponding ancestral protein should have been present in the common ancestor of these two domains of life. Some partial experimental data are known from the *Saccharomyces cerevisiae* COG2042 homolog. Deletion of the *Yor006c* gene was shown to result in a viable phenotype but some apparent moderate growth defects were noticed on a fermentable carbon source [8,9]. Two putative protein partners for *Yor006c* were identified through a high-throughput two-hybrid study [10]: Ydl017w, a serine/threonine kinase also known as the cell division control protein 7 (Cdc7), and Yil025c, a hypothetical ORF. However, the cellular function of COG2042 proteins remains unknown.

A polar region, named RLI, is conserved at the N terminus of COG2042 proteins as well as at the N terminus of another cluster of orthologous proteins, namely COG1245. The latter, exemplified by SSO0287 in *Sulfolobus solfataricus* [11], are large proteins (about 600 residues) that encompass four different domains: a RLI domain, a [4Fe-4S] ferredoxin domain, and two ATPase domains, usually found in ABC transporter. Their putative function is currently subjected to discussion [12,13] but could be related to rRNA metabolism. Indeed, four of the eleven proteins shown to interact with the yeast COG1245 homolog (Ydr091c) were identified as involved in rRNA metabolism (Ymr047c, Ydl213c, Ylr340w, Ylr192c). Experimental data on the human homolog of Ydr091c indicated that this protein reversibly associates with RNaseL, and thus COG1245 proteins were named RNase L inhibitor [14].

Because knowledge of protein structure is of high importance to understand protein function, huge efforts have been recently invested in high-throughput protein structure determination programs [15]. Recent reports indicate that only a relatively small percentage of expressed and purified proteins are amenable to full 3D structure by NMR or crystallography and X-ray diffraction [16,17]. In silico modeling (homology modeling, fold recognition, ab initio and de novo modeling) is the alternative to quickly gain the fold of a protein. However, such approach sometimes remains ambiguous in reliably identifying correct structures for protein sequences remotely-related to those found in PDB database. A promising strategy is the use of experimental data (if possible easily obtained) for model discrimination or refinement [18-20]. For example, the tertiary structure of the bovine basic fibroblast growth factor (FGF)-2 was probed with a lysine-specific cross-linking agent and subjected to tryptic peptide mapping by mass spectrometry to identify the sites of cross-linking [21]. The low resolution interatomic distance information obtained experimentally allowed the authors to distinguish among

threading models in spite of a relatively low sequence similarity (13 % of identical residues). Interestingly, the constant development of novel cross-linking reagents suitable for mass spectrometry [22] enables enrichment of cross-linked peptides facilitating such strategy. A chemical modification approach [23-26], in combination with limited proteolysis procedures [27,28], can also provide useful structural constraints [29] for model refinement.

A step further is to attempt such approaches with proteins having no detectable homologs. In order to get insight into the topology of COG2042 members and if possible to use these experimental data to discriminate among structural protein templates, we combined limited proteolysis, lysine labeling and cross-linking strategies. The protein SSO0551 from the hyperthermophilic archaea *Sulfolobus solfataricus* was chosen as a prototype because of its thermostability and the probable absence of post-translational modifications when produced as a recombinant form in *Escherichia coli*. The SSO0551 protein is monomeric with a low molecular mass (19 kDa). This size is easily amenable to characterization by mass spectrometry. Our results reveal that the polar RLI motif at the N terminus is probably structured and solvent exposed, pointing at a common trait between COG2042 and COG1245 proteins, this latter group being also conserved in Eukaryota and Archaea but absent in Bacteria. The accessible and flexible regions defined by limited proteolysis combined with lysine accessibility assessed by NHS-biotin labeling and DTSSP cross-linking allowed us to discriminate among ten top ranking *de novo* three-dimensional (3D) models.

## Results

### COG2042 comprises members exclusively from Eukaryota and Archaea

The sequence of SSO0551 from *S. solfataricus* was used as query in a PSI-BLAST database search to identify homologous proteins. A constant cutoff expectation value of 10⁻¹⁵ resulted after three iterations in selection of 40 sequences (15 from Archaea and 25 from Eukaryota) that were all aligned over their full length. No close homologs (E-value below 10⁻¹⁰ in the third iteration) with full-length sequence matching to SSO0551 were found among Bacteria. Remarkably, all completely sequenced Eukaryal organisms were found to have one SSO0551 homolog. Fig. 1 shows an unrooted phylogram of the updated COG2042 family (Fig. 1A) and an alignment of a selection of six representative sequences (Fig. 1B), selected on the basis of their phylogenetic distribution. When experimental evidences concerning the protein are unfortunately lacking for ORF description genome annotators usually take into consideration the most upstream initiation codon. For this reason, the most probable start codons of several open reading frames should be
<table><caption>Table 1: Fingerprint identification of recombinant products from pSBTN-AB31 and pSBTN-AB30 constructs.</caption><thead><tr><th colspan="2">pSBTN-AB31 construct<sup>a</sup></th><th colspan="2">pSBTN-AB30 construct<sup>a</sup></th><th colspan="2">Theoretical peptides</th></tr><tr><th>[MH]<sup>+</sup> observed (in amu)</th><th>Δmass (in ppm)</th><th>[MH]<sup>+</sup> observed (in amu)</th><th>Δmass (in ppm)</th><th>[MH]<sup>+</sup> expected (in amu)</th><th>Position<sup>b</sup> (Start-End)</th><th>Peptide assignment</th></tr></thead><tbody><tr><td>1590.64</td><td>-31</td><td>nd<sup>c</sup></td><td>-</td><td>1590.69</td><td>3–16</td><td>GSHHHHHHGMASMK</td></tr><tr><td>nd<sup>c</sup></td><td>-</td><td>1050.54</td><td>-19</td><td>1050.56</td><td>17–24</td><td>VYIIDYHK</td></tr><tr><td>1505.71</td><td>-33</td><td>nd<sup>c</sup></td><td>-</td><td>1505.76</td><td>17–28</td><td>VYIIDYHKDDPK</td></tr><tr><td>1661.83</td><td>-18</td><td>1661.85</td><td>-6</td><td>1661.86</td><td>17–29</td><td>VYIIDYHKDDPKR</td></tr><tr><td>1601.84</td><td>-31</td><td>1601.83</td><td>-37</td><td>1601.89</td><td>49–63</td><td>GVVLDPFAQITLSNK</td></tr><tr><td>1844.98</td><td>-16</td><td>1845.00</td><td>-5</td><td>1845.01</td><td>49–65</td><td>GVVLDPFAQITLSNKD</td></tr><tr><td>2328.28</td><td>-4</td><td>2328.28</td><td>-4</td><td>2328.29</td><td>49–69</td><td>GVVLDPFA...KDKDIVR</td></tr><tr><td>2484.35</td><td>-16</td><td>2484.34</td><td>-20</td><td>2484.39</td><td>49–70</td><td>GVVLDPFA...KDKDIVRR</td></tr><tr><td>2140.06</td><td>0</td><td>2140.06</td><td>0</td><td>2140.06</td><td>71–89</td><td>IGITIVDTSWNNTSQSEFK</td></tr><tr><td>2296.16</td><td>0</td><td>2296.16</td><td>0</td><td>2296.16</td><td>70–89</td><td>RIGITIVDTSWNNTSQSEFK</td></tr><tr><td>2043.17</td><td>10</td><td>2043.18</td><td>15</td><td>2043.15</td><td>97–114</td><td>RIPILFAGNPIHYGIAYK</td></tr><tr><td>1887.05</td><td>0</td><td>1887.06</td><td>5</td><td>1887.05</td><td>98–114</td><td>IPILFAGNPIHYGIAYK</td></tr><tr><td>1244.62</td><td>-16</td><td>1244.59</td><td>-40</td><td>1244.64</td><td>143–152</td><td>WGHTFIELNK</td></tr><tr><td>2091.09</td><td>0</td><td>nd<sup>c</sup></td><td>-</td><td>2091.09</td><td>143–159</td><td>WGHTFIELNKELLEAYK</td></tr><tr><td>nd<sup>c</sup></td><td>-</td><td>865.53</td><td>69</td><td>865.47</td><td>153–159</td><td>ELLEAYK</td></tr></tbody></table>

reconsidered after exhaustive alignment (Fig. 1). For example, atg codon starting at nucleotide 484790 on the Crick strand for SSO0551 from S. solfataricus (NC_002754) should be a more appropriate start codon than atg starting at nucleotide 484916 and mentioned erroneously in current database. From the unrooted phylogram (Fig. 1A), two main lineages (archaeal and eukaryal) can be defined based on organism origin. This suggests that occurrence of these proteins is at least as ancient as divergence of these two phyla. No paralogs, sign of a possible evolution of a new derived function, have been evidenced in entirely sequenced organisms currently available. Although these proteins are of ancient origin, the core sequence appears well conserved as observed in Fig. 1B. Thirty-three residues (38%) are found identical in the core central segment (out of 88 amino acids) between the most distant COG2042 orthologs, namely gi48852409 from Ferroplasma acidarmanus and gi6324579 from Saccharomyces cerevisiae (Fig. 1B). From the alignment, several conserved motifs that may be functionally crucial (cofactor or substrate binding, catalysis, or partner interactions) were detected. A conserved hexapeptide sequence, Val-Val/Ile/Leu-Asp/Glu-Cys-Ser-Trp (motif I in Fig. 1B), is found distant of 14–17 amino acids from another conserved motif of 25 amino acids containing 4 polar, 18 hydrophobic and 3 aromatic residues (motif II). Database searches with these motifs as queries did not allow identification of remotely-related proteins.

All sequences from COG2042 encompass a stretch of 35 conserved amino acids upstream of the core common sequence. This motif, called RLI, is extremely polar (11 basic and 4 acidic residues) and is also found at the N terminus of another group of orthologous sequences, namely COG1245.

**Expression in E. coli of two engineered SSO0551 constructs**

From multiple sequence alignments, SSO0551 should encode a 166 amino acid polypeptide. An N-terminal 6His tagged recombinant construct (pSBTN-AB31) was engineered. As we could not exclude that the 42 amino acids extension at the N terminus was not an annotation artifact, we intended to check experimentally whether this putative extension could have some influence on SSO0551. A second construct (pSBTN-AB30) was simultaneously engineered supposedly allowing production of a 26 kDa N-terminal 6His variant. Unexpectedly, no major difference in expression was detected between the two cellular extracts when they were resolved on SDS-PAGE. Two overexpressed products with both an apparent molecular weight of approximately 20 kDa were obtained upon addition of IPTG (data not shown). Fingerprint identification of these two products was carried out by trypsin proteolysis and mass spectrometry. Table 1 shows the MALDI-TOF mass measurements recorded for the two samples. The tryptic peptides that were detected revealed that both products correspond to native SSO0551

<footnote><sup>a</sup>The two products were isolated from SDS-PAGE and treated with trypsin. The resulting peptides were analyzed by MALDI-TOF spectrometry in reflectron mode. The expected length for the two products are: 180 amino acids for SSO0551 with the 6His-tag at the N terminus (pSBTN-AB31 construct) and 222 amino acids for SSO0551 extended version (42 additional amino acids at the N terminus and a 6His-tag; pSBTN-AB30 construct).</footnote>

<footnote><sup>b</sup>Peptide positions refer to SSO0551.</footnote>

<footnote><sup>c</sup>«nd» denote the peptides not detected in one of the samples.</footnote>
<figure><img src="image_2.png" /></figure>

B

I

II

M. barkeri - 48838408 68-GI:VLDCSWE:EVER:VF:PD:LAK:NL-K:HRA:LY:LL:G:NP:VN:F:GR:F:PK:F:KL:N:S:A:FA:A:FA:A:Y:IL:G:H:KE:A:EK:VLS:F:K:F:N:G:H:S:F:LE:N:K:F:PE:EY-156

M. kandleri - 20094076 61-GL:VA:DCSWE:H:V:HRY-F:G:P:LRR:RC-R:HRI:L:P:YL:A:A:NP:VNY:G:P:C:KL:ST:VE:A:A:A:Y:IL:G:F:R:EA:E:F:IS:F:K:W:G:P:A:F:LE:N:R:R:EAY-148

F. acidarmanus - 48852409 62-G:IV:VID:G:SW:N:R:INT---I:EN:IK:G:K:N:G:KK:L:P:LL:P:LV:N:PV:NY:G:P:K:G:L:S:LE:A:A:A:Y:IL:G:Y:SE:A:LG:D:V:SI:K:Y:SW:A:Q:N:F:K:T:N:IN:P:LN:DY-147

S. solfataricus - 38605619* 58-G:IT:V:DT:SW:N:T:SQ:S--EF:KN:RG-E:HRI:P:IL:F:A:G:N:PI:H:Y:G:IA:Y:KL:S:SI:E:A:LI:AT:LY:I:V:DE:E:A:IK:LS:N:VVK:W:G:T:H:F:IEL:N:K:F:EL:EY-144

G. lamblia - 29250359 179-GL:A:V:D:CSW:N:E:IL:AD:Q:V:P:T:K:L:C:N:H:R:L:P:LL:F:LL:A:SN:VNY:G:P:F:K:L:N:C:A:FA:A:G:L:W:IC:GY:R:D:Q:A:D:T:I:Q:F:F:Y:A:K:A:F:L:D:N:E:LL:D:RY-268

S. cerevisiae - 6324579 106-G:AS:V:E:CSWAR:LEE--VP:F:N:K:IG:K:G:H:ER:L:P:YL:A:VA:A:NV:NY:G:P:RW:L:N:C:VE:A:A:C:F:A:V:G:R:M:D:W:A:ELL:S:H:F:SW:G:M:F:LE:N:K:F:EL:EY-193

H. sapiens - 42661007 57-G:V:A:VID:CSWAR:LEE--TP:F:K:M:RG:SH:L:L:P:YL:A:VA:A:NP:VNY:G:P:RY:L:SC:VE:A:F:A:T:F:C:V:GF:P:DL:A:V:I:L:R:K:F:KW:G:K:F:LD:N:R:QL:L:D:K:Y-144

A. thaliana - 42573329 104-GL:V:V:D:CSWAR:LTD--VP:F:A:LR:CT:A:P:R:LL:P:WL:A:VA:A:NP:VNY:G:P:C:EL:SC:VE:A:SA:A:AL:I:L:C:GE:K:T:A:E:LL:IG:K:F:KW:G:H:A:F:LS:N:KD:I:L:KEY-191

**Figure 1**

**COG2042 sequences comparison.** A – Phylogenetic analysis of archaeal and eukaryal COG2042 sequences. Archaeal and eukaryal homologs were obtained from public databases [54] by PSI-BLAST searches. To get the most conserved alignment between COG2042 polypeptide sequences, most probable start codons should be considered as atg at nucleotide 500978 on Crick strand for MTH554 from *Methanothermobacter thermoautotrophicus* str. DeltaH (NC_000916), gtg at nucleotide 1526308 on Watson strand for Vng2075c from *Halobacterium* sp. NRC-1 (NC_002607), atg at nucleotide 8398 on Watson strand for Mbur141901 from *Methanococcoides burtonii* DSM6242 (NZ_AADH01000008) and atg at nucleotide 484790 on Crick strand for SSO0551 from *S. solfataricus* (NC_002754). Asterisks indicate modified protein sequences according to this new proposed annotation. Multiple alignments were performed by ClustalW. Following removal of a few ambiguously aligned regions, a data set was assembled comprising 40 sequences over 162 amino acid positions. An unrooted evolutionary distance tree was constructed based on Kimura distances and neighbor joining tree reconstruction algorithm. Bootstrap confidence levels at nodes were computed by the Phylips package with 400 replicates. Scale bar represents unit of amino acid substitutions per position. Accession numbers (gi) are indicated beside the organism. B – Conserved sequence blocks in the alignment of COG2042 members. Based on the phylogenetic analysis (Panel A), eight representative sequences were selected out of 40 COG2042 sequences. Four first sequences are from Archaea while last four sequences are from Eukaryota. SSO0551 sequence (gi38605619) from *S. solfataricus*, labeled with an asterisk, has been numbered according to the new annotation proposed in RESULTS section. Invariant residues in the eight sequences are shown in red and conserved residues in blue. Two conserved motifs commented in the results section are indicated with grey boxes.
sequence. From the 6His-SSO0551 product (pSBTN-AB31 construct), thirteen peptides map with the theoretical sequence (57 % sequence coverage). Noteworthy, a peptide (1590.64 amu) was attributed to part of the 6His-modified N terminus (Table 1). The twelve peptides recorded from the 6His-SSO0551 extended version (pSBTN-AB30 construct) fit only to the C terminus of the theoretical construct (43 % sequence coverage). These results along with low molecular weight observation on SDS-PAGE indicate that probably a truncated protein was obtained during expression of the ORF comprising the 126 nt 5'-extension (42 additional amino acids at the N-terminus). This product, corresponding in fact to untagged SSO0551 as confirmed hereafter with purified product, showed no binding on Ni-NTA chromatography. This observation is in agreement with absence of 6His tag at the N terminus.

**Recombinant SSO0551 is structured, thermostable and monomeric**

Crude extract containing native untagged SSO0551 polypeptide from *E. coli* Rosetta(DE3)(pLysS)(pSBTN-AB30) cells was heated at different temperatures. Proteins that remained soluble were analyzed on SDS-PAGE. Most of *E. coli* contaminants were removed by such treatment. SSO0551 polypeptide remained soluble even when cell extract was heated to 80°C and therefore this protein was considered as thermostable. This protein was purified to homogeneity by a three-step purification protocol. A 20 min heat treatment at 70°C (Fig. 2A, lane 3), followed by a Resource-S ion exchange chromatography (Fig. 2A, lane 4) and a Superdex75 gel filtration (Fig. 2A, lane 5), yielded approximately 1.6 mg of pure protein per L of culture. Purified protein was subjected to MALDI-TOF mass analysis. Fig. 2B shows the spectrum recorded. The experimental *m/z* of 19,198 measured for the monocharged polypeptide matches perfectly with theoretical average mass of native untagged SSO0551 protein (average mass of 19,197 Da). This measurement unequivocally confirmed that a truncated protein is produced using *E. coli* Rosetta(DE3)pLysS transformed with pSBTN-AB30. Both SDS-PAGE (Fig. 2A, lane 5) and MALDI-TOF spectrum (Fig. 2B) testify for homogeneity of the sample.

Content of secondary structure elements in SSO0551 was estimated by far-UV circular dichroism. Fig. 3 shows the spectrum recorded at 20°C. Purified protein presents negative ellipticity in the near-UV with minima at 208 (-14.7 10³ deg cm² dmol⁻¹) and 222 nm (-12.7 10³ deg cm² dmol⁻¹). Deconvolution of the CD spectrum leads to an estimation of secondary structural element content of about 28-29 % of α-helices and 14-16 % of β-sheets using K2D neural-software. Predictions of SSO0551 secondary structures by PSIPRED and Jpred web servers gave values of 10-11 % of β-sheets in relative agreement with the circu-

lar dichroïsm data, but overestimated the α-helix average content (54 %). PSIPRED and Jpred predictions are based on neural networks trained on known folds. The overestimation of the α-helix content may be due to the novel fold of these COG2042 proteins as discussed here below.

Native molecular mass of SSO0551 was determined by size-exclusion chromatography on a Superdex 200 HR10/30 calibrated column. Pure protein eluted as a peak centered at 39.1 mL in the assay conditions corresponding to an apparent molecular mass lower than 20 kDa. This elution profile indicates that this structured protein behaves as a compact monomer.

**Limited proteolysis defines Glu<sup>73</sup>-Arg<sup>78</sup> as a hyper-sensitive region**

Purified SSO0551 protein was subjected to limited proteolysis with various endopeptidases (trypsin, chymotrypsin, ArgC and GluC). MALDI-TOF mass spectrometry was used to determine cleavage sites by following the time course generation of peptides. Several protease/substrate ratios were assessed to confirm which preferential sites on entire protein were first attacked (earliest cleavage), thus corresponding to a native state of the protein. The two fragments generated by such cleavage may be more vulnerable to subsequent attacks than native protein and therefore late proteolytic sites are considered less informative. Both small and large peptides generated during proteolysis were evaluated. Partial proteolyzed products obtained with trypsin were first resolved by reverse-phase chromatography and analyzed by MALDI-TOF mass spectrometry. Results recorded from direct analysis of the digestions without prior separation were almost similar to those obtained with separation. Therefore, the latter cost-effective strategy was used for analyzing the numerous conditions tested. Figure 4 shows the MALDI mass spectrum of the main large products obtained from a tryptic digest of SSO0551 (enzyme/protein ratio of 1:20) after 60 sec of reaction. In these conditions, the signal of intact protein was still visible at *m/z* 19198.4, but mixed with signals corresponding to 8 different large fragments. Among these, 7 peptides arose from an N-terminal proteolysis: [32-166], [35-166], [56-166], [57-166], [76-166], [79-166], [101-166] (Fig. 4). Such peptidic profile indicates that SSO0551 N terminus is rather solvent exposed in comparison to C terminus.

During the earliest events of the trypsin proteolysis analyzed in various conditions for detection of large products but also smaller peptides, monocharged cations with following *m/z*: 8614.6 amu, 10603.4 amu, 6489.8 amu, and 12724.1 amu, were attributed to fragments [1-75] (Δmass: -178 ppm), [76-166] (+89 ppm), [1-56] (+145 ppm), [57-166] (+200 ppm), respectively (data not shown). These data clearly indicate that Lys<sup>75</sup> and Arg<sup>56</sup>
<figure><img src="image_3.png" /><figcaption>**Figure 2**
Analysis of untagged SSO0551 purification. A – Various purification steps followed by SDS-PAGE. SDS-PAGE was performed on a 4–12% gradient Novex bis-tris acrylamide gel (Invitrogen) and stained with BlueSafe Coomassie solution (Biorad). Lane M, molecular weight markers; Lane 1, cell-free extract of E. coli Rosetta(DE3)(pLysS)(pSBTN-AB30) ; Lane 2, soluble proteins from cell-free extract; Lane 3, soluble proteins after the 60°C heat treatment; Lane 4, Resource-S eluate ; Lane 5, Superdex75 eluate. Bands corresponding to recombinant SSO0551 polypeptide are indicated with the arrow. B – SSO0551 mass measurement by MALDI-TOF mass spectrometry. Measurements from 250 laser shots of purified protein cristallized with sinapinic acid were averaged.</figcaption></figure>

are two sites of early cleavage by trypsin. Identification of
peptides Val³²-Lys¹⁶⁶ (15478.7 amu, -53 ppm) and Gly³⁵-
Lys¹⁶⁶ (15191.2 amu, +153 ppm) also indicates that Arg³¹
and Lys³⁴ could be two other initial nick-sites.

Similar experiments with endoproteinase Arg-C resulted
in observation of two pairs of complementary peptides
with *m/z* of 1920.9 amu ([1-15] +70 ppm) and 17296.8
amu ([16-166], -95 ppm) on one hand, 8998.1 amu ([1-
78], -176 ppm) and 10217.5 amu ([79-166], +332 ppm)
on the other hand. These data indicated that Arg⁷⁸ and
Arg¹⁵ are the main proteolyzed sites when ArgC enzyme
was used. Chymotrypsin attacks SSO0551 native protein
mainly at Phe⁷⁴ because two complementary peptides

with *m/z* of 8487.0 amu ([1-74], -249 ppm) and 10734.2
amu ([75-166], -157 ppm) were clearly evidenced. Glu⁷³
is the main proteolyzed site when GluC protease was
used, as peptides with *m/z* of 8338.7 amu ([1-73], -118
ppm) and 10880.0 amu ([74-166], +28 ppm) were
detected. For all these analysis, smaller peptidic fragments
that accumulated over time could be attributed from fur-
ther proteolysis of the products arising from initial attacks
(data not shown). All these results are concordant in
pointing out that Glu⁷³-Arg⁷⁸ and Glu²⁸-Arg³¹ are two
accessible solvent-exposed regions of the protein as they
can be proteolyzed by several endopeptidases, the first
cited being definitively hyper-sensitive. Local unfolding
not just surface exposure is necessary for efficient *in vitro*
<figure><img src="image_4.png" /><figcaption>**Figure 3**
Evaluation of SSO0551 secondary structure elements by circular dichroïsm. The molar ellipticity was calculated on the basis of exact amino acid composition of SSO0551 recombinant product.</figcaption></figure>

proteolysis because the polypeptide segment being cleaved must form a specific structure with the associated protease [30]. For this reason, Glu<sup>73</sup>-Arg<sup>78</sup> region should also correspond to a flexible region, i.e. a protruding loop.

**Lysine labeling with NHS-biotin and DTSSP cross-linking confirm that the N terminus is rather solvent-exposed**

The SSO0551 protein contains 21 lysine residues (12%) distributed along the whole polypeptide sequence. Under mild conditions that should keep the native conformation of the protein, specific labeling of these residues with NHS-biotin may give further details about their respective surface accessibility and/or their interactions with other residues [31]. After reaction with various amount of chemical reagent (molar ratio NHS-Biotin/total lysines of 1:40, 1:20, 1:10, 1:2, 1:1, 2:1), protein labeling was monitored by determining the mass of undigested samples. Figure 5 shows the signals measured by MALDI-TOF mass spectrometry for four of these ratios. The fact that some unmodified protein is still present at ratio below 1:20 testifies for mild conditions that should allow modification of protein still in a native state. As expected with NHS-biotin, each peak exhibits the predictable mass increment

(average mass of 226.3 amu per label). Figure 5 shows that at molar ratio of 1:40 a simple modification is obtained, while a more heterogeneous population was detected for higher ratio. For examples, 1 to 3 modifications are detected at ratio 1:20, 2 to 5 modifications at ratio 1:2. However, a limited number of modifications (8–10) are recorded for higher ratio, indicating that among the 21 lysine residues only a fraction is accessible to the chemical.

To localize all labeled residues, NHS-biotin treated samples were subsequently subjected to proteolysis with various endoproteases (trypsin, Arg-C, or Glu-C) and compared to untreated samples. SSO0551 sequence coverage was estimated to be 92% with all 21 lysine residues included in this coverage. Peptides (Δmass below 120 ppm) detected with NHS-biotin treated samples but not detected with untreated samples are listed in Table 2. Using limiting amount of NHS-biotin (molar ratios of 1:10, 1:20 or 1:40), nine reactive residues are unequivocally identified: Lys<sup>10</sup>, Lys<sup>14</sup>, Lys<sup>20</sup>, Lys<sup>23</sup>, Lys<sup>25</sup>, Lys<sup>51</sup>, Lys<sup>75</sup>, Lys<sup>128</sup>, and Lys<sup>154</sup>, assuming that proteases do not cleave after a modified residue. Other residues, such as
<figure><img src="image_5.png" /><figcaption>**Figure 4**</figcaption></figure>

**MALDI-TOF** mass spectrum of large peptide fragments produced by SSO0551 partial trypsin proteolysis. A partial enzymatic proteolysis was carried out during 60 sec at 20°C with a trypsin/SSO0551 protein ratio of 1/20 (w/w). The products were then resolved onto a C8 reverse phase chromatographic column and the different UV absorbing fractions were analyzed by MALDI-TOF. The m/z spectrum obtained with the fraction eluting at 40–50 % acetonitrile is shown. The asterisk labels a peak arising from trypsin autolysis. Peaks that could be assigned are identified with numbers (experimental m/z, residues, Δmass in ppm compared to theoretical [M+H]$^+$ average mass): **M1** (19198.4, [1–166], -1), **M2** (15478.7, [32–166], -53), **M3** (15191.2, [35–166], +153), **M4** (12883.4, [52–162] or [56–166], -46 or -43), **M5** (12724.2, [57–166], +193), **M6** (10603.8, [76–166], +51), **M7** (10220.1, [79–166], +77), **M8** (9257.2, [83–162], +63) and **M9** (7717.0, [101–166], -134). Peptides **C1** and **C2**, complementary of **M5** and **M6**, were not observed in this spectrum but in another fraction from the C8 reverse-phase chromatography corresponding to smaller peptides.
A

<figure><img src="image_6.png" /></figure>

B

<figure><img src="image_7.png" /></figure>

**Figure 5**
MALDI-TOF mass spectra of proteins labeled with NHS-biotin. Protein samples were mixed with various amount of NHS-biotin for 15 min at room temperature. Unlabeled protein mass spectrum is shown in grey solid line. Reagent/total lysine ratios were in Fig. 5A: 1:40 (black dotted line) and 1:20 (black dashed line), and in Fig. 5B: 1:2 (black solid line) and 2:1 (grey dotted line). The number of labeled lysines is indicated above each peak.

Lys³⁴ and Lys⁴⁹ might be also labeled (Table 2). The number of labeled lysines is in agreement with the limited number of modifications recorded at higher ratio. Remarkably, spectra of whole peptide mixture were informative enough to give assignment of all modified peptides without the need of a purification step. Therefore, other amine reactive reagent that creates a mass shift could have been used.

Using a lysine cross-linking reagent, DTSSP, it is possible to assess intra- or inter-molecular protein contacts [21,32]. DTSSP enables cross-linking of amino groups up to 12 Å apart. As SSO0551 was shown to be monomeric and its concentration used in the assays was low (2.5 pM), intramolecular cross-links should be favored over inter-molecular cross-links. In addition, the low reagent concentration used should avoid unwanted conformational changes that may be induced by multiple intramolecular cross-linking. After reaction with DTSSP, products were subjected to trypsin proteolysis and peptides were identified by MALDI-TOF. As the protein is relatively small, mass signals could be attributed with a good confidence (tolerance < 120 ppm). In addition, peak attribution was always confirmed upon reduction of products and sometimes through redundancy due to miss-cleavage.

SSO0551 sequence coverage was 89%. The monoisotopic cations at *m/z* 1169.55, 2077.04, 2715.11 and 2871.14, detected for SSO0551 treated by DTSSP (ratio DTSSP/total polypeptide of 20:1) correspond in mass to addition of a DTSSP moiety on peptides [24-31] (Δmass tolerance: +16 ppm), [123-138] (-7 ppm), [57-78] (+61 ppm) and [56-78] (+82 ppm). Since trypsin does not cleave after a modified lysine, we conclude that Lys²⁵, Lys⁷⁵ and Lys¹²⁸ were modified. After DTT treatment, peaks corresponding to the expected products (-103.993 amu theoretically) were detected at *m/z* 1065.46 (-109 ppm), 1973.03 (0 ppm), 2611.21 (-28 ppm) and 2767.26 (-45 ppm), respectively. Fig. 6 shows two monoisotopic [MH]⁺ ions at *m/z* 1491.76 and 1835.84 corresponding to intrapeptide cross-linked peptides: [21-31] (+31 ppm) and [3-15] (+3 ppm), respectively. These peptides contain two proximal lysine residues (Lys²³-Lys²⁵ and Lys¹⁰-Lys¹⁴). As shown in Fig. 6, these two peaks were absent in mass spectrum following DTT reduction but new peaks at *m/z* 1493.69 and 1837.81 appears at the expected increment (+2.016 amu theoretically). An additional peak at *m/z* 2502.22 could be attributed to peptide [35-55] (+37 ppm) with an intrapeptide cross-link between Lys⁴⁹ and Lys⁵¹. However, the corresponding reduced peak was not detected. Strikingly, every lysines that were reactive with DTSSP were
<table><caption>Table 2: Monoisotopic [M+H]$^+$ peptides generated by various proteases after NHS-biotin labeling of SSO0551.</caption><thead><tr><th>Protease</th><th>[M+H]<sup>+</sup><br>(observed)</th><th>[M+H]<sup>+</sup><br>(expected)</th><th>Δmass (ppm)</th><th>Position<br>(Start-End)</th><th>Peptide assignment<sup>a</sup></th><th>NHS-biotin<br>label<sup>b</sup></th><th>Modified residues<sup>a</sup></th></tr></thead><tbody><tr><td>Trypsin</td><td>1731.94</td><td>1731.84</td><td>-57</td><td>3-14</td><td>VYIIDYHKDDPK</td><td>1</td><td>K<sup>10</sup></td></tr><tr><td>Trypsin</td><td>2114.27</td><td>2114.02</td><td>118</td><td>3-15</td><td>VYIIDYHKDDPKR</td><td>2</td><td>K<sup>10</sup>andK<sup>14</sup></td></tr><tr><td>Trypsin</td><td>1888.15</td><td>1887.94</td><td>111</td><td>3-15</td><td>VYIIDYHKDDPKR</td><td>1</td><td>K<sup>10</sup>orK<sup>14</sup></td></tr><tr><td>Trypsin</td><td>1991.21</td><td>1990.98</td><td>-119</td><td>1-14</td><td>MKVYIIDYHKDDPK</td><td>1</td><td>K<sup>2</sup>orK<sup>10</sup></td></tr><tr><td>Trypsin</td><td>2147.33</td><td>2147.08</td><td>-120</td><td>1-15</td><td>MKVYIIDYHKDDPKR</td><td>1</td><td>K<sup>2</sup>orK<sup>10</sup>orK<sup>14</sup></td></tr><tr><td>Trypsin</td><td>2163.33</td><td>2163.07</td><td>-120</td><td>1-15</td><td>M<sup>ox</sup>KVYIIDYHKDDPKR</td><td>1</td><td>K<sup>2</sup>orK<sup>10</sup>orK<sup>14</sup></td></tr><tr><td>Trypsin</td><td>2373.41</td><td>2373.16</td><td>-107</td><td>1-15</td><td>MKVYIIDYHKDDPKR</td><td>2</td><td>K<sup>2</sup>or/andK<sup>10</sup>or/andK<sup>14</sup></td></tr><tr><td>Trypsin</td><td>1180.77</td><td>1180.69</td><td>-64</td><td>20-25</td><td>KLVKLK</td><td>2</td><td>K<sup>20</sup>andK<sup>23</sup></td></tr><tr><td>Trypsin</td><td>1203.76</td><td>1203.66</td><td>-87</td><td>24-31</td><td>LKIAEFTTR</td><td>1</td><td>K<sup>25</sup></td></tr><tr><td>Trypsin</td><td>1770.16</td><td>1769.98</td><td>102</td><td>21-31</td><td>LVVKLKIAEFTTR</td><td>2</td><td>K<sup>23</sup>andK<sup>25</sup></td></tr><tr><td>Trypsin</td><td>1898.15</td><td>1898.08</td><td>-39</td><td>20-31</td><td>KLVKLKIAEFTTR</td><td>2</td><td>K<sup>20</sup>orK<sup>23</sup>orK<sup>25</sup></td></tr><tr><td>Trypsin</td><td>1127.68</td><td>1127.60</td><td>-71</td><td>50-56</td><td>DKDIVRR</td><td>1</td><td>K<sup>51</sup></td></tr><tr><td>Trypsin</td><td>3064.63</td><td>3064.63</td><td>1</td><td>32-55</td><td>VGKGVVL...NKDKDIVR</td><td>2</td><td>K<sup>34</sup>or/andK<sup>49</sup>or/andK<sup>51</sup></td></tr><tr><td>Trypsin</td><td>2554.57</td><td>2554.37</td><td>-78</td><td>35-55</td><td>GVVL...NKDKDIVR</td><td>1</td><td>K<sup>49</sup>orK<sup>51</sup></td></tr><tr><td>Trypsin</td><td>2710.61</td><td>2710.47</td><td>-51</td><td>35-56</td><td>GVVL...NKDKDIVRR</td><td>1</td><td>K<sup>49</sup>orK<sup>51</sup></td></tr><tr><td>Trypsin</td><td>2838.65</td><td>2838.55</td><td>-34</td><td>32-55</td><td>VGKGVVL...NKDKDIVR</td><td>1</td><td>K<sup>34</sup>orK<sup>49</sup>orK<sup>51</sup></td></tr><tr><td>Arg-C</td><td>2838.75</td><td>2838.55</td><td>-69</td><td>32-55</td><td>VGKGVVL...NKDKDIVR</td><td>1</td><td>K<sup>34</sup>orK<sup>49</sup>orK<sup>51</sup></td></tr><tr><td>Trypsin</td><td>2749.36</td><td>2749.36</td><td>0</td><td>57-78</td><td>IGIT...SEFKNIR</td><td>1</td><td>K<sup>75</sup></td></tr><tr><td>Arg-C</td><td>2749.60</td><td>2749.36</td><td>-87</td><td>57-78</td><td>IGIT...SEFKNIR</td><td>1</td><td>K<sup>75</sup></td></tr><tr><td>Glu-C</td><td>1089.56</td><td>1089.55</td><td>-8</td><td>74-80</td><td>FKNIRGE</td><td>1</td><td>K<sup>75</sup></td></tr><tr><td>Trypsin</td><td>2111.22</td><td>2111.11</td><td>-52</td><td>123-138</td><td>LSNVVKWGHTFIELNK</td><td>1</td><td>K<sup>128</sup></td></tr><tr><td>Glu-C</td><td>2068.19</td><td>2068.10</td><td>-41</td><td>120-135</td><td>AIKLSNVVKWGHTFIE</td><td>1</td><td>K<sup>122</sup>orK<sup>128</sup></td></tr><tr><td>Trypsin</td><td>1330.74</td><td>1330.67</td><td>-55</td><td>146-154</td><td>NKTEEDIKK</td><td>1</td><td>K<sup>147</sup>orK<sup>153</sup></td></tr><tr><td>Glu-C</td><td>1256.65</td><td>1256.67</td><td>13</td><td>151-158</td><td>DIKKIERE</td><td>1</td><td>K<sup>153</sup>orK<sup>154</sup></td></tr><tr><td>Trypsin</td><td>1486.85</td><td>1486.76</td><td>-60</td><td>148-157</td><td>TEEDIKKIEREIIEK</td><td>1</td><td>K<sup>153</sup>orK<sup>154</sup></td></tr><tr><td>Trypsin</td><td>1383.76</td><td>1383.77</td><td>5</td><td>154-162</td><td>KIEREIIIEK</td><td>1</td><td>K<sup>154</sup></td></tr><tr><td>Trypsin</td><td>1867.14</td><td>1867.07</td><td>-36</td><td>154-166</td><td>KIEREIIIEKILEK</td><td>1</td><td>K<sup>162</sup>orK<sup>154</sup></td></tr><tr><td>Trypsin</td><td>1340.71</td><td>1340.75</td><td>29</td><td>158-166</td><td>EIIIEKILEK</td><td>1</td><td>K<sup>162</sup>orK<sup>166</sup></td></tr><tr><td>Trypsin</td><td>1739.01</td><td>1738.98</td><td>-19</td><td>155-166</td><td>IEREIIIEKILEK</td><td>1</td><td>K<sup>162</sup>orK<sup>166</sup></td></tr></tbody></table>

also detected by NHS-biotin labeling. Most modifications are located at the N terminus of the protein where five modified residues belong to the RLI motif.

Discussion

Although COG2042 proteins are distributed among a large number of organisms, no experimental evidences have yet been reported concerning their biochemical characterization and function. As they are not related, even remotely, to any other family of proteins, COG2042 members can be phylogenetically considered as orphans. Figure 7 (Panel A) summarizes the structural information obtained with chemical modification approach, in combination with limited proteolysis procedures. Using MALDI-TOF mass spectrometry to identify protease-accessible sites, we have shown that the most exposed regions

are located at the first half of the protein, the Glu<sup>73</sup>-Arg<sup>78</sup> region being revealed hyper-sensitive to various proteases (Fig. 7A). It probably indicates a protruding loop out of the globular protein. This charged region is relatively conserved among COG2042 orthologs and lies between two highly conserved segments of COG2042 (motif I and II as shown on Fig. 7). Chemical modification agrees with limited proteolysis in that the RLI motif is solvent exposed while the C terminus appeared rather inaccessible (Fig. 7A). The length of the RLI motif, first defined by conserved domain search [33], matches perfectly with two sensitive proteolytic sites (Arg<sup>31</sup> and Lys<sup>34</sup>). The RLI domain is also present at the N terminus of another group of orthologous proteins, namely COG1245. Remarkably, COG1245 proteins only occur in two domains of life (Archaea and Eukarya) similarly to COG2042 proteins.

<footnote>ªUnequivocal modified residues are indicated in bold characters.</footnote>

<footnote>ᵇNHS-biotin label (+226.08 amu).</footnote>
<figure><img src="image_8.png" /><figcaption>**Figure 6**

**DTSSP cross-linked peptides identified by MALDI-TOF mass spectrometry.** Protein samples were mixed with DTSSP for 30 min at room temperature at a molar ratio 1:20 (protein:DTSSP). They were then subjected to trypsin proteolysis. Masses of peptides treated with DTT (grey spectrum) or untreated (black spectrum) were measured by MALDI-TOF mass spectrometry. Peptide sequences are indicated and DTSSP cross-links or DTSSP moieties arising from DTT treatment are depicted schematically. Monoisotopic masses of protonated peptides [MH⁺] are theoretically: 1491.807 amu and 1493.822 amu for LVKLKIAEFTR [21-31] cross-linked with one DTSSP molecule (untreated (C₆₈H₁₁₅N₁₆O₁₇S₂) and DTT-reduced (C₆₈H₁₁₇N₁₆O₁₇S₂), respectively), 1835.846 amu and 1837.862 amu for VYIIDYHKDDPKR [3-15] cross-linked with one DTSSP molecule (untreated (C₈₂H₁₂₃N₂₀O₂₄S₂) and DTT-reduced (C₈₂H₁₂₅N₂₀O₂₄S₂), respectively).</figcaption></figure>

Although co-occurrence of protein members is not strictly identical (for example, pyrococci encompass the information for COG1245 but not for COG2042 polypeptides), such occurrence pattern may reflect a functional link between the two protein families.

Our initial objective was to obtain about SSO0551 as much low-resolution structural information as possible in order to discriminate among putative three-dimensional models representing COG2042 protein structure. However, currently available threading tools applied on SSO0551 failed to detect any structurally related-proteins. Alternatively, we obtained ten different *ab initio* models of SSO0551 using the fully-automated ROBETTA server based on ROSETTA procedures [34]. On these ten models, we applied all the low-resolution structural information

gathered in this work. We predicted for every model location of preferential proteolytic sites using the Nick-Pred software [35]. Models M1, M2 and M6 on one hand, and M9 and M10 on the other, show hypersensitive regions in the RLI motif or C terminus, respectively. These features do not correspond to our experimental data. Only models M4, M7 and M8 predict that the loop Glu<sup>73</sup>-Arg<sup>78</sup> is solvent exposed (data not shown). Among these three models, M4 and M8 respect the ranking of preferential nick-sites for trypsin, chymotrypsin, ArgC and GluC proteases. Solvent accessibility for lysine side chain was evaluated for models M4, M7 and M8 and compared with experimental data (data not shown). All the lysine residues labeled with NHS-biotin are found solvent-exposed in model M8. Manual inspection of cross-linked lysines (Lys<sup>10</sup>-Lys<sup>14</sup> and Lys<sup>23</sup>-Lys<sup>25</sup>) revealed that model M4 is
<figure><img src="image_9.png" /><figcaption>**Figure 7**
Experimental topology information obtained on SSO0551 and theoretical 3D model M8 shown to be compatible. A – Schematic representation of experimental low-resolution structural information. SSO0551 sequence is displayed with its predicted secondary structural elements (α-helices and β-strands are represented by helices and arrows, respectively). The two conserved COG2042 motifs described in this study (I and II) and the 35 amino acids RLI motif are specified. Vertical arrows indicate sites of limited proteolysis that have been experimentally determined. Lysines labeled with NHS-biotin are indicated with asterisks and lysines involved in DTSSP cross-links are schematically bridged. Sequence numbering refers to the native version (untagged) of SSO0551 experimentally characterized. B – Schematic drawing of ab initio model M8. Secondary structure elements are shown as coil for α-helices, arrows for β-strands, and lines for loops. The color code follows a rainbow style where N terminus is in blue and C terminus in red. C – Molecular surface of ab initio model M8. The localization of solvent-exposed lysines, as determined by NHS-biotin labeling experiments, is colored in red whereas other lysines are colored in blue. This view shows the protein in the same orientation as in B.</figcaption></figure>

not valid because of the opposite orientation of Lys¹⁰ and Lys¹⁴. Figure 7 (Panels B & C) shows cartoon views of the M8 model that fulfills all our experimental constraints. For this model, the distance between the two reactive amine groups of Lys¹⁰-Lys¹⁴ and Lys²³-Lys²⁵ pairs are 12.7 Å and 13.3 Å, respectively. Search with DALI for structural homologs using model M8 did not result in significant

scores with any known PDB structures. This is consistent with the PSI-BLAST results and may indicate that COG2042 proteins share a novel fold. COG2042 proteins are thus a target of choice for genomic structural studies.

In conclusion, we have presented a strategy consisting in obtaining low-resolution structural information (deter-
mination of nick-sites, solvent exposed residues, and residue-residue distances) that can be used to distinguish among a large set of theoretical molecular models. Lack of remotely-related structural templates or lack of adequacy between experimental data and most theoretical models indicates that such family of proteins should become a priority in structural genomic projects.

## Methods

### Chemical and biological reagents

Most chemicals used in this study were obtained from Sigma and were of analytical grade. Oligonucleotide primers were purchased from Genset. N-hydroxysuccinimide-biotin (NHS-biotin) and 3,3'-dithio-bis [sulfosuccinimidyl-propionate] (DTSSP) were obtained from Pierce. Matrices for Matrix-assisted Laser Desorption Ionization-Time of Flight (MALDI-TOF) mass spectrometry and calibration standards were purchased from Bruker Daltonics. Sequencing grade proteolytic enzymes were from Roche Applied Science.

### Cloning and overexpression of SSO0551

Two constructs were designed in order to get overexpression of the SSO0551 ORF (starting with an ATG codon at nucleotide 484790 on the Crick strand of S. solfataricus P2 genome (NC_002754)) and an N-terminal extended version of SSO0551 (starting with an ATG codon at nucleotide 484916). For both proteins, an N-terminal 6His tag was added to render the purification of the recombinant products easier. For this purpose, synthetic oligonucleotide primers were oAB22 (5'-gctagcATGAAGCCCAAACCC-3') and oAB49 (5'-gctagcATGAAGGTATATATTATAGAC-3') that both contain an engineered NheI site, oAC34 (5'-cggatcctacTATTTTCAAGTATTTTC-3') and oAE62 (5'-ggatcctcaTCATTTTCA AGTATTTCTC-3') that both contain an engineered BamHI site (restriction sites underlined in the primer sequences and nucleotides not present in the original sequence shown by lower case). Oligonucleotide pairs oAB22/oAC34 and oAB49/oAC34 were used for two distinct PCR amplifications of SSO0551 with S. sulfolobus total DNA as template. A 643-bp fragment (N-ter 6His-tag extended version of SSO0551) and a 517-bp fragment (N-ter 6His-tag SSO0551) were obtained, respectively. They were cloned into pCRScript-cam (Stratagene), resulting in plasmids pSBTN-AB36 and pSBTN-AB37, respectively. The two inserts were removed by digestion with NheI and BamHI and ligated with T4 DNA ligase into plasmid pSBTN-AB23 (Armengaud J. & Chaumont V., unpublished data), a derivative of pCR T7/NT-topo (Invitrogen) containing a T7 promoter and 6 His-tag, previously digested with the same endonucleases. The resulting plasmids pSBTN-AB30 and pSBTN-AB31, respectively, were verified by DNA sequencing in order to ascertain the integrity of the nucleotide sequence. Hyper-

expression of the recombinant SSO0551 constructs was achieved with E. coli Rosetta(DE3)pLysS strain (Novagen), freshly transformed with the plasmids described above. Cultures were carried out at 30°C as described earlier [6].

### Purification of recombinant SSO0551 protein

The purification of recombinant SSO0551 was performed from 44 g (wet material) packed cells. Buffer A consisted of 50 mM K₂HPO₄/KH₂PO₄ buffer (pH 7.2) containing 400 mM K-glutamate. The pellet was thawed on ice and resuspended in 120 mL of buffer A. The cells were disrupted by sonication with a total energy delivered of 71 kJ. The cell-extract was then centrifuged at 30,000 g for 20 min at 4°C to remove cellular debris and aggregated proteins. The supernatant was subjected to a 20 min heat treatment using a water bath maintained at 70°C, and immediately centrifuged a second time at 30,000 g for 20 min at 4°C. Chromatographic steps were performed at room temperature using an Äkta Purifier FPLC system (Amersham Biosciences). The 135 mL supernatant was applied at a flow rate of 2.8 mL/min onto a XK 26 × 20 column (Amersham Biosciences) containing 50 mL of Chelating Sepharose Fast Flow (Amersham Biosciences) and previously loaded with 200 mM NiSO₄, washed with milliQ water and equilibrated with Buffer A containing 50 mM imidazole. The fraction collected during the IMAC loading was shown to contain the SSO0551 protein. This 222 mL fraction was concentrated to a volume of 56 mL by means of Centricon Plus-20 filtration units (Millipore) and then dialyzed overnight at 4°C against 20 mM K₂HPO₄/KH₂PO₄ buffer (pH 7.2) containing 20 mM NaCl (buffer B). The 78 mL supernatant obtained after centrifugation at 30,000 g for 10 min at 4°C was divided and applied in two separate runs onto a 6 mL Resource-S ion-exchange column (30 mm × 16 mm, 15 µm) from Amersham Biosciences, previously equilibrated with buffer B and operated at a flow rate of 3 mL/min. After a 10 column volume wash with buffer B, proteins were resolved with a 25 column volume linear gradient from 20 to 500 mM NaCl in buffer B. Recombinant SSO0551 was eluted at approximately 250 mM NaCl and desalted by overnight dialysis against Buffer B. The resulting 20 mL protein solution was concentrated to a volume of 8 mL by means of Centricon Plus-20 filtration units (Millipore). The sample was again divided and applied in two separate runs onto a superdex75 gel filtration packed into a HR 16/50 column at a flow rate of 1.5 mL/min in 20 mM K₂HPO₄/KH₂PO₄ buffer (pH 7.2) containing 100 mM NaCl. The fractions obtained with the two runs were pooled and dialyzed overnight at 4°C against 10 mM HEPES buffer (pH 7.2). After dialysis, the fraction was centrifuged at 26,000 g for 20 min at 4°C and the protein concentration was measured by spectrophotometry using a molar absorption coefficient of 19060 M⁻¹ cm⁻¹ at 280
nm. The purified protein was flash frozen in liquid nitrogen and stored at -80°C at a concentration of 0.48 mg/ mL.

### Circular dichroism

Far- and near-UV circular dichroism spectra were recorded at 20°C between 200 and 300 nm on a J-810 Jasco spectropolarimeter equipped with a PTC-424S Jasco Peltier, using a quartz cuvette of 1 mm path length, with a 20 nm/ min scanning speed and a band-width of 1 nm. Three spectra of purified SSO0551 at 1.92 µM in 10 mM HEPES buffer (pH 7.2) were averaged and corrected from the baseline for buffer solvent contribution. Experimental data were analyzed using the program K2D [36] described by Andrade et al. [37].

### Determination of native molecular mass by gel filtration

The native molecular mass of SSO0551 was estimated by gel filtration chromatography on a Superdex 200 gel packed into a HR10/30 column (Amersham Biosciences) with a final bed volume of 24 mL. The column was equilibrated at room temperature at a flow rate of 0.5 mL/ min with 50 mM Tris/HCl buffer, pH 8.3, containing 50 mM NaCl and eluted with the same buffer. Protein standards used to calibrate the column were ribonuclease A (15.8 kDa), chymotrypsinogen A (21.2 kDa), ovalbumin (49.4 kDa), albumin (69.8 kDa), aldolase (191 kDa) and catalase (215 kDa), all from Amersham Biosciences. Exclusion limit was evaluated with dextran blue 2000 (Amersham Biosciences). A sample consisting of 90 µL of SSO0551 at 25.2 µM was injected and specific absorptions at 280 and 266 nm were followed.

### Mass spectrometry

Matrix-Assisted Laser Desorption/Ionization Time-Of-Flight (MALDI-TOF) mass measurements were performed using a Biflex IV instrument (Bruker Daltonics) in positive ionization mode. Protein samples and large peptidic fragments (>3500 Da) were applied to the target using sinapinic acid prepared as saturated solution in 30% acetonitrile, 70% milli-Q water and 0.1% TFA as matrix. Samples were prepared using the dried droplet method and measured in linear mode. Small peptide samples were measured in reflectron mode using α-cyano-4-hydroxycinnamic acid in 30% acetonitrile containing 0.1% trifluoroacetic acid as matrix. Mass spectra were obtained by summation of 100–210 laser shots. The instrument was calibrated for determination of entire protein masses using either a mixture of chymotrypsin and bovine serum albumine, or apomyoglobin and aldolase. For peptides, the instrument was calibrated using a pepmix calibration kit (Bruker Daltonics). When necessary, the mass spectrometer was also internally calibrated using some of the theoretical peptide masses.

### Limited protease digestion

For in-solution partial digestion, 0.2 nmol of pure SSO0551 were diluted into buffer D1 (20 mM TRIS/HCl, pH 7.8), buffer D2 (20 mM NH₄HCO₃, pH 7.8) or buffer D3 (20 mM TRIS/HCl, pH 7.8, containing 10 mM CaCl₂ and 5 mM DTT). Trypsin or chymotrypsin was added to SSO0551 diluted into buffer D1, whereas Glu-C or Arg-C was added to the protein diluted into buffer D2 or D3, respectively. Several enzyme/protein ratios (1:50 (w/w), 1:20 (w/w) and 1:2 (w/w)) were tested for each endoprotease. The digestions were performed at room temperature and aliquots were analyzed from 30 sec to 10–240 min. Digested samples were desalted using ZipTip<sub>C18</sub> or ZipTip<sub>C4</sub> pipette tips (Millipore) according to the protocol specified by the manufacturer and their mass directly evaluated by MALDI-TOF. Eventually, partially proteolyzed mixtures of larger quantities (10 nmol of SSO0551) were fractionated by reverse-phase HPLC using an Aquapore RP-300 column (PerkinElmer; 100 × 1.0 mm, 7 µm, 300 Å pore size) developed at 200 µL/min with a linear gradient from 5 to 90% of acetonitrile in TFA 0.1% over 45 min. The elution was monitored at 220 nm with an Agilent 1100 Series HPLC system equipped with a G1315 diode array detector. Individual fractions were concentrated by evaporation in a SpeedVac (Savant) and directly analyzed by MALDI-TOF.

### Lysine labeling by NHS-biotin

N-hydroxysuccinimide-biotin (NHS-biotin) was used to label ε-amino groups of SSO0551 lysines. After reaction the biotin labels resulted coupled to the lysines through a stable amide bond. The increase in mass for each label (C₁₀H₁₄N₂O₂S₁) should be 226.293 amu if average mass is considered or 226.078 amu in monoisotopic mode. Modification of lysine residues was carried out by incubating 1.25 nmol of SSO0551 in 20 mM HEPES, pH 7.2, with various amount of freshly prepared NHS-biotin reagent dissolved in anhydrous dimethylsulfoxide. After 30 min of incubation at room temperature, the reagent in excess was removed by a 30 min micro-dialysis against 20 mM HEPES, pH 7.2. Samples were directly desalted by using ZipTip<sub>C4</sub> (Millipore) prior MALDI-TOF analysis. They were eventually digested overnight with an endoprotease (trypsin, GluC or ArgC) and desalted by using ZipTip<sub>C18</sub> pipette tips (Millipore) prior mass analysis.

### Lysine cross-linking with DTSSP

3,3'-Dithio-bis [sulfosuccinimidyl-propionate] (DTSSP) was used to cross-link two ε-amino groups of SSO0551 lysines, essentially as described in [32]. The mass increase (in monoisotopic mode) for each label should be 191.991 amu (C₆H₈O₃S₂) or 87.998 amu (C₃H₄O₁S₁) when DTT treated. The increase in mass for an intramolecular cross-link between two lysines should be 173.981 amu (C₆H₆O₂S₂) or 175.997 amu (2 × C₃H₄O₁S₁) when
DTT treated. Therefore after reduction of the disulfide bridge by DTT, an additional increase of 2.016 amu should be measured. Reaction was carried out by incubating 0.25 nmol of SSO0551 in 20 mM NaH₂PO₄/Na₂HPO₄, pH 7.5 containing 150 mM NaCl, with various amount of DTSSP reagent (molar ratio of 20, 35, and 50 mol of DTSSP per mol of polypeptide). After 30 min of incubation at room temperature, the reagent in excess was removed by a 30 min micro-dialysis against 20 mM NaH₂PO₄/Na₂HPO₄, pH 7.5 containing 150 mM NaCl. Prior overnight trypsin proteolysis, urea (330 mM final concentration) was added to each sample. Before being desalted by using ZipTip<sub>C18</sub> pipette tips (Millipore), the digested peptide mixture was eventually reduced with 50 mM DTT for 30 minutes at 37°C to reduce the thiol linker.

**In silico analysis**

Sequence searching was performed using PSI-BLAST with default parameters. Multiple sequence alignments were performed using VectorNTI software package (Informatix Inc). Secondary structure predictions were obtained through the PSIPRED v2.4 web-interfaced facilities [38] described by McGuffin et al. [39]. The molar absorption coefficient at 280 nm for SSO0551 was obtained from calculation of the amino acid composition of the recombinant protein [40,41]. Isotopic and average mass of both DTSSP cross-linker and NHS-biotin were calculated using a web-interfaced molecular weight calculator [42]. The peptide assignment and the first attempt for identifying the labeled products and cross-linking products were performed using the FindMod package at ExPaSy [43]. If no match was found, a more detailed search for multiple labels or combinatorial cross-linkable peptide pairs was carried out. Partially proteolyzed products were assigned using the FindPept tool [44]. Tertiary structure predictions were carried out using publicly available online services, including 3D-PSSM [45], FUGUE [46] and PSIPRED [39]. *Ab initio* modeling was performed using the ROBETTA server [34,47]. Each model was analyzed in terms of proteolytic sensitivity using the NICKPRED software [35,48,49]. Residues accessibility have been calculated using a modified version of Connolly's MS program ([50]; Pellequer JL, unpublished results). Structural homologs were searched using DALI web server from the European Bioinformatics Institute [51]. Model views were obtained with the MOLSCRIPT program [52] and rendered using RASTER3D [53].

**List of abbreviations**

amu, atomic mass unit; COG, Cluster of Orthologous Group; DTSSP, 3,3'-dithio-bis [sulfosuccinimidyl-propionate]; IPTG, isopropyl-γ-D-thiogalactopyranoside; HPLC, high performance liquid chromatography; EDTA, Ethylenediaminetetraacetic acid; HEPES, 4-(2-hydroxyethyl)piperazine-1-ethanesulfonic acid; HEPPS, N-(2-hydroxyethyl)piperazine-N'-(3-propanesulfonic acid); IMAC, immobilized metal ion adsorption chromatography; MALDI-TOF, Matrix-assisted Laser Desorption/Ionization Time-of-Flight; NHS-biotin, N-hydroxysuccinimide-biotin; PSI-BLAST, Position-Specific Iterated Blast; Tris, 2-[bis(2-hydroxyethyl)amino]-2-(hydroxymethyl)-propane-1,3-diol.

**Authors' contributions**

JA conceived, coordinated the study and participated in all its experimental aspects. He analyzed the genomic distribution of this family of proteins, designed and engineered the recombinant SSO0551 molecule, conceived the mass spectrometry strategies and interpreted the data. He proposed *ab initio* modeling of SSO0551 and drafted the original manuscript. AD actively participated in conception of the mass spectrometry strategies, advised OS on execution and interpretation of mass spectrometry experiments and assisted in figure design. OS performed and interpreted all mass spectrometry experiments. JLP contributed its experience for the modeling aspects of the project. EQ contributed its experience in mass spectrometry-based topology. All authors participated in manuscript preparation, read and approved the final manuscript.

**Acknowledgements**

We gratefully acknowledge Yvan Zivanovic (CNRS-IGM, Orsay, France) for kind gift of S. sulfolobus total genomic DNA and Patrick Forterre (Université d'Orsay, Orsay, France) for initial discussions of the interest of characterizing SSO0551 protein. We thank our enthusiast technical assistants (CEA-VALRHO): Valérie Chaumont for performing the cloning and overexpression experiments, Charles Marchetti for operating the fermenter facilities, Bernard Fernandez for assistance with chromatography and recording circular dichroism signal, Isabelle Dany for initial fingerprint mass characterization of overproduced SSO0551, and Pascale Richard for technical support.

**References**

1. Woese CR, Kandler O, Wheelis ML: Towards a natural system of organisms: proposal for the domains Archaea, Bacteria, and Eucarya. Proc Natl Acad Sci U S A 1990, 87(12):4576-4579.

2. Dennis PP: Ancient ciphers: translation in Archaea. Cell 1997, 89(7):1007-1010.

3. Olsen GJ, Woese CR: Archaeal genomics: an overview. Cell 1997, 89(7):991-994.

4. Makarova KS, Koonin EV: Comparative genomics of Archaea: how much have we learned in six years, and what's next? Genome Biol 2003, 4(8):115.

5. Matte-Tailliez O, Zivanovic Y, Forterre P: Mining archaeal proteomes for eukaryotic proteins with novel functions: the PACE case. Trends Genet 2000, 16(12):533-536.

6. Armengaud J, Fernandez B, Chaumont V, Rollin-Genetet F, Finet S, Marchetti C, Myllykallio H, Vidaud C, Pellequer JL, Gribaldo S, et al.: Identification, purification, and characterization of an eukaryotic-like phosphopantetheine adenylyltransferase (coenzyme A biosynthetic pathway) in the hyperthermophilic archaeon Pyrococcus abyssi. J Biol Chem 2003, 278(33):31078-31087.

7. Armengaud J, Urbonavicius J, Fernandez B, Chaussinand G, Bujnicki JM, Grosjean H: N2-methylation of guanosine at position 10 in tRNA is catalyzed by a THUMP domain-containing, S-adenosylmethionine-dependent methyltransferase, conserved
in Archaea and Eukaryota. *J Biol Chem* 2004, 279(35):37142-37152.

8. Giaever G, Chu AM, Ni L, Connelly C, Riles L, Veronneau S, Dow S, Luca-Danila A, Anderson K, Andre B, et al.: Functional profiling of the Saccharomyces cerevisiae genome. *Nature* 2002, 418(6896):387-391.

9. Steinmetz LM, Scharfe C, Deutschbauer AM, Mokranjac D, Herman ZS, Jones T, Chu AM, Giaever G, Prokisch H, Oefner PJ, et al.: *Systematic screen for human disease genes in yeast. Nat Genet* 2002, 31(4):400-404.

10. Uetz P, Giot L, Cagney G, Mansfield TA, Judson RS, Knight JR, Lockshon D, Narayan V, Srinivasan M, Pochart P, et al.: A comprehensive analysis of protein-protein interactions in Saccharomyces cerevisiae. *Nature* 2000, 403(6770):623-627.

11. She Q, Singh RK, Confalonieri F, Zivanovic Y, Allard G, Awayez MJ, Chan-Weiher CC, Clausen IG, Curtis BA, De Moors A, et al.: The complete genome of the crenarchaeon Sulfolobus solfataricus P2. *Proc Natl Acad Sci U S A* 2001, 98(14):7835-7840.

12. Gabaldon T, Huynen MA: Prediction of protein function and pathways in the genome era. *Cell Mol Life Sci* 2004, 61(7-8):930-944.

13. Kerr ID: Sequence analysis of twin ATP binding cassette proteins involved in translational control, antibiotic resistance, and ribonuclease L inhibition. *Biochem Biophys Res Commun* 2004, 315(1):166-173.

14. Bisbal C, Martinand C, Silhol M, Lebleu B, Salehzada T: Cloning and characterization of a RNAse L inhibitor. A new component of the interferon-regulated 2-5A pathway. *J Biol Chem* 1995, 270(22):13308-13317.

15. Marx V: Quest: solve elusive. *Genomics & Proteomics* 2004:22-28.

16. Vincentelli R, Bignon C, Gruez A, Canaan S, Sulzenbacher G, Tegoni M, Campanacci V, Cambillau C: Medium-scale structural genomics: strategies for protein expression and crystallization. *Acc Chem Res* 2003, 36(3):165-172.

17. Goulding CW, Perry LJ: Protein production in Escherichia coli for structural studies by X-ray crystallography. *J Struct Biol* 2003, 142(1):133-143.

18. Fancy DA: Elucidation of protein-protein interactions using chemical cross-linking or label transfer techniques. *Curr Opin Chem Biol* 2000, 4(1):28-33.

19. Back JW, de Jong L, Muijsers AO, de Koster CG: Chemical cross-linking and mass spectrometry for protein structural modeling. *J Mol Biol* 2003, 331(2):303-313.

20. Sinz A: Chemical cross-linking and mass spectrometry for mapping three-dimensional structures of proteins and protein complexes. *J Mass Spectrom* 2003, 38(12):1225-1237.

21. Young MM, Tang N, Hempel JC, Oshiro CM, Taylor EW, Kuntz ID, Gibson BW, Dollinger G: High throughput protein fold identification by using experimental constraints derived from intramolecular cross-links and mass spectrometry. *Proc Natl Acad Sci U S A* 2000, 97(11):5802-5806.

22. Trester-Zedlitz M, Kamada K, Burley SK, Fenyo D, Chait BT, Muir TW: A modular cross-linking approach for exploring protein interactions. *J Am Chem Soc* 2003, 125(9):2416-2425.

23. Zappacosta F, Ingallinella P, Scaloni A, Pessi A, Bianchi E, Sollazzo M, Tramontano A, Marino G, Pucci P: Surface topology of Minibody by selective chemical modifications and mass spectrometry. *Protein Sci* 1997, 6(9):1901-1909.

24. Leite JF, Cascio M: Probing the topology of the glycine receptor by chemical modification coupled to mass spectrometry. *Biochemistry* 2002, 41(19):6140-6148.

25. Back JW, Sanz MA, De Jong L, De Koning LJ, Nijtmans LG, De Koster CG, Grivell LA, Van Der Spek H, Muijsers AO: A structure for the yeast prohibitin complex: Structure prediction and evidence from chemical crosslinking and mass spectrometry. *Protein Sci* 2002, 11(10):2471-2478.

26. Schulz DM, Ihling C, Clore GM, Sinz A: Mapping the topology and determination of a low-resolution three-dimensional structure of the calmodulin-melittin complex by chemical cross-linking and high-resolution FTICRMS: direct demonstration of multiple binding modes. *Biochemistry* 2004, 43(16):4703-4715.

27. Zappacosta F, Pessi A, Bianchi E, Venturini S, Sollazzo M, Tramontano A, Marino G, Pucci P: Probing the tertiary structure of proteins by limited proteolysis and mass spectrometry: the case of Minibody. *Protein Sci* 1996, 5(5):802-813.

28. Leite JF, Amoscato AA, Cascio M: Coupled proteolytic and mass spectrometry studies indicate a novel topology for the glycine receptor. *J Biol Chem* 2000, 275(18):13683-13689.

29. D'Ambrosio C, Talamo F, Vitale RM, Amodeo P, Tell G, Ferrara L, Scaloni A: Probing the dimeric structure of porcine aminoacylase I by mass spectrometric and modeling procedures. *Biochemistry* 2003, 42(15):4430-4443.

30. Hubbard SJ, Eisenmenger F, Thornton JM: Modeling studies of the change in conformation required for cleavage of limited proteolytic sites. *Protein Sci* 1994, 3(5):757-768.

31. Glocker MO, Borchers C, Fiedler W, Suckau D, Przybylski M: Molecular characterization of surface topology in protein tertiary structures by amino-acylation and mass spectrometric peptide mapping. *Bioconjug Chem* 1994, 5(6):583-590.

32. Bennett KL, Kussmann M, Bjork P, Godzwon M, Mikkelsen M, Sorensen P, Roeptorff P: Chemical cross-linking with thiol-cleavable reagents combined with differential mass spectrometric peptide mapping – a novel approach to assess intermolecular protein contacts. *Protein Sci* 2000, 9(8):1503-1518.

33. [http://www.ncbi.nlm.nih.gov/Structure/cdd/cdd.shtml].

34. Kim DE, Chivian D, Baker D: Protein structure prediction and analysis using the Robetta server. *Nucleic Acids Res* 2004, 32(Web Server):W526-531.

35. Hubbard SJ, Beynon RJ, Thornton JM: Assessment of conformational parameters as predictors of limited proteolytic sites in native protein structures. *Protein Eng* 1998, 11(5):349-359.

36. [http://www.embl-heidelberg.de/~andrade/k2d/].

37. Andrade MA, Chacon P, Merelo JJ, Moran F: Evaluation of secondary structure of proteins from UV circular dichroism spectra using an unsupervised learning neural network. *Protein Eng* 1993, 6(4):383-390.

38. [http://bioinf.cs.ucl.ac.uk/psipred/].

39. McGuffin LJ, Bryson K, Jones DT: The PSIPRED protein structure prediction server. *Bioinformatics* 2000, 16(4):404-405.

40. Gill SC, von Hippel PH: Calculation of protein extinction coefficients from amino acid sequence data. *Anal Biochem* 1989, 182(2):319-326.

41. [http://www.expasy.org/tools/propparam.html].

42. [http://www.chemie.fu-berlin.de/cgi-bin/molform/].

43. [http://www.expasy.org/tools/findmod/].

44. [http://www.expasy.org/tools/findpept.html].

45. Kelley LA, MacCallum RM, Sternberg MJ: Enhanced genome annotation using structural profiles in the program 3D-PSSM. *J Mol Biol* 2000, 299(2):499-520.

46. Shi J, Blundell TL, Mizuguchi K: FUGUE: sequence-structure homology recognition using environment-specific substitution tables and structure-dependent gap penalties. *J Mol Biol* 2001, 310(1):243-257.

47. [http://robetta.bakerlab.org/].

48. Hubbard SJ: The structural aspects of limited proteolysis of native proteins. *Biochim Biophys Acta* 1998, 1382(2):191-206.

49. [http://wolf.bms.umist.ac.uk/npred/nickpred.html].

50. Connolly ML: Solvent-accessible surfaces of proteins and nucleic acids. *Science* 1983, 221:709-713.

51. [http://www.ebi.ac.uk/dali/].

52. Kraulis PJ: MOLSCRIPT: a program to produce both detailed and schematic plots of protein structures. *J Appl Cryst* 1991, 24:946-950.

53. Merritt EA, Bacon DJ: Raster3D: Photorealistic molecular graphics. *Meth Enzymol* 1997, 277:505-524.

54. [http://www.ncbi.nlm.nih.gov/].