Why Does Trypsin Cleave BPTI so Slowly?

Prebound State Discovered in the Unbinding Pathway of Fluorinated Variants of the Trypsin-BPTI Complex Using Random Acceleration Molecular Dynamics Simulations

The serine protease trypsin forms a tightly bound inhibitor complex with the bovine pancreatic trypsin inhibitor (BPTI). The complex is stabilized by the P1 residue Lys15, which interacts with negatively charged amino acids at the bottom of the S1 pocket. Truncating the P1 residue of wildtype BPTI to α-aminobutyric acid (Abu) leaves a complex with moderate inhibitor strength, which is held in place by additional hydrogen bonds at the protein-protein interface. Fluorination of the Abu residue partially restores the inhibitor strength. The mechanism with which fluorination can restore the inhibitor strength is unknown, and accurate computational investigation requires knowledge of the binding and unbinding pathways. The preferred unbinding pathway is likely to be complex, as encounter states have been described before, and unrestrained umbrella sampling simulations of these complexes suggest additional energetic minima. Here, we use random acceleration molecular dynamics to find a new metastable state in the unbinding pathway of Abu-BPTI variants and wildtype BPTI from trypsin, which we call the prebound state. The prebound state and the fully bound state differ by a substantial shift in the position, a slight shift in the orientation of the BPTI variants, and changes in the interaction pattern. Particularly important is the breaking of three hydrogen bonds around Arg17. Fluorination of the P1 residue lowers the energy barrier of the transition between the fully bound state and prebound state and also lowers the energy minimum of the prebound state. While the effect of fluorination is in general difficult to quantify, here, it is in part caused by favorable stabilization of a hydrogen bond between Gln194 and Cys14. The interaction pattern of the prebound state offers insights into the inhibitory mechanism of BPTI and might add valuable information for the design of serine protease inhibitors.

The monodomain Kunitz protein EgKU-7 from the dog tapeworm Echinococcus granulosus is a high-affinity trypsin inhibitor with two interaction sites

Typical Kunitz proteins (I2 family of the MEROPS database, Kunitz-A family) are metazoan competitive inhibitors of serine peptidases that form tight complexes of 1:1 stoichiometry, mimicking substrates. The cestode Echinococcus granulosus, the dog tapeworm causing cystic echinococcosis in humans and livestock, encodes an expanded family of monodomain Kunitz proteins, some of which are secreted to the dog host interface. The Kunitz protein EgKU-7 contains, in addition to the Kunitz domain with the anti-peptidase loop comprising a critical arginine, a C-terminal extension of ∼20 amino acids. Kinetic, electrophoretic, and mass spectrometry studies using EgKU-7, a C-terminally truncated variant, and a mutant in which the critical arginine was substituted by alanine, show that EgKU-7 is a tight inhibitor of bovine and canine trypsins with the unusual property of possessing two instead of one site of interaction with the peptidases. One site resides in the anti-peptidase loop and is partially hydrolyzed by bovine but not canine trypsins, suggesting specificity for the target enzymes. The other site is located in the C-terminal extension. This extension can be hydrolyzed in a particular arginine by cationic bovine and canine trypsins but not by anionic canine trypsin. This is the first time to our knowledge that a monodomain Kunitz-A protein is reported to have two interaction sites with its target. Considering that putative orthologs of EgKU-7 are present in other cestodes, our finding unveils a novel piece in the repertoire of peptidase-inhibitor interactions and adds new notes to the evolutionary host-parasite concerto.

Structure-function study of a Ca2+-independent metacaspase involved in lateral root emergence

Metacaspases are part of an evolutionarily broad family of multifunctional cysteine proteases, involved in disease and normal development. As the structure-function relationship of metacaspases remains poorly understood, we solved the X-ray crystal structure of an Arabidopsis thaliana type II metacaspase (AtMCA-IIf) belonging to a particular subgroup not requiring calcium ions for activation. To study metacaspase activity in plants, we developed an in vitro chemical screen to identify small molecule metacaspase inhibitors and found several hits with a minimal thioxodihydropyrimidine-dione structure, of which some are specific AtMCA-IIf inhibitors. We provide mechanistic insight into the basis of inhibition by the TDP-containing compounds through molecular docking onto the AtMCA-IIf crystal structure. Finally, a TDP-containing compound (TDP6) effectively hampered lateral root emergence in vivo, probably through inhibition of metacaspases specifically expressed in the endodermal cells overlying developing lateral root primordia. In the future, the small compound inhibitors and crystal structure of AtMCA-IIf can be used to study metacaspases in other species, such as important human pathogens, including those causing neglected diseases.

Protein 3D Hydration: A Case of Bovine Pancreatic Trypsin Inhibitor

Characterization of the hydrated state of a protein is crucial for understanding its structural stability and function. In the present study, we have investigated the 3D hydration structure of the protein BPTI (bovine pancreatic trypsin inhibitor) by molecular dynamics (MD) and the integral equation method in the three-dimensional reference interaction site model (3D-RISM) approach. Both methods have found a well-defined hydration layer around the protein and revealed the localization of BPTI buried water molecules corresponding to the X-ray crystallography data. Moreover, under 3D-RISM calculations, the obtained positions of waters bound firmly to the BPTI sites are in reasonable agreement with the experimental results mentioned above for the BPTI crystal form. The analysis of the 3D hydration structure (thickness of hydration shell and hydration numbers) was performed for the entire protein and its polar and non-polar parts using various cut-off distances taken from the literature as well as by a straightforward procedure proposed here for determining the thickness of the hydration layer. Using the thickness of the hydration shell from this procedure allows for calculating the total hydration number of biomolecules properly under both methods. Following this approach, we have obtained the thickness of the BPTI hydration layer of 3.6 Å with 369 water molecules in the case of MD simulation and 3.9 Å with 333 water molecules in the case of the 3D-RISM approach. The above procedure was also applied for a more detailed description of the BPTI hydration structure near the polar charged and uncharged radicals as well as non-polar radicals. The results presented for the BPTI as an example bring new knowledge to the understanding of protein hydration.

Structural and Biophysical Insights into SPINK1 Bound to Human Cationic Trypsin

(1) The serine protease inhibitor Kazal type 1 (SPINK1) inhibits trypsin activity in zymogen granules of pancreatic acinar cells. Several mutations in the SPINK1 gene are associated with acute recurrent pancreatitis (ARP) and chronic pancreatitis (CP). The most common variant is SPINK1 p.N34S. Although this mutation was identified two decades ago, the mechanism of action has remained elusive. (2) SPINK1 and human cationic trypsin (TRY1) were expressed in E. coli, and inhibitory activities were determined. Crystals of SPINK1-TRY1 complexes were grown by using the hanging-drop method, and phases were solved by molecular replacement. (3) Both SPINK1 variants show similar inhibitory behavior toward TRY1. The crystal structures are almost identical, with minor differences in the mutated loop. Both complexes show an unexpected rotamer conformation of the His63 residue in TRY1, which is a member of the catalytic triad. (4) The SPINK1 p.N34S mutation does not affect the inhibitory behavior or the overall structure of the protein. Therefore, the pathophysiological mechanism of action of the p.N34S variant cannot be explained mechanistically or structurally at the protein level. The observed histidine conformation is part of a mechanism for SPINK1 that can explain the exceptional proteolytic stability of this inhibitor.

KLK4 Inhibition by Cyclic and Acyclic Peptides: Structural and Dynamical Insights into Standard-Mechanism Protease Inhibitors

Sunflower trypsin inhibitor (SFTI-1) is a 14 amino acid serine protease inhibitor. The dual antiparallel β-sheet arrangement of SFTI-1 is stabilized by an N-terminal-C-terminal backbone cyclization and a further disulfide bridge to form a final bicyclic structure. This constrained structure is further rigidified by an extensive network of internal hydrogen bonds. Thus, the structure of SFTI-1 in solution resembles the protease-bound structure, reducing the entropic penalty upon protease binding. When cleaved at the scissile bond, it is thought that the rigidifying features of SFTI-1 maintain its structure, allowing the scissile bond to be reformed. The lack of structural plasticity for SFTI-1 is proposed to favor initial protease binding and continued occupancy in the protease active site, resulting in an equilibrium between the cleaved and uncleaved inhibitor in the presence of a protease. We have determined, at 1.15 Å resolution, the X-ray crystal structures of complexes between human kallikrein-related peptidase 4 (KLK4) and SFTI-FCQR(Asn14) and between KLK4 and an acyclic form of the same inhibitor, SFTI-FCQR(Asn14)[1,14], with the latter displaying a cleaved scissile bond. Structural analysis and MD simulations together reveal the roles of the altered contact sequence, intramolecular hydrogen bonding network, and backbone cyclization in altering the state of SFTI's scissile bond ligation at the protease active site. Taken together, the data presented reveal insights into the role of dynamics in the standard-mechanism inhibition and suggest that modifications on the non-contact strand may be a useful, underexplored approach for generating further potent or selective SFTI-based inhibitors against members of the serine protease family.

Molecular insight into chymotrypsin inhibitor 2 resisting proteolytic degradation

Chymotrypsin inhibitor 2 (CI2) is a special serine protease inhibitor which can resist hydrolysis for several days with a rapid equilibrium between the Michaelis complex and acyl-enzyme intermediate. The energies and conformational changes for subtilisin-catalyzed proteolysis of CI2 were examined in this paper for the first time by employing pseudo bond ab initio QM/MM MD simulations. In the acylation reaction, a low-barrier hydrogen bond between His64 and Asp32 in the transition state together with the lack of covalent backbone constraints makes the peptide bonds of CI2 break more easily than in other serine protease inhibitors. After acyl-enzyme formation, molecular dynamics simulations showed that the access of hydrolytic water to the active site requires partial dissociation of the leaving group. However, retention of the leaving group mainly by the intra- and inter-molecular H-bonding networks hinders the access of water and retards the deacylation reaction. Instead of the dissociation constant of inhibitors, we suggest employing the free energy at the acyl-enzyme state to predict the relative hydrolysis rates of CI2 mutants, which are testified by the experimental relative hydrolysis rates.

Analysis of Kunitz inhibitors from plants for comprehensive structural and functional insights

Legume Kunitz type trypsin inhibitor (KTI) family is one of the most versatile families of proteins. A typical KTI features a single peptide folded in β-trefoil manner, with the molecular weight about 20-22kDa and two disulphide bonds. The members are known to inhibit a wide range of serpins proteases at the same time many of them possess unique features. Copaifera langsdorffii Trypsin inhibitor (CTI) has a β-trefoil fold made up of two non-covalently bound polypeptide chains with only a single disulfide bridge. Delonix regia Trypsin inhibitor (DrTI) has one amino acid insertion between P1 and P2 of the reactive site distorting its conformation. Bauhinia bauhinioides Cruzipain inhibitor (BbCI) has a conservative β-trefoil fold but lacks disulfide bonds. Such subtle differences in structures make Kunitz inhibitors different from other inhibitor families. Most of the studies on these inhibitors are focused towards their proposed role in defense from insect pests and wounding but their exact physiological role in nature is still uncharted. Thus, it would be very interesting to closely analyze the structural details of these inhibitors in order to ascertain their biological role and other fascinating applications.

Soybean P34 Probable Thiol Protease Probably Has Proteolytic Activity on Oleosins

P34 probable thiol protease (P34) and Gly m Bd 30K (30K) show high relationship with the protease of 24 kDa oleosin of soybean oil bodies. In this study, 9 day germinated soybean was used to separate bioprocessed P34 (P32) from bioprocessed 30K (28K). Interestingly, P32 existed as dimer, whereas 28K existed as monomer; a P32-rich sample had proteolytic activity and high cleavage site specificity (Lys-Thr of 24 kDa oleosin), whereas a 28K-rich sample showed low proteolytic activity; the P32-rich sample contained one thiol protease. After mixing with purified oil bodies, all P32 dimers were dissociated and bound to 24 kDa oleosins to form P32-24 kDa oleosin complexes. By incubation, 24 kDa oleosin was preferentially hydrolyzed, and two hydrolyzed products (HPs; 17 and 7 kDa) were confirmed. After most of 24 kDa oleosin was hydrolyzed, some P32 existed as dimer, and the other as P32-17 kDa HP. It was suggested that P32 was the protease.

An Acrobatic Substrate Metamorphosis Reveals a Requirement for Substrate Conformational Dynamics in Trypsin Proteolysis

The molecular basis of enzyme catalytic power and specificity derives from dynamic interactions between enzyme and substrate during catalysis. Although considerable effort has been devoted to understanding how conformational dynamics within enzymes affect catalysis, the role of conformational dynamics within protein substrates has not been addressed. Here, we examine the importance of substrate dynamics in the cleavage of Kunitz-bovine pancreatic trypsin inhibitor protease inhibitors by mesotrypsin, finding that the varied conformational dynamics of structurally similar substrates can profoundly impact the rate of catalysis. A 1.4-Å crystal structure of a mesotrypsin-product complex formed with a rapidly cleaved substrate reveals a dramatic conformational change in the substrate upon proteolysis. By using long all-atom molecular dynamics simulations of acyl-enzyme intermediates with proteolysis rates spanning 3 orders of magnitude, we identify global and local dynamic features of substrates on the nanosecond-microsecond time scale that correlate with enzymatic rates and explain differential susceptibility to proteolysis. By integrating multiple enhanced sampling methods for molecular dynamics, we model a viable conformational pathway between substrate-like and product-like states, linking substrate dynamics on the nanosecond-microsecond time scale with large collective substrate motions on the much slower time scale of catalysis. Our findings implicate substrate flexibility as a critical determinant of catalysis.

Effect of ligand binding on the dynamics of trypsin. Comparison of different approaches

The intramolecular signal transduction induced by the binding of ligands to trypsin was investigated by molecular dynamics simulations. Ligand binding changes the residue-residue interaction energies and suppresses the mobility of loops that are in direct contact with the ligand. The reduced mobility of these loops results in the altered flexibility of the nearby loops and thereby transmits the information from ligand binding site to the remote sites. The analysis of the flexibility of all residues confirmed the coupling between loops L1 (185-188) and L2 (221-224) and the residues in the active center. The significance of S1 pocket residues for the signal transduction from the active center to the substrate-binding site was confirmed by the dynamical network and covariance matrix analyses. Gaussian network model and principal component analysis demonstrated that the active center residues had zero amplitude in the slowest fluctuations acting as hinges or anchors. Overall, our results provide a new insight into protein-ligand interactions and show how the allosteric signaling may occur.

Serine protease acylation proceeds with a subtle re-orientation of the histidine ring at the tetrahedral intermediate

The acylation mechanism of a prototypical serine protease trypsin and its complete free energy reaction profile have been determined by Born-Oppenheimer ab initio QM/MM molecular dynamics simulations with umbrella sampling.

Combined QM/MM mechanistic study of the acylation process in furin complexed with the H5N1 avian influenza virus hemagglutinin's cleavage site

Combined quantum mechanical/molecular mechanical (QM/MM) techniques have been applied to investigate the detailed reaction mechanism of the first step of the acylation process by furin in which the cleavage site of the highly pathogenic avian influenza virus subtype H5N1 (HPH5) acts as its substrate. The energy profile shows a simultaneous mechanism, known as a concerted reaction, of the two subprocesses: the proton transfer from Ser368 to His194 and the nucleophilic attack on the carbonyl carbon of the scissile peptide of the HPH5 cleavage site with a formation of tetrahedral intermediate (INT). The calculated energy barrier for this reaction is 16.2 kcal.mol(-1) at QM/MM B3LYP/6-31+G*//PM3-CHARMM22 level of theory. Once the reaction proceeds, the ordering of the electrostatic stabilization by protein environment is of the enzyme-substrate < transition state < INT complexes. Asp153 was found to play the most important role in the enzymatic reaction by providing the highest degree of intermediate complex stabilization. In addition, the negatively charged carbonyl oxygen of INT is well stabilized by the oxyanion hole constructed by Asn295's carboxamide and Ser368's backbone.

Structure of a serine protease poised to resynthesize a peptide bond

The serine proteases are among the most thoroughly studied enzymes, and numerous crystal structures representing the enzyme-substrate complex and intermediates in the hydrolysis reactions have been reported. Some aspects of the catalytic mechanism remain controversial, however, especially the role of conformational changes in the reaction. We describe here a high-resolution (1.46 A) crystal structure of a complex formed between a cleaved form of bovine pancreatic trypsin inhibitor (BPTI) and a catalytically inactive trypsin variant with the BPTI cleavage site ideally positioned in the active site for resynthesis of the peptide bond. This structure defines the positions of the newly generated amino and carboxyl groups following the 2 steps in the hydrolytic reaction. Comparison of this structure with those representing other intermediates in the reaction demonstrates that the residues of the catalytic triad are positioned to promote each step of both the forward and reverse reaction with remarkably little motion and with conservation of hydrogen-bonding interactions. The results also provide insights into the mechanism by which inhibitors like BPTI normally resist hydrolysis when bound to their target proteases.

Functional and structural roles of the Cys14-Cys38 disulfide of bovine pancreatic trypsin inhibitor

The disulfide bond between Cys14 and Cys38 of bovine pancreatic trypsin inhibitor lies on the surface of the inhibitor and forms part of the protease-binding region. The functional properties of three variants lacking this disulfide, with one or both of the Cys residues replaced with Ser, were examined, and X-ray crystal structures of the complexes with bovine trypsin were determined and refined to the 1.58-A resolution limit. The crystal structure of the complex formed with the mutant with both Cys residues replaced was nearly identical with that of the complex containing the wild-type protein, with the Ser oxygen atoms positioned to replace the disulfide bond with a hydrogen bond. The two structures of the complexes with single replacements displayed small local perturbations with alternate conformations of the Ser side chains. Despite the absence of the disulfide bond, the crystallographic temperature factors show no evidence of increased flexibility in the complexes with the mutant inhibitors. All three of the variants were cleaved by trypsin more rapidly than the wild-type inhibitor, by as much as 10,000-fold, indicating that the covalent constraint normally imposed by the disulfide contributes to the remarkable resistance to hydrolysis displayed by the wild-type protein. The rates of hydrolysis display an unusual dependence on pH over the range of 3.5-8.0, decreasing at the more alkaline values, as compared with the increased hydrolysis rates for normal substrates under these conditions. These observations can be accounted for by a model for inhibition in which an acyl-enzyme intermediate forms at a significant rate but is rapidly converted back to the enzyme-inhibitor complex by nucleophilic attack by the newly created amino group. The model suggests that a lack of flexibility in the acyl-enzyme intermediate, rather than the enzyme-inhibitor complex, may be a key factor in the ability of bovine pancreatic trypsin inhibitor and similar inhibitors to resist hydrolysis.

Analysis of catalytic mechanism of serine proteases. Viability of the ring-flip hypothesis

Quantum calculations are applied to the active site of serine proteases, including four specific residues and a water molecule, as well as a substrate and proton donors in the oxyanion hole. Residues are tethered to the protein backbone of an X-ray structure but otherwise allowed to move freely to their lowest energy positions. The viability of the ring-flip hypothesis, which proposes that a 180 degrees rotation of the His-57 imidazole ring facilitates the catalysis, is assessed by comparison of energies of configurations both before and after such a flip. Specifically considered is the contribution to catalysis of the Ser-214 residue and a water molecule that is observed in the active site. The calculations provide detailed information concerning the nature, geometry, and strength of hydrogen bonds that are formed within the active site at each stage of the enzymatic mechanism.

Rigidification of a flexible protease inhibitor variant upon binding to trypsin

The Tyr35-->Gly replacement in bovine pancreatic trypsin inhibitor (BPTI) has previously been shown to dramatically enhance the flexibility of the trypsin-binding region of the free inhibitor and to destabilize the interaction with the protease by about 3 kcal/mol. The effects of this replacement on the enzyme-inhibitor interaction were further studied here by X-ray crystallography and isothermal titration calorimetry (ITC). The co-crystal structure of Y35G BPTI bound to trypsin was determined using 1.65 A resolution X-ray diffraction data collected from cryopreserved crystals, and a new structure of the complex with wild-type BPTI under the same conditions was determined using 1.62 A data. These structures reveal that, in contrast to the free protein, Y35G BPTI adopts a conformation nearly identical with that of the wild-type protein, with a water-filled cavity in place of the missing Tyr side-chain. The crystallographic temperature factors for the two complexes indicate that the mutant inhibitor is nearly as rigid as the wild-type protein when bound to trypsin. Calorimetric measurements show that the change in enthalpy upon dissociation of the complex is 2.5 kcal/mol less favorable for the complex containing Y35G BPTI than for the complex with the wild-type inhibitor. Thus, the destabilization of the complex resulting from the Y35G replacement is due to a more favorable change in entropy upon dissociation. The heat capacity changes for dissociation of the mutant and wild-type complexes were very similar, suggesting that the entropic effects probably do not arise from solvation effects, but are more likely due to an increase in protein conformational entropy upon dissociation of the mutant inhibitor. These results define the biophysical role of a highly conserved core residue located outside of a protein-binding interface, demonstrating that Tyr35 has little impact on the trypsin-bound BPTI structure and acts primarily to define the structure of the free protein so as to maximize binding affinity.

Long dynamics simulations of proteins using atomistic force fields and a continuum representation of solvent effects: calculation of structural and dynamic properties

Long dynamics simulations were carried out on the B1 immunoglobulin-binding domain of streptococcal protein G (ProtG) and bovine pancreatic trypsin inhibitor (BPTI) using atomistic descriptions of the proteins and a continuum representation of solvent effects. To mimic frictional and random collision effects, Langevin dynamics (LD) were used. The main goal of the calculations was to explore the stability of tens-of-nanosecond trajectories as generated by this molecular mechanics approximation and to analyze in detail structural and dynamical properties. Conformational fluctuations, order parameters, cross correlation matrices, residue solvent accessibilities, pKa values of titratable groups, and hydrogen-bonding (HB) patterns were calculated from all of the trajectories and compared with available experimental data. The simulations comprised over 40 ns per trajectory for ProtG and over 30 ns per trajectory for BPTI. For comparison, explicit water molecular dynamics simulations (EW/MD) of 3 ns and 4 ns, respectively, were also carried out. Two continuum simulations were performed on each protein using the CHARMM program, one with the all-atom PAR22 representation of the protein force field (here referred to as PAR22/LD simulations) and the other with the modifications introduced by the recently developed CMAP potential (CMAP/LD simulations). The explicit solvent simulations were performed with PAR22 only. Solvent effects are described by a continuum model based on screened Coulomb potentials (SCP) reported earlier, i.e., the SCP-based implicit solvent model (SCP-ISM). For ProtG, both the PAR22/LD and the CMAP/LD 40-ns trajectories were stable, yielding C(alpha) root mean square deviations (RMSD) of about 1.0 and 0.8 A respectively along the entire simulation time, compared to 0.8 A for the EW/MD simulation. For BPTI, only the CMAP/LD trajectory was stable for the entire 30-ns simulation, with a C(alpha) RMSD of approximately 1.4 A, while the PAR22/LD trajectory became unstable early in the simulation, reaching a C(alpha) RMSD of about 2.7 A and remaining at this value until the end of the simulation; the C(alpha) RMSD of the EW/MD simulation was about 1.5 A. The source of the instabilities of the BPTI trajectories in the PAR22/LD simulations was explored by an analysis of the backbone torsion angles. To further validate the findings from this analysis of BPTI, a 35-ns SCP-ISM simulation of Ubiquitin (Ubq) was carried out. For this protein, the CMAP/LD simulation was stable for the entire simulation time (C(alpha) RMSD of approximately 1.0 A), while the PAR22/LD trajectory showed a trend similar to that in BPTI, reaching a C(alpha) RMSD of approximately 1.5 A at 7 ns. All the calculated properties were found to be in agreement with the corresponding experimental values, although local deviations were also observed. HB patterns were also well reproduced by all the continuum solvent simulations with the exception of solvent-exposed side chain-side chain (sc-sc) HB in ProtG, where several of the HB interactions observed in the crystal structure and in the EW/MD simulation were lost. The overall analysis reported in this work suggests that the combination of an atomistic representation of a protein with a CMAP/CHARMM force field and a continuum representation of solvent effects such as the SCP-ISM provides a good description of structural and dynamic properties obtained from long computer simulations. Although the SCP-ISM simulations (CMAP/LD) reported here were shown to be stable and the properties well reproduced, further refinement is needed to attain a level of accuracy suitable for more challenging biological applications, particularly the study of protein-protein interactions.

A comparison of staphostatin B with standard mechanism serine protease inhibitors

Staphostatins are the endogenous, highly specific inhibitors of staphopains, the major secreted cysteine proteases from Staphylococcus aureus. We have previously shown that staphostatins A and B are competitive, active site-directed inhibitors that span the active site clefts of their target proteases in the same orientation as substrates. We now report the crystal structure of staphostatin B in complex with wild-type staphopain B at 1.9 A resolution. In the complex structure, the catalytic residues are found in exactly the positions that would be expected for uncomplexed papain-type proteases. There is robust, continuous density for the staphostatin B binding loop and no indication for cleavage of the peptide bond that comes closest to the active site cysteine of staphopain B. The carbonyl carbon atom C of this peptide bond is 4.1 A away from the active site cysteine sulfur Sgamma atom. The carbonyl oxygen atom O of this peptide bond points away from the putative oxyanion hole and lies almost on a line from the Sgamma atom to the C atom. The arrangement is strikingly similar to the "ionmolecule" arrangement for the complex of papain-type enzymes with their substrates but differs significantly from the arrangement conventionally assumed for the Michaelis complex of papain-type enzymes with their substrates and also from the arrangement that is crystallographically observed for complexes of standard mechanism inhibitors and their target serine proteases.

Role of Asp102 in the catalytic relay system of serine proteases: a theoretical study

The role of Asp102 in the catalytic relay system of serine proteases is studied theoretically by calculating the free energy profiles of the single proton-transfer reaction by the Asn102 mutant trypsin and the concerted double proton-transfer reaction (so-called the charge-relay mechanism) of the wild-type trypsin. For each reaction, the reaction free energy profile of the rate-determining step (the tetrahedral intermediate formation step) is calculated by using ab initio QM/MM electronic structure calculations combined with molecular dynamics-free energy perturbation method. In the mutant reaction, the free energy monotonically increases along the reaction path. The rate-determining step of the mutant reaction is the formation of tetrahedral intermediate complex, not the base (His57) abstraction of the proton from Ser195. In contrast to the single proton-transfer reaction of the wild-type, MD simulations of the enzyme-substrate complex show that the catalytically favorable alignment of the relay system (the hydrogen bonding network between the mutant triad, His57, Asn102, and Ser195) is rarely observed even in the presence of a substrate at the active site. In the double proton-transfer reaction, the energy barrier is observed at the proton abstraction step, which corresponds to the rate-determining step of the single proton-transfer reaction of the wild-type. Although both reaction profiles show an increase of the activation barrier by several kcals/mol, these increases have different energetic origins: a large energetic loss of the electrostatic stabilization between His57 and Asn102 in the mutant reaction, while the lack of stabilization by the protein environment in the double proton-transfer reaction. Comparing the present results with the single proton transfer of the wild-type, Asp102 is proven to play two important roles in the catalytic process. One is to stabilize the protonated His57, or ionic intermediate, formed during the acylation, and the other is to fix the configuration around the active site, which is favorable to promote the catalytic process. These two factors are closely related to each other and are indispensable for the efficient catalysis. Also the present calculations suggest the importance of the remote site interaction between His57 and Val213-Ser214 at the catalytic transition state.

Theoretical perspectives on the reaction mechanism of serine proteases: the reaction free energy profiles of the acylation process

The reaction mechanism of serine proteases (trypsin), which catalyze peptide hydrolysis, is studied theoretically by ab initio QM/MM electronic structure calculations combined with Molecular Dynamics-Free Energy Perturbation calculations. We have calculated the entire reaction free energy profiles of the first reaction step of this enzyme (acylation process). The present calculations show that the rate-determining step of the acylation is the formation of the tetrahedral intermediate, and the breakdown of this intermediate has a small energy barrier. The calculated activation free energy for the acylation is approximately 17.8 kcal/mol at QM/MM MP2/(aug)-cc-pVDZ//HF/6-31(+)G/AMBER level, and this reaction is an exothermic process. MD simulations of the enzyme-substrate (ES) complex and the free enzyme in aqueous phase show that the substrate binding induces slight conformational changes around the active site, which favor the alignment of the reactive fragments (His57, Asp102, and Ser195) together in a reactive orientation. It is also shown that the proton transfer from Ser195 to His57 and the nucleophilic attack of Ser195 to the carbonyl carbon of the scissile bond of the substrate occur in a concerted manner. In this reaction, protein environment plays a crucial role to lowering the activation free energy by stabilizing the tetrahedral intermediate compared to the ES complex. The polarization energy calculations show that the enzyme active site is in a very polar environment because of the polar main chain contributions of protein. Also, the ground-state destabilization effect (steric strain) is not a major catalytic factor. The most important catalytic factor of stabilizing the tetrahedral intermediate is the electrostatic interaction between the active site and particular regions of protein: the main chain NH groups in Gly193 and Ser195 (so-called oxyanion hole region) stabilize negative charge generated on the carbonyl oxygen of the scissile bond, and the main chain carbonyl groups in Ile212 approximately Ser214 stabilize a positive charge generated on the imidazole ring of His57.

Identification of a Residue Critical for Maintaining the Functional Conformation of BPTI

The effects of amino acid replacements on the backbone dynamics of bovine pancreatic trypsin inhibitor (BPTI) were examined using 15N NMR relaxation experiments. Previous studies have shown that backbone amide groups within the trypsin-binding region of the wild-type protein undergo conformational exchange processes on the micros time scale, and that replacement of Tyr35 with Gly greatly increases the number of backbone atoms involved in such motions. In order to determine whether these mutational effects are specific to the replacement of this residue with Gly, six additional replacements were examined in the present study. In two of these, Tyr35 was replaced with either Ala or Leu, and the other four were single replacements of Tyr23, Phe33, Asn43 or Asn44, all of which are highly buried in the native structure and conserved in homologous proteins. The Y35A and Y35L mutants displayed dynamic properties very similar to those of the Y35G mutant, with the backbone segments including residues 10-19 and 32-44 undergoing motions revealed by enhanced 15N transverse relaxation rates. On the other hand, the Y23L, N43G and N44A substitutions caused almost no detectable changes in backbone dynamics, on either the ns-ps or ms-micros time scales, even though each of these replacements significantly destabilizes the native conformation. Replacement of Phe33 with Leu caused intermediate effects, with several residues that have previously been implicated in motions in the wild-type protein displaying enhanced transverse relaxation rates. These results demonstrate that destabilizing amino acid replacements can be accommodated in a native protein with dramatically different effects on conformational dynamics and that Tyr35 plays a particularly important role in defining the conformation of the trypsin-binding site of BPTI.

A substrate variant as a high-affinity, reversible inhibitor: insight from the X-ray structure of cilastatin bound to membrane dipeptidase

An analysis of the X-ray structure of cilastatin bound to membrane dipeptidase, together with docking studies, is presented here to reveal how a simple amide may act as a high-affinity, reversible, amidase inhibitor. Cilastatin binds as a normal substrate and is orientated in a perfect near-attack conformer for formation of a tetrahedral intermediate with the zinc-bound water/hydroxide. This intermediate is fated, however, only to revert to its starting components as scission of the amide bond is prevented by the precise fit of cilastatin within the active site. The cilastatin alkyl end groups that are tightly buttressed against amino acid residues on opposite sides of the active site, are aligned along the C-N reaction coordinate axis thereby preventing collapse of the intermediate via rupture of the C-N bond. Such a feature could have more general applicability in the explicit design of substrate variants as selective, tight-binding, and reversible inhibitors.

A clogged gutter mechanism for protease inhibitors

A classical peptide inhibitor of serine proteases that is hydrolyzed approximately 10(7) times more slowly than a good substrate is shown to form an acyl-enzyme intermediate rapidly. Despite this quick first step, further reaction is slowed dramatically because of tight and oriented binding of the cleaved peptide, preventing acyl-enzyme hydrolysis and favoring the reverse reaction. Moreover, this mechanism appears to be common to a large class of tight-binding serine protease inhibitors that mimic good substrates. The arrest of enzymatic reaction at the intermediate stage allowed us to determine that the consensus nucleophilic attack angle is close to 90 degrees in the reactive Michaelis complexes.

Quantum mechanical methods for enzyme kinetics

This review discusses methods for the incorporation of quantum mechanical effects into enzyme kinetics simulations in which the enzyme is an explicit part of the model. We emphasize three aspects: (a) use of quantum mechanical electronic structure methods such as molecular orbital theory and density functional theory, usually in conjunction with molecular mechanics; (b) treating vibrational motions quantum mechanically, either in an instantaneous harmonic approximation, or by path integrals, or by a three-dimensional wave function coupled to classical nuclear motion; (c) incorporation of multidimensional tunneling approximations into reaction rate calculations.

Correlative motions and memory effects in molecular dynamics simulations of molecules: principal components and rescaled range analysis suggest that the motions of native BPTI are more correlated than those of its mutants

In this work MD simulations of the native bovine pancreatic trypsin inhibitor (BPTI) and 16 mutants were done in vacuum in order to study memory effects in the mutants using principal component analysis (PCA) and the rescaled range analysis (Hurst exponents). Both PCA and the rescaled range analysis support our previous proposition, based on PCA of lysozyme, that the motions of a native protein are more correlated than those of mutants. The methods are compared, the nature and applications of the rule and the role of the long-range correlations in MD time series (i.e. memory) are discussed in the context of collective motions.

Simulating enzyme reactions: challenges and perspectives

Elucidating how enzymes enhance the rates of the reactions that they catalyze is a major goal of contemporary biochemistry, and it is an area in which computational and theoretical techniques can make a major contribution. This article outlines some of the processes that need to be investigated if enzyme catalysis is to be understood, reviews the current state-of-the-art in enzyme simulation work, and highlights challenges for the future.

Elucidating the nature of enzyme catalysis utilizing a new twist on an old methodology: quantum mechanical-free energy calculations on chemical reactions in enzymes and in aqueous solution

How do enzymes achieve very large rate enhancements compared to corresponding uncatalyzed reactions in solution? We present a computational approach which combines high-level ab initio quantum mechanical calculations with classical free energy calculations to address this question. Our calculations lead to accurate estimates of DeltaG for both trypsin and catechol O-methyltransferase-catalyzed and reference uncatalyzed reactions and give new insights into the nature of enzyme catalysis. The same methodology applied to steps in the catalytic mechanism of citrate synthase further supports the conclusion that one need not invoke special concepts such as "low-barrier hydrogen bonds" or "pK(a) matching" to explain enzyme catalysis.