A Fourier transform NMR study of the thermal denaturation of ribonuclease A at low pH

Combining conformational flexibility and continuum electrostatics for calculating pK(a)s in proteins

Protein stability and function relies on residues being in their appropriate ionization states at physiological pH. In situ residue pK(a)s also provides a sensitive measure of the local protein environment. Multiconformation continuum electrostatics (MCCE) combines continuum electrostatics and molecular mechanics force fields in Monte Carlo sampling to simultaneously calculate side chain ionization and conformation. The response of protein to charges is incorporated both in the protein dielectric constant (epsilon(prot)) of four and by explicit conformational changes. The pK(a) of 166 residues in 12 proteins was determined. The root mean square error is 0.83 pH units, and >90% have errors of <1 pH units whereas only 3% have errors >2 pH units. Similar results are found with crystal and solution structures, showing that the method's explicit conformational sampling reduces sensitivity to the initial structure. The outcome also changes little with protein dielectric constant (epsilon(prot) 4-20). Multiconformation continuum electrostatics titrations show coupling of conformational flexibility and changes in ionization state. Examples are provided where ionizable side chain position (protein G), Asn orientation (lysozyme), His tautomer distribution (RNase A), and phosphate ion binding (RNase A and H) change with pH. Disallowing these motions changes the calculated pK(a).

Determination of the solvent accessibility of specific aromatic residues in gamma-crystallin by photo-CIDNP NMR measurements

The surface or solvent accessibility of certain individual aromatic residues of calf-gamma II crystallin in solution (1 mM) were measured by the dramatic intensity enhancements of NMR lines generated by the interactions of cyclic radical pair formation of the 3-N-carboxymethyl lumiflavin (flavin I) dye excited (488nm) by an argon laser (5 watts) with the protein. This effect is called photo-chemically induced dynamic nuclear polarization: photo-CIDNP. The "light" and "dark" NMR spectra were taken in alternating scans in the pulsed Fourier transform mode on a Bruker 360 MHz instrument. Subtraction results in the photo-CIDNP difference spectrum containing lines of the polarized residues. With flavin dyes only tyrosine, histidine, and tryptophan can be polarized. The respective theoretical static accessibility of these residues based upon van der Waal's contact radii have been calculated from the atomic coordinates and provide a basis for evaluating the dynamic NMR photo-CIDNP results and for assigning the resonances. These results suggest that while the four tryptophan residues are completely buried, His-113 and His-14 of the five histidines; and Tyr-165 and Tyr-62 of the fifteen tyrosines are sufficiently exposed to elicit a photo-CIDNP effect. These results confirm and extend the observations previously obtained with theoretical electrostatic programs and FT-NMR measurements.

Exchange behavior of the H-bonded amide protons in the 3 to 13 helix of ribonuclease S

The preceding article shows that there are eight highly protected amide protons in the S-peptide moiety of RNAase S at pH 5, 0 degrees C. The residues with protected NH protons are 7 to 13, whose amide protons are H-bonded in the 3 to 13 alpha-helix, and Asp 14, whose NH proton is H-bonded to the CO group of Val47. We describe here the exchange behavior of these eight protected protons as a function of pH. Exchange rates of the individual NH protons are measured by 1H nuclear magnetic resonance in D2O. A procedure is used for specifically labeling with 1H only these eight NH protons. The resonance assignments of the eight protons are made chiefly by partial exchange, through correlating the resonance intensities in spectra taken when the peptide is bound and when it is dissociated from S-protein in 3.5 M-urea-d4, in D2O, pH 2.3, -4 degrees C. The two remaining assignments are made and some other assignments are checked by measurements of the nuclear Overhauser effect between adjacent NH protons of the alpha-helix. There is a transition in exchange behavior between pH 3, where the helix is weakly protected against exchange, and pH 5 where the helix is much more stable. At pH 3.1, 20 degrees C, exchange rates are uniform within the helix within a factor of two, after correction for different intrinsic exchange rates. The degree of protection within the helix is only 10 to 20-fold at this pH. At pH 5.1, 20 degrees C, the helix is more stable by two orders of magnitude and exchange occurs preferentially from the N-terminal end. At both pH values the NH proton of Asp 14, which is just outside the helix, is less protected by an order of magnitude than the adjacent NH protons inside the helix. Opening of the helix can be observed below pH 3.7 by changes in chemical shifts of the NH protons in the helix. At pH 2.4 the changes are 25% of those expected for complete opening. Helix opening is a fast reaction on the n.m.r. time scale (tau much less than 1 ms) unlike the generalized unfolding of RNAase S which is a slow reaction. Dissociation of S-peptide from S-protein in native RNAase S at pH 3.0 also is a slow reaction. Opening of the helix below pH 3.7 is a two-state reaction, as judged by comparing chemical shifts with exchange rates. The exchange rates at pH 3.1 are predicted correctly from the changes in chemical shift by assuming that helix opening is a two-state reaction. At pH values above 3.7, the nature of the helix opening reaction changes. These results indicate that at least one partially unfolded state of RNAase S is populated in the low pH unfolding transition.

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.

The thermal unfolding of ribonuclease A. A 13C NMR study

The 13C nuclear magnetic resonance (NMR) spectra of ribonuclease A over the pH range 1-7 and between 6 and 70 degrees C reveal many of the details of its reversible unfolding. Although the unfolding may loosely be described as 'two-state', evidence is presented for intermediate unfolding stages at least 10 degrees C on either side of the main unfolding transition, particularly at low pH. The first residues to unfold are 17-24, in agreement with other results. The C-terminal region shows a steeper temperature dependence of its unfolding than does the main transition, which itself is shown to lead at all pH values to a semi-structured but internally flexible state which is far from being truly random-coil. This is confirmed by measurements of T1 and of nuclear Overhauser enhancement. Indeed, even at pH 1.1 and 70 degrees C there is evidence for considerable motional restriction of cysteine and proline residues, amongst others. The native protein has more variability of structure at low pH than at neutral pH, and also interchanges more rapidly with the semi-structured, denatured state.