Tetrahedral Intermediate Research Articles

Structure- and Substrate-based Inhibitor Design for Clostridium botulinum Neurotoxin Serotype A

The seven antigenically distinct serotypes of Clostridium botulinum neurotoxins cleave specific soluble N-ethylmaleimide-sensitive factor attachment protein receptor complex proteins and block the release of neurotransmitters that cause flaccid paralysis and are considered potential bioweapons. Botulinum neurotoxin type A is the most potent among the clostridial neurotoxins, and to date there is no post-exposure therapeutic intervention available. To develop inhibitors leading to drug design, it is imperative that critical interactions between the enzyme and the substrate near the active site are known. Although enzyme-substrate interactions at exosites away from the active site are mapped in detail for botulinum neurotoxin type A, information about the active site interactions is lacking. Here, we present the crystal structures of botulinum neurotoxin type A catalytic domain in complex with four inhibitory substrate analog tetrapeptides, viz. RRGC, RRGL, RRGI, and RRGM at resolutions of 1.6-1.8 A. These structures show for the first time the interactions between the substrate and enzyme at the active site and delineate residues important for substrate stabilization and catalytic activity. We show that OH of Tyr(366) and NH(2) of Arg(363) are hydrogen-bonded to carbonyl oxygens of P1 and P1' of the substrate analog and position it for catalytic activity. Most importantly, the nucleophilic water is replaced by the amino group of the N-terminal residue of the tetrapeptide. Furthermore, the S1' site is formed by Phe(194), Thr(215), Thr(220), Asp(370), and Arg(363). The K(i) of the best inhibitory tetrapeptide is 157 nm.

Force‐field parameters for the simulation of tetrahedral intermediates of serine hydrolases

CHARMM force-field parameters are reported for the tetrahedral intermediate of serine hydrolases. The fitting follows the standard protocol proposed for CHARMM22. The reference data include ab initio (RHF/6-31G*) interaction energies of complexes between water and the model compound 1,1-dimethoxyethoxide, torsional profiles of related model compounds from correlated ab initio (MP2/6-311+G*//B3LYP/6-31+G*) calculations, as well as molecular geometries and vibrational frequencies from density functional theory (B3LYP/6-31+G*). The optimized parameters reproduce the target data well. Their utility is demonstrated by a QM/MM study of the tetrahedral intermediate in Bacillus subtilis lipase A, and by classical molecular modeling of enantioselectivity in Pseudomonas aeruginosa lipase and its mutants.

Monomer versus Alcohol Activation in the 4‐Dimethylaminopyridine‐Catalyzed Ring‐Opening Polymerization of Lactide and Lactic O‐Carboxylic Anhydride

Model reactions for the 4-dimethylaminopyridine (DMAP)-catalyzed ring-opening polymerization of lactide and the corresponding lactic O-carboxylic anhydride (lacOCA) have been studied computationally at the B3LYP/6-31G(d) level of theory. The solvent effect of dichloromethane was taken into account through PCM/SCRF single-point calculations at the B3LYP/6-31G(d) level of theory. In marked contrast with that predicted for the reaction of alcohols with acetic anhydride, the mechanism in which nucleophilic activation of the monomer involving acylpyridinium intermediates was found to be energetically less favorable than the base activation of the alcohol through hydrogen bonding. The concerted pathway for the ring-opening of lactide and lacOCA was shown to compete with the traditional stepwise mechanism involving tetrahedral intermediates. Furthermore, DMAP is proposed to act as a bifunctional catalyst through its basic nitrogen center and an acidic ortho-hydrogen atom.

Roles of Histidines 154 and 189 and Aspartate 139 in the Active Site of Serine Acetyltransferase from Haemophilus influenzae

A crystal structure of serine acetyltransferase (SAT) with cysteine bound in the serine subsite of the active site shows that both H154 and H189 are within hydrogen-bonding distance to the cysteine thiol [Olsen, L. R., Huang, B., Vetting, M. W., and Roderick, S. L. (2004) Biochemistry 43, 6013 -6019]. In addition, H154 is in an apparent dyad linkage with D139. The structure suggests that H154 is the most likely catalytic general base and that H189 and D139 may also play important roles during the catalytic reaction. Site-directed mutagenesis was performed to mutate each of these three residues to Asn, one at a time. The V1/Et value of all of the single mutant enzymes decreased, with the largest decrease (approximately 1240-fold) exhibited by the H154N mutant enzyme. Mutation of both histidines, H154N/H189N, gave a V1/Et approximately 23700-fold lower than that of the wild-type enzyme. An increase in K Ser was observed for the H189N, D139N, and H154N/H189N mutant enzymes, while the H154N mutant enzyme gave an 8-fold decrease in K Ser. For all three single mutant enzymes, V1/Et and V1/K Ser Et decrease at low pH and give a pKa of about 7, while the V1/Et of the double mutant enzyme was pH independent. The solvent deuterium kinetic isotope effects on V 1 and V1/K Ser decreased compared to wild type for the H154N mutant enzyme and increased for the H189N mutant enzyme but was about the same as that of wild type for D139N and H154N/H189N. Data suggest that H154, H189, and D139 play different catalytic roles for SAT. H154 likely serves as a general base, accepting a proton from the beta-hydroxyl of serine as the tetrahedral intermediate is formed upon nucleophilic attack on the thioester carbonyl of acetyl-CoA. However, activity is not completely lost upon elimination of H154, and thus, H189 may be able to serve as a backup general base at a lower efficiency compared to H154; it also aids in binding and orienting the serine substrate. Aspartate 139, in dyad linkage with H154, likely facilitates catalysis by increasing the basicity of H154.

Catalytic Features of the Botulinum Neurotoxin A Light Chain Revealed by High Resolution Structure of an Inhibitory Peptide Complex

The Clostridium botulinum neurotoxin serotype A light chain (BoNT/A-LC) is a Zn(II)-dependent metalloprotease that blocks the release of acetylcholine at the neuromuscular junction by cleaving SNAP-25, one of the SNARE proteins required for exocytosis. Because of the potential for use of the toxin in bioterrorism and the increasingly widespread application of the toxin in the medical field, there is significant interest in the development of small-molecule inhibitors of the metalloprotease. Efforts to design such inhibitors have not benefited from knowledge of how peptides bind to the active site since the enzyme-peptide structures available previously either were not occupied in the vicinity of the catalytic Zn(II) ion or did not represent the product of SNAP-25 substrate cleavage. Herein we report the 1.4 A-resolution X-ray crystal structure of a complex between the BoNT/A-LC and the inhibitory peptide N-Ac-CRATKML, the first structure of the light chain with an inhibitory peptide bound at the catalytic Zn(II) ion. The peptide is bound with the Cys S gamma atom coordinating the metal ion. Surprisingly, the cysteine sulfur is oxidized to the sulfenic acid form. Given the unstable nature of this species in solution, is it likely that oxidation occurs on the enzyme. In addition to the peptide-bound structure, we report two structures of the unliganded light chain with and without the Zn(II) cofactor bound at 1.25 and 1.20 A resolution, respectively. The two structures are nearly identical, confirming that the Zn(II) ion plays a purely catalytic role. Additionally, the structure of the Zn(II)-bound uncomplexed enzyme allows identification of the catalytic water molecule and a second water molecule that occupies the same position as the peptidic oxygen in the tetrahedral intermediate. This observation suggests that the enzyme active site is prearranged to stabilize the tetrahedral intermediate of the protease reaction.

Computational Study on the Aminolysis of β-Hydroxy-α,β-Unsaturated Ester via the Favorable Path Including the Formation of α-Oxo Ketene Intermediate

The possible mechanisms of the aminolysis of N-methyl-3-(methoxycarbonyl)-4-hydroxy-2-pyridone (beta-hydroxy-alpha,beta-unsaturated ester) with dimethylamine are investigated at the hybrid density functional theory B3LYP/6-31G(d,p) level in the gas phase. Single-point computations at the B3LYP/6-311++G(d,p) and the Becke88-Becke95 1-parameter model BB1K/6-311++G(d,p) levels are performed for more precise energy predictions. Solvent effects are also assessed by single-point calculations at the integral equation formalism polarized continuum model IEFPCM-B3LYP/6-311++G(d,p) and IEFPCM-BB1K/6-311++G(d,p) levels on the gas-phase optimized geometries. Three possible pathways, the concerted pathway (path A), the stepwise pathway involving tetrahedral intermediates (path B), and the stepwise pathway via alpha-oxo ketene intermediate due to the participation of beta-hydroxy (path C), are taken into account for the title reaction. Moreover, path C includes two sequential processes. The first process is to generate alpha-oxo ketene intermediate via the decomposition of N-methyl-3-(methoxycarbonyl)-4-hydroxy-2-pyridone; the second process is the addition of dimethylamine to alpha-oxo ketene intermediate. Our results indicate that path C is more favorable than paths A and B both in the gas phase and in solvent (heptane). In path C, the first process is the rate-determining step, and the second process is revealed to be a [4+2] pseudopericyclic reaction without the energy barrier. Being independent of the concentration of amine, the first process obeys the first-order rate law.

The Electrostatic Driving Force for Nucleophilic Catalysis inl-Arginine Deiminase: A Combined Experimental and Theoretical Study

L-arginine deiminase (ADI) catalyzes the hydrolysis of L-arginine to form L-citrulline and ammonia via two partial reactions. A working model of the ADI catalytic mechanism assumes nucleophilic catalysis by a stringently conserved active site Cys and general acid-general base catalysis by a stringently conserved active site His. Accordingly, in the first partial reaction, the Cys attacks the substrate guanidino C zeta atom to form a tetrahedral covalent adduct, which is protonated by the His at the departing ammonia group to facilitate the formation of the Cys- S-alkylthiouronium intermediate. In the second partial reaction, the His activates a water molecule for nucleophilic addition at the thiouronium C zeta atom to form the second tetrahedral intermediate, which eliminates the Cys in formation of the L-citrulline product. The absence of a basic residue near the Cys thiol suggested that the electrostatic environment of the Cys thiol, in the enzyme-substrate complex, stabilizes the Cys thiolate anion. The studies described in this paper explore the mechanism of stabilization of the Cys thiolate. First, the log(k(cat)/K(m)) and log k(cat) pH rate profiles were measured for several structurally divergent ADIs to establish the pH range for ADI catalysis. All ADIs were optimally active at pH 5, which suggested that the Cys pKa is strongly perturbed by the prevailing electrostatics of the ADI active site. The p K a of the Bacillus cereus ADI (BcADI) was determined by UV-pH titration to be 9.6. In contrast, the pKa determined by iodoacetamide Cys alkylation is 6.9. These results suggest that the negative electrostatic field from the two opposing Asp carboxylates perturbs the Cys pKa upward in the apoenzyme and that the binding of the iodoacetamide (a truncated analogue of the citrulline product) between the Cys thiol and the two Asp carboxylates shields the Cys thiol, thereby reducing its pKa. It is hypothesized that the bound positively charged guanidinium group of the L-arginine substrate further stabilizes the Cys thiolate. The so-called "substrate-assisted" Cys ionization, first reported by Fast and co-workers to operate in the related enzyme dimethylarginine dimethylaminohydrolase [Stone, E. M., Costello, A. L., Tierney, D. L., and Fast, W. (2006) Biochemistry 45, 5618-5630], was further explored computationally in ADI by using an ab initio quantum mechanics/molecular mechanics method. The energy profiles for formation of the tetrahedral intermediate in the first partial reaction were calculated for three different reaction scenarios. From these results, we conclude that catalytic turnover commences from the active configuration of the ADI(L-arginine) complex which consists of the Cys thiolate (nucleophile) and His imidazolium ion (general acid) and that the energy barriers for the nucleophilic addition of Cys thiolate to the thiouronium C zeta atom and His imidazolium ion-assisted elimination from the tetrahedral intermediate are small.

A Thermal Decarbonylation of Penam β-Lactams

Penam acids 6-8 [i.e., (2S,5R,6R)-, (2S,5S,6R)-, and (2S,5R,6S)-isomers of 6-(1,3-dihydro-1,3-dioxo-2H-isoindol-2-yl)-3,3-dimethyl-7-oxo-4-thia-1-azabicyclo[3.2.0]heptane-2-carboxylic acid] were prepared by deesterification of the corresponding methyl esters 2-4. The same methodology applied to ester 1 did not lead to the (2S,5S,6S)-isomer 5 but rather a 72% yield of the thiazoline derivative 9. High-resolution mass spectrometry analysis of the reaction headspace gases indicated that a stoichiometric amount of carbon monoxide is produced during the deesterification of 1. A mechanism for this decarbonylation reaction is proposed. This appears to represent a new type of fragmentation reaction for a penam carboxylic acid. The free energies of various reaction species along viable decarbonylation reaction coordinates for acids 5 and 7 were computed by using the density functional theory method IEFPCM/M06/6-31+G(d). Anionic and zwitterionic (neutral) variants of the proposed mechanism were considered, but each produced computed activation free energies deemed to be too high (>45 kcal/mol) to be experimentally relevant. The computed activation free energies for the protonated (cationic) variant of the mechanism were 17.3 kcal/mol for 7 vs 8.8 kcal/mol for 5. The value of this difference in energies of activation (DeltaDeltaG++) is quite consistent with experimental observations and supports the proposed mechanism. For a portion of the computed reaction coordinate that involves ring opening of the lactam ring by an internal carboxylic acid group to form a cyclic anhydride, the expected tetrahedral intermediate was circumvented by a direct (concerted) and facile N- to O-acyl migration event. Additional thermal gas-phase reaction products produced during gas chromatographic analysis of the penams 1-8 were characterized with high-resolution mass spectrometry, and possible mechanisms for their formation are presented.

5-Enolpyruvylshikimate-3-phosphate synthase: Determination of the protonation state of active site residues by the semiempirical method

EPSP synthase (EPSPS) catalyzes the addition of shikimate-3-phosphate (S3P) and phosphoenolpyruvate (PEP) to form a tetrahedral intermediate (TI) that is converted to 5-enolpyruvylshikimate-3-phosphate (EPSP) and inorganic phosphate. A semiempirical molecular modeling study of the EPSPS active site containing the TI was implemented for the assignment of the protonation states of four basic residues, Lys22, Lys340, His385, and Lys411, based on the evaluation of 16 different protonation states and comparison of the resulting energy minimized heavy atoms coordinates with available X-ray crystallographic data of the D313A mutant of EPSPS. The results, employing both gas phase and continuum solvent models, are indicative that after the TI formation the histidine residue is most probably in neutral form (N ε-protonated) and the lysine residues are in protonated form, which suggests that none of the presently proposed assignments of aminoacid residues involved in the reaction mechanism could be completely correct. The protonated state of Lys22 in the presence of the TI supports the proposal that this residue is a general acid catalyst for TI breakdown. Modeling of the native enzyme active site suggests that Asp313 residue has only minor effects on the definition of the TI position inside the active site. Hydrogen-bonds distances suggest that, in order to act as a base, Asp313 needs the intermediacy of a hydroxyl group of the TI for effecting the attack on the TI methyl group in the elimination step leading to EPSP, as suggested previously in the literature.

Amino Acids Allosterically Regulate the Thiamine Diphosphate-dependent α-Keto Acid Decarboxylase from Mycobacterium tuberculosis

The gene rv0853c from Mycobacterium tuberculosis strain H37Rv codes for a thiamine diphosphate-dependent alpha-keto acid decarboxylase (MtKDC), an enzyme involved in the amino acid degradation via the Ehrlich pathway. Steady state kinetic experiments were performed to determine the substrate specificity of MtKDC. The mycobacterial enzyme was found to convert a broad spectrum of branched-chain and aromatic alpha-keto acids. Stopped-flow kinetics showed that MtKDC is allosterically activated by alpha-keto acids. Even more, we demonstrate that also amino acids are potent activators of this thiamine diphosphate-dependent enzyme. Thus, metabolic flow through the Ehrlich pathway can be directly regulated at the decarboxylation step. The influence of amino acids on MtKDC catalysis was investigated, and implications for other thiamine diphosphate-dependent enzymes are discussed.

Nucleophilic Attack on the Carbon - Nitrogen Double Bond Leading to Tetrahedral Intermediates with Conformationally Restricted Stereochemistry

Methoxide ion substitution on 4-chloro-1H-2,3-benzoxazine 4 has been investigated. The rates of cyclization have been measured for (Z)- and (E)-O-(2-hydroxyethyl)benzohydroximoyl halides (PhC(X)=NOCH2CH2OH; X = Cl or Br) using potassium t-butoxide to generate alkoxide ions from the alcohols. The element effect for the Z and E isomers are kZBr/kZCl = 3.08 and kEBr/kECl = 2.38. The Z-to-E rate ratios are 0.25 (chlorides) and 0.33 (bromides). These reactions are proceeding by rate-determining nucleophilic attack by the alkoxide ion on the carbon–nitrogen double bond (AN# + DN). The restricted stereochemistry of the tetrahedral intermediates suggests that the E-hydroximoyl bromide and chloride may be undergoing unassisted ionization to form a zwitterionic intermediate that leads to product. Alternatively, the chair conformation of the tetrahedral intermediate may be undergoing a ring flip to a boat conformation in which there are antiperiplanar electron pairs that assist in the ionization of the bromide or chloride ion.

Structural Requirements for the Activation of Escherichia coli CTP Synthase by the Allosteric Effector GTP Are Stringent, but Requirements for Inhibition Are Lax

Cytidine 5'-triphosphate synthase catalyzes the ATP-dependent formation of CTP from UTP using either NH(3) or l-glutamine (Gln) as the source of nitrogen. GTP acts as an allosteric effector promoting Gln hydrolysis but inhibiting Gln-dependent CTP formation at concentrations of >0.15 mM and NH(3)-dependent CTP formation at all concentrations. A structure-activity study using a variety of GTP and guanosine analogues revealed that only a few GTP analogues were capable of activating Gln-dependent CTP formation to varying degrees: GTP approximately 6-thio-GTP > ITP approximately guanosine 5'-tetraphosphate > O(6)-methyl-GTP > 2'-deoxy-GTP. No activation was observed with guanosine, GMP, GDP, 2',3'-dideoxy-GTP, acycloguanosine, and acycloguanosine monophosphate, indicating that the 5'-triphosphate, 2'-OH, and 3'-OH are required for full activation. The 2-NH(2) group plays an important role in binding recognition, whereas substituents at the 6-position play an important role in activation. The presence of a 6-NH(2) group obviates activation, consistent with the inability of ATP to substitute for GTP. Nucleotide and nucleoside analogues of GTP and guanosine, respectively, all inhibited NH(3)- and Gln-dependent CTP formation (often in a cooperative manner) to a similar extent (IC(50) approximately 0.2-0.5 mM). This inhibition appeared to be due solely to the purine base and was relatively insensitive to the identity of the purine with the exception of inosine, ITP, and adenosine (IC(50) approximately 4-12 mM). 8-Oxoguanosine was the best inhibitor identified (IC(50) = 80 microM). Our findings suggest that modifying 2-aminopurine or 2-aminopurine riboside may serve as an effective strategy for developing cytidine 5'-triphosphate synthase inhibitors.

Theoretical study on the hydrolysis mechanism of N,N‐dimethyl‐N′‐(2‐oxo‐1, 2‐dihydro‐pyrimidinyl)formamidine: Water‐assisted mechanism and cluster‐continuum model

The hydrolysis reaction of N,N-dimethyl-N'-(2-oxo-1, 2-dihydro-pyrimidinyl)formamidine (DMPFA), a model compound of the antivirus drug amidine-3TC (3TC = 2', 3'-dideoxy-3'-thiacytidine), is investigated by the hybrid density functional theory B3LYP/6-31+G (d,p) method. The hydrolysis reaction of the title compound is predicted to undergo via two pathways, each of which is a stepwise process. Path A is the addition of H2O to the C=N double bond in the amidine group to form a tetrahedral structure in its first step, and then the transfer of the H atom of hydroxyl leads to the corresponding products via four possible channels. Path B simultaneously involves the nucleophilic attack of H2O to the C atom of the C=N bond and the proton transfer to the N atom of amino group leading to the cleavage of the C-N single bond in the amidine group. The results indicate that path A is more favorable than path B in the gas phase. Moreover, to simulate the title reaction in aqueous solution, water-assisted mechanism and the cluster-continuum model, based on the SCRF/CPCM model, are taken into account in our work. The results indicate that it is rational for two water molecules served as a bridge to assist in the first step of path A and that cytosine rather than the cytosine-substituted formamide should be released from the tetrahedral intermediate via s six-membered cycle transition state (channel 2). Our calculations exhibit that the process toward the tetrahedral intermediate is the rate-determining step both in the gas phase and in aqueous solution.

Reaction of poly‐L‐lysine with aryl acetates and aryl methyl carbonates. A mechanistic study

AbstractThe pH profile (log k vs. pH) of the reactions of poly‐L‐lysine (PL) with a series of aryl acetates and aryl methyl carbonates in aqueous solution show the same conformational changes as those determined by potentiometric titrations. When PL is a random coil, the most probable mechanism for the reactions studied is through the formation of a tetrahedral intermediate and its breakdown to products as the rate‐determining step. The tetrahedral intermediate is stabilized by a hydrogen bond interaction between the nitro groups in the substrate and the NH group of the principal chain or some NH2 groups of the lateral chains. Copyright © 2007 John Wiley & Sons, Ltd.

Transglutaminase 2 Undergoes a Large Conformational Change upon Activation

Human transglutaminase 2 (TG2), a member of a large family of enzymes that catalyze protein crosslinking, plays an important role in the extracellular matrix biology of many tissues and is implicated in the gluten-induced pathogenesis of celiac sprue. Although vertebrate transglutaminases have been studied extensively, thus far all structurally characterized members of this family have been crystallized in conformations with inaccessible active sites. We have trapped human TG2 in complex with an inhibitor that mimics inflammatory gluten peptide substrates and have solved, at 2-Å resolution, its x-ray crystal structure. The inhibitor stabilizes TG2 in an extended conformation that is dramatically different from earlier transglutaminase structures. The active site is exposed, revealing that catalysis takes place in a tunnel, bridged by two tryptophan residues that separate acyl-donor from acyl-acceptor and stabilize the tetrahedral reaction intermediates. Site-directed mutagenesis was used to investigate the acyl-acceptor side of the tunnel, yielding mutants with a marked increase in preference for hydrolysis over transamidation. By providing the ability to visualize this activated conformer, our results create a foundation for understanding the catalytic as well as the non-catalytic roles of TG2 in biology, and for dissecting the process by which the autoantibody response to TG2 is induced in celiac sprue patients.

1.2Å-resolution crystal structures reveal the second tetrahedral intermediates of streptogrisin B (SGPB)

Streptogrisin B (SGPB) has served as one of the models for studying the catalytic activities of serine peptidases. Here we report its native crystal structures at pH 4.2 at a resolution of 1.18Å, and at pH 7.3 at a resolution of 1.23Å. Unexpectedly, outstanding electron density peaks occurred in the active site and the substrate-binding region of SGPB in the computed maps at both pHs. The densities at pH 4.2 were assigned as a tetrapeptide, Asp-Ala-Ile-Tyr, whereas those at pH 7.3 were assigned as a tyrosine molecule and a leucine molecule existing at equal occupancies in both of the SGPB molecules in the asymmetric unit. Refinement with relaxed geometric restraints resulted in molecular structures representing mixtures of the second tetrahedral intermediates and the enzyme-product complexes of SGPB existing in a pH-dependent equilibrium. Structural comparisons with the complexes of SGPB with turkey ovomucoid third domain (OMTKY3) and its variants have shown that, upon the formation of the tetrahedral intermediate, residues Glu192A to Gly193 of SGPB move towards the α-carboxylate O of residue P1 of the bound species, and adjustments in the side-chain conformational angles of His57 and Ser195 of SGPB favor the progression of the catalytic mechanism of SGPB.

Cavitands with Introverted Functionality Stabilize Tetrahedral Intermediates

The synthesis and characterization of two deepened cavitand hosts with introverted functionality--functional groups directed into the cavity--is described. Two functions can be introverted, alcohol and aldehyde, and they show the formation of hemiacetals and hemiketals on binding small guests with complementary functional groups. The structures of the bound hemiacetals are determined by 1D and 2D NMR studies. The arrangements of the guests in the cavitands enhance the equilibrium constants of carbonyl additions, K/K(ctrl), between 13- and 10(5)-fold, compared to their counterparts in solution. The stabilization of the addition products is due to the prior complexation of the guests and the organized solvation provided to the tetrahedral intermediates by a network of secondary amides at the cavitand rim.

Detection of Reactive Tetrahedral Intermediates in a Deep Cavitand with an Introverted Functionality

Labile hemiaminal intermediates are stabilized by binding in a deep cavitand with an introverted aldehyde functionality. The aldehyde is attached to the cavitand via an anthracene spacer that rotates rapidly about the cavitand rim. The half-lives of these hemiaminals vary from 30 min to over 100 h at ambient temperature, due to hydrogen bonding with the organized peptide-like framework at the cavitand rim. The intermediates are sufficiently long-lived to allow study by 2D NMR techniques requiring many hours of acquisition time. Mechanistic analysis of the dehydration step shows first-order kinetics. The analogous "extroverted" reaction was also performed, where the addition took place outside the cavitand, displaying standard steady-state kinetics; no hemiaminal was observed. The cavitand shows strong selectivity based not on binding affinity but upon the rate of the product-forming step. A 10:1 ratio of product imines was obtained, while the initial binding ratio was 1:1. The cavitand acts as a mimic of enzymes in that it uses weak binding forces to stabilize reactive intermediates and isolates them from the medium. The synthetic environment allows direct detection and analysis of the intermediates, as opposed to natural systems that must be analyzed indirectly.

Molecular Dynamics Simulations of the Catalytic Pathway of a Cysteine Protease: A Combined QM/MM Study of Human Cathepsin K

Molecular dynamics simulations using a combined QM/MM potential have been performed to study the catalytic mechanism of human cathepsin K, a member of the papain family of cysteine proteases. We have determined the two-dimensional free energy surfaces of both acylation and deacylation steps to characterize the reaction mechanism. These free energy profiles show that the acylation step is rate limiting with a barrier height of 19.8 kcal/mol in human cathepsin K and of 29.3 kcal/mol in aqueous solution. The free energy of activation for the deacylation step is 16.7 kcal/mol in cathepsin K and 17.8 kcal/mol in aqueous solution. The reduction of free energy barrier is achieved by stabilization of the oxyanion in the transition state. Interestingly, although the "oxyanion hole" has been formed in the Michaelis complex, the amide units do not donate hydrogen bonds directly to the carbonyl oxygen of the substrate, but they stabilize the thiolate anion nucleophile. Hydrogen-bonding interactions are induced as the substrate amide group approaches the nucleophile, moving more than 2 A and placing the oxyanion in contact with Gln19 and the backbone amide of Cys25. The hydrolysis of peptide substrate shares a common mechanism both for the catalyzed reaction in human cathepsin K and for the uncatalyzed reaction in water. Overall, the nucleophilic attack by Cys25 thiolate and the proton-transfer reaction from His162 to the amide nitrogen are highly coupled, whereas a tetrahedral intermediate is formed along the nucleophilic reaction pathway.

NMR Study of the Inhibition of Pepsin by Glyoxal Inhibitors: Mechanism of Tetrahedral Intermediate Stabilization by the Aspartyl Proteases

Z-Ala-Ala-Phe-glyoxal (where Z is benzyloxycarbonyl) has been shown to be a competitive inhibitor of pepsin with a Ki = 89 +/- 24 nM at pH 2.0 and 25 degrees C. Both the ketone carbon (R13COCHO) and the aldehyde carbon (RCO13CHO) of the glyoxal group of Z-Ala-Ala-Phe-glyoxal have been 13C-enriched. Using 13C NMR, it has been shown that when the inhibitor is bound to pepsin, the glyoxal keto and aldehyde carbons give signals at 98.8 and 90.9 ppm, respectively. This demonstrates that pepsin binds and preferentially stabilizes the fully hydrated form of the glyoxal inhibitor Z-Ala-Ala-Phe-glyoxal. From 13C NMR pH studies with glyoxal inhibitor, we obtain no evidence for its hemiketal or hemiacetal hydroxyl groups ionizing to give oxyanions. We conclude that if an oxyanion is formed its pKa must be >8.0. Using 1H NMR, we observe four hydrogen bonds in free pepsin and in pepsin/Z-Ala-Ala-Phe-glyoxal complexes. In the pepsin/pepstatin complex an additional hydrogen bond is formed. We examine the effect of pH on hydrogen bond formation, but we do not find any evidence for low-barrier hydrogen bond formation in the inhibitor complexes. We conclude that the primary role of hydrogen bonding to catalytic tetrahedral intermediates in the aspartyl proteases is to correctly orientate the tetrahedral intermediate for catalysis.