Azepanone-based inhibitors of human and rat cathepsin K.

The synthesis, in vitro activities, and pharmacokinetics of a series of azepanone-based inhibitors of the cysteine protease cathepsin K (EC 3.4.22.38) are described. These compounds show improved configurational stability of the C-4 diastereomeric center relative to the previously published five- and six-membered ring ketone-based inhibitor series. Studies in this series have led to the identification of 20, a potent, selective inhibitor of human cathepsin K (K(i) = 0.16 nM) as well as 24, a potent inhibitor of both human (K(i) = 0.0048 nM) and rat (K(i,app) = 4.8 nM) cathepsin K. Small-molecule X-ray crystallographic analysis of 20 established the C-4 S stereochemistry as being critical for potent inhibition and that unbound 20 adopted the expected equatorial conformation for the C-4 substituent. Molecular modeling studies predicted the higher energy axial orientation at C-4 of 20 when bound within the active site of cathepsin K, a feature subsequently confirmed by X-ray crystallography. Pharmacokinetic studies in the rat show 20 to be 42% orally bioavailable. Comparison of the transport of the cyclic and acyclic analogues through CaCo-2 cells suggests that oral bioavailability of the acyclic derivatives is limited by a P-glycoprotein-mediated efflux mechanism. It is concluded that the introduction of a conformational constraint has served the dual purpose of increasing inhibitor potency by locking in a bioactive conformation as well as locking out available conformations which may serve as substrates for enzyme systems that limit oral bioavailability.

Fragment-based modeling of NAD binding to the catalytic subunits of diphtheria and pertussis toxins.

We describe a novel application of a fragment-based ligand docking technique; similar methods are commonly applied to the de novo design of ligands for target protein binding sites. We have used several new flexible docking and superposition tools, as well as a more conventional rigid-body (fragment) docking method, to examine NAD binding to the catalytic subunits of diphtheria (DT) and pertussis (PT) toxins, and to propose a model of the NAD-PT complex. Docking simulations with the rigid NAD fragments adenine and nicotinamide revealed that the low-energy dockings clustered in three distinct sites on the two proteins. Two of the sites were common to both fragments and were related to the structure of NAD bound to DT in an obvious way; however, the adenine subsite of PT was shifted relative to that of DT. We chose adenine/nicotinamide pairs of PT dockings from these clusters and flexibly superimposed NAD onto these pairs. A Monte Carlo-based flexible docking procedure and energy minimization were used to refine the modeled NAD-PT complexes. The modeled complex accounts for the sequence and structural similarities between PT and DT and is consistent with many results that suggest the catalytic importance of certain residues. A possible functional role for the structural difference between the two complexes is discussed.

Structure of the shiga-like toxin I B-pentamer complexed with an analogue of its receptor Gb3.

Shiga-like toxin I (SLT-I) is a virulence factor of Escherichia coli strains that cause disease in humans. Like other members of the Shiga toxin family, it consists of an enzymatic (A) subunit and five copies of a binding subunit (the B-pentamer). The B-pentamer binds to a specific glycolipid, globotriaosylceramide (Gb3), on the surface of target cells and thereby plays a crucial role in the entry of the toxin. Here we present the crystal structure at 2.8 A resolution of the SLT-I B-pentamer complexed with an analogue of the Gb3 trisaccharide. The structure reveals a surprising density of binding sites, with three trisaccharide molecules bound to each B-subunit monomer of 69 residues. All 15 trisaccharides bind to one side of the B-pentamer, providing further evidence that this side faces the cell membrane. The structural model is consistent with data from site-directed mutagenesis and binding of carbohydrate analogues, and allows the rational design of therapeutic Gb3 analogues that block the attachment of toxin to cells.

Modeling the carbohydrate-binding specificity of pig edema toxin.

The wild-type binding pentamer of Shiga-like toxin IIe (SLT-IIe) binds both the globotriaosylceramide (Gb3) and globotetraosylceramide (Gb4) cell surface glycolipids, whereas the double mutant GT3 (Q65E/K67Q) exhibits a marked preference for Gb3 [Tyrrell, G. J., et al. (1992) Proc. Natl. Acad. Sci. U.S.A. 89, 524-528]. We modeled three unique sites (sites 1-3) for binding of the carbohydrate moiety of Gb3 to GT3 and SLT-IIe, on the basis of the three sites observed for the SLT-I pentamer [Ling, H., et al. (1998) Biochemistry 37, 1777-1788]. Examination of the three sites in light of various mutation and binding data strongly suggested that one of the binding sites plays a role in the change of specificity observed for the GT3 mutant. We applied several modeling techniques, and developed a model for binding of the carbohydrate moiety of Gb4 to this site of the SLT-IIe binding pentamer. This model is consistent with a wide variety of mutation and binding data and clearly shows the importance of the terminal GalNAc residue of Gb4, as well as that of the two mutated residues of GT3, to the intermolecular interaction.

Atomic solvation parameters in the analysis of protein-protein docking results.

Several sets of amino acid surface areas and transfer free energies were used to derive a total of nine sets of atomic solvation parameters (ASPs). We tested the accuracy of each of these sets of parameters in predicting the experimentally determined transfer free energies of the amino acid derivatives from which the parameters were derived. In all cases, the calculated and experimental values correlated well. We then chose three parameter sets and examined the effect of adding an energetic correction for desolvation based on these three parameter sets to the simple potential function used in our multiple start Monte Carlo docking method. A variety of protein-protein interactions and docking results were examined. In the docking simulations studied, the desolvation correction was only applied during the final energy calculation of each simulation. For most of the docking results we analyzed, the use of an octanol-water-based ASP set marginally improved the energetic ranking of the low-energy dockings, whereas the other ASP sets we tested disturbed the ranking of the low-energy dockings in many of the same systems. We also examined the correlation between the experimental free energies of association and our calculated interaction energies for a series of proteinase-inhibitor complexes. Again, the octanol-water-based ASP set was compatible with our standard potential function, whereas ASP sets derived from other solvent systems were not.

Monte Carlo docking with ubiquitin.

The development of general strategies for the performance of docking simulations is prerequisite to the exploitation of this powerful computational method. Comprehensive strategies can only be derived from docking experiences with a diverse array of biological systems, and we have chosen the ubiquitin/diubiquitin system as a learning tool for this process. Using our multiple-start Monte Carlo docking method, we have reconstructed the known structure of diubiquitin from its two halves as well as from two copies of the uncomplexed monomer. For both of these cases, our relatively simple potential function ranked the correct solution among the lowest energy configurations. In the experiments involving the ubiquitin monomer, various structural modifications were made to compensate for the lack of flexibility and for the lack of a covalent bond in the modeled interaction. Potentially flexible regions could be identified using available biochemical and structural information. A systematic conformational search ruled out the possibility that the required covalent bond could be formed in one family of low-energy configurations, which was distant from the observed dimer configuration. A variety of analyses was performed on the low-energy dockings obtained in the experiment involving structurally modified ubiquitin. Characterization of the size and chemical nature of the interface surfaces was a powerful adjunct to our potential function, enabling us to distinguish more accurately between correct and incorrect dockings. Calculations with the structure of tetraubiquitin indicated that the dimer configuration in this molecule is much less favorable than that observed in the diubiquitin structure, for a simple monomer-monomer pair. Based on the analysis of our results, we draw conclusions regarding some of the approximations involved in our simulations, the use of diverse chemical and biochemical information in experimental design and the analysis of docking results, as well as possible modifications to our docking protocol.

Photoinactivation of acetylcholinesterase by erythrosin B and related compounds.

Acetylcholinesterase was rapidly inactivated when exposed to light in the presence of xanthene dyes. Photosensitizing efficiency paralleled the dye triplet state quantum yields, increasing in the order fluorescein less than eosin B less than eosin Y less than erythrosin B less than rose bengal. The observed first-order rate constants of photoinactivation increased hyperbolically with dye concentration. Evidence for the formation of a dye-enzyme complex prior to inactivation was obtained from spectrophotometric and protein fluorescence quenching methods. The latter technique allowed estimates of the dye-enzyme dissociation constants for rose bengal (20 microM) and erythrosin B (30 microM). After photoinactivation, a portion of the dye became covalently bound to the enzyme. The photoinactivation reaction occurs in both aerobic (air saturated) and anaerobic (argon saturated) solution, with the rates of photoinactivation being about three to five times greater under the latter conditions. The aerobic reaction exhibits a large deuterium isotope enhancement effect and is largely (but not completely) quenched by 10(-2) M azide. The anaerobic reaction is unaffected by azide and exhibits only a small deuterium isotope effect. These results indicate that the photoinactivation reaction proceeds mainly by a type II (singlet oxygen mediated) pathway under aerobic conditions and by a type I (radical) pathway under anaerobic conditions. The enzyme was protected from inactivation by edrophonium, a competitive inhibitor, but not by d-tubocurarine, a peripheral-site ligand, indicating that destruction of a crucial residue at or near the catalytic site is an important component of the inactivation process. Extensive destruction of tryptophan undoubtedly occurs, at least under aerobic conditions.(ABSTRACT TRUNCATED AT 250 WORDS)