The pKa Values of Acidic and Basic Residues Buried at the Same Internal Location in a Protein Are Governed by Different Factors

The pK_a values of internal ionizable groups are usually very different from the normal pK_a values of ionizable groups in water. To examine the molecular determinants of pK_a values of internal groups, we compared the properties of Lys, Asp, and Glu at internal position 38 in staphylococcal nuclease. Lys38 titrates with a normal or elevated pK_a, whereas Asp38 and Glu38 titrate with elevated pK_a values of 7.0 and 7.2, respectively. In the structure of the L38K variant, the buried amino group of the Lys38 side chain makes an ion pair with Glu122, whereas in the structure of the L38E variant, the buried carboxyl group of Glu38 interacts with two backbone amides and has several nearby carboxyl oxygen atoms. Previously, we showed that the pK_a of Lys38 is normal owing to structural reorganization and water penetration concomitant with ionization of the Lys side chain. In contrast, the pK_a values of Asp38 and Glu38 are perturbed significantly owing to an imbalance between favorable polar interactions and unfavorable contributions from dehydration and from Coulomb interactions with surface carboxylic groups. Their ionization is also coupled to subtle structural reorganization. These results illustrate the complex interplay between local polarity, Coulomb interactions, and structural reorganization as determinants of pK_a values of internal groups in proteins. This study suggests that improvements to computational methods for pK_a calculations will require explicit treatment of the conformational reorganization that can occur when internal groups ionize.     

Protonation of kanamycin A: Detailing of thermodynamics and protonation sites assignment

Protonation of an aminoglycoside antibiotic kanamycin A sulfate was studied by potentiometric titrations at variable ionic strength, sulfate concentration and temperature. From these results the association constants of differently protonated forms of kanamycin A with sulfate and enthalpy changes for protonation of each amino group were determined. The protonation of all amino groups of kanamycin A is exothermic, but the protonation enthalpy does not correlate with basicity as in a case of simple polyamines. The sites of stepwise protonation of kanamycin A have been assigned by analysis of ¹H-¹³C-HSQC spectra at variable pH in D₂O. Plots of chemical shifts for each H and C atom of kanamycin A vs. pH were fitted to the theoretical equation relating them to pK_a values of ionogenic groups and it was observed that changes in chemical shifts of all atoms in ring C were controlled by ionization of a single amino group with pK_a 7.98, in ring B by ionization of two amino groups with pK_a 6.61 and 8.54, but in ring A all atoms felt ionization of one group with pK_a 9.19 and some atoms felt ionization of a second group with pK_a 6.51, which therefore should belong to amino group at C3 in ring B positioned closer to the ring A while higher pK_a 8.54 can be assigned to the group at C1. This resolves the previously existed uncertainty in assignment of protonation sites in rings B and C.     

Structural origins of high apparent dielectric constants experienced by ionizable groups in the hydrophobic core of a protein

The side chains of Lys66, Asp66, and Glu66 in staphylococcal nuclease are fully buried and surrounded mainly by hydrophobic matter, except for internal water molecules associated with carboxylic oxygen atoms. These ionizable side chains titrate with pK_a values of 5.7, 8.8, and 8.9, respectively. To reproduce these pK_a values with continuum electrostatics calculations, we treated the protein with high dielectric constants. We have examined the structural origins of these high apparent dielectric constants by using NMR spectroscopy to characterize the structural response to the ionization of these internal side chains. Substitution of Val66 with Lys66 and Asp66 led to increased conformational fluctuations of the microenvironments surrounding these groups, even under pH conditions where Lys66 and Asp66 are neutral. When Lys66, Asp66, and Glu66 are charged, the proteins remain almost fully folded, but resonances for a few backbone amides adjacent to the internal ionizable residues are broadened. This suggests that the ionization of the internal groups promotes a local increase in dynamics on the intermediate timescale, consistent with either partial unfolding or increased backbone fluctuations of helix 1 near residue 66, or, less likely, with increased fluctuations of the charged side chains at position 66. These experiments confirm that the high apparent dielectric constants reported by internal Lys66, Asp66, and Glu66 reflect localized changes in conformational fluctuations without incurring detectable global structural reorganization. To improve structure-based pK_a calculations in proteins, we will need to learn how to treat this coupling between ionization of internal groups and local changes in conformational fluctuations explicitly.     

Binding of the bovine pancreatic secretory trypsin inhibitor (Kazal) to bovine serine (pro)enzymes

The effect of temperature and pH on the association equilibrium constant (K_a) for the binding of the bovine pancreatic secretory trypsin inhibitor (bovine PSTI, type I; Kazal inhibitor) to bovine β-trypsin, bovine α-chymotrypsin and bovine trypsinogen has been investigated. The results suggest that serine (pro)enzyme inhibitor interaction involves both rigorous spatial configuration and molecular flexibility.     

Fluorescence polarization studies of the interactions between protamine and enzymes in aqueous solutions

Fluorescence polarization has been used to study the interaction of dansylated protamine with the enzymes: pepsin, α-chymotrypsin, alkaline phosphatase and invertase. These interactions have been compared with those between dansylated protamine and polyacrylate, or polyvinylsulphate. Each of the various complexes was found to be dissociated by the addition of sodium nitrate and a critical electrolyte concentration (CEC) was determined for each system, to allow assessment of the relative order of binding to the dansylated protamine. This order was: polyvinylsulphate >pepsin >polyacrylate >alkaline phosphatase >α-chymotrypsin. The strength of binding was also assessed by determination of a binding constant, K_a. The values of K_a showed the same relative order of binding as the CEC values. Invertase behaved similarly to the other enzymes, but it was not possible to obtain an unambiguous assessment of the comparative strength of binding. In each case, the stoichiometry of the complex was also determined.     

Introduction of non-natural amino acid residues into the substrate-specific P1 position of trypsin inhibitor SFTI-1 yields potent chymotrypsin and cathepsin G inhibitors

A series of trypsin inhibitor SFTI-1compounds modified in substrate-specific P₁ position was synthesized by the solid-phase method. Lys5 present in the wild inhibitor was replaced by Phe derivatives substituted in para position of the phenyl ring, l-pyridylalanine and N-4-nitrobenzylgycine. Their inhibitory activities with bovine α-chymotrypsin and cathepsin G were estimated by determination of association equilibrium constants (K_a). All analogues inhibited bovine α-chymotrypsin. The highest inihbitory activity displayed peptides with the fluorine, nitro and methyl substituents. They were 13–15-fold more active than [Phe⁵]SFTI-1 used as a reference. They are the most potent chymotrypsin inhibitors of this size. Substitution of Lys5 by Phe did not change the cathepsin G inhibitory activity. Introduction of Phe(p-F), Phe(p-NH₂) and Phe(p-CH₃) in this position retained the affinity towards this proteinase, whereas Phe(p-guanidine) gave an inhibitor more than twice as active, which appeared to be stable in human serum. On the other hand, a peptomeric analogue with N-4-nitrobenzylglycine failed to inhibit cathepsin G. Despite the fact the introduced amino acids were non-coded, the peptide bonds formed by them were hydrolyzed by chymotrypsin. We postulate that additional interaction of para-substitutents with the enzyme are responsible for the enhanced inhibitory activity of the analogues.     

Binding changes and apparent pK a shifts of bromthymol blue as tools for mitochondrial reactions

The passive interaction of bromthymol blue and other anionic and neutral dyes with mitochondria and particles is accompanied by an increase of the apparent pK_a, is enhanced at higher ionic strengths, is slightly inhibited by neutral dyes, is high only for the acidic, un-ionized form, and is independent of the presence of the sulfonic side chain. The pH dependence for the binding of the dyes follows the ionization of the chromophoric group: the pK of binding, defined as the pH for 50% binding, is identical with the apparent pK_a of the dye.     

Kinetics of interaction of some α- and β-D-monosaccharides with concanavalin A

The rates of formation and dissociation of concanavalin A with some 4-methylumbelliferyl and p-nitrophenyl derivatives of α- and β-D-mannopyranosides and glucopyranosides were measured by fluorescence and spectral stopped-flow methods. All process examined were uniphasic. The second-order formation rate constants varied only from 6.8 · 10⁴ to 12.8 · 10⁴ M^?. s^?1, whereas the first-order dissociation rate constants ranged from 4.1. to 220 s^?1, all at ph 5.0, I = 0.3 M, and 25°C. Dissociation rates thus controlled the value of binding constant. The effect of temperature on these reactions was examined, from which enthalpies and entropies of activation and of reaction could be calculated. The effects of pH at 25°C on the reaction rates of 4-methylumbelliferyl α-D-mannopyranoside and 4-methylumbelliferyl α-D-glucopyranoside with concanavalin A were examined. The value of the binding constant K_ap (derived from the kinetics) at any pH could be related to the intrinsic binding constant K by the expression K_ap = K_aK(K_a + [H⁺])^?1. The values of K_a, the ionization constant of the protein segment responsive to sugar binding, were 3 · 10^?4 M and 1 · 10^?4 M for 4-methylumbelliferyl α-D-mannopyranoside and 4-methylumbelliferyl α-D-glucopyranoside, respectively. The binding constant of p-nitrophenyl α-D-mannopyranoside is surprisingly much less sensitive to a pH change from 5.0 to 2.7. Ionic strength had little effect on the binding characteristics of 4-methylumbelliferyl α-D-mannopyranoside to concanavalin A at pH 5.2 and 25°C.     

Selective oxidation of Met-192 in bovine α-chymotrypsin. Effect on catalytic and inhibitor binding properties

Catalytic and inhibitor binding properties of bovine α-chymotrypsin, in which the Met-192 residue has been converted by treatment with chloramine T to the sulfoxide derivative (Met(O)192 α-chymotrypsin), have been examined relative to the native enzyme (α-chymotrypsin), between pH 4.5 and 8.0 (μ = 0.1), and/or 5.0°C and 40.0°C. Values of k_cat, k₊₂ and/or k₊₃ for the hydrolysis of all the substrates examined (i.e., tMetAcONp, ZAlaONp, ZLeuONp, ZLysONp and ZTyrONp) catalyzed by native and Met(O)192 α-chymotrypsin are similar, as well as values of K_m for the hydrolysis of ZLeuONp, ZLysONp and ZTyrONp. On the other hand, K_s and K_m values for the hydrolysis of ZAlaONp and tMetAcONp are decreased by about 5-fold. Met-192 oxidation does not affect the kinetic and thermodynamic parameters for the (de)stabilization of the complex formed between the proteinase and the bovine basic pancreatic trypsin inhibitor. On the other hand, the recognition process between between α-chymotrypsin and the recombinant proteinase inhibitor eglin c from the leech Hirudo medicinalis is influenced by the oxidation event. Considering known molecular models, the observed catalytic and inhibitor binding properties of native and Met(O)192 α-chymotrypsin were related to the inferred stereochemistry of the proteinase-substrate and proteinase-inhibitor contact region(s).     

Energy cost for the proper ionization of active site residues in 6-phosphogluconate dehydrogenase from T. brucei

The catalytic mechanism of 6-phosphogluconate dehydrogenase requires the inversion of a Lys/Glu couple from its natural ionization state. The pK_a of these residues in free and substrate bound enzymes has been determined measuring by ITC the proton release/uptake induced by substrate binding at different pH values. Wt 6-phosphogluconate dehydrogenase from Trypanosoma brucei and two active site enzyme mutants, K185H and E192Q were investigated. Substrate binding was accompanied by proton release and was dependent on the ionization of a group with pK_a 7.07 which was absent in the E192Q mutant. Kinetic data highlighted two pK_a, 7.17 and 9.64, in the enzyme–substrate complex, the latter being absent in the E192Q mutant, suggesting that the substrate binding shifts Glu192 pK_a from 7.07 to 9.64. A comparison of wt and E192Q mutant appears to show that the substrate binding shifts Lys185 pK_a from 9.9 to 7.17. By comparing differences in proton release and the binding enthalpy of wt and mutant enzymes, the enthalpic cost of the change in the protonation state of Lys185 and Glu192 was estimated at ≈ 6.1 kcal/mol. The change in protonation state of Lys185 and Glu192 has little effect on Gibbs free energy, 240–325 cal/mol. However proton balance evidences the dissociation of other group(s) that can be collectively described by a single pK_a shift from 9.1 to 7.54. This further change in ionization state of the enzyme causes an increase of free energy with a total cost of 1.2–2.3 kcal/mol to set the enzyme into a catalytically competent form.     

Molecular Asymmetry in Pigeonpea Urease: pH Inactivation Studies

Pigeonpea (Cajanus cajan) urease was inactivated by incubating it in buffer of low pHs i.e., 4.8 and 4.5. The pattern of inactivation at both pHs was found to be biphasic, in which half of the activity was destroyed more rapidly than the remaining half. This distribution of active site into two categories is suggestive of site-site heterogeneity, or more specifically, the half-site reactivity of the enzyme moiety. Our pH studies on the rate of reaction showed the presence of two ionizable groups of pK _a values 6.2 ± 0.1 and 8.8 ± 0.1, respectively (Srivastava PK & Kayastha AM, J Mol Catal B: Enz, 16 (2001) 81-89). The later group corresponds to the pK _a value of cysteine group. Here we correlate the loss of urease activity by low pH treatment is due to the effect on essential thiol residues.     

Mechanism of action of α-galactosidase

The effect ofpH on K_m and V_max values of coconut α-galactosidase indicates the involvement of two ionizing groups with pK_a values of 3.5 and 6.5 in catalysis. Chemical modification has indicated the presence of two carboxyl groups, a tryptophan and a tyrosine, at or near the active site of α-galactosidase. Based on these facts a new mechanism of action for α-galactosidase is proposed in which the ionizing group with a pK_a of 3.5 is a carboxyl group involved in stabilizing a carbonium ion intermediate and the ionizing group with a pK_a of 6.5 is a carboxyl group perturbed due to the presence of a hydrophobic residues in its vicinity which donates a H⁺ ion in catalysis.     

Proton magnetic resonance studies of cytochrome c peroxidase: pD dependence of the isotropically shifted resonances

Proton nuclear magnetic resonance studies of the isotropically shifted resonances of native cytochrome c peroxidase have been carried out at 360 MHz. Proton resonances extend to 84 ppm downfield and 12 ppm upfleld from 2,2-dimethyl-2-silapentane-5-sulfonate and are characteristic of high-spin iron +3 heme proteins. Between pH 8 and 4.1 the isotropic resonances exhibit dramatic pH-dependent behavior which demonstrates the presence of two acid-base interconversions. One process, with a pK_a of 5.8, is slow on the NMR time scale and probably represents a protein conformation change. This process correlates with an apparent pK_a observed in the kinetic properties of the enzymes, with the alkaline form being the enzymatically active species. A second ionization with a pK of 4.9 is fast on the NMR time scale and probably represents a true ionization.     

Electrostatic analysis of the hepatitis C virus NS3 helicase reveals both active and allosteric site locations

Multi-conformation continuum electrostatics (MCCE) was used to analyze various structures of the NS3 RNA helicase from the hepatitis C virus in order to determine the ionization state of amino acid side chains and their pK_as. In MCCE analyses of HCV helicase structures that lacked ligands, several active site residues were identified to have perturbed pK_as in both the nucleic acid binding site and in the distant ATP-binding site, which regulates helicase movement. In all HCV helicase structures, Glu493 was unusually basic and His369 was abnormally acidic. Both these residues are part of the HCV helicase nucleic acid binding site, and their roles were analyzed by examining the pH profiles of site-directed mutants. Data support the accuracy of MCCE predicted pK_a values, and reveal that Glu493 is critical for low pH enzyme activation. Several key residues, which were previously shown to be involved in helicase-catalyzed ATP hydrolysis, were also identified to have perturbed pK_as including Lys210 in the Walker-A motif and the DExD/H-box motif residues Asp290 and His293. When DNA was present in the structure, the calculated pK_as shifted for both Lys210 and Asp290, demonstrating how DNA binding might lead to electrostatic changes that stimulate ATP hydrolysis.     

Chloroplast thylakoid membrane proteins having buried amine buffering groups

Some chloroplast thylakoid membrane proteins have anomalously low pK_a (near 7.8) amine groups, indicating that the buffering groups may be buried in hydrophobic regions and/or close to other positive charges. Other work has shown that the low pK_a amine group array is not in ready equilibrium with either the inner or outer bulk aqueous phases (Laszlo, J.A., Baker, G.M. and Dilley, R.A. (1984) J. Bioenerg. Biomembranes, 16, 37–51). Acetic anhydride reacts with the neutral amine and has been used as a probe for labeling the low pK_a amines. The buried array of buffering groups can be labeled with [³H]acetic anhydride in the dark only after the membranes were made leaky to protons with uncoupler addition. Sodium dodecyl sulfate/urea-polyacrylamide gel electrophoresis was used to separate the polypeptides and identify those that show the uncoupler-dependent labeling increase. Included in that group are polypeptides known to be associated with Photosystem II having M_r 17000, 22000 and 31000, some of the light-harvesting pigment proteins with M_r 24000–28000, the CF₀ component with M_r 8000 and some polypeptides associated with Photosystem I. A protein with M_r 15000 showed very large changes in labeling, but the identity of this polypeptide is unknown. The arrays of buried amine buffering groups are diversely distributed among membrane proteins and it is not clear what role, if any, they play in membrane function.     

Carbon-13 nmr spectra of vitamin B6 compounds: Aqueous solution equilibria and determination of microscopic acid dissociation constants

The aqueous solution equilibria and solute structure of vitamin B₆ compounds and several model compounds have been investigated using ¹³C-nmr spectroscopy. The unsubstituted α-carbon of these compounds is a very good probe for data which permits assignment of the ionization steps to indicidual groups. While the ionizations of the pyridinium and phenolic groups take place simultaneously in 3-hydroxypyridine, they take place in well-separated steps in pyridoxamine (PM), pyridoxamine phosphate (PMP), and pyridoxal phosphate. It has been established that the ionization with a pK_a value of 3.7 is predominantly phenolic in origin in PM and PMP. A zwitterionic structure consistent with the earlier spectroscopic investigations is proposed for the vitamin B₆ compounds in neutral aqueous solution.     

Specific dissociation of αB subunits from α-crystallin

Exposure of bovine α-crystallin to 0.1 M glycine at pH 7 decreases the average molar mass of the protein from 700 to 420 kDa. When the pH is lowered to 2.5, in the same buffer, the αB chains specifically dissociate from the aggregates, leaving a particle of 290 kDa containing only αA chains. The decrease in the molar mass corresponds to the mass of the αB chains in the original aggregate. The pH-dependent dissociation is fully reversible. Similar changes were observed with rat and kangaroo α-crystallins but the dogfish protein was not affected. Sedimentation velocity analyses and fluorescence spectroscopy yielded a pK, for the dissociation, of 3.7 for α-crystallin and 4.0 for a homopolymer constructed from purified αB₂ polypeptides. An αA₂ homopolymer was virtually unaffected by the lowering of pH. The products from the dissociation were isolated and their properties studied by sedimentation analysis and acrylamide quenching of tryptophan fluorescence. The αB chains were found to be completely denatured, whereas the structure of the αA chains, in the 290 kDa, particle, were only slightly altered. Comparisons of the sequences of the various proteins examined suggested that decreased ionization of aspartic acid 127 in the αB chain was responsible for the specific dissociation of this polypeptide.     

Ionic interactions for substituted MCH1R inhibitors studied by pKa values

MCH1R inhibitors with the quinoline moiety having the aromatic amine and aliphatic amine chain were selected, and then the effect of substituents of the quinoline ring on the ionic interaction were studied by calculating pK_a values for these amines at the B3LYP/6-311++G(d,p)//B3LYP/6-31+G(d) level in the gas phase and in water. For substituent with C, N, O, and S atoms next to the quinoline ring, respectively, the pK_a values of aromatic amines are estimated to be 8.98, 12.19, 4.64, and 4.33 and those of the aliphatic amines are 12.65, 10.82, 9.94, and 11.55, respectively.     

Study of the interaction of trypsin inhibitor from the sea anemone Radianthus macrodactylus with proteases

The interaction of the inhibitor VJ (InhVJ), isolated from sea anemone R. macrodactylus, with different proteases was investigated using the method of biosensor analysis. The following enzymes were tested: serine proteases (trypsin, α-chymotrypsin, plasmin, thrombin, kallikrein), cysteina protease (papain) and aspartic protease (pepsin). In the rage of the concentrations studied (10–400 nM) inhibitor VJ interacted only with trypsin and α-chymotrypsin. The intermolecular complexes formation between inhibitor VJ and each of these enzymes was characterized by the following kinetic and thermodynamics parameters: K_D = 7.38 × 10^?8 M and 9.93 × 10^?7 M for pairs InhVJ/trypsin and InhVJ/α-chymotrypsin, respectively.     

Origin of the pKa shift of the catalytic lysine in acetoacetate decarboxylase

The pK_a value of Lys115, the catalytic residue in acetoacetate decarboxylate, was calculated using atomic coordinates of the X-ray crystal structure with consideration of the protonation states of all titratable sites in the protein. The calculated pK_a value of Lys115 (pK_a(Lys115)) was unusually low (≈6) in agreement with the experimentally measured value. Although charged residues impact pK_a(Lys115) considerably in the native protein, the significant pK_a(Lys115) downshift in the protein with respect to aqueous solution was mainly due to loss of the solvation energy in the catalytic active site relative to bulk water.