Conformational changes in HIV-1 proteinase: effect of protonation of the active center on conformation of HIV-1 proteinase in water

At weak acidic pH, where HIV-1 proteinase is most stable and active, its catalytic Asp 25/25' dyad shares one proton. At a physiological pH the dyad is deprotonated, however, 2 ns molecular dynamics simulations of the HIV-1 protease with monoprotonated and deprotonated Asp25/25' dyad is performed, in order to investigate the influence of Asp25/25' protonation state on the proteinase dynamics. For net charge neutralization the 4 Cl- ions were included. In case of deprotonated active site the significant tertiary structure deviation of HIV-1 PR structure from crystal structure is observed, while in the monoprotonated one the tertiary structure fluctuates near starting structure. Possible mechanism of the influence of the Asp25/25' protonation state on proteinase dynamics is discussed.     

Population density analysis for determining the protonation state of the catalytic dyad in BACE1-tertiary carbinamine-based inhibitor complex

BACE1 is an aspartyl protease with a very relevant role in medicinal chemistry related to Alzheimer Disease since it has demonstrated to be a promising therapeutic target for inhibition and possible control for the progress of the peptide accumulation characteristic of this pathology. The enzymatic activity of this protein is given by the aspartic dyad, Asp93 and Asp289, which can adopt several protonation states depending on the chemical nature of its inhibitors, this is, monoprotonated, diprotonated and di-deprotonated states. In the present study, the analysis of the population density, for a series of protein-inhibitor molecular dynamics simulations, was carried out to identify the most feasible protonation state adopted by the catalytic dyad in the presence of tertiary carbinamine (TC) transition state analog inhibitors. The results revealed that the monoprotonated Asp289i state, in which the Asp93 and Asp289 residue side chains are deprotonated and protonated on the inner oxygen, respectively, is the most preferred in the presence of TC family inhibitors. This result was obtained after evaluating, for all 9 possible protonation state configurations, the individual and combined population densities of a set of parameters sensitive to protonation state of the Aspartic dyad, using an X-ray experimental BACE1/TC crystallographic structure as reference. This case study demonstrates again the usefulness of the concept of population density as a quantitative tool to establish the most stable system settings, among all possible, by measuring the level of occurrence of simultaneous events obtained from a sampling over time. These results will help to clear the phenomena related to the TCs inhibitory pathway, as well as assist in the design of better TC inhibitors against Alzheimer’s protease.     

Insights into conformational regulation of PfMATE transporter from Pyrococcus furiosus induced by alternating protonation state of Asp41 residue: A molecular dynamics simulation study

BackgroundMultidrug and toxic compound extrusion (MATE) family transporters induce multiple-drug resistance (MDR) of bacterial pathogens and cancer cells, thus causing critical reductions in the therapeutic efficacies of antibiotics and anti-cancer drugs. Unfortunately, to date, the details and intrinsic reason about conformational regulation mechanism of MATE transporters remain elusive.MethodIn this work, molecular dynamics (MD) simulations were conducted to explore the conformational regulation mechanism of PfMATE transporter from Pyrococcus furiosus based on different protonation state of Asp41. Two (MD) simulation systems were investigated: a system with protonation of Asp41 and a system without protonation of Asp41, which were named by D184(H)D41(H) system and D184(H) system, respectively.Results and conclusionsFirstly, MD simulation results indicate that conformational changes mainly happen in extracellular regions of PfMATE protein. Further analysis reveals that PfMATE protein experiences different motion mode and forms different conformation based on different protonation state of Asp41. In the D184(H)D41(H) system, PfMATE experiences an opening motion and forms a more outward-open conformation. As for the D184(H) system, the protein has an anticlockwise rotational motion with the channel axis of protein and the more outward-open conformation does not appear. It can be inferred that protonation of Asp41 is essential for conformational regulation of PfMATE during transporting substrates.General significanceThese findings provide intrinsic information for understanding the conformational regulation mechanism of PfMATE and will be very meaningful to explore the MDR mechanism of PfMATE further.     

Insights into the functional role of protonation states in the HIV-1 protease-BEA369 complex: molecular dynamics simulations and free energy calculations

The molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) method combined with molecular dynamics (MD) simulations were used to investigate the functional role of protonation in human immunodeficiency virus type 1 (HIV-1) protease complexed with the inhibitor BEA369. Our results demonstrate that protonation of two aspartic acids (Asp25/Asp25′) has a strong influence on the dynamics behavior of the complex, the binding free energy of BEA369, and inhibitor–residue interactions. Relative binding free energies calculated using the MM-PBSA method show that protonation of Asp25 results in the strongest binding of BEA369 to HIV-1 protease. Inhibitor–residue interactions computed by the theory of free energy decomposition also indicate that protonation of Asp25 has the most favorable effect on binding of BEA369. In addition, hydrogen-bond analysis based on the trajectories of the MD simulations shows that protonation of Asp25 strongly influences the water-mediated link of a conserved water molecule, Wat301. We expect that the results of this study will contribute significantly to binding calculations for BEA369, and to the design of high affinity inhibitors.     

Substrate-dependent transport mechanism in AcrB of multidrug resistant bacteria

The multidrug resistance (MDR) system effectively expels antibiotics out of bacteria causing serious issues during bacterial infection. In addition to drug, indole, a common metabolic waste of bacteria, is expelled by MDR system of gram-negative bacteria for their survival. Experimental results suggest that AcrB, one of the key components of MDR system, undergoes large scale conformation changes during the pumping due to proton-motive process. However, due to extremely short time scale, it is difficult to observe (experimentally) those changes in the AcrB, which might facilitate the pumping process. Molecular simulations can shed light to understand the conformational changes for transport of indole in AcrB. Examination of conformational changes using all-atom simulation is, however, impractical. Here, we develop a hybrid coarse-grained force field to study the conformational changes of AcrB in presence of indole in the porter domain of monomer II. Using the coarse-grained force field, we investigated the conformational changes of AcrB for a number of model systems considering the effect of protonation in aspartic acid (Asp) residues Asp407 and Asp408 in the transmembrane domain of monomer II. Our results show that in the presence of indole, protonation of Asp408 or Asp407 residue causes conformational changes from binding state to extrusion state in monomer II, while remaining two monomers (I and III) approach access state in AcrB protein. We also observed that all three AcrB monomers prefer to go back to access state in the absence of indole. Steered molecular dynamics simulations were performed to demonstrate the feasibility of indole transport mechanism for protonated systems. Identification of indole transport pathway through AcrB can be very helpful in understanding the drug efflux mechanism used by the MDR bacteria.     

Probing the Ser-Ser-Lys catalytic triad mechanism of peptide amidase: computational studies of the ground state, transition state, and intermediate

Peptide amidase (Pam), a hydrolytic enzyme that belongs to the amidase signature (AS) family, selectively catalyzes the hydrolysis of the C-terminal amide bond (CO-NH(2)) of peptides. The recent availability of the X-ray structures of Pam, fatty acid amide hydrolase, and malonamidase E2 has led to the proposal of a novel Ser-Ser-Lys catalytic triad mechanism for the amide hydrolysis by the AS enzymes. The molecular dynamics (MD) simulations using the CHARMM force field were performed to explore the catalytic mechanism of Pam. The 1.8 A X-ray crystal structure of Pam in complex with the amide analogue of chymostatin was chosen for the initial coordinates for the MD simulations. The five systems that were investigated are as follows: (i) enzyme.substrate with Lys123-NH(2), (ii) enzyme.substrate with Lys123-NH(3)(+), (iii) enzyme.substrate with Lys123-NH(3)(+) and Ser226-O(-), (iv) enzyme.transition state, and (v) enzyme.tetrahedral intermediate. Our data support the presence of the hydrogen bonding network among the catalytic triad residues, Ser226, Ser202, and Lys123, where Ser226 acts as the nucleophile and Ser202 bridges Ser226 and Lys123. The MD simulation supports the catalytic role of the crystallographic waters, Wat1 and Wat2. In all the systems that have been studied, the backbone amide nitrogens of Asp224 and Thr223 create an oxyanion hole by hydrogen bonding to the terminal amide oxygen of the substrate, and stabilize the oxyanion tetrahedral intermediate. The results from both our computational investigation and previously published experimental pH profile support two mechanisms. In a mechanism that is relevant at lower pH, the Lys123-NH(3)(+)-Ser202 dyad provides structural support to the catalytic residue Ser226, which in turn carries out a nucleophilic attack at the substrate amide carbonyl in concert with Wat1-mediated deprotonation and stabilization of the tetrahedral transition state by the oxyanion hole. In the mechanism operating at higher pH, the Lys123-NH(2)-Ser202 catalytic dyad acts as a general base to assist addition of Ser226 to the substrate amide carbonyl. The results from the MD simulation of the tetrahedral intermediate state show that both Ser202 and Lys123 are possible candidates for protonation of the leaving group, NH(2), to form the acyl-enzyme intermediate.     

Effects of protonation state of Asp181 and position of active site water molecules on the conformation of PTP1B

In protein tyrosine phosphatase 1B (PTP1B), the flexible WPD loop adopts a closed conformation (WPD_closed) in the active state of PTP1B, bringing the catalytic Asp181 close to the active site pocket, while WPD loop is in an open conformation (WPD_open) in the inactive state. Previous studies showed that Asp181 may be protonated at physiological pH, and ordered water molecules exist in the active site. In the current study, molecular dynamics simulations are employed at different Asp181 protonation states and initial positions of active site water molecules, and compared with the existing crystallographic data of PTP1B. In WPD_closed conformation, the active site is found to maintain its conformation only in the protonated state of Asp181 in both free and liganded states, while Asp181 is likely to be deprotonated in WPD_open conformation. When the active site water molecule network that is a part of the free WPD_closed crystal structure is disrupted, intermediate WPD loop conformations, similar to that in the PTPRR crystal structure, are sampled in the MD simulations. In liganded PTP1B, one active site water molecule is found to be important for facilitating the orientation of Cys215 and the phosphate ion, thus may play a role in the reaction. In conclusion, conformational stability of WPD loop, and possibly catalytic activity of PTP1B, is significantly affected by the protonation state of Asp181 and position of active site water molecules, showing that these aspects should be taken into consideration both in MD simulations and inhibitor design. © Proteins 2013. © 2012 Wiley Periodicals, Inc.     

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.     

The influence of amino acid protonation states on molecular dynamics simulations of the bacterial porin OmpF

Several groups, including our own, have found molecular dynamics (MD) calculations to result in the size of the pore of an outer membrane bacterial porin, OmpF, to be reduced relative to its size in the x-ray crystal structure. At the narrowest portion of its pore, loop L3 was found to move toward the opposite face of the pore, resulting in decreasing the cross-section area by a factor of approximately 2. In an earlier work, we computed the protonation states of titratable residues for this system and obtained values different from those that had been used in previous MD simulations. Here, we show that MD simulations carried out with these recently computed protonation states accurately reproduce the cross-sectional area profile of the channel lumen in agreement with the x-ray structure. Our calculations include the investigation of the effect of assigning different protonation states to the one residue, D(127), whose protonation state could not be modeled in our earlier calculation. We found that both assumptions of charge states for D(127) reproduced the lumen size profile of the x-ray structure. We also found that the charged state of D(127) had a higher degree of hydration and it induced greater mobility of polar side chains in its vicinity, indicating that the apparent polarizability of the D(127) microenvironment is a function of the D(127) protonation state.     

Insights into the structure and dynamics of the dinuclear zinc beta-lactamase site from Bacteroides fragilis

Herein, we report quantum chemical calculations and molecular dynamics (MD) simulations of the dinuclear form of the Bacteroides fragilis zinc beta-lactamase. We studied four different configurations which differ in the protonation state of the Asp103 residue and in the presence or absence of a Zn1-OH-Zn2 bridge. The flexibility of the Zn1-OH-Zn2 bridge was studied by means of quantum mechanical (QM) calculations on cluster models while the relative stabilities of the different configurations were estimated from QM linear scaling calculations on the enzyme. Contacts between important residues (Cys104, Asp69, Lys185, etc.), the solvation of the zinc ions, and the conformation of the active site beta-hairpin loop were characterized by the MD analyses. The influence of the buried sodium ion close to the Zn2 position was investigated by carrying out a secondary simulation where the sodium ion was replaced with an internal water molecule. The comparative structural analyses among the different MD trajectories augmented with energetic calculations have demonstrated that the B. fragilis protein efficiently binds the internal Na(+) ion observed crystallographically. Moreover, we found that when Asp103 is unprotonated, a rigid Zn1-OH-Zn2 bridge results, while for neutral Asp103, a fluctuating Zn1-Zn2 distance was possible via the breaking and formation of the Zn1-OH-Zn2 bridge. The mechanistic implications of these observations are discussed in detail.     

Computational analysis of aspartic protease plasmepsin II complexed with EH58 inhibitor: a QM/MM MD study

Plasmepsin (PM) II is one of four enzymes in the food vacuole of Plasmodium falciparum. It has become an attractive target for combating malaria through research regarding its importance in the P. falciparum metabolism and life cycle, making it the target of choice for structure-based drug design. This paper reports the results of hybrid quantum mechanics / molecular mechanics (QM/MM) molecular dynamics (MD) simulations employed to study the details of the interactions established between PM II and N-(3-{(2-benzo[1, 3]dioxol-5-yl-ethyl)[3-(1-methyl-3-oxo-1,3-dihydro-isoindol-2-yl) propionyl]-amino}-1-benzyl-2-(hydroxyl-propyl)-4-benzyloxy-3,5dimethoxy-benzamide (EH58), a well-known potent inhibitor for this enzyme. Electrostatic binding free energy and energy terms decomposition have been computed for PM II complexed with the EH58 inhibitor. The results reveal that there is a strong interaction between Asp34, Val78, Ser79, Tyr192 and Asp214 residues and the EH58 inhibitor. In addition, we have computed the potential of the mean force (PMF) profile in order to assign the protonation state of the two catalytic aspartates in PM II-EH58 complex. The results indicate that the protonation of Asp214 favors a stable active site structure, which is consistent with our electrostatic binding free energy calculation and with previous published works.     

Protonation state of Asp30 exerts crucial influence over surface loop rearrangements responsible for NO release in nitrophorin 4

pK(a) values of ionizable residues were calculated for the crystal structures describing the pH and NO binding dependant conformations of nitrophorin 4, a pH sensitive NO carrier heme protein. Comparison of resultant H-bonding patterns allowed the identification of the amino acids that take part in signaling pH change. We carried out MD simulations to show that the protonation state of Asp30, buried in the closed conformation, is crucial for maintaining the tight packed conformation of the closed form of the complex - presenting a model for the functional decrease of NO binding affinity of nitrophorins at physiological pH.     

Protonation states of important acidic residues in the central Ca²⁺ ion binding sites of the Ca²⁺-ATPase: a molecular modeling study

The P-type ATPases are responsible for the transport of cations across cell membranes. The sarco(endo)plasmic reticulum Ca2?-ATPase (SERCA) transports two Ca2? ions from the cytoplasm to the lumen of the sarco(endo)plasmic reticulum and countertransports two or three protons per catalytic cycle. Two binding sites for Ca2? ions have been located via protein crystallography, including four acidic amino acid residues that are essential to the ion coordination. In this study, we present molecular dynamics (MD) simulations examining the protonation states of these amino acid residues in a Ca2?-free conformation of SERCA. Such knowledge will be important for an improved understanding of atomistic details of the transport mechanism of protons and Ca2? ions. Eight combinations of the protonation of four central acidic residues, Glu309, Glu771, Asp800, and Glu908, are tested from 10 ns MD simulations with respect to protein stability and ability to maintain a structure similar to the crystal structure. The trajectories for the most prospective combinations of protonation states were elongated to 50 ns and subjected to more detailed analysis, including prediction of pK(a) values of the four acidic residues over the trajectories. From the simulations we find that the combination leaving only Asp800 as charged is most likely. The results are compared to available experimental data and explain the observed destabilization upon full deprotonation, resulting in the entry of cytoplasmic K? ions into the Ca2? binding sites during the simulation in which Ca2? ions are absent. Furthermore, a hypothesis for the exchange of protons from the central binding cavity is proposed.     

A molecular dynamics study of the structural stability of HIV-1 protease under physiological conditions: the role of Na+ ions in stabilizing the active site

HIV-1 protease is most active under weakly acidic conditions (pH 3.5-6.5), when the catalytic Asp25 and Asp25' residues share 1 proton. At neutral pH, this proton is lost and the stability of the structure is reduced. Here we present an investigation of the effect of pH on the dynamics of HIV-1 protease using MD simulation techniques. MD simulations of the solvated HIV-1 protease with the Asp25/25' residues monoprotonated and deprotonated have been performed. In addition we investigated the effect of the inclusion of Na(+) and Cl(-) ions to mimic physiological salt conditions. The simulations of the monoprotonated form and deprotonated form including Na(+) show very similar behavior. In both cases the protein remained stable in the compact, "self-blocked" conformation in which the active site is blocked by the tips of the flaps. In the deprotonated system a Na(+) ion binds tightly to the catalytic dyad shielding the repulsion between the COO(-) groups. Ab initio calculations also suggest the geometry of the active site with the Na(+) bound closely resembles that of the monoprotonated case. In the simulations of the deprotonated form (without Na(+) ions), a water molecule bound between the Asp25 Asp25' side-chains. This disrupted the dimerization interface and eventually led to a fully open conformation.     

Conformational sampling and polarization of Asp26 in pKa calculations of thioredoxin

Thioredoxin is a protein that has been used as model system by various computational methods to predict the pK_a of aspartate residue Asp26 which is 3.5 units higher than a solvent exposed one (eg, Asp20). Here, we use extensive atomistic molecular dynamics simulations of two different protonation states of Asp26 in combination with conformational analysis based on RMSD clustering and principle component analysis to identify representative conformations of the protein in solution. For each conformation, the Gibbs free energy of proton transfer between Asp26 and Asp20, which is fully solvated in a loop region of the protein, is calculated with the Amber99sb force field in alchemical transformations. The varying polarization of the two residues in different molecular environments and protonation states is described by Hirshfeld-I (HI) atomic charges obtained from the averaged polarized electron density. Our results show that the Gibbs free energy of proton transfer is dependent on the protein conformation, the proper sampling of the neighboring Lys57 residue orientations and on water molecules entering the hydrophobic cavity upon deprotonating Asp26. The inclusion of the polarization of both aspartate residues in the free energy cycle by HI atomic charges corrects the results from the non-polarizable force field and reproduces the experimental ΔpK_a value of Asp26.     

His ... Asp catalytic dyad of ribonuclease A: histidine pKa values in the wild-type, D121N, and D121A enzymes

Bovine pancreatic ribonuclease A (RNase A) has a conserved His ... Asp catalytic dyad in its active site. Structural analyses had indicated that Asp121 forms a hydrogen bond with His119, which serves as an acid during catalysis of RNA cleavage. The enzyme contains three other histidine residues including His12, which is also in the active site. Here, 1H-NMR spectra of wild-type RNase A and the D121N and D121A variants were analyzed thoroughly as a function of pH. The effect of replacing Asp121 on the microscopic pKa values of the histidine residues is modest: none change by more than 0.2 units. There is no evidence for the formation of a low-barrier hydrogen bond between His119 and either an aspartate or an asparagine residue at position 121. In the presence of the reaction product, uridine 3'-phosphate (3'-UMP), protonation of one active-site histidine residue favors protonation of the other. This finding is consistent with the phosphoryl group of 3'-UMP interacting more strongly with the two active-site histidine residues when both are protonated. Comparison of the titration curves of the unliganded enzyme with that obtained in the presence of different concentrations of 3'-UMP shows that a second molecule of 3'-UMP can bind to the enzyme. Together, the data indicate that the aspartate residue in the His ... Asp catalytic dyad of RNase A has a measurable but modest effect on the ionization of the adjacent histidine residue.     

Determinants of reactivity and selectivity in soluble epoxide hydrolase from quantum mechanics/molecular mechanics modeling

Soluble epoxide hydrolase (sEH) is an enzyme involved in drug metabolism that catalyzes the hydrolysis of epoxides to form their corresponding diols. sEH has a broad substrate range and shows high regio- and enantioselectivity for nucleophilic ring opening by Asp333. Epoxide hydrolases therefore have potential synthetic applications. We have used combined quantum mechanics/molecular mechanics (QM/MM) umbrella sampling molecular dynamics (MD) simulations (at the AM1/CHARMM22 level) and high-level ab initio (SCS-MP2) QM/MM calculations to analyze the reactions, and determinants of selectivity, for two substrates: trans-stilbene oxide (t-SO) and trans-diphenylpropene oxide (t-DPPO). The calculated free energy barriers from the QM/MM (AM1/CHARMM22) umbrella sampling MD simulations show a lower barrier for phenyl attack in t-DPPO, compared with that for benzylic attack, in agreement with experiment. Activation barriers in agreement with experimental rate constants are obtained only with the highest level of QM theory (SCS-MP2) used. Our results show that the selectivity of the ring-opening reaction is influenced by several factors, including proximity to the nucleophile, electronic stabilization of the transition state, and hydrogen bonding to two active site tyrosine residues. The protonation state of His523 during nucleophilic attack has also been investigated, and our results show that the protonated form is most consistent with experimental findings. The work presented here illustrates how determinants of selectivity can be identified from QM/MM simulations. These insights may also provide useful information for the design of novel catalysts for use in the synthesis of enantiopure compounds.     

Ensemble MD simulations restrained via crystallographic data: Accurate structure leads to accurate dynamics

Currently, the best existing molecular dynamics (MD) force fields cannot accurately reproduce the global free‐energy minimum which realizes the experimental protein structure. As a result, long MD trajectories tend to drift away from the starting coordinates (e.g., crystallographic structures). To address this problem, we have devised a new simulation strategy aimed at protein crystals. An MD simulation of protein crystal is essentially an ensemble simulation involving multiple protein molecules in a crystal unit cell (or a block of unit cells). To ensure that average protein coordinates remain correct during the simulation, we introduced crystallography‐based restraints into the MD protocol. Because these restraints are aimed at the ensemble‐average structure, they have only minimal impact on conformational dynamics of the individual protein molecules. So long as the average structure remains reasonable, the proteins move in a native‐like fashion as dictated by the original force field. To validate this approach, we have used the data from solid‐state NMR spectroscopy, which is the orthogonal experimental technique uniquely sensitive to protein local dynamics. The new method has been tested on the well‐established model protein, ubiquitin. The ensemble‐restrained MD simulations produced lower crystallographic R factors than conventional simulations; they also led to more accurate predictions for crystallographic temperature factors, solid‐state chemical shifts, and backbone order parameters. The predictions for ¹⁵N relaxation rates are at least as accurate as those obtained from conventional simulations. Taken together, these results suggest that the presented trajectories may be among the most realistic protein MD simulations ever reported. In this context, the ensemble restraints based on high‐resolution crystallographic data can be viewed as protein‐specific empirical corrections to the standard force fields.     

Caught in the Act: the 1.5 A resolution crystal structures of the HIV-1 protease and the I54V mutant reveal a tetrahedral reaction intermediate

HIV-1 protease (PR) is the target for several important antiviral drugs used in AIDS therapy. The drugs bind inside the active site cavity of PR where normally the viral polyprotein substrate is bound and hydrolyzed. We report two high-resolution crystal structures of wild-type PR (PRWT) and the multi-drug-resistant variant with the I54V mutation (PRI54V) in complex with a peptide at 1.46 and 1.50 A resolution, respectively. The peptide forms a gem-diol tetrahedral reaction intermediate (TI) in the crystal structures. Distinctive interactions are observed for the TI binding in the active site cavity of PRWT and PRI54V. The mutant PRI54V/TI complex has lost water-mediated hydrogen bond interactions with the amides of Ile50 and Ile50' in the flap. Hence, the structures provide insight into the mechanism of drug resistance arising from this mutation. The structures also illustrate an intermediate state in the hydrolysis reaction. One of the gem-diol hydroxide groups in the PRWT complex forms a very short (2.3 A) hydrogen bond with the outer carboxylate oxygen of Asp25. Quantum chemical calculations based on this TI structure are consistent with protonation of the inner carboxylate oxygen of Asp25', in contrast to several theoretical studies. These TI complexes and quantum calculations are discussed in relation to the chemical mechanism of the peptide bond hydrolysis catalyzed by PR.