Discrete charge calculations of potentiometric titrations for globular proteins: sperm whale myoglobin, hemoglobin alpha chain, cytochrome c

Overview of the purification of recombinant proteins

When the first version of this unit was written in 1995, protein purification of recombinant proteins was based on a variety of standard chromatographic methods and approaches, many of which were described and mentioned throughout Current Protocols in Protein Science. In the interim, there has been a shift toward an almost universal usage of the affinity or fusion tag. This may not be the case for biotechnology manufacture where affinity tags can complicate producing proteins under regulatory conditions. Regardless of the protein expression system, questions are asked as to which and how many affinity tags to use, where to attach them in the protein, and whether to engineer a self-cleavage system or simply leave them on. We will briefly address some of these issues. Also, although this overview focuses on E.coli, protein expression and purification, other commonly used expression systems are mentioned and, apart from cell-breakage methods, protein purification methods and strategies are essentially the same.

Overview of the purification of recombinant proteins produced in Escherichia coli

The updated version of this unit presents an overview of recombinant protein purification with special emphasis on proteins expressed in E. coli. The first section deals with information pertinent to protein purification that can be derived from translation of the cDNA sequence. This is followed by a discussion of common problems associated with bacterial protein expression. A flow chart summarizes approaches for establishing solubility and localization of bacterially produced proteins. Purification strategies for both soluble and insoluble proteins are also reviewed. A section on glycoproteins produced in bacteria in the nonglycosylated state is included to emphasize that, although they may not be useful for in vivo studies, such proteins are well suited for structural studies. Finally, protein handling, scale and aims of purification, and specialized equipment needed for recombinant protein purification and characterization are discussed. The methodologies and approaches described here are essentially suitable for laboratory-scale operations.

Characteristics of amino acids

The physical properties of the amino acids determine the structure and function of the proteins in which they are found. This appendix presents some useful details and relevant physical characteristics of the natural amino acids and modified amino acids (i.e., chemical structures, explanations of the role these structures play in enzymes, relative hydrophobicities, shapes and volumes, and values for some common mass changes due to post-translational modifications).

Are the electrospray mass spectra of proteins related to their aqueous solution chemistry?

Institute for Environmental Chemistry, National Research Council of Canada, Ottawa, Ontario, Canada The shape of the profile described by the relative abundances of multiply charged ions of proteins in the electrospray mass spectrum can be described by means of one or more Gaussian functions. An aqueous solution equilibrium model of the distribution of multiply charged ions of equine cytochrome c and myoglobin has been shown to match qualitatively the shape of the distribution of these ions in an electrospray mass spectrum. Monotonic functions such as the quadrupole mass spectrometric transmission efficiency may alter the centroid of the profile, but the shape of the ion abundance pattern appears to be controlled by the aqueous solution chemistry of the proteins. NRCC No. 32938.

Charge effects on folded and unfolded proteins

We develop a theory for the effects of charge on the stabilization of globular proteins. The folding process is modeled as occurring through a fictitious intermediate state along a two-part thermodynamic pathway in which the molecule (i) increases its density and then (ii) rearranges its ionic groups to the protein surface. The equilibrium for the binding of protons in salt solutions is assumed to be driven by the electrical potential due to the charge distribution, in addition to the intrinsic binding affinity and bulk proton concentration. The potential is calculated for inside and outside a porous sphere model of the protein using the Poisson-Boltzmann relation, wherein the interior dielectric constant is taken to be a linear function of the chain density. The model predicts the slope of the titration curves for native myoglobin in agreement with experiments by Breslow and Gurd (1962). From the similar experiments on the unfolded state, and from the experiments of Privalov et al. (1986) on the intrinsic viscosity of the unfolded molecules, the theory shows that the unfolded state has a much higher density than a chain in a theta solvent and that the density increases with ionic strength. In addition, from the free energy of proton binding to the protein, we also calculate the electrostatic contributions to protein stability, a major contribution deriving from changes in ionization. We consider the example of the stability of myoglobin as a function of pH, ionic strength, and ionic groups buried in the native protein structure. We show that although maximum stability of most proteins should occur at their isoelectric point, the burial of nontitratable groups should lead to maximum stabilities at pH values other than the isoelectric point.

Properties of a cleaved two-chain form of recombinant human growth hormone

Escherichia coli cells transformed with plasmids engineered for the expression of recombinant human growth hormone as a secreted product also produced a proteolytically cleaved form of rhGH. This variant is isolated at a high resolution anion exchange chromatography stage during the manufacturing process. The higher isoelectric point of this form is demonstrated by isoelectric focusing and chromatofocusing and the two-chain nature by tryptic mapping, N- and C-terminal sequence analyses, and sodium dodecyl sulfate polyacrylamide gel electrophoresis. These data indicate that the single site of cleavage is between Thr-142 and Tyr-143, in contrast to the two-chain variant isolated from human pituitary glands, which has a clip after residue Phe-139. The recombinant two-chain form was further characterized by reversed-phase high performance liquid chromatography at both acidic and basic pHs. The assay utilizing bicarbonate-containing mobile phases was determined to be the most efficient and sensitive method. The bioactivity of this two-chain form was measured by the in vivo rat weight gain assay and by the in vitro Nb2 cell bioassay. Its immunological similarity to intact one-chain rhGH was demonstrated with an enzyme-linked immunosorbent assay.

Direct determination of macromolecular charge by equilibrium electrophoresis

Charge is a fundamental property of macromolecules in solution. However, estimation of the apparent charge on polyions has confounded science for decades. Presented here is a general method to determine directly the apparent charge on a polyion, regardless of its size or shape. This new method uses equilibrium electrophoresis, a procedure in which opposing solute flows from electrophoresis and from diffusion balance everywhere as the system reaches a steady-state distribution. The method uses only small quantities of materials, is nondestructive, and requires only simple, inexpensive instrumentation. Here we describe a prototype apparatus, demonstrate the phenomenon, and present experimental examples of the procedure.

Purification and characterization of recombinant human interleukin 4. Biological activities, receptor binding and the generation of monoclonal antibodies

A synthetic gene coding for human interleukin 4 (IL-4) was cloned and expressed in Saccharomyces cerevisiae (baker's yeast) as a C-terminal fusion protein with the yeast prepro alpha-mating factor sequence, resulting in secretion of mature IL-4 into the culture medium (0.6-0.8 micrograms/ml). A protocol was developed for purification of this protein. Crude cell-free conditioned medium was passed over a concanavalin A-Sepharose affinity column; bound proteins were eluted and further purified by S-Sepharose Fast Flow cation exchange and C18 reverse-phase h.p.l.c. Highly purified IL-4 was obtained by this method (0.3-0.4 mg per litre of culture) with a recovery of 51%. Thermospray liquid chromatography-mass spectrometry showed the C-terminal N-glycosylation site to be largely unmodified, and also showed that the N-terminus of the purified recombinant IL-4 (rIL-4) was authentic. Thiol titration revealed no free cysteine residues, implying that there are three disulphide groups, the positions of which remain to be determined. We have characterized the biological activities of the purified rIL-4. This material is active in B-cell co-stimulator assays, T-cell proliferation assays and in the induction of cell-surface expression of CD23 (the low-affinity receptor for IgE) on tonsillar B-cells. Half-maximal biological activity of the rIL-4 was achieved at a concentration of 120 pM. We have radioiodinated rIL-4 without loss of biological activity and performed equilibrium binding studies on Raji cells, a human B-cell line. The 125I-rIL-4 bound specifically to a single class of binding studies on Raji cells, a human B-cell line. The 125I-rIL-4 bound specifically to a single class of binding site with high affinity (Kd = 100 pM) and revealed 1100 receptors per cell. Receptor-ligand cross-linking studies demonstrated a single cell-surface receptor with an apparent molecular mass of 124 kDa. Two monoclonal antibodies have been raised to the human rIL-4, one of which blocks both the biological activity of rIL-4 and binding to its receptor.

Molecular dynamics effects on protein electrostatics

Electrostatic calculations have been carried out on a number of structural conformers of tuna cytochrome c. Conformers were generated using molecular dynamics simulations with a range of solvent simulating, macroscopic dielectric formalisms, and one solvent model that explicitly included solvent water molecules. Structures generated using the lowest dielectric models were relatively tight, with side chains collapsed on the surface, while those from the higher dielectric models had more internal and external fluidity, with surface side chains exploring a fuller range of conformational space. The average structure generated with the explicitly solvated model corresponded most closely with the crystal structure. Individual pK values, overall titration curves, and electrostatic potential surfaces were calculated for average structures and structures along each simulation. Differences between structural conformers within each simulation give rise to substantial changes in calculated local electrostatic interactions, resulting in pK value fluctuations for individual sites in the protein that vary by 0.3-2.0 pK units from the calculated time average. These variations are due to the thermal side chain reorientations that produce fluctuations in charge site separations. Properties like overall titration curves and pH dependent stability are not as sensitive to side chain fluctuations within a simulation, but there are substantial effects between simulations due to marked differences in average side chain behavior. These findings underscore the importance of proper dielectric formalism in molecular dynamics simulations when used to generate alternate solution structures from a crystal structure, and suggest that conformers significantly removed from the average structure have altered electrostatic properties that may prove important in episodic protein properties such as catalysis.

Enzyme-assisted semisynthesis of polypeptide active esters and their use

A method is described for the preparation of polypeptides activated uniquely at the C-terminus. The polypeptide is incubated in a concentrated solution of an amino acid active ester, the latter having its amino group free but adequately protected by protonation. The amino acid ester is coupled via its amino group to the C-terminus of the polypeptide by enzymic catalysis (reverse proteolysis). The resulting polypeptide C-terminal active ester is then isolated and coupled to a suitable amino component (generally a polypeptide) in a subsequent chemical coupling. The method appears to be generally applicable; fragments of horse heart cytochrome c, and porcine insulin, are used as examples. Two new analogues of cytochrome c have been prepared by using this method, with yields of up to 60% in the final coupling. Scope and limitations of the method are discussed.

Preparation and characterization of bovine growth hormones produced in recombinant Escherichia coli

Two analogues of bovine growth hormone (BGH) have been produced in Escherichia coli by recombinant DNA techniques. In analogue Delta-1, the N-terminal alanine residue of the full-length bovine sequence is replaced by methionine. In analogue Delta-9, which is expressed at much higher levels than is Delta-1, the full-length bovine sequence is truncated at the N-terminus by eight residues and there is a serine-for-glycine substitution in the first position of the truncated protein. Both analogues, which were characterized by isoelectric focusing (i.e.f.), polyacrylamide-gel electrophoresis in the presence of SDS (SDS/PAGE), amino acid analysis and N-terminal amino acid sequence determination using combined g.l.c.-m.s., are compared with BGH isolated from pituitaries. In contrast with pituitary-derived BGH, the recombinant-derived proteins are homogeneous on SDS/PAGE and on i.e.f. In a radioimmunoassay, a radioreceptor assay and a bioassay in vivo (rat tibia), Delta-9 BGH showed very similar characteristics to the pituitary-derived hormone. Similar results have also been obtained with the Delta-1 analogue.

Chromatofocusing of N-terminally processed forms of proteins. Isolation and characterization of two forms of interleukin-1 beta and of bovine growth hormone

Using chromatofocusing, two fractions have been obtained from recombinant-derived interleukin-1 beta (IL-1 beta) and from pituitary-derived bovine growth hormone (BGH). The forms of both proteins responsible for these fractions have been characterized by N-terminal and C-terminal amino acid sequence determination. Recombinant IL-1 beta, as a mixture of correctly processed polypeptide and an N-terminally methionylated form, was resolved rapidly by chromatofocusing. BGH was resolved into the full-length polypeptide commencing Ala-Phe-Pro-Ala-Met-Ser-Leu- and a form truncated at the N-terminus by four amino acid residues, which thus commences Met-Ser-Leu-; the fraction containing the truncated form also contains a species having N-terminal Phe-Pro-Ala-. These results, and the possible generality of the separation, are discussed.

Electrostatic analysis of the interaction of cytochrome c with native and dimethyl ester heme substituted cytochrome b5

The stability of the complex formed between cytochrome c and dimethyl ester heme substituted cytochrome b5 (DME-cytochrome b5) has been determined under a variety of experimental conditions to evaluate the role of the cytochrome b5 heme propionate groups in the interaction of the two native proteins. Interaction between cytochrome c and the modified cytochrome b5 was found to produce a difference spectrum in the visible range that is very similar to that generated by the interaction of the native proteins and that can be used to monitor complex formation between the two proteins. At pH 8 [25 degrees C (HEPPS), I = 5 mM], DME-cytochrome b5 and cytochrome c form a 1:1 complex with an association constant KA of 3 (1) X 10(6) M-1. This pH is the optimal pH for complex formation between these two proteins and is significantly higher than that observed for the interaction between the two native proteins. The stability of the complex formed between DME-cytochrome b5 and cytochrome c is strongly dependent on ionic strength with KA ranging from 2.4 X 10(7) M-1 at I = 1 mM to 8.2 X 10(4) M-1 at I = 13 mM [pH 8.0 (HEPPS), 25 degrees C]. Calculations for the native, trypsin-solubilized form of cytochrome b5 and cytochrome c confirm that the intermolecular complex proposed by Salemme [Salemme, F. R. (1976) J. Mol. Biol. 102, 563] describes the protein-protein orientation that is electrostatically favored at neutral pH.(ABSTRACT TRUNCATED AT 250 WORDS)

pH-dependent processes in proteins

Recent improvements in the understanding of electrostatic interactions in proteins serve as a focus for the general topic of pH-dependent processes in proteins. The general importance of pH-dependent processes is first set out in terms of hydrogen ion equilibria, stability, ligand interactions, assembly, dynamics, and events in related molecular systems. The development of various theoretical treatments includes various formalisms in addition to the solvent interface model developed by Shire et al. as an extension of the Tanford-Kirkwood treatment. A number of detailed applications of the model are presented and future potentialities are sketched.

Anion binding and pH-dependent electrostatic effects in ribonuclease

The solvent-accessibility-modified, Tanford-Kirkwood, discrete charge model for electrostatic effects is applied to both ribonuclease A and ribonuclease S. The behavior of individual titratable sites and the pH-dependent free energy of denaturation are correctly predicted. The use of the solvent-accessibility factor in reducing charge-site interactions introduces a higher Coulombic shielding for solvent-exposed sites. This shielding is interpreted as a higher local strength or alternatively a higher effective dielectric constant. Specific anion binding sites are determined by locating areas of high positive electrostatic potential at the protein solvent interface. The potential and thus the anion affinity of a given site are calculated and shown to vary with the pH-dependent charge array. pH-dependent anion binding constants are calculated for the ribonuclease S active site. These binding constants and the predicted response of the active-site histidine pK1/2 values to anion binding are shown to agree with experimental determinations.

Effect on protein stability of reversing the charge on amino groups

The amino groups of beta-lactoglobulins A and B, cytochrome c and ribonuclease were progressively converted to acidic groups by reaction with succinic anhydride. The mixtures of modified proteins generated in this way were analyzed by urea-gradient electrophoresis, which separates the molecules on the basis of their net charge and demonstrates visually their urea-induced unfolding transitions. Molecules succinylated to varying extents were resolved by the electrophoresis, so purification of the many modified species was not required. It is demonstrated that accurate estimates of the stability of the folded state of an individual species may be estimated very easily from its urea-gradient electrophoretic pattern. Changes in ionization of the protein upon unfolding may also be detected. The general electrostatic effect of varying the net charge on these proteins was small. Converting the normally basic ribonuclease and cytochrome c to neutral and then to acidic proteins caused the net stabilities of their folded states to vary by no more than a few kJ/mol. However, specific interactions between a few ionized groups appear to be more important in some instances. Succinylation of the 19th, and final, lysine residue of cytochrome c produced unfolding even in the absence of urea, whereas reaction of the first 18 had very little effect. Reaction of the initial amino groups of beta-lactoglobulins A and B produced a small increase in stability in a few instances, a decrease in others; modification of more than about ten groups abruptly caused unfolding in the absence of urea.

Increase in apparent compressibility of cytochrome c upon oxidation

The apparent molal adiabatic compressibilities of ferri- and ferrocytochrome c have been determined from measurements of density and sound velocity. The values found were +2.99 X 10(-8) and -2.40 X 10(-8) cm5 mol-1 dyne-1 for the ferri and ferro forms, respectively. Experiments were performed on identical solutions containing either the oxidized or reduced form of protein. Solutions of ferricytochrome c were found to have significantly greater adiabatic compressibility than equivalent solutions of ferrocytochrome c at 25 degrees C and pH 7.15. The remarkable similarity of the three-dimensional structures of the ferri and ferro proteins [Takano, T. & Dickerson, R.E. (1980) Proc. Natl. Acad. Sci. USA 77, 6371-6375] strongly suggests that this difference in compressibility is due to an increase in volume fluctuations within ferricytochrome c relative to the ferro form rather than a change in equilibrium structure or hydration. Such a difference in the dynamic properties of the structures is consistent with both the crystallographic thermal B factors and the observed increase in amide hydrogen exchange kinetics when ferrocytochrome c is oxidized. The relative magnitude of the root mean square volume fluctuations is approximated from an ideal solution treatment of the compressibility data and yields a ratio of delta Vrms (ferri cyt c)/ delta Vrms (ferro cyt c) = 1.3.

Electrostatic contributions to the energetics of dimer-tetramer assembly in human hemoglobin: pH dependence and effect of specifically bound chloride ions

The pH dependence and effects of specifically bound chloride ions on the electrostatic contribution to the energetics of human hemoglobin dimer-tetramer assembly were computed for deoxy- and liganded hemoglobin. In the absence of bound chloride, the electrostatic contribution models the observed contrasting pH dependence of dimer-tetramer assembly for deoxy- and oxyhemoglobin. The effect of specifically bound chloride on the computations depends on the number and placement of the anions. Deoxy assembly shows a greater sensitivity to anion binding, with effects propagating as far as 32 A from the binding site. This sensitivity suggests a mechanism for electronic communication with the heme. At pH 7.4, 24-34% of the experimental value for deoxy and 73-85% for oxy dimer-tetramer assembly stabilization are predicted. Together with the findings of Chu and Ackers [Chu, A. H., & Ackers, G. K. (1981) J. Biol. Chem. 256, 1199] and other recent work, these results suggest that salt bridge formation is not the dominant energetic factor favoring deoxyhemoglobin dimer-tetramer assembly. Results of this work suggest that the marked electrostatic stabilization favoring oxy dimer-tetramer assembly may be a significant contributor to the quaternary enhancement observed in assembly reactions whereas the nonelectrostatic factors favoring deoxy dimer-tetramer assembly may be largely responsible for quaternary constraint.

Molecular interpretation of kinetic-ionic strength effects

A recent and important approach to investigating electron transfer mechanisms of redox proteins has been through kinetic-ionic strength studies. There is, however, significant controversy as to whether such studies (1) yield information regarding the charge (or location) of the electron transfer site or (2) more simply reflect the influence of net or overall protein charge on the electrostatic interactions. A critical analysis using different theoretical approaches is made of our recent work and of the bulk of the published non-physiological small molecule-protein and protein-protein kinetic ionic strength studies; it is concluded that (1) the approximated Bronsted-Debye-Huckel equation can not be used at all for protein redox reactions, (2) irrespective of the theoretical approaches discussed, such studies do not provide information regarding the charge of the electron transfer site, (3) it is the net charge of the reactants that control the electrostatic interactions, (4) both the equation derived by Wherland and Gray and the full Bronsted-Debye-Huckel equation provide reasonably good approximations of net protein charge, (5) pH changes quantitatively modulate net protein charge, and (6) thus, protein redox rates need to be electrostatically corrected if relevant interpretations of kinetic-ionic strength experiments are to be made.

Reassignment of residue 122 in the myoglobin from the killer whale, Orcinus orca

The complete amino acid sequence of the major component myoglobin from killer whale, Orcinus orca, was determined by automated Edman degradation. In this study residue 122 was found to be glutamic acid instead of glutamine as was originally reported (Castillo et al. 1977). This reassignment affects the phylogenetic relationship of killer whale myoglobin with the myoglobins from other closely related cetacean species and also affects studies concerned with the physical parameters of the protein.

Electrostatic stabilization in sperm whale and harbor seal myoglobins. Identification of groups primarily responsible for changes in anchoring of the A helix

The compact, largely helical structure of sperm whale and harbor seal myoglobins undergoes an abrupt one-step transition between pH 4.5 and 3.5 as monitored by changes in either the heme Soret band absorbance or circular dichroism probes of secondary structure, for which a modified Tanford-Kirkwood theory provides identification of certain dominant electrostatic interactions responsible for the loss of stability. A similar treatment permits identification of the electrostatic interactions primarily responsible for a process in which the anchoring of the A helix to other parts of the molecule is weakened. This process is detected with both myoglobins, in a pH range approximately 1 unit higher than the onset of the overall unfolding process, through changes in the circular dichroic spectra near 295 nm which correspond to the L1 O-O band of the only two tryptophan residues in these proteins, residues 7 and 14. In each case protonation of certain sites in neighboring parts of the molecule can be identified as producing destabilizing interactions with components of the A helix, particularly with lysine 6.

Evolution of the amino acid substitution in the mammalian myoglobin gene

Multivariate statistical analyses were applied to 16 physical and chemical properties of amino acids. Four of these properties; volume, polarity, isoelectric point (charge), and hydrophobicity were found to explain adequately 96% of the total variance of amino acid attributes. Using these four quantitative measures of amino acid properties, a structural discriminate function in the form of a weighted difference sum of squares equation was developed. The discriminate function is weighted by the location of each particular residue within a given tertiary structure and yields a numerical discriminate or difference value for the replacement of these residues by different amino acids. This resulting discriminate value represents an expression of the perturbation in the local positional environment of a protein when an amino acid substitution occurs. With the use of this structural discriminate function, a residue by residue comparison of the known mammalian myoglobin sequences was carried out in an attempt to elucidate the positions of possible deviations from the known tertiary structure of sperm whale myoglobin. Only 11 of the 153 residue positions in myoglobin demonstrated possible structural deviations. From this analysis, indices of difference were calculated for all amino acid exchanges between the various myoglobins. All comparisons yielded indices of difference that were considerably lower than would be expected if mutations had been fixed at random, even if the organization of the genetic code is taken into consideration. On the basis of these results, it is inferred that some form of selection has acted in the evolution of mammalian myoglobins to favor amino acid substitutions that are compatible with the retention of the original conformation of the protein.

Effect of a molecular dipole on the ionic strength dependence of a biomolecular rate constant. Identification of the site of reaction

A theory is proposed for determining the location of a reaction site on a protein of known tertiary structure with an asymmetric charge distribution by an analysis of the effect of ionic strength on the rate of reaction of the protein with a small ion, using equations of Brønsted (J. N. Brønsted, 1922, Z. Phys, Chem. 102:169-207), Debye and Hückel (P. Debye and E. Hückel, 1923, Phys. Z. 24:185-206), and Kirkwood (J. G. Kirkwood, 1934, J. Chem. Phys. 2:351-361). The theory is based on the fact that the dipole moment of the transition complex differs from that of the protein, which will be reflected in the ionic strength dependence of the reaction. The location of the small ion with respect to the dipole axis of the protein can be calculated from this difference. For protein-protein reactions, an a priori assumption has to be made about the orientation of one of the reaction partners, since many different orientations of the reactants with respect to each other result in dipole moments of the same magnitude.

Partitioning of electrostatic and conformational contributions in the redox reactions of modified cytochromes c

The reduction of acetylated, fully succinylated and dicarboxymethyl horse cytochromes c by the radicals CH3CH(OH), CO2.-, O2.-, and e-aq' and the oxidation of the reduced cytochrome c derivatives by Fe(CN)3-6 were studied using the pulse radiolysis technique. Many of the reactions were also examined as a function of ionic strength. By obtaining rate constants for the reactions of differently charged small molecules redox agents with the differently charged cytochrome c derivatives at both zero ionic strength and infinite ionic strength, electrostatic and conformational contributions to the electron transfer mechanism were effectively partioned from each other in some cases. In regard to cytochrome c electron transfer mechanism, the results, especially those for which conformational influences predominate, are supportive of the electron being transferred in the heme edge region.

Electrostatic stabilization in myoglobin. Interactive free energies between individual sites

The pattern of electrostatic interactions between pairs of charge sites in sperm whale ferrimyoglobin was examined as a function of pH in terms of proton site occupancy, static solvent accessibility, and distance of separation. By grouping all examples of the most stabilizing interactions and all examples of the most destabilizing interactions, we can easily show that at pH 7.50 the former is much stronger; that is, the negative contributions to electrostatic free energy far outweigh the positive contributions. Much of the electrostatic energy of stabilization in native myoglobin is provided by specific charge-pair partners that are very highly conserved among 53 mammalian myoglobin species and is invariant substantially from pH 8.5 to 3.5. Destablizing interactions that become most significant, but not actually dominant, near the acid unfolding pH range can be recognized in emerging clusters of uncompensated positive charges. Binding of azide ion by the heme iron effectively reduces the most prominent destabilizing set of such interactions. In general, thoe charged residues that experience the largest summed stabilizing interactions with other groups are the most conserved between species. The histidine residues, however, show their best correlation of conservation with low values of static accessibility. Although histidine residue 64 has an effective pK corresponding to the midpoint of the unfolding transition near pH 4.2 at an ionic strength of 0.10 M and so might be called a "trigger group", its interactions contribute only a modest fraction of the overall pH-dependent free energy change. An examination of the primary stabilizing interactions represented by the charge-pair partners indicates a probably major role of electrostatic interactions in the nucleation and docking stages of the condensation of the polypeptide chain into the compact native structure.

Complete amino acid sequence of the myoglobin from the Pacific spotted dolphin, Stenella attenuata graffmani

The complete amino acid sequence of the major component myoglobin from the Pacific spotted dolphin, Stenella attenuata graffmani, was determined by the automated Edman degradation of several large peptides obtained by specific cleavage of the protein. The acetimidated apomyoglobin was selectively cleaved at its two methionyl residues with cyanogen bromide and at its three arginyl residues by trypsin. By subjecting four of these peptides and the apomyoglobin to automated Edman degradation, over 80% of the primary structure of the protein was obtained. The remainder of the covalent structure was determined by the sequence analysis of peptides that resulted from further digestion of the central cyanogen bromide fragment. This fragment was cleaved at its glutamyl residues with staphylococcal protease and its lysyl residues with trypsin. The action of trypsin was restricted to the lysyl residues by chemical modification of the single arginyl residue of the fragment with 1,2-cyclohexanedione. The primary structure of this myoglobin proved to be identical with that from the Atlantic bottlenosed dolphin and Pacific common dolphin but differs from the myoglobins of the killer whale and pilot whale at two positions. The above sequence identities and differences reflect the close taxonomic relationship of these five species of Cetacea.