A link between protein structure and enzyme catalyzed hydrogen tunneling

We present evidence that the size of an active site side chain may modulate the degree of hydrogen tunneling in an enzyme-catalyzed reaction. Primary and secondary kH/kT and kD/kT kinetic isotope effects have been measured for the oxidation of benzyl alcohol catalyzed by horse liver alcohol dehydrogenase at 25°C. As reported in earlier studies, the relationship between secondary kH/kT and kD/kT isotope effects provides a sensitive probe for deviations from classical behavior. In the present work, catalytic efficiency and the extent of hydrogen tunneling have been correlated for the alcohol dehydrogenase-catalyzed hydride transfer among a group of site-directed mutants at position 203. Val-203 interacts with the opposite face of the cofactor NAD+ from the alcohol substrate. The reduction in size of this residue is correlated with diminished tunneling and a two orders of magnitude decrease in catalytic efficiency. Comparison of the x-ray crystal structures of a ternary complex of a high-tunneling (Phe-93 → Trp) and a low-tunneling (Val-203 → Ala) mutant provides a structural basis for the observed effects, demonstrating an increase in the hydrogen transfer distance for the low-tunneling mutant. The Val-203 → Ala ternary complex crystal structure also shows a hyperclosed interdomain geometry relative to the wild-type and the Phe-93 → Trp mutant ternary complex structures. This demonstrates a flexibility in interdomain movement that could potentially narrow the distance between the donor and acceptor carbons in the native enzyme and may enhance the role of tunneling in the hydride transfer reaction.

Structural basis for efficient phosphorylation of 3′-azidothymidine monophosphate by Escherichia coli thymidylate kinase

The crystal structures of Escherichia coli thymidylate kinase (TmpK) in complex with P1-(5′-adenosyl)-P5-(5′-thymidyl)pentaphosphate and P1-(5′-adenosyl)P5-[5′-(3′-azido-3′-deoxythymidine)] pentaphosphate have been solved to 2.0-Å and 2.2-Å resolution, respectively. The overall structure of the bacterial TmpK is very similar to that of yeast TmpK. In contrast to the human and yeast TmpKs, which phosphorylate 3′-azido-3′-deoxythymidine 5′-monophosphate (AZT-MP) at a 200-fold reduced turnover number (kcat) in comparison to the physiological substrate dTMP, reduction of kcat is only 2-fold for the bacterial enzyme. The different kinetic properties toward AZT-MP between the eukaryotic TmpKs and E. coli TmpK can be rationalized by the different ways in which these enzymes stabilize the presumed transition state and the different manner in which a carboxylic acid side chain in the P loop interacts with the deoxyribose of the monophosphate. Yeast TmpK interacts with the 3′-hydroxyl of dTMP through Asp-14 of the P loop in a bidentate manner: binding of AZT-MP results in a shift of the P loop to accommodate the larger substituent. In E. coli TmpK, the corresponding residue is Glu-12, and it interacts in a side-on fashion with the 3′-hydroxyl of dTMP. This different mode of interaction between the P loop carboxylic acid with the 3′ substituent of the monophosphate deoxyribose allows the accommodation of an azido group in the case of the E. coli enzyme without significant P loop movement. In addition, although the yeast enzyme uses Arg-15 (a glycine in E. coli) to stabilize the transition state, E. coli seems to use Arg-153 from a region termed Lid instead. Thus, the binding of AZT-MP to the yeast TmpK results in the shift of a catalytic residue, which is not the case for the bacterial kinase.

Structural origins of the exonuclease resistance of a zwitterionic RNA

Nuclease resistance and RNA affinity are key criteria in the search for optimal antisense nucleic acid modifications, but the origins of the various levels of resistance to nuclease degradation conferred by chemical modification of DNA and RNA are currently not understood. The 2′-O-aminopropyl (AP)-RNA modification displays the highest nuclease resistance among all phosphodiester-based analogues and its RNA binding affinity surpasses that of phosphorothioate DNA by 1°C per modified residue. We found that oligodeoxynucleotides containing AP-RNA residues at their 3′ ends competitively inhibit the degradation of single-stranded DNA by the Escherichia coli Klenow fragment (KF) 3′-5′ exonuclease and snake venom phosphodiesterase. To shed light on the origins of nuclease resistance brought about by the AP modification, we determined the crystal structure of an A-form DNA duplex with AP-RNA modifications at 1.6-Å resolution. In addition, the crystal structures of complexes between short DNA fragments carrying AP-RNA modifications and wild-type KF were determined at resolutions between 2.2 and 3.0 Å and compared with the structure of the complex between oligo(dT) and the D355A/E357A KF mutant. The structural models suggest that interference of the positively charged 2′-O-substituent with the metal ion binding site B of the exonuclease allows AP-RNA to effectively slow down degradation.

Ordered water molecules as key allosteric mediators in a cooperative dimeric hemoglobin

One of the most remarkable structural aspects of Scapharca dimeric hemoglobin is the disruption of a very well-ordered water cluster at the subunit interface upon ligand binding. We have explored the role of these crystallographically observed water molecules by site-directed mutagenesis and osmotic stress techniques. The isosteric mutation of Thr-72 → Val in the interface increases oxygen affinity more than 40-fold with a surprising enhancement of cooperativity. The only significant structural effect of this mutation is to destabilize two ordered water molecules in the deoxy interface. Wild-type Scapharca hemoglobin is strongly sensitive to osmotic conditions. Upon addition of glycerol, striking changes in Raman spectrum of the deoxy form are observed that indicate a transition toward the liganded form. Increased osmotic pressure, which lowers the oxygen affinity in human hemoglobin, raises the oxygen affinity of Scapharca hemoglobin regardless of whether the solute is glycerol, glucose, or sucrose. Analysis of these results provides an estimate of six water molecules lost upon oxygen binding to the dimer, in good agreement with eight predicted from crystal structures. These experiments suggest that the observed cluster of interfacial water molecules plays a crucial role in communication between subunits.

Solid-state synthesis and mechanical unfolding of polymers of T4 lysozyme

Recent advances in single molecule manipulation methods offer a novel approach to investigating the protein folding problem. These studies usually are done on molecules that are naturally organized as linear arrays of globular domains. To extend these techniques to study proteins that normally exist as monomers, we have developed a method of synthesizing polymers of protein molecules in the solid state. By introducing cysteines at locations where bacteriophage T4 lysozyme molecules contact each other in a crystal and taking advantage of the alignment provided by the lattice, we have obtained polymers of defined polarity up to 25 molecules long that retain enzymatic activity. These polymers then were manipulated mechanically by using a modified scanning force microscope to characterize the force-induced reversible unfolding of the individual lysozyme molecules. This approach should be general and adaptable to many other proteins with known crystal structures. For T4 lysozyme, the force required to unfold the monomers was 64 ± 16 pN at the pulling speed used. Refolding occurred within 1 sec of relaxation with an efficiency close to 100%. Analysis of the force versus extension curves suggests that the mechanical unfolding transition follows a two-state model. The unfolding forces determined in 1 M guanidine hydrochloride indicate that in these conditions the activation barrier for unfolding is reduced by 2 kcal/mol.

Unexpected crucial role of residue 225 in serine proteases

Residue 225 in serine proteases of the chymotrypsin family is Pro or Tyr in more than 95% of nearly 300 available sequences. Proteases with Y225 (like some blood coagulation and complement factors) are almost exclusively found in vertebrates, whereas proteases with P225 (like degradative enzymes) are present from bacteria to human. Saturation mutagenesis of Y225 in thrombin shows that residue 225 affects ligand recognition up to 60,000-fold. With the exception of Tyr and Phe, all residues are associated with comparable or greatly reduced catalytic activity relative to Pro. The crystal structures of three mutants that differ widely in catalytic activity (Y225F, Y225P, and Y225I) show that although residue 225 makes no contact with substrate, it drastically influences the shape of the water channel around the primary specificity site. The activity profiles obtained for thrombin also suggest that the conversion of Pro to Tyr or Phe documented in the vertebrates occurred through Ser and was driven by a significant gain (up to 50-fold) in catalytic activity. In fact, Ser and Phe are documented in 4% of serine proteases, which together with Pro and Tyr account for almost the entire distribution of residues at position 225. The unexpected crucial role of residue 225 in serine proteases explains the evolutionary selection of residues at this position and shows that the structural determinants of protease activity and specificity are more complex than currently believed. These findings have broad implications in the rational design of enzymes with enhanced catalytic properties.

Fibroblast growth factors 1 and 2 are distinct in oligomerization in the presence of heparin-like glycosaminoglycans

Fibroblast growth factor (FGF) 1 and FGF-2 are prototypic members of the FGF family, which to date comprises at least 18 members. Surprisingly, even though FGF-1 and FGF-2 share more than 80% sequence similarity and an identical structural fold, these two growth factors are biologically very different. FGF-1 and FGF-2 differ in their ability to bind isoforms of the FGF receptor family as well as the heparin-like glycosaminoglycan (HLGAG) component of proteoglycans on the cell surface to initiate signaling in different cell types. Herein, we provide evidence for one mechanism by which these two proteins could differ biologically. Previously, it has been noted that FGF-1 and FGF-2 can oligomerize in the presence of HLGAGs. Therefore, we investigated whether FGF-1 and FGF-2 oligomerize by the same mechanism or by a different one. Through a combination of matrix-assisted laser desorption ionization mass spectrometry and chemical crosslinking, we show here that, under identical conditions, FGF-1 and FGF-2 differ in the degree and kind of oligomerization. Furthermore, an extensive analysis of FGF-1 and FGF-2 uncomplexed and HLGAG complexed crystal structures enables us to readily explain why FGF-2 forms sequential oligomers whereas FGF-1 forms only dimers. FGF-2, which possesses an interface capable of protein association, forms a translationally related oligomer, whereas FGF-1, which does not have this interface, forms only a symmetrically related dimer. Taken together, these data show that FGF-1 and FGF-2, despite their sequence homology, differ in their mechanism of oligomerization.

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.

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.

Inhibition of HIV type 1 infectivity by constrained α-helical peptides: Implications for the viral fusion mechanism

Linear peptides derived from the membrane proximal region of the gp41 ectodomain are effective inhibitors of HIV type 1 (HIV-1)-mediated fusion events. These inhibitory peptides lack structure in solution, rendering mechanistic interpretation of their activity difficult. Using structurally constrained analogs of these molecules, we demonstrate that the peptides inhibit infectivity by adopting a helical conformation. Moreover, we show that a specific face of the helix must be exposed to block viral infectivity. Recent crystal structures show that the region of gp41 corresponding to the inhibitory peptides is helical and uses the analogous face to pack against a groove formed by an N-terminal coiled-coil trimer. Our results provide a direct link between the inhibition of HIV-1 infectivity by these peptides and the x-ray structures, and suggest that the conformation of gp41 observed by crystallography represents the fusogenic state. Other agents that block HIV-1 infectivity by binding to this groove may hold promise for the treatment of AIDS.

Mutagenesis of the regulatory domain of phenylalanine hydroxylase

The regulatory domain of phenylalanine hydroxylase (PAH, EC 1.14.16.1) consists of more than 100 amino acids at the N terminus, the removal of which significantly activates the enzyme. To study the regulatory properties controlled by the N terminus, a series of truncations and site-specific mutations were made in this region of rat PAH. These enzymes were expressed highly in Escherichia coli and purified through a pterin-conjugated Sepharose affinity column. The removal of the first 26 amino acids of the N terminus increased the activity by about 20-fold, but removal of the first 15 amino acids increased the activity by only 2-fold. Replacing serine-29 of rat PAH with cysteine from the same site of human PAH increased the activity by more than 4-fold. Mutation of serine to other amino acids with varying side chains: alanine, methionine, leucine, aspartic acid, asparagine, and arginine also resulted in significant activation, indicating a serine-specific inhibitory effect. But these site-specific mutants showed 30–40% lower activity when assayed with 6-methyl-5,6,7,8-tetrahydropterin. Stimulation of hydroxylase activity by preincubation of the enzyme with phenylalanine was inversely proportional to the activation state of all these mutants. Combined with recent crystal structures of PAH [Kobe, B. et al. (1999) Nat. Struct. Biol. 6, 442–448; and Erlandsen, H., Bjorgo, E., Flatmark, T. & Stevens, R. C. (2000) Biochemistry 39, 2208–2217], these data suggest that residues 16–26 have a controlling regulatory effect on the activity by interaction with the dihydroxypropyl side chain of (6R)-5,6,7,8-tetrahydrobiopterin. The serine/cysteine switch explains the difference in regulatory properties between human and rat PAH. The N terminus as a whole is important for maintaining rat PAH in an optimum catalytic conformation.

One ring or two? Determination of ring number in carotenoids by lycopene ɛ-cyclases

Carotenoids in the photosynthetic membranes of plants typically contain two β-rings (e.g., β-carotene and zeaxanthin) or one ɛ- and one β-ring (e.g., lutein). Carotenoids with two ɛ-rings are uncommon. We reported earlier that the Arabidopsis thaliana lycopene ɛ-cyclase (LCYe) adds one ɛ-ring to the symmetrical linear substrate lycopene, whereas the structurally related lycopene β-cyclase (LCYb) adds two β-rings. Here we describe a cDNA encoding LCYe in romaine lettuce (Lactuca sativa var. romaine), one of the few plant species known to accumulate substantial quantities of a carotenoid with two ɛ-rings: lactucaxanthin. The product of the lettuce cDNA, similar in sequence to the Arabidopsis LCYe (77% amino acid identity), efficiently converted lycopene into the bicyclic ɛ-carotene in a heterologous Escherichia coli system. Regions of the lettuce and Arabidopsis ɛ-cyclases involved in the determination of ring number were mapped by analysis of chimeric ɛ-cyclases constructed by using an inverse PCR approach. A single amino acid was found to act as a molecular switch: lettuce LCYe mutant H457L added only one ɛ-ring to lycopene, whereas the complementary Arabidopsis LCYe mutant, L448H, added two ɛ-rings. An R residue in this position also yields a bi-ɛ-cyclase for both the lettuce and Arabidopsis enzymes. Construction and analysis of chimera of related enzymes with differing catalytic activities provide an informative approach that may be of particular utility for studying membrane-associated enzymes that cannot easily be crystallized or modeled to existing crystal structures.

Solution structure of the A loop of 23S ribosomal RNA

The A loop is an essential RNA component of the ribosome peptidyltransferase center that directly interacts with aminoacyl (A)-site tRNA. The A loop is highly conserved and contains a ubiquitous 2′-O-methyl ribose modification at position U2552. Here, we present the solution structure of a modified and unmodified A-loop RNA to define both the A-loop fold and the structural impact of the U2552 modification. Solution data reveal that the A-loop RNA has a compact structure that includes a noncanonical base pair between C2556 and U2552. NMR evidence is presented that the N3 position of C2556 has a shifted pKa and that protonation at C2556-N3 changes the C-U pair geometry. Our data indicate that U2552 methylation modifies the A-loop fold, in particular the dynamics and position of residues C2556 and U2555. We compare our structural data with the structure of the A loop observed in a recent 50S crystal structure [Ban, N., Nissen, P., Hansen, J., Moore, P. B. & Steitz, T. A. (2000) Science 289, 905–920; Nissen, P., Hansen, J., Ban, N., Moore, P. B. & Steitz, T. A. (2000) Science 289, 920–930]. The solution and crystal structures of the A loop are dramatically different, suggesting that a structural rearrangement of the A loop must occur on docking into the peptidyltransferase center. Possible roles of this docking event, the shifted pKa of C2556 and the U2552 2′-O-methylation in the mechanism of translation, are discussed.

Structure-assisted design of mechanism-based irreversible inhibitors of human rhinovirus 3C protease with potent antiviral activity against multiple rhinovirus serotypes

Human rhinoviruses, the most important etiologic agents of the common cold, are messenger-active single-stranded monocistronic RNA viruses that have evolved a highly complex cascade of proteolytic processing events to control viral gene expression and replication. Most maturation cleavages within the precursor polyprotein are mediated by rhinovirus 3C protease (or its immediate precursor, 3CD), a cysteine protease with a trypsin-like polypeptide fold. High-resolution crystal structures of the enzyme from three viral serotypes have been used for the design and elaboration of 3C protease inhibitors representing different structural and chemical classes. Inhibitors having α,β-unsaturated carbonyl groups combined with peptidyl-binding elements specific for 3C protease undergo a Michael reaction mediated by nucleophilic addition of the enzyme’s catalytic Cys-147, resulting in covalent-bond formation and irreversible inactivation of the viral protease. Direct inhibition of 3C proteolytic activity in virally infected cells treated with these compounds can be inferred from dose-dependent accumulations of viral precursor polyproteins as determined by SDS/PAGE analysis of radiolabeled proteins. Cocrystal-structure-assisted optimization of 3C-protease-directed Michael acceptors has yielded molecules having extremely rapid in vitro inactivation of the viral protease, potent antiviral activity against multiple rhinovirus serotypes and low cellular toxicity. Recently, one compound in this series, AG7088, has entered clinical trials.

Synthetic zeolites and other microporous oxide molecular sieves

Use of synthetic zeolites and other microporous oxides since 1950 has improved insulated windows, automobile air-conditioning, refrigerators, air brakes on trucks, laundry detergents, etc. Their large internal pore volumes, molecular-size pores, regularity of crystal structures, and the diverse framework chemical compositions allow “tailoring” of structure and properties. Thus, highly active and selective catalysts as well as adsorbents and ion exchangers with high capacities and selectivities were developed. In the petroleum refining and petrochemical industries, zeolites have made possible cheaper and lead-free gasoline, higher performance and lower-cost synthetic fibers and plastics, and many improvements in process efficiency and quality and in performance. Zeolites also help protect the environment by improving energy efficiency, reducing automobile exhaust and other emissions, cleaning up hazardous wastes (including the Three Mile Island nuclear power plant and other radioactive wastes), and, as specially tailored desiccants, facilitating the substitution of new refrigerants for the ozone-depleting chlorofluorocarbons banned by the Montreal Protocol.