Histidine hydrogen-deuterium exchange mass spectrometry for probing the microenvironment of histidine residues in dihydrofolate reductase

Background

Histidine Hydrogen-Deuterium Exchange Mass Spectrometry (His-HDX-MS) determines the HDX rates at the imidazole C₂-hydrogen of histidine residues. This method provides not only the HDX rates but also the pK _a values of histidine imidazole rings. His-HDX-MS was used to probe the microenvironment of histidine residues of E. coli dihydrofolate reductase (DHFR), an enzyme proposed to undergo multiple conformational changes during catalysis.

Methodology/Principal Findings

Using His-HDX-MS, the pK _a values and the half-lives (t _1/2) of HDX reactions of five histidine residues of apo-DHFR, DHFR in complex with methotrexate (DHFR-MTX), DHFR in complex with MTX and NADPH (DHFR-MTX-NADPH), and DHFR in complex with folate and NADP⁺ (DHFR-folate-NADP⁺) were determined. The results showed that the two parameters (pK _a and t _1/2) are sensitive to the changes of the microenvironment around the histidine residues. Although four of the five histidine residues are located far from the active site, ligand binding affected their pK _a, t _1/2 or both. This is consistent with previous observations of ligand binding-induced distal conformational changes on DHFR. Most of the observed pK _a and t _1/2 changes could be rationalized using the X-ray structures of apo-DHFR, DHFR-MTX-NADPH, and DHFR-folate-NADP⁺. The availability of the neutron diffraction structure of DHFR-MTX enabled us to compare the protonation states of histidine imidazole rings.

Conclusions/Significance

Our results demonstrate the usefulness of His-HDX-MS in probing the microenvironments of histidine residues within proteins.     

The pKa values of two histidine residues in human haemoglobin, the Bohr effect, and the dipole moments of alpha-helices

Studies of abnormal and chemically modified haemoglobins indicate that in 0.1 m-NaCl about 40% of the alkaline Bohr effect of human haemoglobin is contributed by the C-terminal histidine HC3(146)β. In deoxyhaemoglobin, the imidazole of this histidine forms a salt bridge with aspartate FG1(94)β, in oxyhaemoglobin or carbonmonoxyhaemoglobin it accepts a hydrogen bond from its own NH group instead. Kilmartin et al. (1973) showed that in 0.2 m-NaCl + 0.2 m-phosphate this change of ligation lowered the pK_a of the histidine from 8.0 in Hb³ to 7.1 in HbCO, but Russu et al. (1980) claimed that in bis-Tris buffer without added NaCl its pK_a in HbCO dropped no lower than 7.85, and that in this medium the C-terminal histidine made only a negligible contribution to the alkaline Bohr effect.We have compared the histidine resonances of HbCO A with those of three abnormal haemoglobins: HbCO Cowtown (His HC3(146)β → Leu), HbCO Wood (His FG4(97)β → Leu) and HbCO Malmø (His FG4(97)β → Gln). Our results show that the resonance assigned by Russu et al. to His HC3(146)β in fact belongs to His FG4(97)β. Although in Hb the pK_a of His HC3(146)β is 8.05 ± 0.05 independent of ionic strength, in HbCO its pK_a drops sharply with diminishing ionic strength, so that in the buffer employed by Russu et al. it has a pK_a of 6.2 and makes a contribution to the alkaline Bohr effect that is 57% larger than in the phosphate buffer employed by Kilmartin et al. (1973).In HbCO A, His FG4(97)β does not contribute to the Bohr effect, but in HbCO from which His HC3(146)β has been cleaved (HbCO des-His), His FG4(97)β is in equilibrium between two conformations with different pK_a values. This equilibrium varies with ionic strength and pH, and presumably also with degree of ligation of the haem moiety.In HbCO A, His FG4(97)β has a pK_a of 7.8 compared to the pK_a value of about 6.6 characteristic of free histidines at the surface of proteins. This high pK_a is accounted for by its interaction with the negative pole at the C terminus of helices F and FG. It corresponds to a free energy change of the same order as that observed in the interaction of histidines with carboxylate ions and confirms the strongly dipolar character of α-helices, which manifests itself even when they lie on the surface of the protein.     

Ionization constants and thermodynamic parameters of histidine and derivatives

The pK_a values for the proton dissociation of carboxyl, imidazolium, and ammonium groups for histidine and ten of its derivatives were determined electrometrically at seven temperatures in the range 10–40°C. The ΔH and ΔS values were estimated from the temperature dependence of the dissociation constants of histidine and its derivatives. These results and the pK_a values compared in terms of inductive effect suggest an ion-dipole interaction between the protonated amino group and the unprotonated imidazole ring. The charge and the solvation effects of the neighboring groups are the main factors that determine the imidazole group pK_a in histidine and its studied derivatives. The N^τ-H tautomer is favored over the N^π-H by 1.6 kcal/mol, indicating that the inductive substituent effect at position 4 of the imidazole ring is the major component in determining this tautomeric preference.     

Proton nuclear magnetic resonance studies of histidine residues in rat and other rodent pancreatic ribonucleases. Effects of pH and inhibitors

Proton nuclear magnetic resonance spectra of the histidine residues in bovine and rat ribonuclease have been compared. The changes in chemical shift on titration and on binding of cytidine-3′-monophosphate and cytidine-2′-monophosphate have been followed. In the presence of the cytidine derivatives the spectra of both enzymes resemble each other more than those of the free enzymes. With these inhibitors, two histidines in rat ribonuclease exhibit the same pK values and shifts as the active site residues histidine 12 and 119 in the bovine enzyme. Their pK values in the inhibitor-free rat enzyme are about 0.4 higher than in the beef enzyme, which can be explained by the substitution at the entrance of the active site cleft of arginine 39 in the beef enzyme by serine in the rat enzyme. Rat ribonuclease contains one histidine with a rather high pK value of 7.6. The cytidine derivatives affect its chemical shift in exactly the same way as the shift of histidine 48 in bovine ribonuclease. The high pK value of this residue in rat ribonuclease can be explained by assuming a strong hydrogen bridge with glutamic acid 16. The other two histidines in rat ribonuclease have rather low pK values of 6.1 and 6.3. The histidine with a pK value of 6.3 has been assigned to position 105 and that with a pK value of 6.1 to position 73.The closer resemblance of the active sites of bovine and rat ribonuclease in the presence of inhibitors than in the inhibitor-free enzymes makes the concept of induced fit interesting from an evolutionary point of view.The characteristic downfield shift of the protonated form of histidine 119 in the complexes of bovine and rat ribonuclease with cytidine-3′-monophosphate is not observed with uridine-3′-monophosphate, suggesting non-identical binding of these pyrimidine nucleotides.Some preliminary results on the nuclear magnetic resonance properties of the histidine residues in coypu and chinchilla pancreatic ribonuclease have been obtained.     

1H nuclear magnetic resonance studies of histidine-containing di- and tripeptides. Estimation of the effects of charged groups on the pKa value of the imidazole ring.

The nmr titration curves of chemical shifts versus pH were observed for the protons of various histidine-containing di- and tripeptides. With these results, the macroscopic pK_a values and the chemical shifts intrinsic to each ionic species were determined by a computer curve-fitting based on a simple acid dissociation sequence. The pK_a value of the imidazole ring in N-acetyl-L -histidine methylamide was assumed to represent the intrinsic (or unperturbed) pK_a of the imidazole rings of histidine having peptide linkages at both the CO and NH sides. The pK_a values of the imidazole rings observed for most di- and tripeptides were reasonably reproduced by simple calculations using the intrinsic value and the perturbations due to the CO₂^? and NH₃⁺ groups located at various positions. Some other factors affecting the pK_a value of the imidazole ring are also discussed.     

Functional tuning of the catalytic residue pKa in a de novo designed esterase

AlleyCatE is a de novo designed esterase that can be allosterically regulated by calcium ions. This artificial enzyme has been shown to hydrolyze p‐nitrophenyl acetate (pNPA) and 4‐nitrophenyl‐(2‐phenyl)‐propanoate (pNPP) with high catalytic efficiency. AlleyCatE was created by introducing a single‐histidine residue (His¹⁴⁴) into a hydrophobic pocket of calmodulin. In this work, we explore the determinants of catalytic properties of AlleyCatE. We obtained the pK_a value of the catalytic histidine using experimental measurements by NMR and pH rate profile and compared these values to those predicted from electrostatics pK_a calculations (from both empirical and continuum electrostatics calculations). Surprisingly, the pK_a value of the catalytic histidine inside the hydrophobic pocket of calmodulin is elevated as compared to the model compound pK_a value of this residue in water. We determined that a short‐range favorable interaction with Glu¹²⁷ contributes to the elevated pK_a of His¹⁴⁴. We have rationally modulated local electrostatic potential in AlleyCatE to decrease the pK_a of its active nucleophile, His¹⁴⁴, by 0.7 units. As a direct result of the decrease in the His¹⁴⁴ pK_a value, catalytic efficiency of the enzyme increased by 45% at pH 6. This work shows that a series of simple NMR experiments that can be performed using low field spectrometers, combined with straightforward computational analysis, provide rapid and accurate guidance to rationally improve catalytic efficiency of histidine‐promoted catalysis. Proteins 2017; 85:1656–1665. © 2017 Wiley Periodicals, Inc.     

Comparative studies of the interaction of human and bovine platelet factor 4 with heparin using histidine NMR resonances as spectroscopic probes

ThepK _a values of His-38 and His-50 of the heparin-binding protein, bovine platelet factor 4, are 5.6 and 6.5, respectively, as determined by¹H NMR spectroscopy. The¹H NMR resonance of His-38 of bovine platelet factor 4 which exhibits the lowerpK _a value is perturbed upon heparin binding to a greater degree than the resonance of His-50. Human platelet factor 4 contains the homologous residues His-23 and His-35. ThepK _a values of the two histidine residues of human platelet factor 4 are 5.3 and 6.4. The¹H NMR resonance of the histidine of human platelet factor 4 exhibiting the lowerpK _a value also is perturbed upon heparin binding to a greater degree than the histidine resonance exhibiting the higherpK _a , thereby suggesting comparable heparin-protein interactions in bovine and human platelet factor 4.     

Mutation of non-conserved amino acids surrounding catalytic site to shift pH optimum of Bacillus circulans xylanase

Improvement of enzyme function by engineering pH dependence of enzymatic activity is of importance for industrial application of Bacillus circulans xylanases. Target mutation sites were selected by structural alignment between B. circulans xylanase and other xylanases having different pH optima. We selected non-conserved mutant sites within 8 Å from the catalytic residues, to see whether these residues have some role in modulating pK_as of the catalytic residues. We hypothesized that the non-conserved residues which may not have any role in enzyme catalysis might perturb pK_as of the catalytic residues. Change in pK_a of a titratable group due to change in electrostatic potential of a mutation was calculated and the change in pH optimum was predicted from the change in pK_a of the catalytic residues. Our strategy is proved to be useful in selection of promising mutants to shift the pH optimum of the xylanases towards desired side.     

Application of high-performance liquid chromatography to competitive labeling studies: The chemical properties of functional groups of glucaton

A competitive-labeling study of glucagon was carried out using [³H]- and [¹⁴C]-1-fluoro-2,4-dinitrobenzene to determine simultaneously the chemical properties of the α-amino and imidazole groups of the N-terminal histidine residue, and the lysine and tyrosine residues, under conditions where glucagon is in its physiologically active monomer form. The dinitrophenyl derivatives of these groups were purified by high-performance liquid chromatography which greatly simplified the separation steps of the procedure. The results showed the α-amino and tyrosine groups to have relatively normal behavior, with pK values of 7.98 and 10.22, respectively, while the lysine had a low pK of 8.46. The imidazole function had an apparent pK of 7.84, substantially higher than previous estimates. This difference may be accounted for by the effect of the charged form of the adjacent α-amino group on the nucleophilicity of the imidazole group.     

Characterization and determination of titratable groups of proteins by linearization of titration curves. II. Application to lysozyme

The potentiometric acid-base titration curve of fully protonated lysozyme at ionic strengths of 0.10 and 1.0 m has been performed. The stoichiometry and the pK_a values of each titratable group have been determined through the linearization of titration curves. Two types of carboxylic groups with pK_a values of 3.76 and 5.02, the imidazole group with pK_a 7.37 and the amine group with pK_a 9.63, have been identified at an ionic strength of 0.10 m at 25.0°C. The number of titratable groups found per mole of protein has been 5.12 and 5.60 for the two types of carboxylic groups, 1.13 for the imidazole group, and 3.19 for the amino groups. The endpoint of the titration of the protein obtained by this method accords quite well with the endpoint obtained by the use of Gran function applied to the excess of strong base.     

Histidine proton resonances of carbonmonoxyhaemoglobins A and Cowtown in chloride-free buffer

A re-examination of the C-2 histidine proton resonances of haemoglobins A and Cowtown (His HC3(146)β → Leu) in chloride-free Hepes buffer has shown that all the resonances present in haemoglobin A are present in haemoglobin Cowtown, so that the pK_a of His HC3(146)β cannot be determined by nuclear magnetic resonance in this buffer.     

An Extracellular Ion Pathway Plays a Central Role in the Cooperative Gating of a K2P K+ Channel by Extracellular pH

Proton-gated TASK-3 K⁺ channel belongs to the K_2P family of proteins that underlie the K⁺ leak setting the membrane potential in all cells. TASK-3 is under cooperative gating control by extracellular [H⁺]. Use of recently solved K_2P structures allows us to explore the molecular mechanism of TASK-3 cooperative pH gating. Tunnel-like side portals define an extracellular ion pathway to the selectivity filter. We use a combination of molecular modeling and functional assays to show that pH-sensing histidine residues and K⁺ ions mutually interact electrostatically in the confines of the extracellular ion pathway. K⁺ ions modulate the pK_a of sensing histidine side chains whose charge states in turn determine the open/closed transition of the channel pore. Cooperativity, and therefore steep dependence of TASK-3 K⁺ channel activity on extracellular pH, is dependent on an effect of the permeant ion on the channel pH_o sensors.     

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.     

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.     

Probing the role of active site histidine residues in the catalytic activity of lacrimal gland peroxidase

The role of active site histidine residues in SCN^– oxidation by lacrimal gland peroxidase (LGP) has been probed after modification with diethylpyrocarbonate (DEPC). The enzyme is irreversibly inactivated following pseudo-first order kinetics with a second order rate constant of 0.26 M^–1 sec^–1 at 25°C. The pH dependent rate of inactivation shows an inflection point at 6.6 indicating histidine derivatization. The UV difference spectrum of the modified versus native enzyme shows a peak at 242 nm indicating formation of N-carbethoxyhistidine. Carbethoxyhistidine formation and associated inactivation are reversed by hydroxylamine indicating histidine modification. The stoichiometry of histidine modification and the extent of inactivation show that out of five histidine residues modified, modification of two residues inactivates the enzyme. Substrate protection with SCN^– during modification indicates that although one histidine is protected, it does not prevent inactivation. The spectroscopically detectable compound II formation is lost due to modification and is not evident after SCN^– protection. The data indicate that out of two histidines, one regulates compound I formation while the other one controls SCN^– binding. SCN^– protected enzyme is inactive due to loss of compound I formation. SCN^– binding studies by optical difference spectroscopy indicate that while the native enzyme binds SCN^– with the K_d of 15 mM, the modified enzyme shows very weak binding with the K_d of 660 mM. From the pH dependent binding of SCN^–, a plot of log K_d vs. pH shows a sigmoidal curve from which the involvement of an enzyme ionizable group of pKa 6.6 is ascertained and attributed to the histidine residue controlling SCN^– binding. LGP has thus two distinctly different essential histidine residues – one regulates compound I formation while the other one controls SCN^– binding.     

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.     

Thermodynamic and structural investigation of the specific SDS binding of humicola insolens cutinase

The interaction of lipolytic enzymes with anionic surfactants is of great interest with respect to industrially produced detergents. Here, we report the interaction of cutinase from the thermophilic fungus Humicola insolens with the anionic surfactant SDS, and show the enzyme specifically binds a single SDS molecule under nondenaturing concentrations. Protein interaction with SDS was investigated by NMR, ITC and molecular dynamics simulations. The NMR resonances of the protein were assigned, with large stretches of the protein molecule not showing any detectable resonances. SDS is shown to specifically interact with the loops surrounding the catalytic triad with medium affinity (K_a ≈ 10⁵ M⁻¹). The mode of binding is closely similar to that seen previously for binding of amphiphilic molecules and substrate analogues to cutinases, and hence SDS acts as a substrate mimic. In addition, the structure of the enzyme has been solved by X-ray crystallography in its apo form and after cocrystallization with diethyl p-nitrophenyl phosphate (DNPP) leading to a complex with monoethylphosphate (MEP) esterified to the catalytically active serine. The enzyme has the same fold as reported for other cutinases but, unexpectedly, esterification of the active site serine is accompanied by the ethylation of the active site histidine which flips out from its usual position in the triad.     

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.     

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.     

Histidine Residues in the Na+-coupled Ascorbic Acid Transporter-2 (SVCT2) Are Central Regulators of SVCT2 Function,Modulating pH Sensitivity,Transporter Kinetics,Na+ Cooperativity,Conformational Stability,and Subcellular Localization

Na⁺-coupled ascorbic acid transporter-2 (SVCT2) activity is impaired at acid pH, but little is known about the molecular determinants that define the transporter pH sensitivity. SVCT2 contains six histidine residues in its primary sequence, three of which are exofacial in the transporter secondary structure model. We used site-directed mutagenesis and treatment with diethylpyrocarbonate to identify histidine residues responsible for SVCT2 pH sensitivity. We conclude that five histidine residues, His¹⁰⁹, His²⁰³, His²⁰⁶, His²⁶⁹, and His⁴¹³, are central regulators of SVCT2 function, participating to different degrees in modulating pH sensitivity, transporter kinetics, Na⁺ cooperativity, conformational stability, and subcellular localization. Our results are compatible with a model in which (i) a single exofacial histidine residue, His⁴¹³, localized in the exofacial loop IV that connects transmembrane helices VII-VIII defines the pH sensitivity of SVCT2 through a mechanism involving a marked attenuation of the activation by Na⁺ and loss of Na⁺ cooperativity, which leads to a decreased V_max without altering the transport K_m; (ii) exofacial histidine residues His²⁰³, His²⁰⁶, and His⁴¹³ may be involved in maintaining a functional interaction between exofacial loops II and IV and influence the general folding of the transporter; (iii) histidines 203, 206, 269, and 413 affect the transporter kinetics by modulating the apparent transport K_m; and (iv) histidine 109, localized at the center of transmembrane helix I, might be fundamental for the interaction of SVCT2 with the transported substrate ascorbic acid. Thus, histidine residues are central regulators of SVCT2 function.