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.     

1H NMR studies of a two-iron: two-sulfur ferredoxin from halobacterium of the dead sea

Ferredoxin isolated from Halobacterium of the Dead Sea (HFd) was found to be stable and retain its conformation in 4–0.5 M salt solutions. Reconstitution of the denatured protein to the oxidized form in ²H₂O indicated that the resonances shifted to the 8–10 ppm region, which include 18 protons, are nonexchangeable -NH protons. The C₂H and C₄H resonances of His-119 were assigned in both oxidized and reduced HFd. pH titration curves of these resonances yielded a pK_a for this His of 6.57 ± 0.1 and 6.65 ± 0.1 in oxidized and reduced HFd, respectively. pH titration curves, T₁ relaxation times, and the temperature dependence of the chemical shift were obtained for resonances between 6 and 10 ppm of oxidized HFd. In oxidized HFd a paramagnetically shifted resonance was observed at 15 ppm with 1 H intensity, and an anti-Curie temperature dependence. In reduced HFd eight resonances each with 1 H intensity were shifted downfield by 10–50 ppm and one resonance with 1 H intensity was shifted upfield to ?6.8 ppm. Four of these resonances exhibited an anti-Curie temperature dependence, two exhibited a moderate Curie dependence, and three were temperature independent.     

Carboxyl pKa Values and Acid Denaturation of BBL

The protein BBL undergoes structural transitions and acid denaturation between pH 1.2 and 8.0. Using NMR spectroscopy, we measured the pK_a values of all the carboxylic residues in this pH range. We employed ¹³C direct-detection two-dimensional IPAP (in-phase antiphase) CACO NMR spectroscopy to monitor the ionization state of different carboxylic groups and demonstrated its advantages over other NMR techniques in measuring pK_a values of carboxylic residues. The two residues Glu161 and Asp162 had significantly lowered pK_a values, showing that these residues are involved in a network of stabilizing electrostatic interactions, as is His166. The other carboxylates had unperturbed values. The pH dependence of the free energy of denaturation was described quantitatively by the ionizations of those three residues of perturbed pK_a, and, using thermodynamic cycles, we could calculate their pK_as in the native and denatured states as well as the equilibrium constants for denaturation of the different protonation states. We also measured ¹³C^α chemical shifts of individual residues as a function of pH. These shifts sense structural transitions rather than ionizations, and they titrated with pH consistent with the change in equilibrium constant for denaturation. Kinetic measurements of the folding of BBL E161Q indicated that, at pH 7, the stabilizing interactions with Glu161 are formed mainly in the transition state. We also found that local interactions still exist in the acid-denatured state of BBL, which attenuate somewhat the flexibility of the acid-denatured state.     

pH dependence of conformational fluctuations of the protein backbone

Proton binding equilibria (pK_a values) of ionizable groups in proteins are exquisitely sensitive to their microenvironments. Apparent pK_a values measured for individual ionizable residues with NMR spectroscopy are actually population‐weighted averages of the pK_a in different conformational microstates. NMR spectroscopy experiments with staphylococcal nuclease were used to test the hypothesis that pK_a values of surface Glu and Asp residues are affected by pH‐sensitive fluctuations of the backbone between folded and locally unfolded conformations. ¹⁵N spin relaxation studies showed that as the pH decreases from the neutral into the acidic range the amplitudes of backbone fluctuations in the ps‐ns timescale increase near carboxylic residues. Hydrogen exchange experiments suggested that backbone conformational fluctuations promoted by decreasing pH also reflect slower local or sub‐global unfolding near carboxylic groups. This study has implications for structure‐based pK_a calculations: (1) The timescale of the backbone's response to ionization events in proteins can range from ps to ms, and even longer; (2) pH‐sensitive fluctuations of the backbone can be localized to both the segment the ionizable residue is attached to or the one that occludes the ionizable group; (3) Structural perturbations are not necessarily propagated through Coulomb interactions; instead, local fluctuations appear to be coupled through the co‐operativity inherent to elements of secondary structure and to networks of hydrogen bonds. These results are consistent with the idea that local conformational fluctuations and stabilities are important determinants of apparent pK_a values of ionizable residues in proteins. Proteins 2014; 82:3132–3143. © 2014 Wiley Periodicals, Inc.     

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.     

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.     

Breaking the Carboxyl Rule: LYSINE 96 FACILITATES REPROTONATION OF THE SCHIFF BASE IN THE PHOTOCYCLE OF A RETINAL PROTEIN FROM EXIGUOBACTERIUM SIBIRICUM*?

A lysine instead of the usual carboxyl group is in place of the internal proton donor to the retinal Schiff base in the light-driven proton pump of Exiguobacterium sibiricum (ESR). The involvement of this lysine in proton transfer is indicated by the finding that its substitution with alanine or other residues slows reprotonation of the Schiff base (decay of the M intermediate) by more than 2 orders of magnitude. In these mutants, the rate constant of the M decay linearly decreases with a decrease in proton concentration, as expected if reprotonation is limited by the uptake of a proton from the bulk. In wild type ESR, M decay is biphasic, and the rate constants are nearly pH-independent between pH 6 and 9. Proton uptake occurs after M formation but before M decay, which is especially evident in D₂O and at high pH. Proton uptake is biphasic; the amplitude of the fast phase decreases with a pK_a of 8.5 ± 0.3, which reflects the pK_a of the donor during proton uptake. Similarly, the fraction of the faster component of M decay decreases and the slower one increases, with a pK_a of 8.1 ± 0.2. The data therefore suggest that the reprotonation of the Schiff base in ESR is preceded by transient protonation of an initially unprotonated donor, which is probably the ϵ-amino group of Lys-96 or a water molecule in its vicinity, and it facilitates proton delivery from the bulk to the reaction center of the protein.     

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.     

Elucidating Relayed Proton Transfer through a His–Trp–His Triad of a Transmembrane Proton Channel by Solid-State NMR

Proton transfer through membrane-bound ion channels is mediated by both water and polar residues of proteins, but the detailed molecular mechanism is challenging to determine. The tetrameric influenza A and B virus M2 proteins form canonical proton channels that use an HxxxW motif for proton selectivity and gating. The BM2 channel also contains a second histidine (His), H27, equidistant from the gating tryptophan, which leads to a symmetric H¹⁹xxxW²³xxxH²⁷ motif. The proton-dissociation constants (pK_a's) of H19 in BM2 were found to be much lower than the pK_a's of H37 in AM2. To determine if the lower pK_a's result from H27-facilitated proton dissociation of H19, we have now investigated a H27A mutant of BM2 using solid-state NMR. ¹⁵N NMR spectra indicate that removal of the second histidine converted the protonation and tautomeric equilibria of H19 to be similar to the H37 behavior in AM2, indicating that the peripheral H27 is indeed the origin of the low pK_a's of H19 in wild-type BM2. Measured interhelical distances between W23 sidechains indicate that the pore constriction at W23 increases with the H19 tetrad charge but is independent of the H27A mutation. These results indicate that H27 both accelerates proton dissociation from H19 to increase the inward proton conductance and causes the small reverse conductance of BM2. The proton relay between H19 and H27 is likely mediated by the intervening gating tryptophan through cation–π interactions. This relayed proton transfer may exist in other ion channels and has implications for the design of imidazole-based synthetic proton channels.     

Electrostatic contributions to the stability of the GCN4 leucine zipper structure

Ion pairs are ubiquitous in X-ray structures of coiled coils, and mutagenesis of charged residues can result in large stability losses. By contrast, pK_a values determined by NMR in solution often predict only small contributions to stability from charge interactions. To help reconcile these results we used triple-resonance NMR to determine pK_a values for all groups that ionize between pH 1 and 13 in the 33 residue leucine zipper fragment, GCN4p. In addition to the native state we also determined comprehensive pK_a values for two models of the GCN4p denatured state: the protein in 6 M urea, and unfolded peptide fragments of the protein in water. Only residues that form ion pairs in multiple X-ray structures of GCN4p gave large pK_a differences between the native and denatured states. Moreover, electrostatic contributions to stability were not equivalent for oppositely charged partners in ion pairs, suggesting that the interactions between a charge and its environment are as important as those within the ion pair. The pH dependence of protein stability calculated from NMR-derived pK_a values agreed with the stability profile measured from equilibrium urea-unfolding experiments as a function of pH. The stability profile was also reproduced with structure-based continuum electrostatic calculations, although contributions to stability were overestimated at the extremes of pH. We consider potential sources of errors in the calculations, and how pK_a predictions could be improved. Our results show that although hydrophobic packing and hydrogen bonding have dominant roles, electrostatic interactions also make significant contributions to the stability of the coiled coil.     

Metal ion-biomolecule interactions. II. Methylmercuration of the deprotonated amino groups in adenine,guanine, and cytosine derivatives,and its relationship to amino group acidity

Proton nmr spectroscopic evidence is presented for methylmercury(II) binding to the deprotonated amino groups in adenosine, 9-methyladenine, guanosine, 1-methylguanosine, and cytidine under basic conditions. Except for the guanosine case, ¹H nmr spectra of the products from aqueous or ethanolic 1:1 mixtures of substrate and MeHgOH are consistent with methylmercuration of the deprotonated amino groups. Guanosine undergoes initial binding of MeHg to N₁, and a second equivalent of MeHgOH is necessary to effect amino binding. The nmr spectra of the complexed adenine derivatives suggest that different geometrical isomers exist in (CD₃)₂SO solution, reflecting the partial double bond character of the C₆N bond in these systems. Using a correlation relating the magnitude of the ¹⁹⁹Hg-¹H coupling constant (J) for MeHg-ligand complexes with the ligand pK_a (J = ?3.88 pK_a + 248.5, extending over 13 pK units, based on a variety of N and O donor ligands), estimates (± 0.3 pK unit) of the pK_as of the amino groups of the above substrates have been made. In this way, pK_a values of 15.5 (cytidine), 17.0 (adenosine and 9-methyladenine), 15.1 (guanosine), and 14.9 (1-methylguanosine) are obtained. In the cases where comparisons with literature pK_a data can be made, good agreement is found.     

Determination of the ionization constants and reactivities of the amino-termini of -chymotrypsin

The ionization constants and reactivities of the three amino-termini of α-chymotrypsin have been determined by competitive labelling with acetic anhydride. In 0.1 m-KCl at 10 °C the half cystine-1 amino-terminus of the A chain has a pK_a of 7.9 and is fully exposed on the surface of the protein; the isoleucine-16 amino-terminus of the B chain has a pK_a of 8.9 and is buried; the alanine-149 amino-terminus of the B chain has a pK_a of 8.5 and is slightly buried. The reactivities of the half cystine-1 and isoleucine-16 amino-termini show that above pH 9.8 α-chymotrypsin undergoes a structural change. The reactivity of the alanine-149 amino-terminus, but not the reactivity of the half cystine-1 aminoterminus, is dependent on the ionization state of the isoleucine-16 aminoterminus. It is concluded that the ionization state of the isoleucine-16 aminoterminus controls the activity of α-chymotrypsin but the internal proteolytic cleavages, which generate the alanine-149 amino-terminus, weaken the hydrophobic interactions in the amino-terminal region of the B chain and make the enzyme more sensitive to the effects of deprotonating the isoleucine-16 aminoterminus.     

The deacylation of substituted benzol-alpha-chymotrypsins.

An extensive comparison of deacylation rates of mono- and disubstituted benzoyl-α-chymotrypsins indicates that no steric effects on rate or apparent pK_a of deacylation are detectable within this series. Some anomalous effects on deacylation rate appear to be associated with fluoro- and nitro-substituents in particular positions on the ring and may be attributable to specific interactions at the enzyme active site. The extensive series of structurally similar acyl-enzymes prepared has allowed a thorough analysis of the effect of acyl group pK_a on the apparent pK_a of deacylation. The data indicates that polar effects on the apparent pK_a are probably negligible. Rho for the deacylation reaction is in good agreement with model reactions for an imidazole general base-catalyzed model reaction.     

Ionisations within a subtilisin–glyoxal inhibitor complex

Z-Ala-Pro-Phe-glyoxal (where Z is benzyloxycarbonyl) has been shown to be a competitive inhibitor of subtilisin with a K_i=2.3±0.2 μM at pH 7.0 and 25 °C. Using Z-Ala-Pro-[2-¹³C]Phe-glyoxal we have detected a signal at 107.3 ppm by ¹³C NMR, which we assign to the tetrahedral adduct formed between the hydroxy group of serine-195 and the ¹³C-enriched keto-carbon of the inhibitor. The chemical shift of this signal is pH independent from pH 4.2 to 7.0 and we conclude that the oxyanion pK_a<3. This is the first observation of oxyanion formation in a reversible subtilisin–inhibitor complex. The inhibitor is bound as a hemiketal which is in slow exchange with the free inhibitor. Inhibitor binding depends on a pK_a of ～6.5 in the free enzyme and on a pK_a<3.0 when the inhibitor is bound to subtilisin. Protonation of the oxyanion promotes the disassociation of the inhibitor. We show that oxyanion formation cannot be rate limiting during catalysis and that subtilisin stabilises the oxyanion by at least 45.1 kJ mol^?1. We conclude that if the energy required for oxyanion stabilisation is utilised as binding energy in drug design it should make a significant contribution to inhibitor potency.     

Solid-state NMR and functional studies on proteorhodopsin

Proteins of the proteorhodopsin (PR) family are found abundantly in many marine bacteria in the photic zone of the oceans. They are colour-tuned to their environment. The green absorbing species has been shown to act as a light-driven proton pump and thus could form a potential source of energy. The pK_a of the primary proton acceptor is close to the pH of seawater which could also indicate a regulatory role. Here, we review and summarize our own recent findings in the context of known data and present some new results. Proton transfer in vitro by PR is shown by a fluorescence assay which confirms a pH dependent vectoriality. Previously reported low diffracting 2D crystal preparations of PR are assessed for their use for solid-state NMR by two dimensional ¹³C-¹³C DARR spectra. ¹⁵N-¹H HETCOR MAS NMR experiments show bound water in the vicinity of the protonated Schiff base which could play a role in proton transfer. The effect of highly conserved H75 onto the properties of the chromophore has been investigated by single site mutations. They do show a pronounced effect onto the optical absorption maximum and the pK_a of the proton acceptor but have only a small effect onto the ¹⁵N chemical shifts of the protonated Schiff base.     

Correlation of exchangeable (NH) and nonexchangeable (C2H) histidine resonances in the proton NMR spectrum of ribonuclease A in aqueous solution.

The ionization characteristics of the hydrogen-bonded His 12 N₁ proton observed to titrate between 11 to 13 ppm in the nmr spectrum of ribonuclease A in H₂O solution are compared with the ionization characteristics of the four histidine C₂ protons in the enzyme. Comparison of the pK_a's of the enzyme in H₂O and D₂O in the absence and presence of cytidine monophosphate (?5′, ?3′, and ?2′) inhibitors, line widths in the presence of Cu II at pH 3.6 and 5.6, and chemical shifts in the presence of AgNO₃ permit a correlation of the exchangeable His 12 N₁ proton with the active site histidine C₂ proton exhibiting the lower ionization pK_a. The histidines with pK_a of 5.1 and 5.6 in ribonuclease A in the absence of salt are assigned in this study to His 12 and His 119, respectively.     

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.     

13C magnetic resonance study of the ionization of N-acetyl-dl-proline in aqueous solution

NMR titration curves have been recorded for all the ¹³C resonances of cis and trans$\text{N-}\text{acetyl}\text{-}\text{dl}\text{-}\text{proline}$ in ²H₂O. the measured pK_2H values are 3.4 ± 0.8 and 4.13 ± 0.08 respectively; the free energy of ionization for the trans isomer being (3.8 kJ/mole) greater than for the cis. The ionization shifts of the two isomers differ significantly only at the acetyl carbonyl and C_γ positions. It is suggested that these are related to conformational changes which stabilize the trans form at low p²H.     

Effects of ionization of p-nitrocatechol sulfate on its behavior as a substrate for arylsulfatases

Evidence is presented showing that arylsulfatases reacting at acidic pH are more active towards p-nitrocatechol sulfate substrate, while the reverse is true for the enzymes having alkaline pH_max. An explanation of this phenomenon is suggested based on the ionization of p-nitrocatechol sulfate (pK_a = 6.4) and a much higher reactivity of sulfatases with the acidic form of the substrate.     

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.