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.

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.

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.

The crystal structure of the Rev binding element of HIV-1 reveals novel base pairing and conformational variability

The crystal and molecular structure of an RNA duplex corresponding to the high affinity Rev protein binding element (RBE) has been determined at 2.1-Å resolution. Four unique duplexes are present in the crystal, comprising two structural variants. In each duplex, the RNA double helix consists of an annealed 12-mer and 14-mer that form an asymmetric internal loop consisting of G-G and G-A noncanonical base pairs and a flipped-out uridine. The 12-mer strand has an A-form conformation, whereas the 14-mer strand is distorted to accommodate the bulges and noncanonical base pairing. In contrast to the NMR model of the unbound RBE, an asymmetric G-G pair with N2-N7 and N1-O6 hydrogen bonding, is formed in each helix. The G-A base pairing agrees with the NMR structure in one structural variant, but forms a novel water-mediated pair in the other. A backbone flip and reorientation of the G-G base pair is required to assume the RBE conformation present in the NMR model of the complex between the RBE and the Rev peptide.

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.

An important 2′-OH group for an RNA–protein interaction

We have investigated the role of 2′-OH groups in the specific interaction between the acceptor stem of Escherichia coli tRNACys and cysteine-tRNA synthetase. This interaction provides for the high aminoacylation specificity observed for cysteine-tRNA synthetase. A synthetic RNA microhelix that recapitulates the sequence of the acceptor stem was used as a substrate and variants containing systematic replacement of the 2′-OH by 2′-deoxy or 2′-O-methyl groups were tested. Except for position U73, all substitutions had little effect on aminoacylation. Interestingly, the deoxy substitution at position U73 had no effect on aminoacylation, but the 2′-O-methyl substitution decreased aminoacylation by 10-fold and addition of the even bulkier 2′-O-propyl group decreased aminoacylation by another 2-fold. The lack of an effect by the deoxy substitution suggests that the hydrogen bonding potential of the 2′-OH at position U73 is unimportant for aminoacylation. The decrease in activity upon alkyl substitution suggests that the 2′-OH group instead provides a monitor of the steric environment during the RNA–synthetase interaction. The steric role was confirmed in the context of a reconstituted tRNA and is consistent with the observation that the U73 base is the single most important determinant for aminoacylation and therefore is a site that is likely to be in close contact with cysteine-tRNA synthetase. A steric role is supported by an NMR-based structural model of the acceptor stem, together with biochemical studies of a closely related microhelix. This role suggests that the U73 binding site for cysteine-tRNA synthetase is sterically optimized to accommodate a 2′-OH group in the backbone, but that the hydroxyl group itself is not involved in specific hydrogen bonding interactions.

Evidence for DNA phosphate backbone alkylation and cleavage by pyrrolo[1,2-a]benzimidazoles: small molecules capable of causing base-pair-specific phosphodiester bond hydrolysis.

This report presents evidence that a reduced pyrrolo[1,2-a]benzimidazole (PBI) cleaves DNA as a result of phosphate alkylation followed by hydrolysis of the resulting phosphate triester. The base-pair specificity of the phosphate alkylation results from Hoogsteen-type hydrogen bonding of the reduced PBI in the major groove at only A.T and G.C base pairs. Alkylated phosphates were detected by 31P NMR and the cleavage products were detected by 1H NMR and HPLC. Evidence is also presented that a reduced PBI interacts with DNA in the major groove rather than in the minor groove or by intercalation.

Calicheamicin-DNA complexes: warhead alignment and saccharide recognition of the minor groove.

The solution structures of calicheamicin gamma 1I, its cycloaromatized analog (calicheamicin epsilon), and its aryl tetrasaccharide complexed to a common DNA hairpin duplex have been determined by NMR and distance-refined molecular dynamics computations. Sequence specificity is associated with carbohydrate-DNA recognition that places the aryl tetrasaccharide component of all three ligands in similar orientations in the minor groove at the d(T-C-C-T).d(A-G-G-A) segment. The complementary fit of the ligands and the DNA minor groove binding site creates numerous van der Waals contacts as well as hydrogen bonding interactions. Notable are the iodine and sulfur atoms of calicheamicin that hydrogen bond with the exposed amino proton of the 5'- and 3'-guanines, respectively, of the d(A-G-G-A) segment. The sequence-specific carbohydrate binding orients the enediyne aglycone of calicheamicin gamma 1I such that its C3 and C6 proradical centers are adjacent to the cleavage sites. While the enediyne aglycone of calicheamicin gamma 1I is tilted relative to the helix axis and spans the minor groove, the cycloaromatized aglycone is aligned approximately parallel to the helix axis in the respective complexes. Specific localized conformational perturbations in the DNA have been identified from imino proton complexation shifts and changes in specific sugar pucker patterns on complex formation. The helical parameters for the carbohydrate binding site are comparable with corresponding values in B-DNA fibers while a widening of the groove is observed at the adjacent aglycone binding site.

Crystal and molecular structure of paclitaxel (taxol).

Paclitaxel (formerly called taxol), an important anticancer drug, inhibits cell replication by binding to and stabilizing microtubule polymers. As drug-receptor interactions are governed by the three-dimensional stereochemistries of both participants, we have determined the crystal structure of paclitaxel to identify its conformational preferences that may be related to biological activity. The monoclinic crystals contain two independent paclitaxel molecules in the asymmetric unit plus several water and dioxane solvent molecules. Taxane ring conformation is very similar in both paclitaxel molecules and is similar to the taxane ring conformation found in the crystal structure of the paclitaxel analogue docetaxel (formerly called taxotere). The two paclitaxel molecules have carbon-13 side-chain conformations that differ from each other and from that of the corresponding side chain in the docetaxel crystal structure. The carbon-13 side-chain conformation of one paclitaxel molecule is similar to what was proposed from NMR studies done in polar solvents, while that of the other paclitaxel molecule is different and hitherto unobserved. The paclitaxel molecules interact with each other and with solvent atoms through an extensive network of hydrogen bonds. Analysis of the hydrogen-bonding network together with structure-activity studies may suggest which atoms of paclitaxel are important for binding to microtubule receptors.

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.