NMR evidence for the participation of a low-barrier hydrogen bond in the mechanism of delta 5-3-ketosteroid isomerase.

Delta 5-3-Ketosteroid isomerase (EC 5.3.3.1) promotes an allylic rearrangement involving intramolecular proton transfer via a dienolic intermediate. This enzyme enhances the catalytic rate by a factor of 10(10). Two residues, Tyr-14, the general acid that polarizes the steroid 3-carbonyl group and facilitates enolization, and Asp-38 the general base that abstracts and transfers the 4 beta-proton to the 6 beta-position, contribute 10(4.7) and 10(5.6) to the rate increase, respectively. A major mechanistic enigma is the huge disparity between the pKa values of the catalytic groups and their targets. Upon binding of an analog of the dienolate intermediate to isomerase, proton NMR detects a highly deshielded resonance at 18.15 ppm in proximity to aromatic protons, and with a 3-fold preference for protium over deuterium (fractionation factor, phi = 0.34), consistent with formation of a short, strong (low-barrier) hydrogen bond to Tyr-14. The strength of this hydrogen bond is estimated to be at least 7.1 kcal/mol. This bond is relatively inaccessible to bulk solvent and is pH insensitive. Low-barrier hydrogen bonding of Tyr-14 to the intermediate, in conjunction with the previously demonstrated tunneling contribution to the proton transfer by Asp-38, provide a plausible and quantitative explanation for the high catalytic power of this isomerase.

Kinetics of hydrogen bond breakage in the process of unfolding of ribonuclease A measured by pulsed hydrogen exchange.

A sensitive test for kinetic unfolding intermediates in ribonuclease A (EC 3.1.27.5) is performed under conditions where the enzyme unfolds slowly (10 degrees C, pH 8.0, 4.5 M guanidinium chloride). Exchange of peptide NH protons (2H-1H) is used to monitor structural opening of individual hydrogen bonds during unfolding, and kinetic models are developed for hydrogen exchange during the process of protein unfolding. The analysis indicates that the kinetic process of unfolding can be monitored by EX1 exchange (limited by the rate of opening) for ribonuclease A in these conditions. Of the 49 protons whose unfolding/exchange kinetics was measured, 47 have known hydrogen bond acceptor groups. To test whether exchange during unfolding follows the EX2 (base-catalyzed) or the EX1 (uncatalyzed) mechanism, unfolding/exchange was measured both at pH 8.0 and at pH 9.0. A few faster-exchanging protons were found that undergo exchange by both EX1 and EX2 processes, but the 43 slower-exchanging protons at pH 8 undergo exchange only by the EX1 mechanism, and they have closely similar rates. Thus, it is likely that all 49 protons undergo EX1 exchange at the same rate. The results indicate that a single rate-limiting step in unfolding breaks the entire network of peptide hydrogen bonds and causes the overall unfolding of ribonuclease A. The additional exchange observed for some protons that follows the EX2 mechanism probably results from equilibrium unfolding intermediates and will be discussed elsewhere.

The change in hydrogen bond strength accompanying charge rearrangement: Implications for enzymatic catalysis

The equilibrium for formation of the intramolecular hydrogen bond (KHB) in a series of substituted salicylate monoanions was investigated as a function of ΔpKa, the difference between the pKa values of the hydrogen bond donor and acceptor, in both water and dimethyl sulfoxide. The dependence of log KHB upon ΔpKa is linear in both solvents, but is steeper in dimethyl sulfoxide (slope = 0.73) than in water (slope = 0.05). Thus, hydrogen bond strength can undergo substantially larger increases in nonaqueous media than aqueous solutions as the charge density on the donor or acceptor atom increases. These results support a general mechanism for enzymatic catalysis, in which hydrogen bonding to a substrate is strengthened as charge rearranges in going from the ground state to the transition state; the strengthening of the hydrogen bond would be greater in a nonaqueous enzymatic active site than in water, thus providing a rate enhancement for an enzymatic reaction relative to the solution reaction. We suggest that binding energy of an enzyme is used to fix the substrate in the low-dielectric active site, where the strengthening of the hydrogen bond in the course of a reaction is increased.

Unusual 1H NMR chemical shifts support (His) Cɛ1—H⋅⋅⋅O⩵C H-bond: Proposal for reaction-driven ring flip mechanism in serine protease catalysis

13C-selective NMR, combined with inhibitor perturbation experiments, shows that the Cɛ1—H proton of the catalytic histidine in resting α-lytic protease and subtilisin BPN′ resonates, when protonated, at 9.22 ppm and 9.18 ppm, respectively, which is outside the normal range for such protons and ≈0.6 to 0.8 ppm further downfield than previously reported. They also show that the previous α-lytic protease assignments [Markley, J. L., Neves, D. E., Westler, W. M., Ibanez, I. B., Porubcan, M. A. & Baillargeon, M. W. (1980) Front. Protein Chem. 10, 31–61] were to signals from inactive or denatured protein. Simulations of linewidth vs. pH demonstrate that the true signal is more difficult to detect than corresponding signals from inactive derivatives, owing to higher imidazole pKa values and larger chemical shift differences between protonated and neutral forms. A compilation and analysis of available NMR data indicates that the true Cɛ1—H signals from other serine proteases are similarly displaced downfield, with past assignments to more upfield signals probably in error. The downfield displacement of these proton resonances is shown to be consistent with an H-bond involving the histidine Cɛ1—H as donor, confirming the original hypothesis of Derewenda et al. [Derewenda, Z. S., Derewenda, U. & Kobos, P. M. (1994) J. Mol. Biol. 241, 83–93], which was based on an analysis of literature x-ray crystal structures of serine hydrolases. The invariability of this H-bond among enzymes containing Asp-His-Ser triads indicates functional importance. Here, we propose that it enables a reaction-driven imidazole ring flip mechanism, overcoming a major dilemma inherent in all previous mechanisms, namely how these enzymes catalyze both the formation and productive breakdown of tetrahedral intermediates.

On the electronic nature of low-barrier hydrogen bonds in enzymatic reactions

The electronic nature of low-barrier hydrogen bonds (LBHBs) in enzymatic reactions is discussed based on combined low temperature neutron and x-ray diffraction experiments and on high level ab initio calculations by using the model substrate benzoylacetone. This molecule has a LBHB, as the intramolecular hydrogen bond is described by a double-well potential with a small barrier for hydrogen transfer. From an “atoms in molecules” analysis of the electron density, it is found that the hydrogen atom is stabilized by covalent bonds to both oxygens. Large atomic partial charges on the hydrogen-bonded atoms are found experimentally and theoretically. Therefore, the hydrogen bond gains stabilization from both covalency and from the normal electrostatic interactions found for long, weak hydrogen bonds. Based on comparisons with other systems having short-strong hydrogen bonds or LBHBs, it is proposed that all short-strong and LBHB systems possess similar electronic features of the hydrogen-bonded region, namely polar covalent bonds between the hydrogen atom and both heteroatoms in question.

Fractionation factors and activation energies for exchange of the low barrier hydrogen bonding proton in peptidyl trifluoromethyl ketone complexes of chymotrypsin

NMR investigations have been carried out of complexes between bovine chymotrypsin Aα and a series of four peptidyl trifluoromethyl ketones, listed here in order of increasing affinity for chymotrypsin: N-Acetyl-l-Phe-CF3, N-Acetyl-Gly-l-Phe-CF3, N-Acetyl-l-Val-l-Phe-CF3, and N-Acetyl-l-Leu-l-Phe-CF3. The D/H fractionation factors (φ) for the hydrogen in the H-bond between His 57 and Asp 102 (His 57-Hδ1) in these four complexes at 5°C were in the range φ = 0.32–0.43, expected for a low-barrier hydrogen bond. For this series of complexes, measurements also were made of the chemical shifts of His 57-Hɛ1 (δ2,2-dimethylsilapentane-5-sulfonic acid 8.97–9.18), the exchange rate of the His 57-Hδ1 proton with bulk water protons (284–12.4 s−1), and the activation enthalpies for this hydrogen exchange (14.7–19.4 kcal⋅mol−1). It was found that the previously noted correlations between the inhibition constants (Ki 170–1.2 μM) and the chemical shifts of His 57-Hδ1 (δ2,2-dimethylsilapentane-5-sulfonic acid 18.61–18.95) for this series of peptidyl trifluoromethyl ketones with chymotrypsin [Lin, J., Cassidy, C. S. & Frey, P. A. (1998) Biochemistry 37, 11940–11948] could be extended to include the fractionation factors, hydrogen exchange rates, and hydrogen exchange activation enthalpies. The results support the proposal of low barrier hydrogen bond-facilitated general base catalysis in the addition of Ser 195 to the peptidyl carbonyl group of substrates in the mechanism of chymotrypsin-catalyzed peptide hydrolysis. Trends in the enthalpies for hydrogen exchange and the fractionation factors are consistent with a strong, double-minimum or single-well potential hydrogen bond in the strongest complexes. The lifetimes of His 57-Hδ1, which is solvent shielded in these complexes, track the strength of the hydrogen bond. Because these lifetimes are orders of magnitude shorter than those of the complexes themselves, the enzyme must have a pathway for hydrogen exchange at this site that is independent of dissociation of the complexes.

Interresidue hydrogen bonding in a peptide nucleic acid.RNA heteroduplex.

Previous molecular mechanics calculations suggest that strands of peptide nucleic acids (PNAs) and complementary oligonucleotides form antiparallel duplexes stabilized by interresidue hydrogen bonds. In the computed structures, the amide carbonyl oxygen nearest the nucleobase (O7') forms an interresidue hydrogen bond with the backbone amide proton of the following residue, (n + 1)H1'. Of the 10 published two dimensional 1H NMR structures of a hexameric PNA.RNA heteroduplex. PNA(GAACTC).r(GAGUUC), 9 exhibit two to five potential interresidue hydrogen bonds. In our minimized average structure, created from the coordinates of these 10 NMR structures, three of the five possible interresidue hydrogen bond sites within the PNA backbone display the carbonyl oxygen (O7') and the amide proton (n + 1)H1' distances and N1'-H1'-(n - 1)O7' angles optimal for hydrogen bond formation. The finding of these interresidue hydrogen bonds supports the results of our previous molecular mechanics calculations.

A zinc finger directory for high-affinity DNA recognition

We have used two monovalent phage display libraries containing variants of the Zif268 DNA-binding domain to obtain families of zinc fingers that bind to alterations in the last 4 bp of the DNA sequence of the Zif268 consensus operator, GCG TGGGCG. Affinity selection was performed by altering the Zif268 operator three base pairs at a time, and simultaneously selecting for sets of 16 related DNA sequences. In this way, only four experiments were required to select for all possible 64 combinations of DNA triplet sequences. The results show that (i) for high-affinity DNA binding in the range observed for the Zif268 wild-type complex (Kd = 0.5–5 nM), finger 1 specifically requires the arginine at the carboxy terminus of its recognition helix that forms a bidentate hydrogen-bond with the guanine base (G) in the crystal structure of Zif268 complexed to its DNA operator sequence GCG TGG GCG; (ii) when the guanine base (G) is replaced by A, C, or T, a lower-affinity family (Kd ⩾ 50 nM) can be detected that shows an overall tendency to bind G-rich DNA; (iii) the residues at position 2 on the finger 2 recognition helix do not appear to interact strongly with the complementary 5′ base in the finger 1 binding site; and (iv) unexpected substitutions at the amino terminus of finger 1 can occasionally result in specificity for the 3′ base in the finger 1 binding site. A DNA recognition directory was constructed for high-affinity zinc fingers that recognize all three bases in a DNA triplet for seven sequences of the type GNN. Similar approaches may be applied to other zinc fingers to broaden the scope of the directory.

NMR scalar couplings across Watson–Crick base pair hydrogen bonds in DNA observed by transverse relaxation-optimized spectroscopy

This paper describes the NMR observation of 15N—15N and 1H—15N scalar couplings across the hydrogen bonds in Watson–Crick base pairs in a DNA duplex, hJNN and hJHN. These couplings represent new parameters of interest for both structural studies of DNA and theoretical investigations into the nature of the hydrogen bonds. Two dimensional [15N,1H]-transverse relaxation-optimized spectroscopy (TROSY) with a 15N-labeled 14-mer DNA duplex was used to measure hJNN, which is in the range 6–7 Hz, and the two-dimensional hJNN-correlation-[15N,1H]-TROSY experiment was used to correlate the chemical shifts of pairs of hydrogen bond-related 15N spins and to observe, for the first time, hJHN scalar couplings, with values in the range 2–3.6 Hz. TROSY-based studies of scalar couplings across hydrogen bonds should be applicable for large molecular sizes, including protein-bound nucleic acids.

Three-dimensional structure of human electron transfer flavoprotein to 2.1-Å resolution

Mammalian electron transfer flavoproteins (ETF) are heterodimers containing a single equivalent of flavin adenine dinucleotide (FAD). They function as electron shuttles between primary flavoprotein dehydrogenases involved in mitochondrial fatty acid and amino acid catabolism and the membrane-bound electron transfer flavoprotein ubiquinone oxidoreductase. The structure of human ETF solved to 2.1-Å resolution reveals that the ETF molecule is comprised of three distinct domains: two domains are contributed by the α subunit and the third domain is made up entirely by the β subunit. The N-terminal portion of the α subunit and the majority of the β subunit have identical polypeptide folds, in the absence of any sequence homology. FAD lies in a cleft between the two subunits, with most of the FAD molecule residing in the C-terminal portion of the α subunit. Alignment of all the known sequences for the ETF α subunits together with the putative FixB gene product shows that the residues directly involved in FAD binding are conserved. A hydrogen bond is formed between the N5 of the FAD isoalloxazine ring and the hydroxyl side chain of αT266, suggesting why the pathogenic mutation, αT266M, affects ETF activity in patients with glutaric acidemia type II. Hydrogen bonds between the 4′-hydroxyl of the ribityl chain of FAD and N1 of the isoalloxazine ring, and between αH286 and the C2-carbonyl oxygen of the isoalloxazine ring, may play a role in the stabilization of the anionic semiquinone. With the known structure of medium chain acyl-CoA dehydrogenase, we hypothesize a possible structure for docking the two proteins.

NMR reveals hydrogen bonds between oxygen and distal histidines in oxyhemoglobin

Compared with free heme, the proteins hemoglobin (Hb) and myoglobin (Mb) exhibit greatly enhanced affinity for oxygen relative to carbon monoxide. This physiologically vital property has been attributed to either steric hindrance of CO or stabilization of O2 binding by a hydrogen bond with the distal histidine. We report here the first direct evidence of such a hydrogen bond in both α- and β-chains of oxyhemoglobin, as revealed by heteronuclear NMR spectra of chain-selectively labeled samples. Using these spectra, we have assigned the imidazole ring 1H and 15N chemical shifts of the proximal and distal histidines in both carbonmonoxy- and oxy-Hb. Because of their proximity to the heme, these chemical shifts are extremely sensitive to the heme pocket conformation. Comparison of the measured chemical shifts with values predicted from x-ray structures suggests differences between the solution and crystal structures of oxy-Hb. The chemical shift discrepancies could be accounted for by very small displacements of the proximal and distal histidines. This suggests that NMR could be used to obtain very high-resolution heme pocket structures of Hb, Mb, and other heme proteins.

Tight Asp-85–Thr-89 association during the pump switch of bacteriorhodopsin

Unidirectional proton transport in bacteriorhodopsin is enforced by the switching machinery of the active site. Threonine 89 is located in this region, with its O—H group forming a hydrogen bond with Asp-85, the acceptor for proton transfer from the Schiff base of the retinal chromophore. Previous IR spectroscopy of [3-18O]threonine-labeled bacteriorhodopsin showed that the hydrogen bond of the O—D group of Thr-89 in D2O is strengthened in the K photocycle intermediate. Here, we show that the strength and orientation of this hydrogen bond remains unchanged in the L intermediate and through the M intermediate. Furthermore, a strong interaction between Asp-85 and the O—H (O—D) group of Thr-89 in M is indicated by a shift in the C⩵O stretching vibration of the former because of 18O substitution in the latter. Thus, the strong hydrogen bond between Asp-85 and Thr-89 in K persists through M, contrary to structural models based on x-ray crystallography of the photocycle intermediates. We propose that, upon photoisomerization of the chromophore, Thr-89 forms a tight, persistent complex with one of the side-chain oxygens of Asp-85 and is thereby precluded from participating in the switching process. On the other hand, the loss of hydrogen bonding at the other oxygen of Asp-85 in M may be related to the switching event.

Polar residues drive association of polyleucine transmembrane helices

Although many polar residues are directly involved in transmembrane protein functions, the extent to which they contribute to more general structural features is still unclear. Previous studies have demonstrated that asparagine residues can drive transmembrane helix association through interhelical hydrogen bonding [Choma, C., Gratkowski, H., Lear, J. D. & DeGrado, W. F. (2000) Nat. Struct. Biol. 7, 161–166; and Zhou, F. X., Cocco, M. J., Russ, W. P., Brunger, A. T. & Engelman, D. M. (2000) Nat. Struct. Biol. 7, 154–160]. We have studied the ability of other polar residues to promote helix association in detergent micelles and in biological membranes. Our results show that polyleucine sequences with Asn, Asp, Gln, Glu, and His, residues capable of being simultaneously hydrogen bond donors and acceptors, form homo- or heterooligomers. In contrast, polyleucine sequences with Ser, Thr, and Tyr do not associate more than the polyleucine sequence alone. The results therefore provide experimental evidence that interactions between polar residues in the helices of transmembrane proteins may serve to provide structural stability and oligomerization specificity. Furthermore, such interactions can allow structural flexibility required for the function of some membrane proteins.

Barriers to protein folding: formation of buried polar interactions is a slow step in acquisition of structure.

In the MYL mutant of the Arc repressor dimer, sets of partially buried salt-bridge and hydrogen-bond interactions mediated by Arg-31, Glu-36, and Arg-40 in each subunit are replaced by hydrophobic interactions between Met-31, Tyr-36, and Leu-40. The MYL refolding/dimerization reaction differs from that of wild type in being 10- to 1250-fold faster, having an earlier transition state, and depending upon viscosity but not ionic strength. Formation of the wild-type salt bridges in a hydrophobic environment clearly imposes a kinetic barrier to folding, which can be lowered by high salt concentrations. The changes in the position of the transition state and viscosity dependence can be explained if denatured monomers interact to form a partially folded dimeric intermediate, which then continues folding to form the native dimer. The second step is postulated to be rate limiting for wild type. Replacing the salt bridge with hydrophobic interactions lowers this barrier for MYL. This makes the first kinetic barrier rate limiting for MYL refolding and creates a downhill free-energy landscape in which most molecules which reach the intermediate state continue to form native dimers.