Proteins as random coils. II. Hydrogen ion titration curve of ribonuclease in 6 M guanidine hydrochloride

Probing Denatured State Conformational Bias in a Three-Helix Bundle, UBA(2), Using a Cytochrome c Fusion Protein

Previous work with the four-helix-bundle protein cytochrome c' from Rhodopseudomonas palustris using histidine-heme loop formation methods revealed fold-specific deviations from random coil behavior in its denatured state ensemble. To examine the generality of this finding, we extend this work to a three-helix-bundle polypeptide, the second ubiquitin-associated domain, UBA(2), of the human DNA excision repair protein. We use yeast iso-1-cytochrome c as a scaffold, fusing the UBA(2) domain at the N-terminus of iso-1-cytochrome c. We have engineered histidine into highly solvent accessible positions of UBA(2), creating six single histidine variants. Guanidine hydrochloride denaturation studies show that the UBA(2)-cytochrome c fusion protein unfolds in a three-state process with iso-1-cytochrome c unfolding first. Furthermore, engineered histidine residues in UBA(2) strongly destabilize the iso-1-cytochrome c domain. Equilibrium and kinetic histidine-heme loop formation measurements in the denatured state at 4 and 6 M guanidine hydrochloride show that loop stability decreases as the size of the histidine-heme loop increases, in accord with the Jacobson-Stockmayer equation. However, we observe that the His27-heme loop is both more stable than expected from the Jacobson-Stockmayer relationship and breaks more slowly than expected. These results show that the sequence near His27, which is in the reverse turn between helices 2 and 3 of UBA(2), is prone to persistent interactions in the denatured state. Therefore, consistent with our results for cytochrome c', this reverse turn sequence may help to establish the topology of this fold by biasing the conformational distribution of the denatured state.

Plasmodium falciparum merozoite surface protein 2: epitope mapping and fine specificity of human antibody response against non-polymorphic domains

Background

Two long synthetic peptides representing the dimorphic and constant C-terminal domains of the two allelic families of Plasmodium falciparum merozoite surface proteins 2 are considered promising malaria vaccine candidates. The aim of the current study is to characterize the immune response (epitope mapping) in naturally exposed individuals and relate immune responses to the risk of clinical malaria.

Methods

To optimize their construction, the fine specificity of human serum antibodies from donors of different age, sex and living in four distinct endemic regions was determined in ELISA by using overlapping 20 mer peptides covering the two domains. Immune purified antibodies were used in Western blot and immunofluorescence assay to recognize native parasite derivate proteins.

Results

Immunodominant epitopes were characterized, and their distribution was similar irrespective of geographic origin, age group and gender. Acquisition of a 3D7 family and constant region-specific immune response and antibody avidity maturation occur early in life while a longer period is needed for the corresponding FC27 family response. In addition, the antibody response to individual epitopes within the 3D7 family-specific region contributes to protection from malaria infection with different statistical weight. It is also illustrated that affinity-purified antibodies against the dimorphic or constant regions recognized homologous and heterologous parasites in immunofluorescence and homologous and heterologous MSP2 and other polypeptides in Western blot.

Conclusion

Data from this current study may contribute to a development of MSP2 vaccine candidates based on conserved and dimorphic regions thus bypassing the complexity of vaccine development related to the polymorphism of full-length MSP2.

Electronic supplementary material

The online version of this article (doi:10.1186/1475-2875-13-510) contains supplementary material, which is available to authorized users.

Reaction of bovine and equine growth hormones with tetranitromethane

Bovine and equine growth hormones were chemically modified with tetranitromethane, at pH 7.4 during 5 h and at pH 8.0 in the presence of 8 M urea during 1 h. a) Both hormones have very similar but not identical reactivities. b) The nitration of the reactive tyrosines and tryptophan residues at pH 7.4 produces no detectable changes in their immunological or somatotrophic activities. C) The nitration of all tyrosine residues in both hormones gives rise to a complete loss of somatographic activity with no alteration of the immunological activity.

The secondary structure of histone IV and its stability

The secondary structure within histone IV and its fragments obtained by cyanogen bromide (CNBr) and cleavage at Met 84 has been examined by circular dichroism and spectophotometric pH titration measurements. These studies have confirmed the existence of stable secondary structure within the C-terminal fragment of histone IV (C-peptide which can be perturbed only by 6M urea at pH greater than 8 or 8 M guanidine-HCL. In contrast, the N-terminal fragment (N-peptide) appears to lack significant secondary structure at low ionic strengths but acquires approximately 15% betasheet conformation and 5% alpha-helix upon aggregation at ionic strengths larger than or equal to 0.4. The rates of nitration of the N- and C-peptides by tetranitromethane (TNM) have also been measured as a function of ionic strengths. Under comparable conditions, the rate constant for nitration of the N-peptide was found to be about six times greater than that for the C-peptide, further evidence in support of the presence of stable secondary structure within the C-terminal region of histone IV. After binding these histone IV fragments to DNA, however, the nitration reaction rate constants for the N- and C-peptide in the bound form are found to be 2% and 27% of the corresponding free peptides. Reconstituted nucleohistone IV is about 10% as reactive to TNM as histone IV at comparable ionic strength.

Electrostatic effects in unfolded staphylococcal nuclease

Structure-based calculations of pKa values and electrostatic free energies of proteins assume that electrostatic effects in the unfolded state are negligible. In light of experimental evidence showing that this assumption is invalid for many proteins, and with increasing awareness that the unfolded state is more structured and compact than previously thought, a detailed examination of electrostatic effects in unfolded proteins is warranted. Here we address this issue with structure-based calculations of electrostatic interactions in unfolded staphylococcal nuclease. The approach involves the generation of ensembles of structures representing the unfolded state, and calculation of Coulomb energies to Boltzmann weight the unfolded state ensembles. Four different structural models of the unfolded state were tested. Experimental proton binding data measured with a variant of nuclease that is unfolded under native conditions were used to establish the validity of the calculations. These calculations suggest that weak Coulomb interactions are an unavoidable property of unfolded proteins. At neutral pH, the interactions are too weak to organize the unfolded state; however, at extreme pH values, where the protein has a significant net charge, the combined action of a large number of weak repulsive interactions can lead to the expansion of the unfolded state. The calculated pKa values of ionizable groups in the unfolded state are similar but not identical to the values in small peptides in water. These studies suggest that the accuracy of structure-based calculations of electrostatic contributions to stability cannot be improved unless electrostatic effects in the unfolded state are calculated explicitly.

Ligand effects on the protein ensemble: unifying the descriptions of ligand binding, local conformational fluctuations, and protein stability

Detailed description of the structural and physical basis of allostery, cooperativity, and other manifestations of long-range communication between binding sites in proteins remains elusive. Here we describe an ensemble-based structural-thermodynamic model capable of treating explicitly the coupling between ligand binding reactions, local fluctuations in structure, and global conformational transitions. The H(+) binding reactions of staphylococcal nuclease and the effects of pH on its stability were used to illustrate the properties of proteins that can be described quantitatively with this model. Each microstate in the native ensemble was modeled to have dual structural character; some regions were treated as folded and retained the same atomic geometry as in the crystallographic structure while other regions were treated thermodynamically as if they were unfolded. Two sets of pK(a) values were used to describe the affinity of each H(+) binding site. One set, calculated with a standard continuum electrostatics method, describes H(+) binding to sites in folded parts of the protein. A second set of pK(a) values, obtained from model compounds in water, was used to describe H(+) binding to sites in unfolded regions. An empirical free energy function, parameterized to reproduce folding thermodynamics measured by differential scanning calorimetry, was used to calculate the probability of each microstate. The effects of pH on the distribution of microstates were determined by the H(+) binding properties of each microstate. The validity of the calculations was established by comparison with a number of different experimental observables.

Modeling the electrophoretic mobility and diffusion of weakly charged peptides

A bead model to determine the electrophoretic mobilities and translational diffusion constants of weakly charged peptides is developed that is based on a approximate structural model of peptides and is also grounded in electrohydrodynamic theory. A peptide made up of X amino acids is modeled as N=2X beads with 2 beads representing each amino acid in the chain. For the two beads representing a particular amino acid in a peptide, the radius of one bead is set to one-half the nearest neighbor Calpha-Calpha distance, and the radius of the other bead is chosen on the basis of the diffusion constant of the free amino acid. Peptide conformations, which are defined by a set of psi-phi dihedral angles, are randomly generated by using the transformation matrix approach of Flory (Flory, P. Statistical Mechanics of Chain Molecules; John Wiley: New York, 1969) and rejecting conformations which result in bead overlap. The mobility and diffusion constants are computed for each conformation and at least 100 independent conformations are examined for each peptide. In general, the mobility is found to depend only weakly on peptide conformation. Model and experimental mobilities are compared by examining the data of Janini and co-workers (Janini, G.; et al. J. Chromatogr. 1999, 848, 417-433). A total of 58 peptides consisting of from 2 to 39 amino acids are considered. The average relative error between experimental and model mobilities is found to be 1.0% and the rms relative error 7.7%. In specific cases, the discrepancy can be substantial and possible reasons for this are discussed. It should be emphasized that the input parameters of the peptide model are totally independent of experimental mobilities. It is hoped that the peptide model developed here will be useful in the prediction of peptide mobility as well as in using peptide mobilities to extract information about peptide structure, conformation, and charge. Finally, we show how simultaneous measurements of translational diffusion and mobility can be used to estimate peptide charge.

Free energy of proton binding in proteins

In this article we use literature data on the titration of denatured ribonuclease to test the accuracy of proton-binding distributions obtained using our recent approach employing moments. We find that using only the local slope of the titration curve at a small number of points (five, for example) we can reproduce the detailed proton-binding distribution at all pH values. Our method gives the complete proton-binding polynomial for a given protein and each coefficient in this polynomial in turn yields the free energy for binding a given number of protons in all ways to the protein. Using these net free energies, we can then compute the average proton-binding free energy per proton as a function of the fraction of protons bound. We find that this function is remarkably similar for different proteins, even for proteins that exhibit quite different titration behavior. For the special case of binding to independent sites, we obtain simple relations for the first and last terms in the free energy per-proton function. For this special case we also can calculate the distribution functions giving the probability that a molecule has a given number of positive or negative charges and the joint distribution that a molecule simultaneously has a given number of positive and negative charge.

Steady-state tryptophan fluorescence spectroscopy study to probe tertiary structure of proteins in solid powders

The purpose of this work was to obtain information about protein tertiary structure in solid state by using steady state tryptophan (Trp) fluorescence emission spectroscopy on protein powders. Beta-lactoglobulin (betaLg) and interferon alpha-2a (IFN) powder samples were studied by fluorescence spectroscopy using a front surface sample holder. Two different sets of dried betaLg samples were prepared by vacuum drying of solutions: one containing betaLg, and the other containing a mixture of betaLg and guanidine hydrochloride. Dried IFN samples were prepared by vacuum drying of IFN solutions and by vacuum drying of polyethylene glycol precipitated IFN. The results obtained from solid samples were compared with the emission scans of these proteins in solutions. The emission scans obtained from protein powders were slightly blue-shifted compared to the solution spectra due to the absence of water. The emission scans were red-shifted for betaLg samples dried from solutions containing GuHCl. The magnitude of the shifts in lambda(max) depended on the extent of drying of the samples, which was attributed to the crystallization of GuHCl during the drying process. The shifts in the lambda(max) of the Trp emission spectrum are associated with the changes in the tertiary structure of betaLg. In the case of IFN, the emission scans obtained from PEG-precipitated and dried sample were different compared to the emission scans obtained from IFN in solution and from vacuum dried IFN. The double peaks observed in this sample were attributed to the unfolding of the protein. In the presence of trehalose, the two peaks converged to form a single peak, which was similar to solution emission spectra, whereas no change was observed in the presence of mannitol. We conclude that Trp fluorescence spectroscopy provides a simple and reliable means to characterize Trp microenvironment in protein powders that is related to the tertiary conformation of proteins in the solid state. This study shows that the use of fluorescence spectroscopy of proteins can be extended from simple protein aqueous solutions to protein powders, precipitates, and semidried protein samples to gain understanding of protein tertiary structure in these physical states.

Residual charge interactions in unfolded staphylococcal nuclease can be explained by the Gaussian-chain model

The discrepancy of the pH dependence of the unfolding free energy for staphylococcal nuclease from what is expected from an idealized model for the unfolded state is accounted for by the recently developed Gaussian-chain model. Residual electrostatic effects in the unfolded state are attributed to nonspecific interactions dominated by charges close along the sequence. The dominance of nonspecific local interactions appears to be supported by some experimental evidence.

Polyproline II structure in a sequence of seven alanine residues

A sequence of seven alanine residues-too short to form an alpha-helix and whose side chains do not interact with each other-is a particularly simple model for testing the common description of denatured proteins as structureless random coils. The (3)J(HN alpha) coupling constants of individual alanine residues have been measured from 2 to 56 degrees C by using isotopically labeled samples. The results display a thermal transition between different backbone conformations, which is confirmed by CD spectra. The NMR results suggest that polyproline II is the dominant conformation at 2 degrees C and the content of beta strand is increased by approximately 10% at 55 degrees C relative to that at 2 degrees C. The polyproline II conformation is consistent with recent studies of short alanine peptides, including structure prediction by ab initio quantum mechanics and solution structures for both a blocked alanine dipeptide and an alanine tripeptide. CD and other optical spectroscopies have found structure in longer "random coil" peptides and have implicated polyproline II, which is a major backbone conformation in residues within loop regions of protein structures. Our result suggests that the backbone conformational entropy in alanine peptides is considerably smaller than estimated by the random coil model. New thermodynamic data confirm this suggestion: the entropy loss on alanine helix formation is only 2.2 entropy units per residue.

pH corrections and protein ionization in water/guanidinium chloride

More than 30 years ago, Nozaki and Tanford reported that the pK values for several amino acids and simple substances in 6 M guanidinium chloride differed little from the corresponding values in low salt (Nozaki, Y., and C. Tanford. 1967. J. Am. Chem. Soc. 89:736-742). This puzzling and counter-intuitive result hinders attempts to understand and predict the proton uptake/release behavior of proteins in guanidinium chloride solutions, behavior which may determine whether the DeltaG(N-D) values obtained from guanidinium chloride-induced denaturation data can actually be interpreted as the Gibbs energy difference between the native and denatured states (Bolen, D. W., and M. Yang. 2000. Biochemistry. 39:15208-15216). We show in this work that the Nozaki-Tanford result can be traced back to the fact that glass-electrode pH meter readings in water/guanidinium chloride do not equal true pH values. We determine the correction factors required to convert pH meter readings in water/guanidinium chloride into true pH values and show that, when these corrections are applied, the effect of guanidinium chloride on the pK values of simple substances is found to be significant and similar to that of NaCl. The results reported here allow us to propose plausible guanidinium chloride concentration dependencies for the pK values of carboxylic acids in proteins and, on their basis, to reproduce qualitatively the proton uptake/release behavior for the native and denatured states of several proteins (ribonuclease A, alpha-chymotrypsin, staphylococcal nuclease) in guanidinium chloride solutions. Finally, the implications of the pH correction for the experimental characterization of protein folding energetics are briefly discussed.

Native-state conformational dynamics of GART: a regulatory pH-dependent coil-helix transition examined by electrostatic calculations

Glycinamide ribonucleotide transformylase (GART) undergoes a pH-dependent coil-helix transition with pK(a) approximately 7. An alpha-helix is formed at high pH spanning 8 residues of a 21-residue-long loop, comprising the segment Thr120-His121-Arg122-Gln123-Ala124-Leu125-Glu126-Asn127. To understand the electrostatic nature of this loop-helix, called the activation loop-helix, which leads to the formation and stability of the alpha-helix, pK(a) values of all ionizable residues of GART have been calculated, using Poisson-Boltzmann electrostatic calculations and crystallographic data. Crystallographic structures of high and low pH E70A GART have been used in our analysis. Low pK(a) values of 5.3, 5.3, 3.9, 1.7, and 4.7 have been calculated for five functionally important histidines, His108, His119, His121, His132, and His137, respectively, using the high pH E70A GART structure. Ten theoretical single and double mutants of the high pH E70A structure have been constructed to identify pairwise interactions of ionizable residues, which have aided in elucidating the multiplicity of electrostatic interactions of the activation loop-helix, and the impact of the activation helix on the catalytic site. Based on our pK(a) calculations and structural data, we propose that: (1) His121 forms a molecular switch for the coil-helix transition of the activation helix, depending on its protonation state; (2) a strong electrostatic interaction between His132 and His121 is observed, which can be of stabilizing or destabilizing nature for the activation helix, depending on the relative orientation and protonation states of the rings of His121 and His132; (3) electrostatic interactions involving His119 and Arg122 play a role in the stability of the activation helix; and (4) the activation helix contains the helix-promoting sequence Arg122-Gln123-Ala124-Leu125-Glu126, but its alignment relative to the N and C termini of the helix is not optimal, and is possibly of a destabilizing nature. Finally, we provide electrostatic evidence that the formation and closure of the activation helix create a hydrophobic environment for catalytic-site residue His108, to facilitate catalysis.

pH dependence of stability of staphylococcal nuclease: evidence of substantial electrostatic interactions in the denatured state

The pH dependence of stability of staphylococcal nuclease was studied with two independent equilibrium thermodynamic approaches. First, by measurement of stability in the pH range 9 to 3.5 by fluorescence-monitored denaturation with urea (Delta), GdnHCl (Delta), and heat (Delta). Second, by numerical integration of H(+) titration curves (Delta) measured potentiometrically under native (100 mM KCl) and unfolding (6.0 M GdnHCl) conditions. The pH dependence of stability described by Delta, Delta, and Delta was comparable but significantly different from the one described by Delta. The decrease in Delta between pH 9 and pH 4 was 4 kcal/mol greater than the decrease in Delta, Delta, and Delta in the same pH range. In 6 M GdnHCl, all the ionizable groups titrated with the pK(a) values of model compounds. Therefore, Delta represents the free energy difference between the native state (N) and an ensemble of unstructured, or expanded, and highly screened conformations. In contrast, the shallower pH dependence of stability described by Delta and by Delta between pH 9 and 5 was consistent with the titration of histidines with depressed, nativelike pK(a) values in the denatured state (D). These depressed pK(a) values likely reflect long-range electrostatic interactions with the other 29 basic groups and are a consequence of the compact character of the D state. The steep change in Delta and Delta at pH < 5 suggests that near pH 5 the structural and thermodynamic character of the D state shifts toward a state in which acidic residues titrate with normal pK(a) values, presumably because the electrostatic interactions with basic residues are lost, maybe as a consequence of an expansion.

Realistic modeling of the denatured states of proteins allows accurate calculations of the pH dependence of protein stability

Computational techniques based on continuum electrostatics treatments have been successful in predicting and interpreting the pKa values of ionizable amino acids in folded proteins. Despite this progress, efforts to reproduce the pH-dependence of protein stability have met with only limited success: agreement with experimental results has been only qualitative. It has been argued previously that the most likely reason for discrepancies is the presence of residual electrostatic interactions in the unfolded state, which cause pKa values to be shifted from their model compound values. Here we show that by constructing atomistic models of the unfolded state with a simple molecular mechanics protocol that uses the native state as a starting point, much improved reproduction of pH effects on protein stability can be obtained. In contrast, when a fully extended model of the unfolded state is used, no such improvement is obtained, a result that suggests that local interactions with residues nearby in the sequence are not sufficient to properly account for the pKa shifts in the unfolded state. In comparison to model compound values, the pKa values of acidic residues in "native-like" unfolded states are typically found to be shifted downwards by approximately 0.3 pH unit, in good agreement with the average downward shift deduced from experimental measurements. Given its success in the present situation, the protocol employed here for developing simple models of the unfolded state may prove useful in other computer simulation applications.

Simulations of the structural and dynamical properties of denatured proteins: the "molten coil" state of bovine pancreatic trypsin inhibitor

The dynamic nature of denatured, unfolded proteins makes it difficult to characterize their structures experimentally. To complement experiment and to obtain more detailed information about the structure and dynamic behavior of the denatured state, we have performed eleven 2.5 ns molecular dynamics simulations of reduced bovine pancreatic trypsin inhibitor (BPTI) at high temperature in water and a control simulation at 298 K, for a total of 30 ns of simulation time. In a neutral pH environment (acidic residues ionized), the unfolded protein structures were compact with an average radius of gyration 9% greater than the native state. The compact conformations resulted from the transient formation of non-native hydrophobic clusters, turns and salt bridges. However, when the acidic residues were protonated, the protein periodically expanded to a radius of gyration of 18 to 20 A. The early steps in unfolding were similar in the different simulations until passing through the major transition state of unfolding. Afterwards, unfolding proceeded through one of two general pathways with respect to secondary structure: loss of the C-terminal helix followed by loss of beta-structure or the opposite. To determine whether the protein preferentially sampled particular conformational substates in the denatured state, pairwise Calpha root-mean-square deviations were measured between all structures, but similar structures were found between only two trajectories. Yet, similar composite properties (secondary structure content, side-chain and water contacts, solvent accessible surface area, etc.) were observed for the structures that unfolded through different pathways. Somewhat surprisingly, the unfolded structures are in agreement with both past experiments suggesting that reduced BPTI is a random coil and more recent experiments providing evidence for non-random structure, demonstrating how ensembles of fluctuating structures can give rise to experimental observables that are seemingly at odds.

pKA values of carboxyl groups in the native and denatured states of barnase: the pKA values of the denatured state are on average 0.4 units lower than those of model compounds

We have determined the pKA values of the 12 carboxyl residues in the native and denatured state of barnase by a combination of thermodynamic measurements on mutants of charged residues and NMR titration data. The pKA values of the 11 residues titrating under folding conditions (above pH 2.2) were determined by two-dimensional 1H NMR. The pKA value of the remaining residue, Asp 93 which forms a salt link with Arg 69 and titrates at much lower pH values, was determined by changes in the pH dependence of the stability of the protein upon mutation to Asn: pKAsp93A at low ionic strength (50 mM) and pKAsp93A at high ionic strength (600 mM). The overall titration of the native state is nonideal, and the protein retains fractionally ionized residues other than Asp 93 throughout the experimental pH range of 0.2-6.3. Protonation events taking place at pH values below 2 were further characterized by the pH dependence of the unfolding kinetics of wild-type and charge-mutant proteins. By comparing the observed pH dependence of the protein stability with that calculated from the pKA values for the native protein, we demonstrate that the pKA values of the denatured state are significantly lower than those reported for model compounds: the pKA values of the denatured state appear on average 0.4 units lower than previous estimates in the presence of chemical denaturant. The results have direct implications for calculations of the energetics of proton equilibria and suggest that the acid/thermally denatured state is not an extended coil where the residues are isolated from one another by the intervening solvent but is compact and involves intramolecular charge repulsion.

The nature of the initial step in the conformational folding of disulphide-intact ribonuclease A

Here we investigate conformational folding reaction of disulphide-intact ribonuclease A in the absence of the complicating effects due to non-native interactions (such as cis/trans proline isomerization) in the unfolded state. The conformational folding process is found to be intrinsically very fast occurring on the milliseconds time scale. The kinetic data indicate that the conformational folding of ribonuclease A proceeds through the formation of a hydrophobically collapsed intermediate with properties similar to those of equilibrium molten-globules. Furthermore, the data suggest that the rate-limiting transition states on the unfolding and refolding pathways are substantially different with the refolding transition state having non-native-like properties.

Molecular and biochemical analyses of combining sites of monoclonal anti-morphine antibodies

V region nucleotide sequences were determined by mRNA sequencing for 11 monoclonal anti-morphine antibodies with slightly different specificities for morphine-related opiates. The VH region nucleotide sequences of the antibodies MOR8, MOR33, MOR35, MOR44, MOR83, and MOR76 were classified into the VH-5 (7183) family, while the antibodies MOR39, MOR115, MOR131, MOR158 and MOR180 used VH-1 (J558) family genes. MOR39, MOR115 and MOR131 used the V lambda-1 gene for their L chain V region. MOR158 and MOR180 used the Vk-10 gene. MOR8, MOR33, OR35, MOR44, MOR76 and MOR83 used VK-21D. The antibody sets MOR158 and MOR180; MOR39 and MOR131; and MOR8, MOR33, MOR35, MOR44, MOR76 and MOR83 appeared to be somatic mutants derived from the same clones since they showed the same VH/VL usage and V(D)J recombination pattern. The pH-reactivity profiles for these antibodies revealed that the binding of morphine to the antibodies is highly dependent on the pH value of the assay solution, suggesting the importance of the electrostatic interaction between the positive charge of morphine and the negative charges at or near the combining sites. Direct UV-photoaffinity labeling with 3H-morphine was carried out in order to estimate the orientation of morphine in the combining sites. The H chains were preferentially labeled in MOR8, MOR33, MOR35, MOR76, and MOR83, whereas most of the crosslinked hapten was found in the L chains in MOR39, MOR115, MOR131, MOR158 and MOR180. Thus, these 11 antibodies were classified into two types in terms of reactivity in the photoaffinity labeling.

Biophysical models of protein denaturation. II. Effects of denaturants and of pH

In order to broaden the scope and increase the utility of differential scanning calorimetry, a theoretical model of calorimetric thermograms is presently proposed which facilitates their biophysical interpretation and accounts explicitly for their modifications induced by denaturing agents and/or pH. The model rests mainly on statistical-physical considerations, the denaturation-linked increase of the number of binding sites for denaturants (including H+) serving as the conceptual basis for thermogram modelling. Denaturants were envisioned as contributing indirectly to thermal denaturation by forming complexes preferentially with unfolded protein molecules, shifting thus the equilibrium towards the denatured phase. After postulating the probability of complex formation, mean numbers of the relevant molecular species were computed by ensemble averaging. Finally, an eight-parameter expression has been derived defining protein heat capacity as a function of both temperature and denaturant concentration (or pH), each of the eight parameters having a distinct biophysical meaning. The model has been tested by applying it to the prediction of the pH-dependence of thermograms. Four proteins have been considered (lysozyme, myoglobin, apomyoglobin, and ribonuclease A), each represented by a series of three to four published thermograms recorded under different pH conditions. Model equations, fitted simultaneously to all thermograms in a pH series, reproduced correctly experimental tracings. Parameter values obtained as best-fit requirements (particularly those representing the number of binding sites unmasked by denaturation and the free energy of ion binding) were in close agreement with empirical, mainly potentiometric, data from literature. The empirically established pH-independence of the total enthalpy of denaturation, the phenomenon of cold denaturation, the pH-dependence of the Gibbs free energy of denaturation, of the melting temperature and of the temperature of cold denaturation, were all correctly predicted by the model. Combined effects of multiple denaturants, including the effects of pH in the presence of denaturants other than protons, are also predictable by the model.

Stability of ribonuclease A in solution and the freeze-dried state

The solution stability of bovine pancreatic ribonuclease A (RNase) in phosphate buffer (pH 4.0, 6.4, and 10.0) at 45 degrees C decreased with increasing pH. Soluble aggregates were formed at each pH and corresponded qualitatively to the loss of enzymatic activity in the samples. Freeze drying of RNase resulted in no immediate loss of enzymatic activity in both the presence and absence of phosphate buffer salts. Freeze-dried RNase stored at 45 degrees C lost enzymatic activity and formed nondissociable aggregates at rates described by the following rank order of formulation contents: distilled water less than pH 6.4 phosphate buffer, less than pH 4.0 phosphate buffer less than pH 10.0 phosphate buffer. The amount of residual moisture remaining in the freeze-dried cakes was directly related to the rate of enzymatic activity loss and aggregate formation. The degradation rate was also directly related to the concentration of phosphate buffer salts added to the freeze-dried formulation.

Unfolding free energy changes determined by the linear extrapolation method. 2. Incorporation of delta G degrees N-U values in a thermodynamic cycle

The linear extrapolation method was used to evaluate the unfolding free energy changes (delta G degrees N-U) for phenylmethanesulfonyl chymotrypsin (PMS-Ct) at pH 6.0. The nonlinear least-squares fits of difference spectral data using urea and guanidinium chloride as denaturants gave identical values for delta G degrees N-U and delta epsilon degrees U, the latter being extinction coefficient differences between native and unfolded forms of the protein in the limit of zero concentration of denaturant. The independence of these parameters from the nature of solvent suggests strongly that they are characteristic properties of the protein alone. The delta G degrees N-U data at pH 6.0 and 4.0, which differ by more than 100-fold in stability of the protein, were incorporated into a thermodynamic cycle involving free energy changes for titration of native and unfolded PMS-Ct from pH 4.0 to 6.0. The purpose of the cycle was to test whether delta G degrees N-U obtained by use of the linear extrapolation method exhibits the characteristics required of a thermodynamic function of state. Within error, the thermodynamic cycle was found to accommodate the delta G degrees N-U quantities obtained at pH 4.0 and 6.0 for PMS-Ct.

Cloning, characterization and nucleotide sequences of two cDNAs encoding human pancreatic trypsinogens

Two cDNA clones encoding two major human trypsinogen isozymes were isolated from a human pancreatic cDNA library. The deduced amino acid (aa) sequences of the two trypsinogen precursors are found to have 89% sequence homology, and have the same number of aa (247), including 15 aa for a signal peptide and 8 aa for an activation peptide. Southern blot analysis of human genomic DNA using the cloned cDNA as a probe, revealed that the human trypsinogen genes constitute a multigene family of more than ten genes.

Unfolding pathway of myoglobin. Evidence for a multistate process

The free energy of unfolding of horse myoglobin has been calculated from the denaturation pattern induced by guanidine hydrochloride as well as by acid. The delta GH2O, i.e., the value in the absence of denaturant obtained by using the two-state transition model, was found to be 25% lower than that determined from the acid denaturation pattern, i.e., 12.0 kcal/mol, although the extent of protein denaturation produced by acid was much lower. The amount of helical structure surviving the acid-induced conformational change was estimated to be 50% of that present in the native protein, and it could be destroyed only after exposure of myoglobin samples kept at pH 3.0 to concentrated guanidine. From the guanidine denaturation pattern at acidic pH, a further variation of free energy of unfolding of 5.5 kcal/mol could be calculated, thus indicating that the overall free energy of unfolding determined from the two consecutive processes corresponds to 17.5 kcal/mol. The discrepancy between the two sets of data, i.e., guanidine unfolding at neutral pH and acid unfolding followed by addition of denaturant, has been considered to depend on the general assumption that the guanidine unfolding of myoglobin is a two-state process in the transition region. According to the recent experimental evidence showing the occurrence of at least two molecular events during the guanidine unfolding of apomyoglobin [Colonna, G., Balestrieri, C., Bismuto, E., Servillo, L., & Irace, G. (1982) Biochemistry 21, 212-215], the guanidine denaturation pattern of myoglobin was analyzed in terms of two independent steps.(ABSTRACT TRUNCATED AT 250 WORDS)

Effects of neutral salts on the circular dichroism spectra of ribonuclease A and ribonuclease S

The circular dichroism (CD) spectra of ribonuclease A, ribonuclease S, and N-acetyltyrosineamide were recorded as a function of pH in the presence of various concentrations of inorganic salts. Above pH 9.0 salting-in of tyrosine residues increases their intramolecular associations. This association enhances the contribution from these residues to the CD spectrum leading to an apparent titration curve that is shifted toward lower pH. The data indicate that unfolding of ribonuclease A and S by inorganic salts does not begin with disrupting existing electrostatic interactions. But, as the unfolding process progresses, disruption of electrostatic interactions may take place. This is consistent with our previous calorimetric studies which suggest that unfolding of ribonuclease A by salts proceeds initially by energetically favorable solvation of the folded protein. An increase in ellipiticity at 275 nm of partially unfolded protein in salt was observed as the pH was changed from 7.0 to 4.0. This observation may suggest that the isothermal unfolding of the protein by salts at low pH proceeds through an intermediate step which involves histidine residues and causes a conformational change in the tyrosine's asymmetric environment.

The pH dependence of the reversible unfolding of ovalbumin A1 by guanidine hydrochloride

The pH dependence of the reversible guanidine hydrochloride denaturation of the major fraction of ovalbumin (ovalbumin A1) was studied by a viscometric method in the pH range 1-7, at 25 degrees C and at six different denaturant concentrations (1.5-2.6 M). At any denaturant concentrationa reduction in pH favoured the transition from the native to the denatured state. The latter was essentially 'structureless', as revealed by the fact that the reduced viscosity of the acid and guanidine hydrochloride denatured state of ovalbumin A1 (obtained at different denaturant concentrations in acidic solutions) was measured (at a protein concentration of 3.8 mg/ml) to be 29.2 ml/g which is identical to that found in 6 M guanidine hydrochloride wherein the protein behaves as a cross-linked random coil. A quantitative analysis of the results on the pH dependence of the equilibrium constant for the denaturation process showed that on denaturation the intrinsic pK of two carboxyl groups in ovalbumin A1 went up from 3.1 in the native state to 4.4 in the denatured state of the protein.

Molecular weights and Stokes radii of soluble elastins

Soluble elastin was isolated from lathyritic chick aorta using neutral salt solutions in the presence of beta-amino propionitrile. The effect of a carboxy-methylation step in conjunction with proteolytic inhibitors was investigated. Hydrodynamic (Stokes) radii of soluble elastins were measured by gel filtration and the molecular size and weight distribution in purified fractions are reported.