Genetically-encoded affinity tags constitute an important strategy for purifying proteins. Here, we have designed a novel affinity matrix based on the bis-arsenical fluorescein dye FlAsH, which specifically recognizes short α-helical peptides containing the sequence CCXXCC (Griffin BA, Adams SR, Tsien RY, 1998, Science 281:269–272). We find that kinesin tagged with this cysteine-containing helix binds specifically to FlAsH resin and can be eluted in a fully active form. This affinity tag has several advantages over polyhistidine, the only small affinity tag in common use. The protein obtained with this single chromatographic step from crude Escherichia coli lysates is purer than that obtained with nickel affinity chromatography of 6xHis tagged kinesin. Moreover, unlike nickel affinity chromatography, which requires high concentrations of imidazole or pH changes for elution, protein bound to the FlAsH column can be completely eluted by dithiothreitol. Because of these mild elution conditions, FlAsH affinity chromatography is ideal for recovering fully active protein and for the purification of intact protein complexes.

The crystal structure of ligand-free tryptophanyl-tRNA synthetase (TrpRS) was solved at 2.9 Å using a combination of molecular replacement and maximum-entropy map/phase improvement. The dimeric structure (R = 23.7, Rfree = 26.2) is asymmetric, unlike that of the TrpRS tryptophanyl–5′AMP complex (TAM; Doublié S, Bricogne G, Gilmore CJ, Carter CW Jr, 1995, Structure 3:17–31). In agreement with small-angle solution X-ray scattering experiments, unliganded TrpRS has a conformation in which both monomers open, leaving only the tryptophan-binding regions of their active sites intact. The amino terminal αA-helix, TIGN, and KMSKS signature sequences, and the distal helical domain rotate as a single rigid body away from the dinucleotide-binding fold domain, opening the AMP binding site, seen in the TAM complex, into two halves. Comparison of side-chain packing in ligand-free TrpRS and the TAM complex, using identification of nonpolar nuclei (Ilyin VA, 1994, Protein Eng 7:1189–1195), shows that significant repacking occurs between three relatively stable core regions, one of which acts as a bearing between the other two. These domain rearrangements provide a new structural paradigm that is consistent in detail with the “induced-fit” mechanism proposed for TyrRS by Fersht et al. (Fersht AR, Knill-Jones JW, Beduelle H, Winter G, 1988, Biochemistry 27:1581–1587). Coupling of ATP binding determinants associated with the two catalytic signature sequences to the helical domain containing the presumptive anticodon-binding site provides a mechanism to coordinate active-site chemistry with relocation of the major tRNA binding determinants.

Distant homologies between proteins are often discovered only after three-dimensional structures of both proteins are solved. The sequence divergence for such proteins can be so large that simple comparison of their sequences fails to identify any similarity. New generation of sensitive alignment tools use averaged sequences of entire homologous families (profiles) to detect such homologies. Several algorithms, including the newest generation of BLAST algorithms and BASIC, an algorithm used in our group to assign fold predictions for proteins from several genomes, are compared to each other on the large set of structurally similar proteins with little sequence similarity. Proteins in the benchmark are classified according to the level of their similarity, which allows us to demonstrate that most of the improvement of the new algorithms is achieved for proteins with strong functional similarities, with almost no progress in recognizing distant fold similarities.

It is also shown that details of profile calculation strongly influence its sensitivity in recognizing distant homologies. The most important choice is how to include information from diverging members of the family, avoiding generating false predictions, while accounting for entire sequence divergence within a family. PSI-BLAST takes a conservative approach, deriving a profile from core members of the family, providing a solid improvement without almost any false predictions. BASIC strives for better sensitivity by increasing the weight of divergent family members and paying the price in lower reliability. A new FFAS algorithm introduced here uses a new procedure for profile generation that takes into account all the relations within the family and matches BASIC sensitivity with PSI-BLAST like reliability.

Pigeon liver malic enzyme was inactivated and cleaved at Asp141, Asp194, and Asp464 by the Cu2+-ascorbate system in acidic environment. Site-specific mutagenesis was performed at these putative metal-binding sites. Three point mutants, D141N, D194N, and D464N; three double mutants, D(141,194)N, D(194,464)N, and D(141,464)N; and a triple mutant, D(141,194,464)N; as well as the wild-type malic enzyme (WT) were successfully cloned and expressed in Escherichia coli cells. All recombinant enzymes, except the triple mutant, were purified to apparent homogeneity by successive Q-Sepharose and adenosine-2′,5′-bisphosphate-agarose columns. The mutants showed similar apparent Km,NADP values to that of the WT. The Km,Mal value was increased in the D141N and D194N mutants. The Km,Mn value, on the other hand, was increased only in the D141N mutant by 14-fold, corresponding to ∼1.6 kcal/mol for the Asp141-Mn2+ binding energy. Substrate inhibition by l-malate was only observed in WT, D464N, and D(141,464)N. Initial velocity experiments were performed to derive the various kinetic parameters. The possible interactions between Asp141, Asp194, and Asp464 were analyzed by the double-mutation cycles and triple-mutation box. There are synergistic weakening interactions between Asp141 and Asp194 in the metal binding that impel the D(141,194)N double mutant to an overall specificity constant [kcat/(Kd,MnKm,MalKm,NADP)] at least four orders of magnitude smaller than the WT value. This difference corresponds to an increase of 6.38 kcal/mol energy barrier for the catalytic efficiency. Mutation at Asp464, on the other hand, has partial additivity on the mutations at Asp141 and Asp194. The overall specificity constants for the double mutants D(194,464)N and D(141,464)N or the triple mutant D(141,194,464)N were decreased by only 10- to 100-fold compared to the WT. These results strongly suggest the involvement of Asp141 in the Mn2+-l-malate binding for the pigeon liver malic enzyme. The Asp194 and Asp464, which may be oxidized by nonspecific binding of Cu2+, are involved in the Mn2+-l-malate binding or catalysis indirectly by modulating the binding affinity of Asp141 with the Mn2+.

The relationship between the structure of a free ligand in solution and the structure of its bound form in a complex is of great importance to the understanding of the energetics and mechanism of molecular recognition and complex formation. In this study, we use a structure-based thermodynamic approach to study the dissociation of the complex between the toxin microcystin-LR (MLR) and the catalytic domain of protein phosphatase-1 (PP-1c) for which the crystal structure of the complex is known. We have calculated the thermodynamic parameters (enthalpy, entropy, heat capacity, and free energy) for the dissociation of the complex from its X-ray structure and found the calculated dissociation constant (4.0 × 10−11) to be in excellent agreement with the reported inhibitory constant (3.9 × 10−11). We have also calculated the thermodynamic parameters for the dissociation of 47 PP-1c:MLR complexes generated by docking an ensemble of NMR solution structures of MLR onto the crystal structure of PP-1c. In general, we observe that the lower the root-mean-square deviation (RMSD) of the docked complex (compared to the X-ray complex) the closer its free energy of dissociation (ΔG°d) is to that calculated from the X-ray complex. On the other hand, we note a significant scatter between the ΔG°d and the RMSD of the docked complexes. We have identified a group of seven docked complexes with ΔG°d values very close to the one calculated from the X-ray complex but with significantly dissimilar structures. The analysis of the corresponding enthalpy and entropy of dissociation shows a compensation effect suggesting that MLR molecules with significant structural variability can bind PP-1c and that substantial conformational flexibility in the PP-1c:MLR complex may exist in solution.

The structure of tick anticoagulant peptide (TAP) has been determined by X-ray crystallography at 1.6 Å resolution complexed with bovine pancreatic trypsin inhibitor (BPTI). The TAP–BPTI crystals are tetragonal, a = b = 46.87, c = 50.35 Å, space group P41, four complexes per unit cell. The TAP molecules are highly dipolar and form an intermolecular helical array along the c-axis with a diameter of about 45 Å. Individual TAP units interact in a head-to-tail fashion, the positive end of one molecule associating with the distal negative end of another, and vice versa. The BPTI molecules have a uniformly distributed positively charged surface that interacts extensively through 14 hydrogen bonds and two hydrogen bonded salt bridges with the helical groove around the helical TAP chains. Comparing the structure of TAP in TAP–BPTI with TAP bound to factor Xa(Xa) suggests a massive reorganization in the N-terminal tetrapeptide and the first disulfide loop of TAP (Cys5T –Cys15T) upon binding to Xa. The Tyr1TOH atom of TAP moves 14.2 Å to interact with Asp189 of the S1 specificity site, Arg3TCZ moves 5.0 Å with the guanidinium group forming a cation–π-electron complex in the S4 subsite of Xa, while Lys7TNZ differs in position by 10.6 Å in TAP-BPTI and TAP-Xa, all of which indicates a different pre-Xa–bound conformation for the N-terminal of TAP in its native state. In contrast to TAP, the BPTI structure of TAP–BPTI is practically the same as all those of previously determined structures of BPTI, only arginine and lysine side-chain conformations showing significant differences.

Protein molecules can accommodate a large number of mutations without noticeable effects on their stability and folding kinetics. On the other hand, some mutations can have quite strong effects on protein conformational properties. Such mutations either destabilize secondary structures, e.g., α-helices, are incompatible with close packing of protein hydrophobic cores, or lead to disruption of some specific interactions such as disulfide cross links, salt bridges, hydrogen bonds, or aromatic-aromatic contacts. The Met8 → Leu mutation in CMTI-I results in significant destabilization of the protein structure. This effect could hardly be expected since the mutation is highly conservative, and the side chain of residue 8 is situated on the protein surface. We show that the protein destabilization is caused by rearrangement of a hydrophobic cluster formed by side chains of residues 8, Ile6, and Leu17 that leads to partial breaking of a hydrogen bond formed by the amide group of Leu17 with water and to a reduction of a hydrophobic surface buried within the cluster. The mutation perturbs also the protein folding. In aerobic conditions the reduced wild-type protein folds effectively into its native structure, whereas more then 75% of the mutant molecules are trapped in various misfolded species. The main conclusion of this work is that conservative mutations of hydrophobic residues can destabilize a protein structure even if these residues are situated on the protein surface and partially accessible to water. Structural rearrangement of small hydrophobic clusters formed by such residues can lead to local changes in protein hydration, and consequently, can affect considerably protein stability and folding process.

The distance between Ca2+-binding site III in the C-terminal domain and Cys35 in the N-terminal domain in cardiac muscle troponin C (cTnC) was determined with a single-tryptophan mutant using bound Tb3+ as the energy donor and iodoacetamidotetramethylrhodamine linked to the cysteine residue as energy acceptor. The luminescence of bound Tb3+ was generated through sensitization by the tryptophan located in the 12-residue binding loop of site III upon irradiation at 295 nm, and this sensitized luminescence was the donor signal transferred to the acceptor. In the absence of bound cation at site II, the mean interdomain distance was found to be 48–49 Å regardless of whether the cTnC was unbound or bound to cardiac troponin I, or reconstituted into cardiac troponin. These results suggest that cTnC retains its overall length in the presence of bound target proteins. The distribution of the distances was wide (half-width >9 Å) and suggests considerable interdomain flexibility in isolated cTnC, but the distributions became narrower for cTnC in the complexes with the other subunits. In the presence of bound cation at the regulatory site II, the interdomain distance was shortened by 6 Å for cTnC, but without an effect on the half-width. The decrease in the mean distance was much smaller or negligible when cTnC was complexed with cTnI or cTnI and cTnT under the same conditions. Although free cTnC has considerable interdomain flexibility, this dynamics is slightly reduced in troponin. These results indicate that the transition from the relaxed state to an activated state in cardiac muscle is not accompanied by a gross alteration of the cTnC conformation in cardiac troponin.

Native state hydrogen exchange of cold shock protein A (CspA) has been characterized as a function of the denaturant urea and of the stabilizing agent trimethylamine N-oxide (TMAO). The structure of CspA has five strands of β-sheet. Strands β1–β4 have strongly protected amide protons that, based on experiments as a function of urea, exchange through a simple all-or-none global unfolding mechanism. By contrast, the protection of amide protons from strand β5 is too weak to measure in water. Strand β5 is hydrogen bonded to strands β3 and β4, both of which afford strong protection from solvent exchange. Gaussian network model (GNM) simulations, which assume that the degree of protection depends on tertiary contact density in the native structure, accurately predict the strong protection observed in strands β1–β4 but fail to account for the weak protection in strand β5. The most conspicuous feature of strand β5 is its low sequence hydrophobicity. In the presence of TMAO, there is an increase in the protection of strands β1–β4, and protection extends to amide protons in more hydrophilic segments of the protein, including strand β5 and the loops connecting the β-strands. TMAO stabilizes proteins by raising the free energy of the denatured state, due to highly unfavorable interactions between TMAO and the exposed peptide backbone. As such, the stabilizing effects of TMAO are expected to be relatively independent of sequence hydrophobicity. The present results suggest that the magnitude of solvent exchange protection depends more on solvent accessibility in the ensemble of exchange susceptible conformations than on the strength of hydrogen-bonding interactions in the native structure.

We describe a simple experimental approach for the rapid determination of protein global folds. This strategy utilizes site-directed spin labeling (SDSL) in combination with isotope enrichment to determine long-range distance restraints between amide protons and the unpaired electron of a nitroxide spin label using the paramagnetic effect on relaxation rates. The precision and accuracy of calculating a protein global fold from only paramagnetic effects have been demonstrated on barnase, a well-characterized protein. Two monocysteine derivatives of barnase, (H102C) and (H102A/Q15C), were 15N enriched, and the paramagnetic nitroxide spin label, MTSSL, attached to the single Cys residue of each. Measurement of amide 1H longitudinal relaxation times, in both the oxidized and reduced states, allowed the determination of the paramagnetic contribution to the relaxation processes. Correlation times were obtained from the frequency dependence of these relaxation processes at 800, 600, and 500 MHz. Distances in the range of 8 to 35 Å were calculated from the magnitude of the paramagnetic contribution to the relaxation processes and individual amide 1H correlation times. Distance restraints from the nitroxide spin to amide protons were used as restraints in structure calculations. Using nitroxide to amide 1H distances as long-range restraints and known secondary structure restraints, barnase global folds were calculated having backbone RMSDs <3 Å from the crystal structure. This approach makes it possible to rapidly obtain the overall topology of a protein using a limited number of paramagnetic distance restraints.

The Epidermal Growth Factor (EGF) receptor is a tyrosine kinase that mediates the biological effects of ligands such as EGF and transforming growth factor alpha. An understanding of the molecular basis of its action has been hindered by a lack of structural and mutational data on the receptor. We have constructed comparative models of the four extracellular domains of the EGF receptor that are based on the structure of the first three domains of the insulin-like growth factor-1 (IGF-1) receptor. The first and third domains of the EGF receptor, L1 and L2, are right-handed beta helices. The second and fourth domains of the EGF receptor, S1 and S2, consist of the modules held together by disulfide bonds, which, except for the first module of the S1 domain, form rod-like structures. The arrangement of the L1 and S1 domains of the model are similar to that of the first two domains of the IGF-1 receptor, whereas that of the L2 and S2 domains appear to be significantly different. Using the EGF receptor model and limited information from the literature, we have proposed a number of regions that may be involved in the functioning of the receptor. In particular, the faces containing the large beta sheets in the L1 and L2 domains have been suggested to be involved with ligand binding of EGF to its receptor.

The lipocalin superfamily of proteins functions in the binding and transport of a variety of important hydrophobic molecules. Tear lipocalin is a promiscuous lipid binding member of the family and serves as a paradigm to study the molecular determinants of ligand binding. Conserved regions in the lipocalins, such as the G strand and the F-G loop, may play an important role in ligand binding and delivery. We studied structural changes in the G strand of holo- and apo-tear lipocalin using spectroscopic methods including circular dichroism analysis and site-directed tryptophan fluorescence. Apo-tear lipocalin shows the same general structural characteristics as holo-tear lipocalin including alternating periodicity of a β-strand, orientation of amino acid residues 105, 103, 101, and 99 facing the cavity, and progressive depth in the cavity from residues 105 to 99. For amino acid residues facing the internal aspect of cavity, the presence of a ligand is associated with blue shifted spectra. The collisional rate constants indicate that these residues are not less exposed to solvent in holo-tear lipocalin than in apo-tear lipocalin. Rather the spectral blue shifts may be accounted for by a ligand induced rigidity in holo-TL.

Amino acid residues 94 and 95 are consistent with positions in the F–G loop and show greater exposure to solvent in the holo- than the apo-proteins. These findings are consistent with the general hypothesis that the F–G loop in the holo-proteins of the lipocalin family is available for receptor interactions and delivery of ligands to specific targets.

Site-directed tryptophan fluorescence was used in combination with a nitroxide spin labeled fatty acid analog to elucidate dynamic ligand interactions with specific amino acid residues. Collisional quenching constants of the nitroxide spin label provide evidence that at least three amino acids of the G strand residues interact with the ligand. Stern–Volmer plots are inconsistent with a ligand that is held in a static position in the calyx, but rather suggest that the ligand is in motion. The combination of site-directed tryptophan fluorescence with quenching by nitroxide labeled species has broad applicability in probing specific interactions in the solution structure of proteins and provides dynamic information that is not attainable by X-ray crystallography.

Electrospray ionization mass spectrometry (ESI-MS) was used to measure the binding of Cu2+ ions to synthetic peptides corresponding to sections of the sequence of the mature prion protein (PrP). ESI-MS demonstrates that Cu2+ is unique among divalent metal ions in binding to PrP and defines the location of the major Cu2+ binding site as the octarepeat region in the N-terminal domain, containing multiple copies of the repeat ProHisGlyGlyGlyTrpGlyGln. The stoichiometries of the complexes measured directly by ESI-MS are pH dependent: a peptide containing four octarepeats chelates two Cu2+ ions at pH 6 but four at pH 7.4. At the higher pH, the binding of multiple Cu2+ ions occurs with a high degree of cooperativity for peptides C-terminally extended to incorporate a fifth histidine. Dissociation constants for each Cu2+ ion binding to the octarepeat peptides, reported here for the first time, are mostly in the low micromolar range; for the addition of the third and fourth Cu2+ ions to the extended peptides at pH 7.4, KD's are <100 nM. N-terminal acetylation of the peptides caused some reduction in the stoichiometry of binding at both pH's. Cu2+ also binds to a peptide corresponding to the extreme N-terminus of PrP that precedes the octarepeats, arguing that this region of the sequence may also make a contribution to the Cu2+ complexation. Although the structure of the four-octarepeat peptide is not affected by pH changes in the absence of Cu2+, as judged by circular dichroism, Cu2+ binding induces a modest change at pH 6 and a major structural perturbation at pH 7.4. It is possible that PrP functions as a Cu2+ transporter by binding Cu2+ ions from the extracellular medium under physiologic conditions and then releasing some or all of this metal upon exposure to acidic pH in endosomes or secondary lysosomes.

The in vitro refolding of hen egg-white lysozyme is studied in the presence of various osmolytes. Proline is found to prevent aggregation during protein refolding. However, other osmolytes used in this study fail to exhibit a similar property. Experimental evidence suggests that proline inhibits protein aggregation by binding to folding intermediate(s) and trapping the folding intermediate(s) into enzymatically inactive, “aggregation-insensitive” state(s). However, elimination of proline from the refolded protein mixture results in significant recovery of the bacteriolytic activity. At higher concentrations (>1.5 M), proline is shown to form loose, higher-order molecular aggregate(s). The supramolecular assembly of proline is found to possess an amphipathic character. Formation of higher-order aggregates is believed to be crucial for proline to function as a protein folding aid. In addition to its role in osmoregulation under water stress conditions, the results of this study hint at the possibility of proline behaving as a protein folding chaperone.

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.

The role of the switch II region in the conformational transition of activation of Ha-ras-p21 has been investigated by mutating residues predicted to act as hinges for the conformational transition of this loop (Ala59, Gly60, and Gly75) (Diaz JF, Wroblowski B, Schlitter J, Engelborghs Y, 1997, Proteins 28:434–451), as well as mutating the catalytic residue Gln61. The proposed mutations of the hinge residues decrease the rate of the conformational transition of activation as measured by the binding of BeF3− to the GDP-p21 complex. Also, the thermodynamic parameters of the binding reaction are altered by a factor between three and five, depending on the temperature. (Due to changes in activation and reaction enthalpies, partially compensated by entropy changes.) The control mutation Q61H in which only the catalytic residue is changed has only a limited effect on the kinetic rate constants of the conformational transition and on the thermodynamic parameters of the reaction.

The fact that mutations of the hinge residues of the switch II region affect both the binding of the phosphate analog and the conformational transition of activation indicates that the switch II is implicated both in the early and the late states of the transition.

Mutant human lysozymes (Ile56Thr & Asp67His) have been reported to form amyloid deposits in the viscera. From the standpoint of understanding the mechanism of amyloid formation, we searched for conditions of amyloid formation in vitro using hen egg lysozyme, which has been extensively studied from a physicochemical standpoint. It was found that the circular dichroism spectra in the far-ultraviolet region of the hen egg lysozyme changed to those characteristic of a β-structure from the native α-helix rich spectrum in 90% ethanol solution. When the concentration of protein was increased to 10 mg/mL, the protein solution formed a gel in the presence of 90% ethanol, and precipitated on further addition of 10 mM NaCl. The precipitates were examined by electron microscopy, their ability to bind Congo red, and X-ray diffraction to determine whether amyloid fibrils were formed in the precipitates. Electron micrographs displayed unbranched protofilament with a diameter of ∼70 Å. The peak point of the difference spectrum for the Congo red binding assay was 541 nm, which is characteristic of amyloid fibrils. The X-ray diffraction pattern showed a sharp and intense diffraction ring at 4.7 Å, a reflection that arises from the interstrand spacing in β-sheets. These results indicate that the precipitates of hen egg lysozyme are amyloid protofilament, and that the amyloid protofilament formation of hen egg lysozyme closely follows upon the destruction of the helical and tertiary structures.

The human small GTPase, RhoA, expressed in Saccharomyces cerevisiae is post-translationally processed and, when co-expressed with its cytosolic inhibitory protein, RhoGDI, spontaneously forms a heterodimer in vivo. The RhoA/RhoGDI complex, purified to greater than 98% at high yield from the yeast cytosolic fraction, could be stoichiometrically ADP-ribosylated by Clostridium botulinum C3 exoenzyme, contained stoichiometric GDP, and could be nucleotide exchanged fully with [3H]GDP or partially with GTP in the presence of submicromolar Mg2+. The GTP-RhoA/RhoGDI complex hydrolyzed GTP with a rate constant of 4.5 × 10−5 s−1, considerably slower than free RhoA. Hydrolysis followed pseudo-first-order kinetics indicating that the RhoA hydrolyzing GTP was RhoGDI associated. The constitutively active G14V-RhoA mutant expressed as a complex with RhoGDI and purified without added nucleotide also bound stoichiometric guanine nucleotide: 95% contained GDP and 5% GTP. Microinjection of the GTP-bound G14V-RhoA/RhoGDI complex (but not the GDP form) into serum-starved Swiss 3T3 cells elicited formation of stress fibers and focal adhesions. In vitro, GTP-bound-RhoA spontaneously translocated from its complex with RhoGDI to liposomes, whereas GDP-RhoA did not. These results show that GTP-triggered translocation of RhoA from RhoGDI to a membrane, where it carries out its signaling function, is an intrinsic property of the RhoA/RhoGDI complex that does not require other protein factors or membrane receptors.

The stability of protein is defined not only by the hydrogen bonding, hydrophobic effect, van der Waals interactions, and salt bridges. Additional, much more subtle contributions to protein stability can arise from surface residues that change their properties upon unfolding. The recombinant major cold shock protein of Escherichia coli CspA an all-β protein unfolds reversible in a two-state manner, and behaves in all other respects as typical globular protein. However, the enthalpy of CspA unfolding strongly depends on the pH and buffer composition. Detailed analysis of the unfolding enthalpies as a function of pH and buffers with different heats of ionization shows that CspA unfolding in the pH range 5.5–9.0 is linked to protonation of an amino group. This amino group appears to be the N-terminal α-amino group of the CspA molecule. It undergoes a 1.6 U shift in pKa values between native and unfolded states. Although this shift in pKa is expected to contribute ∼5 kJ/mol to CspA stabilization energy the experimentally observed stabilization is only ∼1 kJ/mol. This discrepancy is related to a strong enthalpy–entropy compensation due, most likely, to the differences in hydration of the protonated and deprotonated forms of the α-amino group.

Electrospray ionization (ESI) mass spectra of both well-characterized and novel metallothioneins (MTs) from various species were recorded to explore their metal-ion-binding modes and stoichiometries. The ESI mass spectra of the zinc- and cadmium-binding MTs showed a single main peak corresponding to metal-to-protein ratios of 4, 6, or 7. These findings combined with data obtained by other methods suggest that these MTs bind zinc or cadmium in a single predominant form and are consistent with the presence of three- and four-metal clusters. An unstable copper-specific MT isoform from Roman snails (Helix pomatia) could be isolated intact and was shown to preferentially bind 12 copper ions. To obtain additional information on the formation and relative stability of metal-thiolate clusters in MTs, a mass spectrometric titration study was conducted. One to seven molar equivalents of zinc or of cadmium were added to metal-free human MT-2 at neutral pH, and the resulting complexes were measured by ESI mass spectrometry. These experiments revealed that the formation of the four-metal cluster and of the thermodynamically less stable three-metal cluster is sequential and largely cooperative for both zinc and cadmium. Minor intermediate forms between metal-free MT, Me4MT, and fully reconstituted Me7MT were also observed. The addition of increasing amounts of cadmium to metal-free blue crab MT-I resulted in prominent peaks whose masses were consistent with apoMT, Cd3MT, and Cd6MT, reflecting the known structure of this MT with two Me3Cys9 centers. In a similar reconstitution experiment performed with Caenorhabditits elegans MT-II, a series of signals corresponding to apoMT and Cd3MT to Cd6MT species were observed.