Use of nitrogen-15 kinetic isotope effects to elucidate details of the chemical mechanism of human immunodeficiency virus 1 protease

Ten-Fold Solvent Kinetic Isotope Effect for the Nonradiative Relaxation of the Aqueous Ferrate(VI) Ion

Hypervalent iron intermediates have been invoked in the catalytic cycles of many metalloproteins, and thus, it is crucial to understand how the coupling between such species and their environment can impact their chemical and physical properties in such contexts. In this work, we take advantage of the solvent kinetic isotope effect (SKIE) to gain insight into the nonradiative deactivation of electronic excited states of the aqueous ferrate(VI) ion. We observe an exceptionally large SKIE of 9.7 for the nanosecond-scale relaxation of the lowest energy triplet ligand field state to the ground state. Proton inventory studies demonstrate that a single solvent O-H bond is coupled to the ion during deactivation, likely due to the sparse vibrational structure of ferrate(VI). Such a mechanism is consistent with that reported for the deactivation of f-f excited states of aqueous trivalent lanthanides, which exhibit comparably large SKIE values. This phenomenon is ascribed entirely to dissipation of energy into a higher overtone of a solvent acceptor mode, as any impact on the apparent relaxation rate due to a change in solvent viscosity is negligible.

Catalytic Mechanism of Human T-Cell Leukemia Virus Type 1 Protease Investigated by Combined QM/MM Molecular Dynamics Simulations

Combined quantum mechanical and molecular mechanical (QM/MM) molecular dynamics simulations were performed to investigate the catalytic mechanism of human T-cell leukemia virus type 1 (HTLV-1) protease, a retroviral aspartic protease that is a potential therapeutic target for curing HTLV-1-associated diseases. To elucidate the proteolytic cleavage mechanism, we determined the two-dimensional free energy surfaces of the HTLV-1 protease-catalyzed reactions through various possible pathways. The free energy simulations suggest that the catalytic reactions of the HTLV-1 protease occur in the following sequential steps: (1) a proton is transferred from the lytic water to Asp32', followed by the nucleophilic addition of the resulting hydroxyl to the carbonyl carbon of the scissile bond, forming a tetrahedral oxyanion intermediate, and (2) a proton is transferred from Asp32 to the peptide nitrogen of the scissile bond, leading to the spontaneous breakage of the scissile bond. The rate-limiting step of this catalytic process is the proton transfer from Asp32 to the peptide nitrogen of the scissile bond, with a free energy of activation of 21.1 kcal/mol. This free energy barrier is close to the experimentally determined free energy of activation (16.3 kcal/mol) calculated from the measured catalytic rate constant (k_cat). This mechanistic study provides detailed dynamic and structural information that will facilitate the design of mechanism-based inhibitors for the treatment of HTLV-1-associated diseases.

Exploring the concerted mechanistic pathway for HIV-1 PR-substrate revealed by umbrella sampling simulation

HIV-1 protease (HIV-1 PR) is an essential enzyme for the replication process of its virus, and therefore considered an important target for the development of drugs against the acquired immunodeficiency syndrome (AIDS). Our previous study shows that the catalytic mechanism of subtype B/C-SA HIV-1 PR follows a one-step concerted acyclic hydrolysis reaction process using a two-layered ONIOM B3LYP/6-31++G(d,p) method. This present work is aimed at exploring the proposed mechanism of the proteolysis catalyzed by HIV-1 PR and to ensure our proposed mechanism is not an artefact of a single theoretical technique. Hence, we present umbrella sampling method that is suitable for calculating potential mean force (PMF) for non-covalent ligand/substrate-enzyme association/dissociation interactions which provide thermodynamic details for molecular recognition. The free activation energy results were computed in terms of PMF analysis within the hybrid QM(DFTB)/MM approach. The theoretical findings suggest that the proposed mechanism corresponds in principle with experimental data. Given our observations, we suggest that the QM/MM MD method can be used as a reliable computational technique to rationalize lead compounds against specific targets such as the HIV-1 protease.

Visualizing Tetrahedral Oxyanion Bound in HIV-1 Protease Using Neutrons: Implications for the Catalytic Mechanism and Drug Design

HIV-1 protease is indispensable for virus propagation and an important therapeutic target for antiviral inhibitors to treat AIDS. As such inhibitors are transition-state mimics, a detailed understanding of the enzyme mechanism is crucial for the development of better anti-HIV drugs. Here, we used room-temperature joint X-ray/neutron crystallography to directly visualize hydrogen atoms and map hydrogen bonding interactions in a protease complex with peptidomimetic inhibitor KVS-1 containing a reactive nonhydrolyzable ketomethylene isostere, which, upon reacting with the catalytic water molecule, is converted into a tetrahedral intermediate state, KVS-1_TI. We unambiguously determined that the resulting tetrahedral intermediate is an oxyanion, rather than the gem-diol, and both catalytic aspartic acid residues are protonated. The oxyanion tetrahedral intermediate appears to be unstable, even though the negative charge on the oxyanion is delocalized through a strong n → π* hyperconjugative interaction into the nearby peptidic carbonyl group of the inhibitor. To better understand the influence of the ketomethylene isostere as a protease inhibitor, we have also examined the protease structure and binding affinity with keto-darunavir (keto-DRV), which similar to KVS-1 includes the ketomethylene isostere. We show that keto-DRV is a significantly less potent protease inhibitor than DRV. These findings shed light on the reaction mechanism of peptide hydrolysis catalyzed by HIV-1 protease and provide valuable insights into further improvements in the design of protease inhibitors.

Novel radial distribution function approach in the study of point mutations: the HIV-1 protease case study

Background: Mutations are one of the engines of evolution. Under constant stress pressure, mutations can lead to the emergence of unwanted, drug-resistant entities. Methodology: The radial distribution function weighted by the number of valence shell electrons is used to design quantitative structure-activity relationship (QSAR) model relating descriptors with the inhibition constant for a series of wild-type HIV-1 protease inhibitor complexes. The residuals of complexes with mutant HIV-1 protease were correlated with the energy of the highest occupied molecular orbitals of the residues introduced to enzyme via point mutations. Conclusion: Successful identification of residues Ile3, Asp25, Val32 and Ile50 as the one whose substitution influences the inhibition constant the most, demonstrates the potential of the proposed methodology for the study of the effects of point mutations.

The Influence of Hydrogen Atoms on the Performance of Radial Distribution Function-Based Descriptors in the Chemoinformatic Studies of HIV-1 Protease Complexes with Inhibitors

AIMS

The aim of this letter is to explore the influence of adding hydrogen atoms to the crystallographic structures of HIV-1 protease complexes with a series of inhibitors on the performance of radial distribution function based descriptors recently introduced in chemoinformatic studies.

BACKGROUND

Quite recently the successful application of molecular descriptors based on a radial distribution function to correlate it with biologically interesting properties of a ligand - enzyme complex was demonstrated. Except its predictive power, the analysis of atoms with dominant contributions to the RDFs can be used to identify relevant atoms and interactions. Since original paper was published on dataset consisting of the X-ray structures of complexes without hydrogen atoms, we wonder weather addition of light atoms can provide us new piece of information.

OBJECTIVE

The primarily objective is to create the model correlating the RDF based descriptors and physicochemical properties of the HIV-1 protease complexes with inhibitors with hydrogen atoms. Then, we will compare the performance of new model with previous one, where the hydrogen atoms were discarded. Information about interactions between the enzyme and the inhibitors will be extracted from the analysis of the RDF.

METHODS

The radial distribution function descriptor weighted by the number of valence shell electrons has proven to be sensitive to the changes in the structure of the enzyme and enzyme-ligand complexes. For each structure in our data set, RDF will be calculated and using multiple linear regression method the mathematical model will be designed correlating RDF based descriptors and the physicochemical properties. Statistical analysis of the atom's contribution to the total RDF will reveal relevant interactions.

RESULTS

The applicability of RDF based descriptor for the correlation with pKi and EC50 values is demonstrated, while simple models containing only two or three parameters are able to explain 78 and 86 % of the variance, respectively. The models with explicitly included hydrogens are of comparable quality with the previous models without hydrogens. The analysis of the atom's dominant contributions highlighted the importance of the hydroxyl groups of the inhibitor near the Asp25 and Asp25' residues when it is bounded to the protease.

CONCLUSION

Models based on the RDF weighted by the number of valence shell electrons for correlating small number of molecular descriptors and physicocehmical properties for structures with and without hydrogens are of comparable quality and both can be used for identification of relevant functional groups and interactions. Other: Our approach can be integrated to the next generation virtual screening methods, because is fast, reliable with high predictability potential.

Using the Relative Energy Gradient Method with Interacting Quantum Atoms to Determine the Reaction Mechanism and Catalytic Effects in the Peptide Hydrolysis in HIV-1 Protease

The reaction mechanism in an active site is of the utmost importance when trying to understand the role that an enzyme plays in biological processes. In a recently published paper [Theor. Chem. Acc. 2017, 136, 86], we formalised the Relative Energy Gradient (REG) method for automating an Interacting Quantum Atoms (IQA) analysis. Here, the REG method is utilised to determine the mechanism of peptide hydrolysis in the aspartic active site of the enzyme HIV-1 Protease. Using the REG method along with the IQA approach we determine the mechanism of peptide hydrolysis without employing any arbitrary parameters and with remarkable ease (albeit at large computational cost: the system contains 133 atoms, which means that there are 17 689 individual IQA terms to be calculated). When REG and IQA work together it is possible to determine a reaction mechanism at atomistic resolution from data directly derived from quantum calculations, without arbitrary parameters. Moreover, the mechanism determined by this novel method gives concrete insight into how the active site residues catalyse peptide hydrolysis.

Dynamic and Electrostatic Effects on the Reaction Catalyzed by HIV-1 Protease

HIV-1 Protease (HIV-1 PR) is one of the three enzymes essential for the replication process of HIV-1 virus, which explains why it has been the main target for design of drugs against acquired immunodeficiency syndrome (AIDS). This work is focused on exploring the proteolysis reaction catalyzed by HIV-1 PR, with special attention to the dynamic and electrostatic effects governing its catalytic power. Free energy surfaces for all possible mechanisms have been computed in terms of potentials of mean force (PMFs) within hybrid QM/MM potentials, with the QM subset of atoms described at semiempirical (AM1) and DFT (M06-2X) level. The results suggest that the most favorable reaction mechanism involves formation of a gem-diol intermediate, whose decomposition into the product complex would correspond to the rate-limiting step. The agreement between the activation free energy of this step with experimental data, as well as kinetic isotope effects (KIEs), supports this prediction. The role of the protein dynamic was studied by protein isotope labeling in the framework of the Variational Transition State Theory. The predicted enzyme KIEs, also very close to the values measured experimentally, reveal a measurable but small dynamic effect. Our calculations show how the contribution of dynamic effects to the effective activation free energy appears to be below 1 kcal·mol^-1. On the contrary, the electric field created by the protein in the active site of the enzyme emerges as being critical for the electronic reorganization required during the reaction. These electrostatic properties of the active site could be used as a mold for future drug design.

Metabolic origin of δ15 N values in nitrogenous compounds from Brassica napus L. leaves

Nitrogen isotope composition (δ(15) N) in plant organic matter is currently used as a natural tracer of nitrogen acquisition efficiency. However, the δ(15) N value of whole leaf material does not properly reflect the way in which N is assimilated because isotope fractionations along metabolic reactions may cause substantial differences among leaf compounds. In other words, any change in metabolic composition or allocation pattern may cause undesirable variability in leaf δ(15) N. Here, we investigated the δ(15) N in different leaf fractions and individual metabolites from rapeseed (Brassica napus) leaves. We show that there were substantial differences in δ(15) N between nitrogenous compounds (up to 30‰) and the content in ((15) N enriched) nitrate had a clear influence on leaf δ(15) N. Using a simple steady-state model of day metabolism, we suggest that the δ(15) N value in major amino acids was mostly explained by isotope fractionation associated with isotope effects on enzyme-catalysed reactions in primary nitrogen metabolism. δ(15) N values were further influenced by light versus dark conditions and the probable occurrence of alternative biosynthetic pathways. We conclude that both biochemical pathways (that fractionate between isotopes) and nitrogen sources (used for amino acid production) should be considered when interpreting the δ(15) N value of leaf nitrogenous compounds.

Transition states of native and drug-resistant HIV-1 protease are the same

HIV-1 protease is an important target for the treatment of HIV/AIDS. However, drug resistance is a persistent problem and new inhibitors are needed. An approach toward understanding enzyme chemistry, the basis of drug resistance, and the design of powerful inhibitors is to establish the structure of enzymatic transition states. Enzymatic transition structures can be established by matching experimental kinetic isotope effects (KIEs) with theoretical predictions. However, the HIV-1 protease transition state has not been previously resolved using these methods. We have measured primary (14)C and (15)N KIEs and secondary (3)H and (18)O KIEs for native and multidrug-resistant HIV-1 protease (I84V). We observed (14)C KIEs ((14)V/K) of 1.029 ± 0.003 and 1.025 ± 0.005, (15)N KIEs ((15)V/K) of 0.987 ± 0.004 and 0.989 ± 0.003, (18)O KIEs ((18)V/K) of 0.999 ± 0.003 and 0.993 ± 0.003, and (3)H KIEs ((3)V/K) KIEs of 0.968 ± 0.001 and 0.976 ± 0.001 for the native and I84V enzyme, respectively. The chemical reaction involves nucleophilic water attack at the carbonyl carbon, proton transfer to the amide nitrogen leaving group, and C-N bond cleavage. A transition structure consistent with the KIE values involves proton transfer from the active site Asp-125 (1.32 Å) with partial hydrogen bond formation to the accepting nitrogen (1.20 Å) and partial bond loss from the carbonyl carbon to the amide leaving group (1.52 Å). The KIEs measured for the native and I84V enzyme indicate nearly identical transition states, implying that a true transition-state analogue should be effective against both enzymes.

Validating the vitality strategy for fighting drug resistance

The current challenge in designing effective drugs against HIV-1 is to find novel candidates with high potency, but with a lower susceptibility to mutations associated with drug resistance. Trying to address this challenge, we developed in our previous study (Ishikita and Warshel, Angew Chem Int Ed Engl 2008; 47:697-700) a novel computational strategy for fighting drug resistance by predicting the likely moves of the virus through constraints on binding and catalysis. This has been based on calculating the ratio between the vitality values ((K(i) k(cat)/K(M))(mutant)/(K(i) k(cat)/K(M))(wild-type)) and using it as a guide for predicting the moves of the virus. The corresponding calculations of the binding affinity, K(i), were carried out using the semi-macroscopic version of the protein dipole Langevin dipole (PDLD/S) in its linear response approximation (LRA) in its β version (PDLD/S-LRA/β). We also calculate the proteolytic efficiency, k(cat)/K(M), by evaluating the transition state (TS) binding free energies using the PDLD/S-LRA/β method. Here we provide an extensive validation of our strategy by calculating the vitality of six existing clinical and experimental drug candidates. It is found that the computationally determined vitalities correlate reasonably well with those derived from the corresponding experimental data. This indicates that the calculated vitality may be used to identify mutations that would be most effective for the survival of the virus. Thus, it should be possible to use our approach in screening for mutations that would provide the most effective resistance to any proposed antiviral drug. This ability should be very useful in guiding the design of drug molecules that will lead to the slowest resistance.

HYBRID QUANTUM MECHANICAL/MOLECULAR MECHANICAL APPROACH TO ENZYMATIC REACTIONS BY UTILIZING THE REAL-SPACE GRID TECHNIQUE

We have developed a novel quantum mechanical/molecular mechanical (QM/MM) code based on the real-space grids in order to realize high parallel efficiency. The details of the methodology and its parallel implementation have been presented. We have computed the electronic state of the QM subsystem using the Kohn–Sham density functional theory, where the one-electron wave functions have been expressed by the real-space grids distributed over a cubic cell. We have performed QM/MM simulations for the peptide hydrolysis in human immunodeficiency virus type-1 aspartyl protease in order to examine the reliability of the present QM/MM approach. The activation energy obtained by the present calculations shows a good agreement with the experimental results and that of the other QM/MM method. Finally, we have parallelized the whole code and found that the grid approach can afford high parallel efficiency (~80%) in such a large scale electronic structure calculation. We conclude that the QM/MM approach utilizing real-space grids is adequate and efficient for the study of the enzymatic reactions.

Computational insights into the development of novel therapeutic strategies for Alzheimer's disease

BACKGROUND

β-amyloidosis and oxidative stress have been implicated as root causes of Alzheimer's disease (AD). Current potential therapeutic strategies for the treatment of AD include inhibition of amyloid β (Aβ) production, stimulation of Aβ degradation and prevention of Aβ oligomerization. However, efforts in this direction are hindered by the lack of understanding of the biochemical processes occurring at the atomic level in AD.

DISCUSSION

A radically different approach to achieve this goal would be the application of comprehensive theoretical and computational techniques such as molecular dynamics, quantum mechanics, hybrid quantum mechanics/molecular mechanics, bioinformatics and rotational spectroscopy to investigate complex chemical and physical processes in β-amyloidosis and the oxidative stress mechanism.

CONCLUSION

Results obtained from these studies will provide an atomic level understanding of biochemical processes occurring in AD and advance efforts to develop effective therapeutic strategies for this disease.

Natural 15N/14N isotope composition in C3 leaves: are enzymatic isotope effects informative for predicting the 15N-abundance in key metabolites?

Although nitrogen isotopes are viewed as important tools for understanding plant N acquisition and allocation, the current interpretation of natural 15N-abundances (δ15N values) is often impaired by substantial variability among individuals or between species. Such variability is likely to stem from the fact that 15N-abundance of assimilated N is not preserved during N metabolism and redistribution within the plant; that is, 14N/15N isotope effects associated with N metabolic reactions are certainly responsible for isotopic shifts between organic-N (amino acids) and absorbed inorganic N (nitrate). Therefore, to gain insights into the metabolic origin of 15N-abundance in plants, the present paper reviews enzymatic isotope effects and integrates them into a metabolic model at the leaf level. Using simple steady-state equations which satisfactorily predict the δ15N value of amino acids, it is shown that the sensitivity of δ15N values to both photorespiratory and N-input (reduction by nitrate reductase) rates is quite high. In other words, the variability in δ15N values observed in nature might originate from subtle changes in metabolic fluxes or environment-driven effects, such as stomatal closure that in turn changes v0, the Rubisco-catalysed oxygenation rate.

Effect of the active site D25N mutation on the structure, stability, and ligand binding of the mature HIV-1 protease

All aspartic proteases, including retroviral proteases, share the triplet DTG critical for the active site geometry and catalytic function. These residues interact closely in the active, dimeric structure of HIV-1 protease (PR). We have systematically assessed the effect of the D25N mutation on the structure and stability of the mature PR monomer and dimer. The D25N mutation (PR(D25N)) increases the equilibrium dimer dissociation constant by a factor >100-fold (1.3 +/- 0.09 microm) relative to PR. In the absence of inhibitor, NMR studies reveal clear structural differences between PR and PR(D25N) in the relatively mobile P1 loop (residues 79-83) and flap regions, and differential scanning calorimetric analyses show that the mutation lowers the stabilities of both the monomer and dimer folds by 5 and 7.3 degrees C, respectively. Only minimal differences are observed in high resolution crystal structures of PR(D25N) complexed to darunavir (DRV), a potent clinical inhibitor, or a non-hydrolyzable substrate analogue, Ac-Thr-Ile-Nle-r-Nle-Gln-Arg-NH(2) (RPB), as compared with PR.DRV and PR.RPB complexes. Although complexation with RPB stabilizes both dimers, the effect on their T(m) is smaller for PR(D25N) (6.2 degrees C) than for PR (8.7 degrees C). The T(m) of PR(D25N).DRV increases by only 3 degrees C relative to free PR(D25N), as compared with a 22 degrees C increase for PR.DRV, and the mutation increases the ligand dissociation constant of PR(D25N).DRV by a factor of approximately 10(6) relative to PR.DRV. These results suggest that interactions mediated by the catalytic Asp residues make a major contribution to the tight binding of DRV to PR.

Crystal structure of HIV-1 protease in situ product complex and observation of a low-barrier hydrogen bond between catalytic aspartates

HIV-1 protease is an effective target for designing drugs against AIDS, and structural information about the true transition state and the correct mechanism can provide important inputs. We present here the three-dimensional structure of a bi-product complex between HIV-1 protease and the two cleavage product peptides AETF and YVDGAA. The structure, refined against synchrotron data to 1.65 A resolution, shows the occurrence of the cleavage reaction in the crystal, with the product peptides still held in the enzyme active site. The separation between the scissile carbon and nitrogen atoms is 2.67 A, which is shorter than a normal van der Waal separation, but it is much longer than a peptide bond length. The substrate is thus in a stage just past the G'Z intermediate described in Northrop's mechanism [Northrop DB (2001) Acc Chem Res 34:790-797]. Because the products are generated in situ, the structure, by extrapolation, can give insight into the mechanism of the cleavage reaction. Both oxygens of the generated carboxyl group form hydrogen bonds with atoms at the catalytic center: one to the OD2 atom of a catalytic aspartate and the other to the scissile nitrogen atom. The latter hydrogen bond may have mediated protonation of scissile nitrogen, triggering peptide bond cleavage. The inner oxygen atoms of the catalytic aspartates in the complex are 2.30 A apart, indicating a low-barrier hydrogen bond between them at this stage of the reaction, an observation not included in Northrop's proposal. This structure forms a template for designing mechanism-based inhibitors.

An artificial aspartic proteinase system

A series of aza crown ether derivatives with or without carboxyl groups in their side arms were synthesized and the former showed deacylation activities toward amino acid p-nitrophenyl ester hydrohalides. Substrate-selective phenomena were also observed. The relationship between the structures and deacylation activities of corresponding compounds suggested a nucleophilic catalytic mechanism. The results partially simulate some aspartic proteinases in the case of catalytic mechanism and are also useful for us to understand the detailed catalytic process of aspartic proteinases.

The role of hydrogen bonding in the enzymatic reaction catalyzed by HIV-1 protease

The hydrogen-bond network in various stages of the enzymatic reaction catalyzed by HIV-1 protease was studied through quantum-classical molecular dynamics simulations. The approximate valence bond method was applied to the active site atoms participating directly in the rearrangement of chemical bonds. The rest of the protein with explicit solvent was treated with a classical molecular mechanics model. Two possible mechanisms were studied, general-acid/general-base (GA/GB) with Asp 25 protonated at the inner oxygen, and a direct nucleophilic attack by Asp 25. Strong hydrogen bonds leading to spontaneous proton transfers were observed in both reaction paths. A single-well hydrogen bond was formed between the peptide nitrogen and outer oxygen of Asp 125. The proton was diffusely distributed with an average central position and transferred back and forth on a picosecond scale. In both mechanisms, this interaction helped change the peptide-bond hybridization, increased the partial charge on peptidyl carbon, and in the GA/GB mechanism, helped deprotonate the water molecule. The inner oxygens of the aspartic dyad formed a low-barrier, but asymmetric hydrogen bond; the proton was not positioned midway and made a slightly elongated covalent bond, transferring from one to the other aspartate. In the GA/GB mechanism both aspartates may help deprotonate the water molecule. We observed the breakage of the peptide bond and found that the protonation of the peptidyl amine group was essential for the peptide-bond cleavage. In studies of the direct nucleophilic mechanism, the peptide carbon of the substrate and oxygen of Asp 25 approached as close as 2.3 A.

Molecular dynamics simulations of the first steps of the reaction catalyzed by HIV-1 protease

The mechanism of the first steps of the reaction catalyzed by HIV-1 protease was studied through molecular dynamics simulations. The potential energy surface in the active site was generated using the approximate valence bond method. The approximate valence bond (AVB) method was parameterized based on density functional calculations. The surrounding protein and explicit water environment was modeled with conventional, classical force field. The calculations were performed based on HIV-1 protease complexed with the MVT-101 inhibitor that was modified to a model substrate. The protonation state of the catalytic aspartates was determined theoretically. Possible reaction mechanisms involving the lytic water molecule are accounted for in this study. The modeled steps include the dissociation of the lytic water molecule and proton transfer onto Asp-125, the nucleophilic attack followed by a proton transfer onto peptide nitrogen. The simulations show that in the active site most preferable energetically are structures consisting of ionized or polarized molecular fragments that are not accounted for in conventional molecular dynamics. The mobility of the lytic water molecule, the dynamics of the hydrogen bond network, and the conformation of the aspartates in the active center were analyzed.

Ab initio molecular dynamics studies on HIV-1 reverse transcriptase triphosphate binding site: implications for nucleoside-analog drug resistance

Quantum-chemical methods are used to shed light on the functional role of residues involved in the resistance of HIV-1 reverse transcriptase against nucleoside-analog drugs. Ab initio molecular dynamics simulations are carried out for models representing the adduct between the triphosphate substrate and the nucleoside binding site. The triphosphate is considered either deprotonated or protonated at the gamma-position. Although the protonated form already experiences large rearrangements in the ps time scale, the fully deprotonated state exhibits a previously unrecognized low-barrier hydrogen bond between Lys65 and gamma-phosphate. Absence of this interaction in Lys65-->Arg HIV-1 RT might play a prominent role in the resistance of this mutant for nucleoside analogs (Gu Z et al., 1994b, Antimicrob Agents Chemother 38:275-281; Zhang D et al., 1994, Antimicrob Agents Chemother 38:282-287). Water molecules present in the active site, not detected in the X-ray structure, form a complex H-bond network. Among these waters, one may be crucial for substrate recognition as it bridges Gln151 and Arg72 with the beta-phosphate. Absence of this stabilizing interaction in Gln151-->Met HIV-1 RT mutant may be a key factor for the known drug resistance of this mutant toward dideoxy-type drugs and AZT (Shirasaka T et al., 1995, Proc Natl Acad Sci USA 92:2398-2402: Iversen AK et al., 1996, J Virol 70:1086-1090).

Theozymes and compuzymes: theoretical models for biological catalysis

A theozyme is a theoretical enzyme constructed by computing the optimal geometry for transition-state stabilization by functional groups. It is created in order to permit quantitative assessment of catalytic function. Theozymes have been used to elucidate the role of transition-state stabilization in the mechanisms underlying enzyme- and antibody-catalyzed hydroxyepoxide cyclizations, eliminations and decarboxylations, peptide and ester hydrolyses, and pericyclic and radical reactions. The enzymes studied include orotodine monophosphate decarboxylase, HIV protease and ribonucleotide reductase.

Synthesis of N-glyoxylyl peptides and their in vitro evaluation as HIV-1 protease inhibitors

A series of novel synthetic peptides containing an N-terminal glyoxylyl function (CHOCO-) have been tested as inhibitors of HIV-1 protease. The N-glyoxylyl peptide CHOCO-Pro-Ile-Val-NH2, which fulfills the specificity requirements of the MA/CA protease cleavage site together with the criteria of transition state analogue of the catalyzed reaction, was found to be a moderate competitive inhibitor although favorable interactions were visualized between its hydrated form and the catalytic aspartates using molecular modeling. Increasing the length of the peptide sequence led to compounds acting only as substrates.

Inhibition and catalytic mechanism of HIV-1 aspartic protease

The structure of the HIV-1 protease in complex with a pseudo-C2 symmetric inhibitor, which contains a central difluoroketone motif, has been determined with X-ray diffraction data extending to 1.7 A resolution. The electron density map clearly indicates that the inhibitor is bound in a symmetric fashion as the hydrated, or gemdiol, form of the difluoroketone. Refinement of the complex reveals a unique, and almost symmetric, set of interactions between the geminal hydroxyl groups, the geminal fluorine atoms, and the active-site aspartate residues. Several hydrogen bonding patterns are consistent with that conformation. The lowest energy hydrogen disposition, as determined by semiempirical energy calculations, shows only one active site aspartate protonated. A comparison between the corresponding dihedral angles of the difluorodiol core and those of a hydrated peptide bond analog, calculated ab-initio, shows that the inhibitor core is a mimic of a hydrated peptide bond in a gauche conformation. The feasibility of an anti-gauche transition for a peptide bond after hydration is verified by extensive molecular dynamics simulations. The simulations suggest that rotation about the C-N scissile bond would readily occur after hydration and would be driven by the optimization of the interactions of peptide side-chains with the enzyme. These results, together with the characterization of a transition state leading to bond breakage via a concerted exchange of two protons, suggest a proteolysis mechanism whereby only one active site aspartate is initially protonated. The steps of this mechanism are: asymmetric binding of the substrate; hydration of the peptidic carbonyl by an active site water; proton translocation between the active site aspartate residues simultaneously with carbonyl hydration; optimization of the binding of the entire substrate facilitated by the flexible structure of the hydrated peptide bond, which, in turn, forces the hydrated peptide bond to assume a gauche conformation; simultaneous proton exchange whereby one hydroxyl donates a proton to the charged aspartate, and, at the same time, the nitrogen lone pair accepts a proton from the other aspartate; and, bond breakage and regeneration of the initial protonation state of the aspartate residues.

Use of steady state kinetic methods to elucidate the kinetic and chemical mechanisms of retroviral proteases

Despite the current plethora of structural data of HIV-1 protease and the availability of potent inhibitors, whose structures are based in part on the presumed mechanism of action of this enzyme, our actual understanding of its chemical mechanism has been until now based largely on the precedents of the mammalian and fungal aspartic proteases and static three-dimensional data. The available steady state kinetic data of the protease, as reviewed here, constitute a first step in a detailed description of the mechanism of the enzyme to complement the structural data.