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.

Recognition by viral and cellular DNA polymerases of nucleosides bearing bases with nonstandard hydrogen bonding patterns.

The ability of DNA polymerases (pols) to catalyze the template-directed synthesis of duplex oligonucleotides containing a nonstandard Watson-Crick base pair between a nucleotide bearing a 5-(2,4-diaminopyrimidine) heterocycle (d kappa) and a nucleotide bearing either deoxyxanthosine (dX) or N1-methyloxoformycin B (pi) has been investigated. The kappa-X and kappa-pi base pairs are jointed by a hydrogen bonding pattern different from and exclusive of those joining the AT and GC base pairs. Reverse transcriptase from human immunodeficiency virus type 1 (HIV-1) incorporates dXTP into an oligonucleotide opposite d kappa in a template with good fidelity. With lower efficiency and fidelity, HIV-1 reverse transcriptase also incorporates d kappa TP opposite dX in the template. With d pi in the template, no incorporation of d kappa TP was observed with HIV reverse transcriptase. The Klenow fragment of DNA pol I from Escherichia coli does not incorporate d kappa TP opposite dX in a template but does incorporate dXTP opposite d kappa. Bovine DNA pols alpha, beta, and epsilon accept neither dXTP opposite d kappa nor d kappa TP opposite d pi. DNA pols alpha and epsilon (but not beta) incorporate d kappa TP opposite dX in a template but discontinue elongation after incorporating a single additional base. These results are discussed in light of the crystal structure for pol beta and general considerations of how polymerases must interact with an incoming base pair to faithfully copy genetic information.

Direct measurement of hydrogen bonding in DNA nucleotide bases by atomic force microscopy.

We have used self-assembled purines and pyrimidines on planar gold surfaces and on gold-coated atomic force microscope (AFM) tips to directly probe intermolecular hydrogen bonds. Electron spectroscopy for chemical analysis (ESCA) and thermal programmed desorption (TPD) measurements of the molecular layers suggested monolayer coverage and a desorption energy of about 25 kcal/mol. Experiments were performed under water, with all four DNA bases immobilized on AFM tips and flat surfaces. Directional hydrogen-bonding interaction between the tip molecules and the surface molecules could be measured only when opposite base-pair coatings were used. The directional interactions were inhibited by excess nucleotide base in solution. Nondirectional van der Waals forces were present in all other cases. Forces as low as two interacting base pairs have been measured. With coated AFM tips, surface chemistry-sensitive recognition atomic force microscopy can be performed.

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.

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.

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.

Peptide hybrids containing α- and β-amino acids: Structure of a decapeptide β-hairpin with two facing β-phenylalanine residues

A β-hairpin conformation has been characterized in crystals of the decapeptide t-butoxycarbonyl-Leu-Val-βPhe-Val-DPro-Gly-Leu-βPhe-Val-Val-methyl ester [βPhe; (S)-β3 homophenylalanine] by x-ray diffraction. The polypeptide chain reversal is nucleated by the centrally positioned DPro-Gly segment, which adopts a type-I′ β-turn conformation. Four intramolecular cross-strand hydrogen bonds stabilize the peptide fold. The βPhe(3) and βPhe(8) residues occupy facing positions on the hairpin, with the side chains projecting on opposite faces of the β-sheet. At the site of insertion of β-residues, the polarity of the peptide units along each strand reverses, as compared with the α-peptide segments. In this analog, a small segment of a polar sheet is observed, where adjacent CO and NH groups line up in opposite directions in each strand. In the crystal, an extended β-sheet is formed by hydrogen bonding between strands of antiparallel pairs of β-hairpins. The crystallographic parameters for C65H102N10O13⋅ 3H2O are: space group P212121; a = 19.059(8) Å, b = 19.470(2) Å, c = 21.077(2) Å; Z = 4; agreement factor R1 = 9.12% for 3,984 data observed >4σ(F) and a resolution of 0.90 Å.

Inter-strand C-H...O hydrogen bonds stabilizing four-stranded intercalated molecules: stereoelectronic effects of O4' in cytosine-rich DNA.

DNA fragments with stretches of cytosine residues can fold into four-stranded structures in which two parallel duplexes, held together by hemiprotonated cytosine.cytosine+ (C.C+) base pairs, intercalate into each other with opposite polarity. The structural details of this intercalated DNA quadruplex have been assessed by solution NMR and single crystal x-ray diffraction studies of cytosine-rich sequences, including those present in metazoan telomeres. A conserved feature of these structures is the absence of stabilizing stacking interactions between the aromatic ring systems of adjacent C.C+ base pairs from intercalated duplexes. Effective stacking involves only the exocyclic keto groups and amino groups of the cytidine bases. The apparent absence of stability provided by stacking interactions between the bases in this intercalated DNA has prompted us to examine the available structures in detail, in particular with regard to unusual features that could compensate for the lack of base stacking. In addition to base-on-deoxyribose stacking and intra-cytidine C-H...O hydrogen bonds, this analysis reveals the presence of a hitherto unobserved, systematic intermolecular C-H...O hydrogen bonding network between the deoxyribose sugar moieties of antiparallel backbones in the four-stranded molecule.

Local densities orthogonal to beta-sheet amide planes: patterns of packing in globular proteins.

We have investigated the efficiency of packing by calculating intramolecular packing density above and below peptide planes of internal beta-pleated sheet residues in five globular proteins. The orientation of interest was chosen to allow study of regions that are approximately perpendicular to the faces of beta-pleated sheets. In these locations, nonbonded van der Waals packing interactions predominate over hydrogen bonding and solvent interactions. We observed considerable variability in packing densities within these regions, confirming that the interior packing of a protein does not result in uniform occupation of the available space. Patterns of fluctuation in packing density suggest that the regular backbone-to-backbone network of hydrogen bonds is not likely to be interrupted to maximize van der Waals interactions. However, high-density packing tends to occur toward the ends of beta-structure strands where hydrogen bonds are more likely to involve nonpolar side-chain groups or solvent molecules. These features result in internal protein folding with a central low-density core surrounded by a higher-density subsurface shell, consistent with our previous calculations regarding overall protein packing density.

Designing conditions for in vitro formation of amyloid protofilaments and fibrils

We have been able to convert a small α/β protein, acylphosphatase, from its soluble and native form into insoluble amyloid fibrils of the type observed in a range of pathological conditions. This was achieved by allowing slow growth in a solution containing moderate concentrations of trifluoroethanol. When analyzed with electron microscopy, the protein aggregate present in the sample after long incubation times consisted of extended, unbranched filaments of 30–50 Å in width that assemble subsequently into higher order structures. This fibrillar material possesses extensive β-sheet structure as revealed by far-UV CD and IR spectroscopy. Furthermore, the fibrils exhibit Congo red birefringence, increased fluorescence with thioflavine T and cause a red-shift of the Congo red absorption spectrum. All of these characteristics are typical of amyloid fibrils. The results indicate that formation of amyloid occurs when the native fold of a protein is destabilized under conditions in which noncovalent interactions, and in particular hydrogen bonding, within the polypeptide chain remain favorable. We suggest that amyloid formation is not restricted to a small number of protein sequences but is a property common to many, if not all, natural polypeptide chains under appropriate conditions.

Role of a critical water in scytalone dehydratase-catalyzed reaction

Scytalone dehydratase (EC 4.2.1.94) catalyzes the dehydration of two important intermediates in the biosynthesis of melanin, and it functions without metal ions or any cofactors. Using molecular orbital theory, we have examined the role of a critical water molecule in the mechanism of scytalone dehydratase. The water, together with an internal hydrogen bonding, contributes significantly to the stabilization of the transition state (or the enolate intermediate). The role of two active site tyrosines (Tyr-50 and Tyr-30) is (i) to hold the critical water in place so that it may stabilize the transition state without much structural rearrangement during the catalytic reaction, and (ii) to polarize the water, making it a better general acid. The stereochemistry of the scytalone dehydratase-catalyzed dehydration is also discussed.

Low frequency vibrational modes in proteins: Changes induced by point-mutations in the protein-cofactor matrix of bacterial reaction centers

As a step toward understanding their functional role, the low frequency vibrational motions (<300 cm−1) that are coupled to optical excitation of the primary donor bacteriochlorophyll cofactors in the reaction center from Rhodobacter sphaeroides were investigated. The pattern of hydrogen-bonding interaction between these bacteriochlorophylls and the surrounding protein was altered in several ways by mutation of single amino acids. The spectrum of low frequency vibrational modes identified by femtosecond coherence spectroscopy varied strongly between the different reaction center complexes, including between different mutants where the pattern of hydrogen bonds was the same. It is argued that these variations are primarily due to changes in the nature of the individual modes, rather than to changes in the charge distribution in the electronic states involved in the optical excitation. Pronounced effects of point mutations on the low frequency vibrational modes active in a protein-cofactor system have not been reported previously. The changes in frequency observed indicate a strong involvement of the protein in these nuclear motions and demonstrate that the protein matrix can increase or decrease the fluctuations of the cofactor along specific directions.

Non-enzymatic and enzymatic hydrolysis of alkyl halides: A haloalkane dehalogenation enzyme evolved to stabilize the gas-phase transition state of an SN2 displacement reaction

The semiempirical PM3 method, calibrated against ab initio HF/6–31+G(d) theory, has been used to elucidate the reaction of 1,2-dichloroethane (DCE) with the carboxylate of Asp-124 at the active site of haloalkane dehalogenase of Xanthobacter autothropicus. Asp-124 and 13 other amino acid side chains that make up the active site cavity (Glu-56, Trp-125, Phe-128, Phe-172, Trp-175, Leu-179, Val-219, Phe-222, Pro-223, Val-226, Leu-262, Leu-263, and His-289) were included in the calculations. The three most significant observations of the present study are that: (i) the DCE substrate and Asp-124 carboxylate, in the reactive ES complex, are present as an ion-molecule complex with a structure similar to that seen in the gas-phase reaction of AcO− with DCE; (ii) the structures of the transition states in the gas-phase and enzymatic reaction are much the same where the structure formed at the active site is somewhat exploded; and (iii) the enthalpies in going from ground states to transition states in the enzymatic and gas-phase reactions differ by only a couple kcal/mol. The dehalogenase derives its catalytic power from: (i) bringing the electrophile and nucleophile together in a low-dielectric environment in an orientation that allows the reaction to occur without much structural reorganization; (ii) desolvation; and (iii) stabilizing the leaving chloride anion by Trp-125 and Trp-175 through hydrogen bonding.