Gelsolin is an actin-binding protein (82 kDa) consisting of six repeated segments (S1–S6), each approximately 120 residues long. It interacts with phospholipids and we previously showed that phosphatidylinositol 4,5-bisphosphate promotes phosphorylation of gelsolin by the tyrosine kinase c-Src. We used a combination of different methods, such as thin-layer chromatography and anti-phosphotyrosine-agarose immunoprecipitation of phosphopeptides combined with matrix assisted laser desorption ionization-mass spectrometry (MALDI-MS) and post source decay (PSD) analysis, to identify the phosphorylation sites in gelsolin. The major phosphorylation site (Tyr438) was located in subdomain 4 (S4). Phosphorylation of gelsolin in the gelsolin-actin2 complex was inhibited by 90%. Gelsolin phosphorylation by c-Src in the presence of lysophosphatidic acid also revealed Tyr438 as the most prominent site. Additional minor sites were found using the anti-phosphotyrosine bead immunoprecipitation method followed by MALDI-MS and PSD analysis. These sites, representing ∼5% of the total phosphate incorporation, were identified as Tyr59, Tyr382, Tyr576, and Tyr624. Based on these results we generated antibodies which specifically recognize Tyr438 phosphorylated gelsolin.

Cellulases are increasingly being used for industrial purposes, particularly in washing powders, yet little is known of the factors governing the stability of proteins in detergent solutions. We present a comparative analysis of the behavior of the cellulase Cel45 from Humicola insolens in the presence of the denaturant guanidinium chloride and the anionic detergent C12-LAS. Although Cel45 unfolds in GdmCl according to a simple two-state model under equilibrium conditions, it accumulates a transient intermediate during refolding. The four disulfide bonds do not contribute detectably to the stability of the native state. Cel45 is unfolded by very low concentrations of C12-LAS (1–4 mM). An analysis of 16 mutants of Cel45 shows a very weak correlation between unfolding rates in denaturant and detergent; mutants that have the same unfolding rate in GdmCl (within a factor of 1.5) vary 1,000-fold in their unfolding rates in C12-LAS. The data support a simple model for unfolding by detergent, in which the introduction of positive charges or removal of negative charges greatly increases detergent sensitivity, while interactions with the hydrophobic detergent tail contribute to a smaller extent. This implies that different detergent-mediated unfolding pathways exist, whose accessibilities depend on individual residues. Double-mutant cycles reveal that mutations in two proximal residues lead to repulsion and a destabilization greater than the sum of the individual mutations as measured by GdmCl denaturation, but they also reduce the affinity for LAS and therefore actually stabilize the protein relative to wild-type. Ligands that interact strongly with the denatured state may therefore alter the unfolding process.

β-Galactosidase (lacZ) from Escherichia coli is a 464 kDa homotetramer. Each subunit consists of five domains, the third being an α/β barrel that contains most of the active site residues. A comparison is made between each of the domains and a large set of proteins representative of all structures from the protein data bank. Many structures include an α/β barrel. Those that are most similar to the α/β barrel of E. coli β-galactosidase have similar catalytic residues and belong to the so-called “4/7 superfamily” of glycosyl hydrolases. The structure comparison suggests that β-amylase should also be included in this family. Of three structure comparison methods tested, the “ProSup” procedure of Zu-Kang and Sippl and the “Superimpose” procedure of Diederichs were slightly superior in discriminating the members of this superfamily, although all procedures were very powerful in identifying related protein structures. Domains 1, 2, and 4 of E. coli β-galactosidase have topologies related to “jelly-roll barrels” and “immunoglobulin constant” domains. This fold also occurs in the cellulose binding domains (CBDs) of a number of glycosyl hydrolases. The fold of domain 1 of E. coli β-galactosidase is closely related to some CBDs, and the domain contributes to substrate binding, but in a manner unrelated to cellulose binding by the CBDs. This is typical of domains 1, 2, 4, and 5, which appear to have been recruited to play roles in β-galactosidase that are unrelated to the functions that such domains provide in other contexts. It is proposed that β-galactosidase arose from a prototypical single domain α/β barrel with an extended active site cleft. The subsequent incorporation of elements from other domains could then have reduced the size of the active site from a cleft to a pocket to better hydrolyze the disaccharide lactose and, at the same time, to facilitate the production of inducer, allolactose.

Apomyoglobin from sperm whale is often used for studies of ligand binding, protein folding, and protein stability. In an effort to describe its conformational properties in solution, homonuclear and heteronuclear (13C and 15N) NMR methods were applied to the protein in its native state. Assignments were confirmed for nuclear Overhauser effects (NOEs) involving side chain and backbone protons in the folded regions of the structure. These NOEs were used to derive distance restraints. The shifts induced by the hydrophobic dye 8-anilino-1-naphthalenesulfonic acid (ANS) were inspected in the regions remote from its binding site and served as an indicator of conformational flexibility. 3JαH-NH values were obtained to assess dihedral angle averaging and to provide additional restraints. A family of structures was calculated with X-PLOR and an ab initio simulated annealing protocol using holomyoglobin as a template. Where the structure appeared well defined by chemical shift, line width, ANS perturbation, and density of NOEs, the low resolution model of apomyoglobin provides a valid approximation for the structure. The new model offers an improved representation of the folded regions of the protein, which encompass the A, B, E, helices as well as parts of the G and H helices. Regions that are less well defined at this stage of calculations include the CD corner and the end of the H-helix. The EF-F-FG segment remains uncharacterized.

The N-terminal 17 residues of ubiquitin have been shown by 1H NMR to fold autonomously into a β-hairpin structure in aqueous solution. This structure has a specific, native-like register, though side-chain contacts differ in detail from those observed in the intact protein. An autonomously folding hairpin has previously been identified in the case of streptococcal protein G, which is structurally homologous with ubiquitin, but remarkably, the two are not in topologically equivalent positions in the fold. This suggests that the organization of folding may be quite different for proteins sharing similar tertiary structures. Two smaller peptides have also been studied, corresponding to the isolated arms of the N-terminal hairpin of ubiquitin, and significant differences from simple random coil predictions observed in the spectra of these subfragments, suggestive of significant limitation of the backbone conformational space sampled, presumably as a consequence of the strongly β-structure favoring composition of the sequences. This illustrates the ability of local sequence elements to express a propensity for β-structure even in the absence of actual sheet formation. Attempts were made to estimate the population of the folded state of the hairpin, in terms of a simple two-state folding model. Using published “random coil” values to model the unfolded state, and values derived from native ubiquitin for the putative unique, folded state, it was found that the apparent population varied widely for different residues and with different NMR parameters. Use of the spectra of the subfragment peptides to provide a more realistic model of the unfolded state led to better agreement in the estimates that could be obtained from chemical shift and coupling constant measurements, while making it clear that some other approaches to population estimation could not give meaningful results, because of the tendency to populate the β-region of conformational space even in the absence of the hairpin structure.

We report here that in vitro exposure of monomeric actin to hydrogen peroxide leads to a conversion of 6 of the 16 methionine residues to methionine sulfoxide residues. Although the initial effect of H2O2 on actin is the oxidation of Cys374, we have found that Met44, Met47, Met176, Met190, Met269, and Met355 are the other sites of the oxidative modification. Met44 and Met47 are the methionyl sites first oxidized. The methionine residues that are oxidized are not simply related to their accessibility to the external medium and are found in all four subdomains of actin. The conformations of subdomain 1, a region critical for the functional binding of different actin-binding proteins, and subdomain 2, which plays important roles in the polymerization process and stabilization of the actin filament, are changed upon oxidation. The conformational changes are deduced from the increased exposure of hydrophobic residues, which correlates with methionine sulfoxide formation, from the perturbations in tryptophan fluorescence, and from the decreased susceptibility to limited proteolysis of oxidized actin.

Mg2+ binds to calmodulin without inducing the changes in secondary structure that are characteristic of Ca2+ binding, or the exposure of hydrophobic surfaces that are involved in typical Ca2+-dependent target interactions. The binding of Mg2+ does, however, produce significant spectroscopic changes in residues located in the Ca2+-binding loops, and the Mg-calmodulin complex is significantly different from apo-calmodulin in loop conformation. Direct measurement of Mg2+ binding constants, and the effects of Mg2+ on Ca2+ binding to calmodulin, are consistent with specific binding of Mg2+, in competition with Ca2+. Mg2+ increases the thermodynamic stability of calmodulin, and we conclude that under resting, nonstimulated conditions, cellular Mg2+ has a direct role in conferring stability on both domains of apo-calmodulin. Apo-calmodulin binds typical target sequences from skeletal muscle myosin light chain kinase and neuromodulin with Kd ∼ 70–90 nM (at low ionic strength). These affinities are virtually unchanged by 5 mM Mg2+, in marked contrast to the strong enhancement of peptide affinity induced by Ca2+. Under conditions of stimulation and increased [Ca2+], Mg2+ has a role in directing the mode of initial target binding preferentially to the C-domain of calmodulin, due to the opposite relative affinities for binding of Ca2+ and Mg2+ to the two domains. Mg2+ thus amplifies the intrinsic differences of the domains, in a target specific manner. It also contributes to setting the Ca2+ threshold for enzyme activation and increases the importance of a partially Ca2+-saturated calmodulin–target complex that can act as a regulatory kinetic and equilibrium intermediate in Ca2+-dependent target interactions.

We have characterized the thermodynamic stability of the SH3 domain from the Saccharomyces cerevisiae Abp1p protein and found it to be relatively low compared to most other SH3 domains, with a Tm of 60 °C and a ΔGu of 3.08 kcal/mol. Analysis of a large alignment of SH3 domains led to the identification of atypical residues at eight positions in the wild-type Abp1p SH3 domain sequence that were subsequently replaced by the residue seen most frequently at that position in the alignment. Three of the eight mutants constructed in this way displayed increases in Tm ranging from 8 to 15 °C with concomitant increases in ΔGu of up to 1.4 kcal/mol. The effects of these substitutions on folding thermodynamics and kinetics were entirely additive, and a mutant containing all three was dramatically stabilized with a Tm greater than 90 °C and a ΔGu more than double that of the wild-type domain. The folding rate of this hyperstable mutant was 10-fold faster than wild-type, while its unfolding rate was fivefold slower. All of the stabilized mutants were still able to bind a target peptide with wild-type affinity. We have analyzed the stabilizing amino acid substitutions isolated in this study and several other similar sequence alignment based studies. In approximately 25% of cases, increased stability can be explained by enhanced propensity of the substituted residue for the local backbone conformation at the mutagenized site.

A technique is described for identifying and locating posttranslational modifications (PTMs) in peptides and proteins of known sequence by interpretation of cn ion signals generated by in-source decay during delayed ion extraction in matrix-assisted laser desorption/ionization time-of-flight mass spectrometry. Sites of phosphorylation in seven synthetic peptides were determined, as was the location of both the heme group and N,N,N-trimethyllysine in yeast cytochrome c. A semi-automated data analysis process facilitates the identification of segments of the sequence on each side of the PTM, permitting its placement at the junction of the segments and definition of the added mass. A graphical display facilitates illustration of both the location and mass of the PTM.

Although the three-dimensional structure of the dimeric class 3 rat aldehyde dehydrogenase has recently been published (Liu ZJ et al., 1997, Nature Struct Biol 4:317–326), few mechanistic studies have been conducted on this isoenzyme. We have characterized the enzymatic properties of recombinant class 3 human stomach aldehyde dehydrogenase, which is very similar in amino acid sequence to the class 3 rat aldehyde dehydrogenase. We have determined that the rate-limiting step for the human class 3 isozyme is hydride transfer rather than deacylation as observed for the human liver class 2 mitochondrial enzyme. No enhancement of NADH fluorescence was observed upon binding to the class 3 enzyme, while fluorescence enhancement of NADH has been previously observed upon binding to the class 2 isoenzyme. It was also observed that binding of the NAD cofactor inhibited the esterase activity of the class 3 enzyme while activating the esterase activity of the class 2 enzyme. Site-directed mutagenesis of two conserved glutamic acid residues (209 and 333) to glutamine residues indicated that, unlike in the class 2 enzyme, Glu333 served as the general base in the catalytic reaction and E209Q had only marginal effects on enzyme activity, thus confirming the proposed mechanism (Hempel J et al., 1999, Adv Exp Med Biol 436:53–59). Together, these data suggest that even though the subunit structures and active site residues of the isozymes are similar, the enzymes have very distinct properties besides their oligomeric state (dimer vs. tetramer) and substrate specificity.

CORA is a suite of programs for multiply aligning and analyzing protein structural families to identify the consensus positions and capture their most conserved structural characteristics (e.g., residue accessibility, torsional angles, and global geometry as described by inter-residue vectors/contacts). Knowledge of these structurally conserved positions, which are mostly in the core of the fold and of their properties, significantly improves the identification and classification of newly-determined relatives. Information is encoded in a consensus three-dimensional (3D) template and relatives found by a sensitive alignment method, which employs a new scoring scheme based on conserved residue contacts. By encapsulating these critical “core” features, templates perform more reliably in recognizing distant structural relatives than searches with representative structures.

Parameters for 3D-template generation and alignment were optimized for each structural class (mainly-α, mainly-β, α-β), using representative superfold families. For all families selected, the templates gave significant improvements in sensitivity and selectivity in recognizing distant structural relatives. Furthermore, since templates contain less than 70% of fold positions and compare fewer positions when aligning structures, scans are at least an order of magnitude faster than scans using selected structures. CORA was subsequently tested on eight other broad structural families from the CATH database.

Diagnostics plots are generated automatically and provide qualitative assistance for classifying newly determined relatives. They are demonstrated here by application to the large globin-like fold family. CORA templates for both homologous superfamilies and fold families will be stored in CATH and used to improve the classification and analysis of newly determined structures.

The partly folded states of α-lactalbumin (α-LA) exposed to acid solution at pH 2.0 (A-state) or at neutral pH upon EDTA-mediated removal of the single protein-bound calcium ion (apo form) have been probed by limited proteolysis experiments. These states are nowadays commonly considered to be molten globules and thus protein-folding intermediates. Pepsin was used for proteolysis at acid pH, while proteinase K and chymotrypsin at neutral pH. The expectations were that these proteolytic probes would detect sites and/or chain regions in the partly folded states of α-LA sufficiently dynamic, or even unfolded, capable of binding and adaptation to the specific stereochemistry of the protease's active site. A time-course analysis of the proteolytic events revealed that the fast, initial proteolytic cuts of the 123-residue chain of α-LA in its A-state or apo form by the three proteases occur at the same chain region 39–54, the actual site(s) of cleavage depending upon the protease employed. This region in native α-LA encompasses the β-sheets of the protein. Subsequent cleavages occur mostly at chain regions 31–35 and 95–105. Four fragment species of α-LA have been isolated by reverse-phase high-performance liquid chromatography, and their conformational properties examined by circular dichroism and fluorescence emission spectroscopy. The single chain fragment 53–103, containing all the binding sites for calcium in native α-LA and cross-linked by two disulfide bridges, maintains in aqueous buffer and in the presence of calcium ions a folded structure characterized by the same content of α-helix of the corresponding chain segment in native α-LA. Evidence for some structure was also obtained for the two-chain species 1–40 and 104–123, as well as 1–31 and 105–123, both systems being covalently linked by two disulfide bonds. In contrast, the protein species given by fragment 1–34 connected to fragment 54–123 or 57–123 via four disulfide bridges adopts in solution a folded structure with the helical content expected for a native-like conformation. Of interest, the proteolytic fragment species herewith isolated correspond to the structural domains and subdomains of α-LA that can be identified by computational analysis of the three-dimensional structure of native α-LA (Siddiqui AS, Barton GI, 1995, Protein Sci 4:872–884). The fast, initial cleavages at the level of the β-sheet region of native α-LA indicate that this region is highly mobile or even unfolded in the α-LA molten globule(s), while the rest of the protein chain maintains sufficient structure and rigidity to prevent extensive proteolysis. The subsequent cleavages at chain segment 95–105 indicate that also this region is somewhat mobile in the A-state or apo form of the protein. It is concluded that the overall domain topology of native α-LA is maintained in acid or at neutral pH upon calcium depletion. Moreover, the molecular properties of the partly folded states of α-LA deduced here from proteolysis experiments do correlate with those derived from previous NMR and other physicochemical measurements.

Prolyl oligopeptidase, an enzyme implicated in memory disorders, is a member of a new serine peptidase family. Crystallographic studies (Fülöp et al., 1998) revealed a novel oxyanion binding site containing a tyrosine residue, Tyr473. To study the importance of Tyr473 OH, we have produced prolyl oligopeptidase and its Tyr473Phe variant in Escherichia coli. The specificity rate constant, kcat/Km, for the modified enzyme decreased by a factor of 8–40 with highly specific substrates, Z-Gly-Pro-Nap, and a fluorogenic octapeptide. With these compounds, the decline in kcat was partly compensated for by reduction in Km, a difference from the extensively studied subtilisin. With the less specific suc-Gly-Pro-Nap, the Km value, which approximates Ks, was not significantly changed, resulting in greater diminution (∼500-fold) in kcat/Km. The second-order rate constant for the reaction with Z-Pro-prolinal, a slow tight-binding transition-state analogue inhibitor, and the Ki values for a slow substrate and two product-like inhibitors were not significantly affected by the Tyr473 OH group. The mechanism of transition-state stabilization was markedly dependent upon the nature of substrate and varied with pH as the enzyme interconverted between its two catalytically competent forms.

Transient absorbance measurements following laser flash photolysis have been used to measure the rate constants for electron transfer (et) from reduced Anabaena ferredoxin (Fd) to wild-type and seven site-specific charge-reversal mutants of Anabaena ferredoxin:NADP+ reductase (FNR). These mutations have been designed to probe the importance of specific positively charged amino acid residues on the surface of the FNR molecule near the exposed edge of the FAD cofactor in the protein–protein interaction during et with Fd. The mutant proteins fall into two groups: overall, the K75E, R16E, and K72E mutants are most severely impaired in et, and the K138E, R264E, K290E, and K294E mutants are impaired to a lesser extent, although the degree of impairment varies with ionic strength. Binding constants for complex formation between the oxidized proteins and for the transient et complexes show that the severity of the alterations in et kinetics for the mutants correlate with decreased stabilities of the protein–protein complexes. Those mutated residues, which show the largest effects, are located in a region of the protein in which positive charge predominates, and charge reversals have large effects on the calculated local surface electrostatic potential. In contrast, K138, R264, K290, and K294 are located within or close to regions of intense negative potential, and therefore the introduction of additional negative charges have considerably smaller effects on the calculated surface potential. We attribute the relative changes in et kinetics and complex binding constants for these mutants to these characteristics of the surface charge distribution in FNR and conclude that the positively charged region of the FNR surface located in the vicinity of K75, R16, and K72 is especially important in the binding and orientation of Fd during electron transfer.

We have determined the three-dimensional structure of the potassium channel inhibitor HsTX1, using nuclear magnetic resonance and molecular modeling. This protein belongs to the scorpion short toxin family, which essentially contains potassium channel blockers of 29 to 39 amino acids and three disulfide bridges. It is highly active on voltage-gated Kv1.3 potassium channels. Furthermore, it has the particularity to possess a fourth disulfide bridge. We show that HsTX1 has a fold similar to that of the three-disulfide–bridged toxins and conserves the hydrophobic core found in the scorpion short toxins. Thus, the fourth bridge has no influence on the global conformation of HsTX1. Most residues spatially analogous to those interacting with voltage-gated potassium channels in the three-disulfide–bridged toxins are conserved in HsTX1. Thus, we propose that Tyr21, Lys23, Met25, and Asn26 are involved in the biological activity of HsTX1. As an additional positively charged residue is always spatially close to the aromatic residue in toxins blocking the voltage-gated potassium channels, and as previous mutagenesis experiments have shown the critical role played by the C-terminus in HsTX1, we suggest that Arg33 is also important for the activity of the four disulfide-bridged toxin. Docking calculations confirm that, if Lys23 and Met25 interact with the GYGDMH motif of Kv1.3, Arg33 can contact Asp386 and, thus, play the role of the additional positively charged residue of the toxin functional site. This original configuration of the binding site of HsTX1 for Kv1.3, if confirmed experimentally, offers new structural possibilities for the construction of a molecule blocking the voltage-gated potassium channels.

The structural organization of the B1 domain of streptococcal protein G (PGA) has been probed using molecular dynamics simulations, with a particular emphasis on the role of the solvent exposed Ile6 residue. In addition to the native protein (WT-PGA), three single-mutants (I6G-PGA, I6F-PGA, and I6T-PGA), one double-mutant (I6T,T53G-PGA), and three isolated peptide fragments (corresponding to the helix and the two β-hairpins) were studied in the presence of explicit water molecules. Comparative analysis of the various systems showed that the level of perturbation was directly related to the hydrophobicity and the size of the side chain of residue 6, the internal rigidity of the proteins decreasing in the order I6T-PGA > I6G-PGA > WT-PGA > I6F-PGA. The results emphasized the importance of residue 6 in controlling both the integrity of the sheet's surface and the orientation of the helix in relation to the sheet by modulation of surface/core interactions. The effects of mutations were delocalized across the structure, and glycine residues, in particular, absorbed most of the introduced strain. A qualitative structural decomposition of the native fold into elementary building-blocks was achieved using principal component analysis and mechanical response matrices. Within this framework, internal motions of the protein were described as coordinated articulations of these structural units, mutations affecting mostly the amplitude of the motions rather than the structure/location of the building-blocks. Analysis of the isolated peptidic fragments suggested that packing did not play a determinant role in defining the elementary building-blocks, but that chain topology was mostly responsible.

A 20-residue peptide E5 containing five glutamates, an analog of the fusion peptide of influenza virus hemagglutinin (HA) exhibiting fusion activity at acidic pH lower than 6.0–6.5 was studied by circular dichroism (CD), Fourier transform infrared, and 1H-NMR spectroscopy in water, water/trifluoroethanol (TFE) mixtures, dodecylphosphocholine (DPC) micelles, and phospholipid vesicles. E5 became structurally ordered at pH ≤6 and the helical content in the peptide increased in the row: water < water/TFE < DPC ∼ phospholipid vesicle while the amount of β-structure was approximately reverse. 1H-NMR data and line-broadening effect of 5-, 16-doxylstearates on proton resonances of DPC bound peptide showed E5 forms amphiphilic α-helix in residues 2–18, which is flexible in 11–18 part. The analysis of the proton chemical shifts of DPC bound and CD intensity at 220 nm of phospholipid bound E5 showed that the pH dependence of helical content is characterized by the same pKa ≈5.6. Only Glu11 and Glu15 in DPC bound peptide showed such elevated pKas, presumably due to transient hydrogen bond(s) Glu11 (Glu15) δCOO−(H+) … HN Glu15 that dispose(s) the side chain of Glu11 (Glu15) residue(s) close to the micelle/water interface. These glutamates are present in the HA-fusion peptide and the experimental half-maximal pH of fusion for HA and E5 peptides is ∼5.6. Therefore, a specific anchorage of these peptides onto membrane necessary for fusion is likely driven by the protonation of the carboxylate group of Glu11 (Glu15) residue(s) participating in transient hydrogen bond(s).

A database of functional sites for proteins with known structures, SITE, is constructed and used in conjunction with a simple pattern matching program SiteMatch to evaluate possible function conservation in a recently constructed database of fold predictions for Escherichia coli proteins (Rychlewski L et al., 1999, Protein Sci 8:614–624). In this and other prediction databases, fold predictions are based on algorithms that can recognize weak sequence similarities and putatively assign new proteins into already characterized protein families. It is not clear whether such sequence similarities arise from distant homologies or general similarity of physicochemical features along the sequence. Leaving aside the important question of nature of relations within fold superfamilies, it is possible to assess possible function conservation by looking at the pattern of conservation of crucial functional residues. SITE consists of a multilevel function description based on structure annotations and structure analyses. In particular, active site residues, ligand binding residues, and patterns of hydrophobic residues on the protein surface are used to describe different functional features. SiteMatch, a simple pattern matching program, is designed to check the conservation of residues involved in protein activity in alignments generated by any alignment method. Here, this procedure is used to study conservation of functional features in alignments between protein sequences from the E. coli genome and their optimal structural templates. The optimal templates were identified and alignments taken from the database of genomic structural predictions was described in a previous publication (Rychlewski L et al., 1999, Protein Sci 8:614–624). An automated assessment of function conservation is used to analyze the relation between fold and function similarity for a large number of fold predictions. For instance, it is shown that identifying low significance predictions with a high level of functional residue conservations can be used to extend the prediction sensitivity for fold prediction methods. Over 100 new fold/function predictions in this class were obtained in the E. coli genome. At the same time, about 30% of our previous fold predictions are not confirmed as function predictions, further highlighting the problem of function divergence in fold superfamilies.

The genome of the unicellular cyanobacterium Synechocystis sp. PCC 6803 contains a gene (slr2097, glbN) encoding a 123 amino-acid product with sequence similarity to globins. Related proteins from cyanobacteria, ciliates, and green algae bind oxygen and have a pronounced tendency to coordinate the heme iron with two protein ligands. To study the structural and functional properties of Synechocystis sp. PCC 6803 hemoglobin, slr2097 was cloned and overexpressed in Escherichia coli. Purification of the hemoglobin was performed after addition of hemin to the clarified cell lysate. Recombinant, heme-reconstituted ferric Synechocystis sp. PCC 6803 hemoglobin was found to be a stable helical protein, soluble to concentrations higher than 500 μM. At neutral pH, it yielded an electronic absorption spectrum typical of a low-spin ferric species, with maxima at 410 and 546 nm. The proton NMR spectrum revealed sharp lines spread over a chemical shift window narrower than 40 ppm, in support of low-spin hexacoordination of the heme iron. Nuclear Overhauser effects demonstrated that the heme is inserted in the protein matrix to produce one major equilibrium form. Addition of dithionite resulted in an absorption spectrum with maxima at 426, 528, and 560 nm. This reduced form appeared capable of carbon monoxide binding. Optical data also suggested that cyanide ions could bind to the heme in the ferric state. The spectral properties of the putative Synechocystis sp. PCC 6803 hemoglobin confirmed that it can be used for further studies of an ancient hemoprotein structure.

The effects of engineered disulfide bonds on protein stability are poorly understood because they can influence the structure, dynamics, and energetics of both the native and denatured states. To explore the effects of two engineered disulfide bonds on the stability of barnase, we have conducted a combined molecular dynamics and NMR study of the denatured state of the two mutants. As expected, the disulfide bonds constrain the denatured state. However, specific extended β-sheet structure can also be detected in one of the mutant proteins. This mutant is also more stable than would be predicted. Our study suggests a possible cause of the very high stability conferred by this disulfide bond: the wild-type denatured ensemble is stabilized by a nonnative hydrophobic cluster, which is constrained from occurring in the mutant due to the formation of secondary structure.