Determination of dissociation constants of (4RS)-[4-carboxy-5,8,11-tris(carboxymethyl)-1-phenyl-2-oxa-5,8,11-triazatridecan-13-oic acid] in water and biological fluids by means of nuclear magnetic resonance spectroscopy and the DISCO software package

The pK(A) values of (4RS)-[4-carboxy-5,8,11-tris(carboxymethyl)-1-phenyl-2-oxa-5,8,11-triazatridecan-13-oic acid] (BOPTA), a polyprotic molecule whose gadolinium complex is an important magnetic resonance imaging contrast agent for clinical use, have been determined in water, in physiologic solution (PS), in serum (S), and in cerebrospinal fluid (CSF), by means of 13C nuclear magnetic resonance spectroscopy data processed by a dedicated software package called DISCO. The aim of this study was to supply the BOPTA pK(A) values in media very similar to the in vivo environment and, consequently, to get a picture of the in vivo behavior of its Gd complex, whose thermodynamic stability is directly linked to the pK(A) values. The pK(A) values appeared to be almost equal both in D(2)O and in PS, while pK(1) and pK(5) values in CSF differ a little. In S, only pK(2) and pK(3) were calculated due to the narrow pH range used for data collection. However, these pK(A) values were found equal to those in the other media. These results represent the first direct spectroscopic evidence of a substantial invariability of BOPTA behavior in different media and they justify the extrapolation to biological fluids of the data obtained in water. The values also confirmed the high-quality performance of DISCO in calculating pK(A) values of polyprotic molecules in complex media.     

Evidence of two-step deprotonation of D-mannitol in aqueous solution

Deprotonation of D-mannitol was studied in aqueous basic solutions by means of potentiometry and (13)C NMR spectroscopy. Two-step dissociation in the pH range from 12 to 13.8 was shown, and successive dissociation constants K(a1) and K(a2) were determined. In a solution with ionic strength I = 1.0 M (NaOH + NaNO(3)) pK(a1) = 13.1 +/- 0.1 and pK(a2) = 13.8 +/- 0.2. With increasing ionic strength from 0.75 to 3.0 M, both pK(a1) and pK(a2) values decrease. Deprotonation-induced chemical shifts in pH-variable (13)C NMR spectra show that the OH-groups next to internal carbon atoms C-3 and C-4 dissociate to a greater extent compared to OH-groups next to external carbon atoms C-1 and C-6.     

Disruption of an intermonomer salt bridge in the p53 tetramerization domain results in an increased propensity to form amyloid fibrils

We describe in molecular detail how disruption of an intermonomer salt bridge (Arg337-Asp352) leads to partial destabilization of the p53 tetramerization domain and a dramatically increased propensity to form amyloid fibrils. At pH 4.0 and 37 degrees C, a p53 tetramerization domain mutant (p53tet-R337H), associated with adrenocortical carcinoma in children, readily formed amyloid fibrils, while the wild-type (p53tet-wt) did not. We characterized these proteins by equilibrium denaturation, 13C(alpha) secondary chemical shifts, (1H)-15N heteronuclear NOEs, and H/D exchange. Although p53tet-R337H was thermodynamically less stable, NMR data indicated that the two proteins had similar secondary structure and molecular dynamics. NMR derived pK(a) values indicated that at low pH the R337H mutation partially disrupted an intermonomer salt bridge. Backbone H/D exchange results showed that for at least a small population of p53tet-R337H molecules disruption of this salt bridge resulted in partial destabilization of the protein. It is proposed that this decrease in p53tet-R337H stability resulted in an increased propensity to form amyloid fibrils.     

Complete (1)H and (13)C NMR chemical shift assignments of mono-, di-, and trisaccharides as basis for NMR chemical shift predictions of polysaccharides using the computer program casper

The computer program casper uses (1)H and (13)C NMR chemical shift data of mono- to trisaccharides for the prediction of chemical shifts of oligo- and polysaccharides. In order to improve the quality of these predictions the (1)H and (13)C, as well as (31)P when applicable, NMR chemical shifts of 30 mono-, di-, and trisaccharides were assigned. The reducing sugars gave two distinct sets of NMR resonances due to the α- and β-anomeric forms. In total 35 (1)H and (13)C NMR chemical shift data sets were obtained from the oligosaccharides. One- and two-dimensional NMR experiments were used for the chemical shift assignments and special techniques were employed in some cases such as 2D (1)H,(13)C-HSQC Hadamard Transform methodology which was acquired approximately 45 times faster than a regular t(1) incremented (1)H,(13)C-HSQC experiment and a 1D (1)H,(1)H-CSSF-TOCSY experiment which was able to distinguish spin-systems in which the target protons were only 3.3Hz apart. The (1)H NMR chemical shifts were subsequently refined using total line-shape analysis with the PERCH NMR software. The acquired NMR data were then utilized in the casper program (http://www.casper.organ.su.se/casper/) for NMR chemical shift predictions of the O-antigen polysaccharides from Klebsiella O5, Shigella flexneri serotype X, and Salmonella arizonae O62. The data were compared to experimental data of the polysaccharides from the two former strains and the lipopolysaccharide of the latter strain showing excellent agreement between predicted and experimental (1)H and (13)C NMR chemical shifts.     

Proton-nuclear-magnetic-resonance study of the conformation of neurotoxin II from Middle-Asian cobra (Naja naja oxiana) venom.

A proton nuclear magnetic resonance (NMR) study at 100 and 300 MHz of neurotoxin II from the venom of Middle-Asian cobra Naja naja oxiana has been performed in 2H2O and H2O solutions. By means of chemical modification and double resonance all the aromatic residue resonances have been assigned. From the NMR titration curves, pK values of histidine 4 and histidine 31 residues have been determined. For one of the two neighbouring tryptophan residues pH dependence (in the 2-8-pH range) of the chemical shifts of indole protons has been revealed. According to the different sensitivity of the linewidth of indole NH resonances to pH in H2O solution, the accessibility of each of the tryptophan residues has been estimated. Temperature dependence has been observed for the linewidth of the aromatic resonances of the tyrosine 24 residue. Deuterium exchange rates have been measured for amide protons as well as for C(2)H histidine resonances. The NMR data obtained have allowed the conclusions to be made that the two histidine residues and one of the tryptophan residues should be localized on the surface of the protein globule, that arginine residues should be present in the environment of histidine 4, that histidine 31 and the buried tryptophan are possibly localized in close spatial proximity and that the side chain of tyrosine 24 is buried within the protein globule.     

Remeasuring HEWL pK(a) values by NMR spectroscopy: methods, analysis, accuracy, and implications for theoretical pK(a) calculations

Site-specific pK(a) values measured by NMR spectroscopy provide essential information on protein electrostatics, the pH-dependence of protein structure, dynamics and function, and constitute an important benchmark for protein pK(a) calculation algorithms. Titration curves can be measured by tracking the NMR chemical shifts of several reporter nuclei versus sample pH. However, careful analysis of these curves is needed to extract residue-specific pK(a) values since pH-dependent chemical shift changes can arise from many sources, including through-bond inductive effects, through-space electric field effects, and conformational changes. We have re-measured titration curves for all carboxylates and His 15 in Hen Egg White Lysozyme (HEWL) by recording the pH-dependent chemical shifts of all backbone amide nitrogens and protons, Asp/Glu side chain protons and carboxyl carbons, and imidazole protonated carbons and protons in this protein. We extracted pK(a) values from the resulting titration curves using standard fitting methods, and compared these values to each other, and with those measured previously by 1H NMR (Bartik et al., Biophys J 1994;66:1180–1184). This analysis gives insights into the true accuracy associated with experimentally measured pK(a) values. We find that apparent pK(a) values frequently differ by 0.5–1.0 units depending upon the nuclei monitored, and that larger differences occasionally can be observed. The variation in measured pK(a) values, which reflects the difficulty in fitting and assigning pH-dependent chemical shifts to specific ionization equilibria, has significant implications for the experimental procedures used for measuring protein pK(a) values, for the benchmarking of protein pK(a) calculation algorithms, and for the understanding of protein electrostatics in general.     

Thermodynamic linkage between the binding of protons and inhibitors to HIV-1 protease

The aspartyl dyad of free HIV-1 protease has apparent pK(a)s of approximately 3 and approximately 6, but recent NMR studies indicate that the aspartyl dyad is fixed in the doubly protonated form over a wide pH range when cyclic urea inhibitors are bound, and in the monoprotonated form when the inhibitor KNI-272 is bound. We present computations and measurements related to these changes in protonation and to the thermodynamic linkage between protonation and inhibition. The Poisson-Boltzmann model of electrostatics is used to compute the apparent pK(a)s of the aspartyl dyad in the free enzyme and in complexes with four different inhibitors. The calculations are done with two parameter sets. One assigns epsilon = 4 to the solute interior and uses a detailed model of ionization; the other uses epsilon = 20 for the solute interior and a simplified representation of ionization. For the free enzyme, both parameter sets agree well with previously measured apparent pK(a)s of approximately 3 and approximately 6. However, the calculations with an internal dielectric constant of 4 reproduce the large pKa shifts upon binding of inhibitors, but the calculations with an internal dielectric constant of 20 do not. This observation has implications for the accurate calculation of pK(a)s in complex protein environments. Because binding of a cyclic urea inhibitor shifts the pK(a)s of the aspartyl dyad, changing the pH is expected to change its apparent binding affinity. However, we find experimentally that the affinity is independent of pH from 5.5 to 7.0. Possible explanations for this discrepancy are discussed.     

The structure and dipole moment of globular proteins in solution and crystalline states: use of NMR and X-ray databases for the numerical calculation of dipole moment

The large dipole moment of globular proteins has been well known because of the detailed studies using dielectric relaxation and electro-optical methods. The search for the origin of these dipolemoments, however, must be based on the detailed knowledge on protein structure with atomic resolutions. At present, we have two sources of information on the structure of protein molecules: (1) x-ray databases obtained in crystalline state; (2) NMR databases obtained in solution state. While x-ray databases consist of only one model, NMR databases, because of the fluctuation of the protein folding in solution, consist of a number of models, thus enabling the computation of dipole moment repeated for all these models. The aim of this work, using these databases, is the detailed investigation on the interdependence between the structure and dipole moment of protein molecules. The dipole moment of protein molecules has roughly two components: one dipole moment is due to surface charges and the other, core dipole moment, is due to polar groups such as N--H and C==O bonds. The computation of surface charge dipole moment consists of two steps: (A) calculation of the pK shifts of charged groups for electrostatic interactions and (B) calculation of the dipole moment using the pK corrected for electrostatic shifts. The dipole moments of several proteins were computed using both NMR and x-ray databases. The dipole moments of these two sets of calculations are, with a few exceptions, in good agreement with one another and also with measured dipole moments.     

Charge equilibria and pK(a) of malvidin-3-glucoside by electrophoresis

Paper electrophoresis has been used over the pH range 1.2 to 10.4 to measure apparent pK(a) values for malvidin-3-O-glucoside of pK(a(1)) 1.76+/-0.07, pK(a(2)) 5.36+/-0.04, and pK(a(3)) 8.39+/-0.07. Using solvent partitioning between buffered aqueous solutions and n-octanol, several micro-pK(a) constants for malvidin-3-O-glucoside were also identified, highlighting the complex nature of malvidin-3-glucoside equilibria. As a nonspectrophotometric procedure, the charge-dependent electrophoretic mobility method provided independent information on the net charge and color of anthocyanin species at wine pH (ca. 3.6). At this pH, the color of malvidin-3-glucoside in red wines is consistent only with the uncharged quinonoidal base as a major colored component of the equilibria.     

Synthesis and 31P NMR characterization of new low toxic highly sensitive pH probes designed for in vivo acidic pH studies

With the aim to provide sensitive 31P NMR probes of intra- and extracellular pH gradients that may reach cellular acidic compartments in biological systems, new alpha-aminophosphonates were designed to meet basic requirements such as a low pK(a)s and a great chemical difference (Deltadelta(ab)) between the limiting 31P NMR chemical shifts in acidic (delta(a)) and basic (delta(b)) media. A series of six phosphorylated pyrrolidines and linear aminophosphonates were synthesized using aminophosphorylation reactions and were screened for cytotoxicity on cultured Müller cells. Among the compounds not being toxic under these conditions, three molecules were selected since they displayed the best in vitro (in several phosphate buffers and in a cytosol-like solution) properties as 31P NMR acidic pH markers, that is 3, 5 and 9, having the pK(a) values of 3.63, 5.89 and 5.66, respectively. The Deltadelta(ab) values of these pH markers were at least 3 times larger than that of standard 31P NMR probes, with a low sensitivity to ionic strength changes. From these data, it was proposed that 3, 5 and 9 could be used as reporting probes of subtle proton movements in acidic compartments, an area that still remains poorly investigated using non invasive 31P NMR methods.     

13C-nuclear magnetic resonance studies of 85% 13C-enriched amino acids and small peptides. pH effects on the chemical shifts, coupling constants, kinetics of cis-trans isomerisation and conformation aspects.

The 13C chemical shifts of several 85% 13C-enriched amino acids and small peptides were studied as a function of pH. The results show that the chemical shifts of carbon atoms of ionizable groups vary significantly within the zone of their pK. Generally with the pH GOING FROM 7 to 1 all the deltaC are shifted more or less upfield with the exception of the carbonyl group carbon of the second last residue which is shifted slightly downfield. This suggests the formation of an hydrogen bond at acid pH involving in a seven-membered ring the C=O in question and the COOH terminal. The percentage of cis and trans conformers of glycyl-L-proline and glycyl-L-prolylglycine were studied as a function of pH. The trans form is always preponderant whatever the pH. The accessibility of the carbonyl group to protonation of the proline residue strongly influences the cis-trans equilibrium. Thus, with the pH varying from 7 to 1, the trans isomer changes from 61 to 85% for glycyl-L-proline and only from 77 to 80% for glycyl-L-prolylglycine. The proton NMR studies underline the important differences existing between the two molecular forms of glycyl-L-proline. The cis conformation is characterized with regard to the trans form by the non-equivalence of the alpha-protons of the glycine residue, by a lower pK(1) and by a larger deltadeltaHalpha of the proline residue as a function of pH. These results could suggest an end-to-end interaction in the cis form of the glycyl-L-proline molecule. The 13C-13C coupling constants were also studied as a function of pH. The results show that J(Co-Calpha) of a C-terminal residue, varying from 5 to 6 Hz and reflecting thhe pK of the carboxylate group, is a linear function of delta(Co) and delta(Calpha) as in the case of the amino acids. The total variation of the electron density of those two carbons in an amino acid is approximately 40% weaker than in a C-terminal residue. The charge distribution along the Calpha-C(o) bond, however, is practically the same in both cases. Finally the ratios of the conversion rate constants of the two isomers cis-trans of glycyl-proline were calculated at different pH values; the relations between the isomer percentages and delta(Co), delta(Calpha) on the one hand and the J(Co-Calpha) on the other were established.     

How do calcium ions induce free radical oxidation of hydroxy-1,4-naphthoquinone? Ca2+ stabilizes the naphthosemiquinone anion-radical of echinochrome A

A surprising effect is the direct action of Ca(2+) on redox reactions of ortho-quinoid compounds. The effect of Ca(2+) on oxidation of the sea urchin pigment 6-ethyl-2,3,5,7,8-pentahydroxy-1,4-naphthoquinone (echinochrome A) has been studied by electron paramagnetic resonance (EPR) spectroscopy, by UV/VIS absorbance spectroscopy, and by measurement of oxygen consumption. Echinochrome A per se reacted with dioxygen only in an alkaline solution; 2,3-semiquinone anion-radical of echinochrome A and superoxide anion-radical were the intermediates of the oxidation. Addition of calcium ions sharply increased the rate of echinochrome A autooxidation at alkaline pH and provoked oxidation at neutral pH. To explain this phenomenon we have focused on changes of the acid-base properties of echinochrome A in the presence of calcium and on stabilization of 2,3-semiquinone anion-radical of echinochrome A by Ca(2+). Dissociation constants (pK(a1), pK(a2), and pK(a3)) of echinochrome A determined by potentiometric titration were 5.20, 6.78, and >10 in calcium-free solution and 5.00, 6.10, and 7.15 in the presence of Ca(2+). We have found that Ca(2+) forms an insoluble adduct with the 2,3-semiquinone anion-radical. Thus, the effect of redox-inert calcium on the free radical reactions could be explained (i) by additional deprotonation of echinochrome A and (ii) by formation of a Ca(2+)-naphtho-2,3-semiquinone complex (calcium semiquinonate). Additionally, we have shown that the dried red spines from Strongylocentrotus intermedius possess paramagnetic properties; the EPR signal of the natural spines was similar to that of calcium semiquinonate obtained in our artificial chemical system.     

Evidence that 31P NMR is a sensitive indicator of small conformational changes in the coenzyme of aspartate aminotransferase

The pH dependence of 31P-NMR spectra of pig cytosolic aspartate aminotransferase, containing either N-(5'-phosphopyridoxyl)-L-aspartate or pyridoxal 5'-deoxymethylenephosphonate in place of the normal coenzyme pyridoxal 5'-phosphate, has been analysed. The chemical shifts of phosphopyridoxylaspartate and of pyridoxal 5'-deoxymethylenephosphonate model Schiff base in free solution show pK values of 6.3 and 7.4, attributable to the second deprotonation step of phosphate and phosphonate, respectively. However, these compounds behave very differently when bound to apoaspartate aminotransferase. 31P-NMR spectra of these enzyme derivatives indicate that the phosph(on)ate group remains dianionic throughout the pH range 4-8.5. A clear correlation between apparent pK values obtained from spectrophotometric titration of the coenzyme chromophore and those obtained by 31P NMR indicates that the same ionisation is being reported by both methods. The data are interpreted, on the basis of available crystallographic structures of chicken mitochondrial aspartate aminotransferase, to indicate that in each case the alteration in 31P chemical shift results from a conformational change in the coenzyme 5' side chain, in which one of the structures involves a near-eclipsed pair of bonds. Such a stressed conformation produces slight alterations in bond angles around the phosphorus atom, which in turn cause the observed change in 31P chemical shift. The evidence is taken to indicate that in this case 31P NMR is a sensitive reporter of stress in enzyme-bound pyridoxal 5'-phosphate and its derivatives.     

Nitrotyrosine chelation of nuclear magnetic resonance shift probes in proteins: application to bovine pancreatic trypsin inhibitor

The interactions of Pr(III) and Eu(III) with specifically nitrated derivatives of the basic bovine pancreatic trypsin inhibitor have been studied using optical spectroscopy and nuclear magnetic resonance (NMR) at 250 and 270 MHz. Stability constants for proton and metal binding to nitrotyrosines 10 and 21 determined optically are in good agreement with those from NMR. Observations of the Eu(III)-induced NMR shifts of the ring protons of nitrotyrosine 21 allowed calibration of the magnetic interactions for this binding site. The Pr(III)-induced shifts for several resolved nonexchangeable backbone proton resonances were compared with calculated shifts using the known x-ray structure. With several simplifying assumptions, the Pr(III)-induced shifts were used to assign one alpha-CH and five NH protons to compatible sets of backbone positions which are consistent with the known pH dependence and resistance to exchange with solvent D2O. Some of the more general aspects of lanthanide-induced shifts are discussed with reference to their use in proteins. Due to the complexities of the analysis of the shift data, the most straightforward use of this technique is in conjunction with the relaxation probe Gd(III) for measurement of intramolecular distances.     

Structural study of the interaction of vanadate with the ligand 1,2-dimethyl-3-hydroxy-4-pyridinone (Hdmpp) in aqueous solution

The interaction of vanadate with the ligand 1,2-dimethyl-3-hydroxy-4-pyridinone (Hdmpp) was studied in aqueous solution using a combination of multinuclear NMR and EPR spectroscopies, as well as potentiometry and cyclic voltammetry. The different species in solution were identified and characterized, and their pKa values and stability constants determined. The vanadium complexes formed in solution are strongly dependent on media composition (ionic strength, presence of buffer), pH and metal-to-ligand ratio (M:L). Two major species--V(V)/dmpp and V(V)/(dmpp)2--are formed in a 140 mM NaCl solution within the pH range 4.5 to 9.0, when M:L = 1:2. In the presence of excess ligand (M:L < or = 1:5), only the 1:2 complex is present, and at pH < 4 paramagnetic species are detected by EPR in solution, thus indicating a reducing capacity of the ligand. Cyclic voltammetry shows that redox processes in solution are not just electron transfer, but are accompanied by chemical reactions. The pK, values and stability constants were determined both by 51V NMR spectroscopy and potentiometry. The present results have a particular interest in the understanding of the aqueous solution chemistry in aerobic conditions of bis(1,2-dimethyl-3-hydroxy-4-pyridinonato) oxovanadium(IV) complex, VO(dmpp)2, a vanadium compound with potential insulin-mimetic properties.     

Tautomerism,acid-base equilibria,and H-bonding of the six histidines in subtilisin BPN' by NMR

We have determined by (15)N, (1)H, and (13)C NMR, the chemical behavior of the six histidines in subtilisin BPN' and their PMSF and peptide boronic acid complexes in aqueous solution as a function of pH in the range of from 5 to 11, and have assigned every (15)N, (1)H, C(epsilon 1), and C(delta2) resonance of all His side chains in resting enzyme. Four of the six histidine residues (17, 39, 67, and 226) are neutrally charged and do not titrate. One histidine (238), located on the protein surface, titrates with pK(a) = 7.30 +/- 0.03 at 25 degrees C, having rapid proton exchange, but restricted mobility. The active site histidine (64) in mutant N155A titrates with a pK(a) value of 7.9 +/- 0.3 and sluggish proton exchange behavior, as shown by two-site exchange computer lineshape simulation. His 64 in resting enzyme contains an extremely high C(epsilon 1)-H proton chemical shift of 9.30 parts per million (ppm) owing to a conserved C(epsilon 1)-H(.)O=C H-bond from the active site imidazole to a backbone carbonyl group, which is found in all known serine proteases representing all four superfamilies. Only His 226, and His 64 at high pH, exist as the rare N(delta1)-H tautomer, exhibiting (13)C(delta1) chemical shifts approximately 9 ppm higher than those for N(epsilon 2)-H tautomers. His 64 in the PMSF complex, unlike that in the resting enzyme, is highly mobile in its low pH form, as shown by (15)N-(1)H NOE effects, and titrates with rapid proton exchange kinetics linked to a pK(a) value of 7.47 +/- 0.02.     

NMR studies of the conformation of the natural sweetener rebaudioside A

Rebaudioside A is a natural sweetener from Stevia rebaudiana in which four β-d-glucopyranose units are attached to the aglycone steviol. Its ¹H and ¹³C NMR spectra in pyridine-d₅ were assigned using 1D and 2D methods. Constrained molecular dynamics of solvated rebaudioside using NMR constraints derived from ROESY cross peaks yielded the orientation of the β-d-glucopyranose units. Hydrogen bonding was examined using the temperature coefficients of the hydroxyl chemical shifts, ROESY and long-range COSY spectra, and proton-proton coupling constants.     

Titration of E. coli transhydrogenase domain III with bound NADP+ or NADPH studied by NMR reveals no pH-dependent conformational change in the physiological pH range

A pH-titration 2D NMR study of Escherichia coli transhydrogenase domain III with bound NADP(+) or NADPH has been carried out, in which the pH was varied between 5.4 and 12. In this analysis, individual amide protons served as reporter groups. The apparent pK(a) values of the amide protons, determined from the pH-dependent chemical shift changes, were attributed to actual pK(a) values for several titrating residues in the protein. The essential Asp392 is shown to be protonated at neutral pH in both the NADP(+) and NADPH forms of domain III, but with a marked difference in pK(a) not only attributable to the charge difference between the substrates. Titrating residues found in loop D/alpha5 point to a conformational difference of these structural elements that is redox-dependent, but not pH dependent. The observed apparent pK(a) values of these residues are discussed in relation to the crystal structure of Rhodospirillum rubrum domain III, the solution structure of E. coli domain III and the mechanism of intact proton-translocating transhydrogenase.     

Electrostatic interactions in ubiquitin: stabilization of carboxylates by lysine amino groups

To explore electrostatic interactions in ubiquitin, pK(a) values have been determined by NMR for all 12 carboxyl groups in wild-type ubiquitin and in variants where single lysines have been replaced by neutral residues. Aspartate pK(a) values in ubiquitin range from 3.1 to 3.8 and are generally less than model compound values. Most aspartate pK(a) values are within 0.2 pH unit of those predicted with a simple Tanford-Kirkwood model. Glutamate pK(a) values range from 3.8 to 4.5, close to model compound values and differing by 0.1-0.8 pH unit from calculated values. To determine the role of positive charges in modulating carboxyl pK(a) values, we mutated lysines at positions 11, 29, and 33 to glutamine and threonine. NMR studies with these six single-site mutants reveal significant interactions of Lys 11 and Lys 29 with Glu 34 and Asp 21, respectively: pK(a) values for Glu 34 and Asp 21 increase by approximately 0.5-0.8 pH unit, similar to predicted values, when the lysines are replaced by neutral residues. In contrast, the predicted interaction between Lys 33 and Glu 34 is not observed experimentally. In some instances, substitution of lysine by glutamine and threonine did not lead to the same changes in carboxyl pK(a) values. These may reflect new short-range interactions between the mutated residues and the carboxyl groups. Carboxyl pK(a) shifts > 0.5 pH unit result from mutations at groups that are <5 A from the carboxyl group. No interactions are observed at >10 A.