Structural basis for chemical inhibition of human blood coagulation factor Xa

Factor Xa, the converting enzyme of prothrombin to thrombin, has emerged as an alternative (to thrombin) target for drug discovery for thromboembolic diseases. An inhibitor has been synthesized and the crystal structure of the complex between Des[1–44] factor Xa and the inhibitor has been determined by crystallographic methods in two different crystal forms to 2.3- and 2.4-Å resolution. The racemic mixture of inhibitor FX-2212, (2RS)-(3′-amidino-3-biphenylyl)-5-(4-pyridylamino)pentanoic acid, inhibits factor Xa activity by 50% at 272 nM in vitro. The S-isomer of FX-2212 (FX-2212a) was found to bind to the active site of factor Xa in both crystal forms. The biphenylamidine of FX-2212a occupies the S1-pocket, and the pyridine ring makes hydrophobic interactions with the factor Xa aryl-binding site. Several water molecules meditate inhibitor binding to residues in the active site. In contrast to the earlier crystal structures of factor Xa, such as those of apo-Des[1–45] factor Xa and Des[1–44] factor Xa in complex with a naphthyl inhibitor DX-9065a, two epidermal growth factor-like domains of factor Xa are well ordered in both our crystal forms as well as the region between the two domains, which recently was found to be the binding site of the effector cell protease receptor-1. This structure provides a basis for designing next generation inhibitors of factor Xa.

The roles of the two proton input channels in cytochrome c oxidase from Rhodobacter sphaeroides probed by the effects of site-directed mutations on time-resolved electrogenic intraprotein proton transfer

The crystal structures of cytochrome c oxidase from both bovine and Paracoccus denitrificans reveal two putative proton input channels that connect the heme-copper center, where dioxygen is reduced, to the internal aqueous phase. In this work we have examined the role of these two channels, looking at the effects of site-directed mutations of residues observed in each of the channels of the cytochrome c oxidase from Rhodobacter sphaeroides. A photoelectric technique was used to monitor the time-resolved electrogenic proton transfer steps associated with the photo-induced reduction of the ferryl-oxo form of heme a3 (Fe4+ = O2−) to the oxidized form (Fe3+OH−). This redox step requires the delivery of a “chemical” H+ to protonate the reduced oxygen atom and is also coupled to proton pumping. It is found that mutations in the K channel (K362M and T359A) have virtually no effect on the ferryl-oxo-to-oxidized (F-to-Ox) transition, although steady-state turnover is severely limited. In contrast, electrogenic proton transfer at this step is strongly suppressed by mutations in the D channel. The results strongly suggest that the functional roles of the two channels are not the separate delivery of chemical or pumped protons, as proposed recently [Iwata, S., Ostermeier, C., Ludwig, B. & Michel, H. (1995) Nature (London) 376, 660–669]. The D channel is likely to be involved in the uptake of both “chemical” and “pumped” protons in the F-to-Ox transition, whereas the K channel is probably idle at this partial reaction and is likely to be used for loading the enzyme with protons at some earlier steps of the catalytic cycle. This conclusion agrees with different redox states of heme a3 in the K362M and E286Q mutants under aerobic steady-state turnover conditions.

Nucleotide-dependent conformational change near the fulcrum region in Dictyostelium myosin II

In skeletal muscle myosin, the reactive thiols (SH1 and SH2) are close to a proposed fulcrum region that is thought to undergo a large conformational change. The reactive thiol region is thought to transmit the conformational changes induced by the actin–myosin–ATP interactions to the lever arm, which amplifies the power stroke. In skeletal muscle myosin, SH1 and SH2 can be chemically cross-linked in the presence of nucleotide, trapping the nucleotide in its pocket. Although the flexibility of the reactive thiol region has been well studied in skeletal muscle myosin, crystal structures of truncated nonmuscle myosin II from Dictyostelium in the presence of various ATP analogs do not show changes at the reactive thiol region that would be consistent with the SH1–SH2 cross-linking observed for muscle myosin. To examine the dynamics of the reactive thiol region in Dictyostelium myosin II, we have examined a modified myosin II that has cysteines at the muscle myosin SH1 and SH2 positions. This myosin is specifically cross-linked at SH1–SH2 by a chemical cross-linker in the presence of ADP, but not in its absence. Furthermore, the cross-linked species traps the nucleotide, as in the case of muscle myosin. Thus, the Dictyostelium myosin II shares the same dynamic behavior in the fulcrum region of the molecule as the skeletal muscle myosin. This result emphasizes the importance of nucleotide-dependent changes in this part of the molecule.

The conducting form of gramicidin A is a right-handed double-stranded double helix

The linear pentadecapeptide antibiotic, gramicidin D, is a naturally occurring product of Bacillus brevis known to form ion channels in synthetic and natural membranes. The x-ray crystal structures of the right-handed double-stranded double-helical dimers (DSDHℛ) reported here agree with 15N-NMR and CD data on the functional gramicidin D channel in lipid bilayers. These structures demonstrate single-file ion transfer through the channels. The results also indicate that previous crystal structure reports of a left-handed double-stranded double-helical dimer in complex with Cs+ and K+ salts may be in error and that our evidence points to the DSDHℛ as the major conformer responsible for ion transport in membranes.

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.